Commits April 2012

commits@lists.einsteintoolkit.org

10 participants
167 discussions

[svn:einsteintoolkit] Workshop_Spring_2012/numerical_relativity/ (Rev. 39)
by diener＠cct.lsu.edu 04 Apr '12

04 Apr '12

User: diener Date: 2012/04/04 09:19 AM Modified: /numerical_relativity/ numerical_relativity.tex Log: A few more slides. File Changes: Directory: /numerical_relativity/ ================================= File [modified]: numerical_relativity.tex Delta lines: +43 -2 =================================================================== --- numerical_relativity/numerical_relativity.tex 2012-04-04 14:08:10 UTC (rev 38) +++ numerical_relativity/numerical_relativity.tex 2012-04-04 14:19:38 UTC (rev 39) @@ -8,6 +8,8 @@ \begin{document} +\newcommand{\code}[1]{\texttt{#1}} + \frame{\titlepage} \frame{\frametitle{Outline} @@ -17,6 +19,7 @@ \item Puncture data. \item Finite differencing. \item The method of lines. + \item McLachlan. \end{itemize} } @@ -96,7 +99,7 @@ } \frame{\frametitle{The BSSN formulation (continued)} -The evolution equation for all the BSSN variables are +The evolution equation for all the BSSN variables are\footnote{Note: the matter terms are left out.} \begin{align} \partial_t \tg_{ij} = {} & -2 \alpha \tA_{ij} + \lbt{\tg}{i}{j}{k}, \\ \partial_t \phi = {} & - \frac{1}{6} \alpha K + \beta^{k}\partial_{k}\phi @@ -298,8 +301,46 @@ integrator.\pause Often Runge-Kutta integration schemes are the scheme of choice. +} +\frame{\frametitle{The method of lines (continued)} +Advantages of using the method of lines: +\begin{itemize} + \item It is easy to change the time integration scheme (e.g.\ going to higher + order).\pause + \item It is easy to couple different evolution equations maintaining + high order coupling.\pause +\end{itemize} +The method of lines is implemented in Cactus in the thorn +\code{CactusNumerical/MoL}.\pause +To use: +\begin{itemize} + \item Schedule a routine to register the evolution variables and RHS + variables with \code{MoL}.\pause + \item Schedule routines to set the RHS variables.\pause + \item Set parameters when launching job to choose the time integration + scheme. +\end{itemize} +} +\frame{\frametitle{McLachlan} +\code{McLachlan} is the Einstein Toolkit implementation of BSSN.\pause + +\code{McLachlan} is named in honor of the Canadian Singer/Song writer Sarah +McLachlan.\pause + +It is based on finite differencing and the method of lines.\pause + +It supports high order finite differencing (8th order).\pause + +Since it is generated from the tensor equations by Kranc it is easy to +maintain and modify.\pause + +It is highly optimized (supports vectorization and OpenMP) and an effort is +ongoing to make it be able to run on GPU's as well.\pause + +In the future it may be necessary to generate different versions optimized +for specific computer architectures. } @@ -428,9 +469,9 @@ } \frame{\frametitle{Conservation Laws} +} - \end{document}

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/example_wavetoy/Cactus/misc/ (Rev. 38)
by roland.haas＠physics.gatech.edu 04 Apr '12

04 Apr '12

User: rhaas Date: 2012/04/04 09:08 AM Added: /example_wavetoy/Cactus/misc/ WaveToy.th Log: add WaveToy thorn list File Changes: Directory: /example_wavetoy/Cactus/misc/ ======================================== File [added]: WaveToy.th Delta lines: +124 -0 =================================================================== --- example_wavetoy/Cactus/misc/WaveToy.th (rev 0) +++ example_wavetoy/Cactus/misc/WaveToy.th 2012-04-04 14:08:10 UTC (rev 38) @@ -0,0 +1,124 @@ +# Component list for the Einstein Toolkit <http://einsteintoolkit.org/> +# $Revision$ +# $Date$ +# $HeadURL$ + +!CRL_VERSION = 1.0 + +!DEFINE ROOT = Cactus +!DEFINE ARR = $ROOT/arrangements +!DEFINE COMPONENTLIST_TARGET = $ROOT/thornlists/ + +!DEFINE ET_RELEASE = ET_2011_10 + +# This thorn list +!TARGET = $ROOT +!TYPE = svn +!AUTH_URL = https://svn.einsteintoolkit.org/$1/branches/$ET_RELEASE +!URL = http://svn.einsteintoolkit.org/$1/branches/$ET_RELEASE +!CHECKOUT = manifest + +# Cactus Flesh +!TARGET = $ROOT +!TYPE = svn +!AUTH_URL = https://svn.cactuscode.org/flesh/branches/$ET_RELEASE +!URL = http://svn.cactuscode.org/flesh/branches/$ET_RELEASE +!CHECKOUT = Cactus +!NAME = . + +# Simulation Factory +!TARGET = $ROOT +!TYPE = svn +!URL = https://svn.cct.lsu.edu/repos/numrel/$1/branches/$ET_RELEASE +!CHECKOUT = simfactory2 +!NAME = simfactory + +# Various Cactus utilities +!TARGET = $ROOT +!TYPE = svn +!AUTH_URL = https://svn.cactuscode.org/Utilities/branches/$ET_RELEASE +!URL = http://svn.cactuscode.org/Utilities/branches/$ET_RELEASE +!CHECKOUT = . +!NAME = utils + +# The GetComponents script +!TARGET = $ROOT/bin +!TYPE = git +!URL = git://github.com/gridaphobe/CRL.git +!AUTH_URL = git@github.com:gridaphobe/CRL.git +!REPO_BRANCH = $ET_RELEASE +!CHECKOUT = GetComponents + + + +# Cactus thorns +!TARGET = $ARR +!TYPE = svn +!AUTH_URL = https://svn.cactuscode.org/arrangements/$1/$2/branches/$ET_RELEASE +!URL = http://svn.cactuscode.org/arrangements/$1/$2/branches/$ET_RELEASE +!CHECKOUT = +CactusBase/Boundary +CactusBase/CartGrid3D +CactusBase/CoordBase +CactusBase/Fortran +CactusBase/InitBase +CactusBase/IOUtil +CactusBase/SymBase +CactusBase/Time + +CactusConnect/HTTPD +CactusConnect/HTTPDExtra +CactusConnect/Socket + +CactusIO/IOJpeg + +CactusNumerical/LocalInterp +CactusNumerical/LocalReduce +CactusNumerical/MoL + +CactusPUGHIO/IOHDF5 +CactusPUGHIO/IOHDF5Util + +# Additional Cactus thorns +!TARGET = $ARR +!TYPE = svn +!AUTH_URL = https://svn.cactuscode.org/projects/$1/$2/branches/$ET_RELEASE +!URL = http://svn.cactuscode.org/projects/$1/$2/branches/$ET_RELEASE +!CHECKOUT = +ExternalLibraries/GSL +ExternalLibraries/HDF5 +ExternalLibraries/libjpeg +ExternalLibraries/OpenSSL +ExternalLibraries/zlib + +# Carpet, the AMR driver +!TARGET = $ARR +!TYPE = hg +!AUTH_URL = ssh://carpetmercurial@carpetcode.org/carpet +!URL = http://www.carpetcode.org/hg/carpet +!REPO_BRANCH = $ET_RELEASE +!CHECKOUT = Carpet/doc +Carpet/Carpet +Carpet/CarpetIOASCII +Carpet/CarpetIOBasic +Carpet/CarpetIOHDF5 +Carpet/CarpetIOScalar +Carpet/CarpetInterp +Carpet/CarpetLib +Carpet/CarpetReduce +Carpet/CarpetSlab +Carpet/LoopControl + + +# WaveToy example +!TARGET = $ARR +!TYPE = svn +!URL = https://svn.einsteintoolkit.org/documents/Workshop_Spring_2012/example_wave… +!CHECKOUT = +Introduction/WaveToy + +# Private thorns +!TARGET = $ARR +!TYPE = ignore +!CHECKOUT = +

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/ (Rev. 37)
by roland.haas＠physics.gatech.edu 04 Apr '12

04 Apr '12

User: rhaas Date: 2012/04/04 08:40 AM Modified: / preamble.tex Log: re-add beamer options for color names File Changes: Directory: / ============ File [modified]: preamble.tex Delta lines: +1 -1 =================================================================== --- preamble.tex 2012-04-04 13:06:33 UTC (rev 36) +++ preamble.tex 2012-04-04 13:40:07 UTC (rev 37) @@ -1,4 +1,4 @@ -\documentclass{beamer} +\documentclass[xcolor=usenames,xcolor=dvipsnames]{beamer} \usepackage[utf8]{inputenc} \usepackage{beamerthemesplit}

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/cactuswriting/ (Rev. 36)
by jfaber＠einsteintoolkit.org 04 Apr '12

04 Apr '12

User: jfaber Date: 2012/04/04 08:06 AM Modified: /cactuswriting/ cactuswriting.aux, cactuswriting.log, cactuswriting.nav, cactuswriting.pdf, cactuswriting.tex Log: Update for toady's talks File Changes: Directory: /cactuswriting/ ========================== File [modified]: cactuswriting.aux Delta lines: +7 -5 =================================================================== --- cactuswriting/cactuswriting.aux 2012-04-04 12:52:23 UTC (rev 35) +++ cactuswriting/cactuswriting.aux 2012-04-04 13:06:33 UTC (rev 36) @@ -45,8 +45,10 @@ \@writefile{nav}{\headcommand {\beamer@framepages {15}{15}}} \@writefile{nav}{\headcommand {\slideentry {0}{0}{16}{16/16}{}{0}}} \@writefile{nav}{\headcommand {\beamer@framepages {16}{16}}} -\@writefile{nav}{\headcommand {\beamer@partpages {1}{16}}} -\@writefile{nav}{\headcommand {\beamer@subsectionpages {1}{16}}} -\@writefile{nav}{\headcommand {\beamer@sectionpages {1}{16}}} -\@writefile{nav}{\headcommand {\beamer@documentpages {16}}} -\@writefile{nav}{\headcommand {\def \inserttotalframenumber {16}}} +\@writefile{nav}{\headcommand {\slideentry {0}{0}{17}{17/17}{}{0}}} +\@writefile{nav}{\headcommand {\beamer@framepages {17}{17}}} +\@writefile{nav}{\headcommand {\beamer@partpages {1}{17}}} +\@writefile{nav}{\headcommand {\beamer@subsectionpages {1}{17}}} +\@writefile{nav}{\headcommand {\beamer@sectionpages {1}{17}}} +\@writefile{nav}{\headcommand {\beamer@documentpages {17}}} +\@writefile{nav}{\headcommand {\def \inserttotalframenumber {17}}} File [modified]: cactuswriting.log Delta lines: +48 -46 =================================================================== --- cactuswriting/cactuswriting.log 2012-04-04 12:52:23 UTC (rev 35) +++ cactuswriting/cactuswriting.log 2012-04-04 13:06:33 UTC (rev 36) @@ -1,4 +1,4 @@ -This is pdfTeX, Version 3.1415926-1.40.10 (TeX Live 2009) (format=pdflatex 2009.11.7) 3 APR 2012 23:41 +This is pdfTeX, Version 3.1415926-1.40.10 (TeX Live 2009) (format=pdflatex 2009.11.7) 4 APR 2012 09:05 entering extended mode \write18 enabled. %&-line parsing enabled. @@ -1146,133 +1146,135 @@ ] [3 -] <config.png, id=104, 773.89125pt x 165.61874pt> +] [4 + +] <config.png, id=130, 773.89125pt x 165.61874pt> File: config.png Graphic file (type png) <use config.png> -Overfull \hbox (91.32469pt too wide) in paragraph at lines 40--40 +Overfull \hbox (91.32469pt too wide) in paragraph at lines 61--61 [][] [] -[4 +[5 - <./config.png>] <interface1.png, id=133, 1108.14pt x 447.6725pt> + <./config.png>] <interface1.png, id=159, 1108.14pt x 447.6725pt> File: interface1.png Graphic file (type png) <use interface1.png> -Overfull \hbox (19.05469pt too wide) in paragraph at lines 48--48 +Overfull \hbox (19.05469pt too wide) in paragraph at lines 69--69 [][] [] -<interface2.png, id=134, 1379.1525pt x 45.16875pt> +<interface2.png, id=160, 1379.1525pt x 45.16875pt> File: interface2.png Graphic file (type png) <use interface2.png> -Overfull \hbox (40.73592pt too wide) in paragraph at lines 48--48 +Overfull \hbox (40.73592pt too wide) in paragraph at lines 69--69 [][] [] -[5 +[6 <./interface1.png> <./interface2.png>] -<interface3.png, id=163, 1233.60875pt x 249.93375pt> +<interface3.png, id=189, 1233.60875pt x 249.93375pt> File: interface3.png Graphic file (type png) <use interface3.png> -Overfull \hbox (55.18968pt too wide) in paragraph at lines 54--54 +Overfull \hbox (55.18968pt too wide) in paragraph at lines 75--75 [][] [] -[6 +[7 - <./interface3.png>] <param1.png, id=191, 1145.27875pt x 513.92pt> + <./interface3.png>] <param1.png, id=218, 1145.27875pt x 513.92pt> File: param1.png Graphic file (type png) <use param1.png> -Overfull \hbox (55.18968pt too wide) in paragraph at lines 60--60 +Overfull \hbox (55.18968pt too wide) in paragraph at lines 81--81 [][] [] -[7 +[8 - <./param1.png>] <param2.png, id=220, 1633.10126pt x 492.84125pt> + <./param1.png>] <param2.png, id=246, 1633.10126pt x 492.84125pt> File: param2.png Graphic file (type png) <use param2.png> -Overfull \hbox (55.18968pt too wide) in paragraph at lines 64--64 +Overfull \hbox (55.18968pt too wide) in paragraph at lines 85--85 [][] [] -[8 +[9 - <./param2.png>] <schedule1.png, id=248, 1020.81375pt x 391.4625pt> + <./param2.png>] <schedule1.png, id=274, 1020.81375pt x 391.4625pt> File: schedule1.png Graphic file (type png) <use schedule1.png> -Overfull \hbox (55.18968pt too wide) in paragraph at lines 71--71 +Overfull \hbox (55.18968pt too wide) in paragraph at lines 92--92 [][] [] -[9 +[10 - <./schedule1.png>] <schedule2.png, id=276, 1149.29375pt x 521.95pt> + <./schedule1.png>] <schedule2.png, id=302, 1149.29375pt x 521.95pt> File: schedule2.png Graphic file (type png) <use schedule2.png> -Overfull \hbox (55.18968pt too wide) in paragraph at lines 78--78 +Overfull \hbox (55.18968pt too wide) in paragraph at lines 99--99 [][] [] -[10 +[11 - <./schedule2.png>] <schedule3.png, id=304, 1102.1175pt x 458.71375pt> + <./schedule2.png>] <schedule3.png, id=330, 1102.1175pt x 458.71375pt> File: schedule3.png Graphic file (type png) <use schedule3.png> -Overfull \hbox (55.18968pt too wide) in paragraph at lines 83--83 +Overfull \hbox (55.18968pt too wide) in paragraph at lines 104--104 [][] [] -[11 +[12 - <./schedule3.png>] <make.png, id=332, 632.3625pt x 148.555pt> + <./schedule3.png>] <make.png, id=358, 632.3625pt x 148.555pt> File: make.png Graphic file (type png) <use make.png> -Overfull \hbox (55.18968pt too wide) in paragraph at lines 88--88 +Overfull \hbox (55.18968pt too wide) in paragraph at lines 109--109 [][] [] -[12 +[13 - <./make.png>] <srccode.png, id=360, 1349.04pt x 740.7675pt> + <./make.png>] <srccode.png, id=387, 1349.04pt x 740.7675pt> File: srccode.png Graphic file (type png) <use srccode.png> -Overfull \hbox (55.18968pt too wide) in paragraph at lines 93--93 +Overfull \hbox (55.18968pt too wide) in paragraph at lines 114--114 [][] [] -[13 +[14 - <./srccode.png>] <srccode2.png, id=389, 1117.17375pt x 221.82875pt> + <./srccode.png>] <srccode2.png, id=415, 1117.17375pt x 221.82875pt> File: srccode2.png Graphic file (type png) <use srccode2.png> -Overfull \hbox (55.18968pt too wide) in paragraph at lines 98--98 +Overfull \hbox (55.18968pt too wide) in paragraph at lines 119--119 [][] [] -[14 +[15 - <./srccode2.png>] <thornlist.png, id=417, 1569.865pt x 181.67876pt> + <./srccode2.png>] <thornlist.png, id=443, 1569.865pt x 181.67876pt> File: thornlist.png Graphic file (type png) <use thornlist.png> -Overfull \hbox (55.18968pt too wide) in paragraph at lines 103--103 +Overfull \hbox (55.18968pt too wide) in paragraph at lines 124--124 [][] [] -[15 +[16 - <./thornlist.png>] [16 + <./thornlist.png>] [17 ] \tf@nav=\write6 @@ -1291,10 +1293,10 @@ ) Here is how much of TeX's memory you used: - 14951 strings out of 493848 - 276028 string characters out of 3149207 + 14955 strings out of 493848 + 276063 string characters out of 3149207 353008 words of memory out of 3000000 - 17773 multiletter control sequences out of 15000+200000 + 17775 multiletter control sequences out of 15000+200000 6821 words of font info for 25 fonts, out of 3000000 for 9000 718 hyphenation exceptions out of 8191 47i,16n,57p,341b,534s stack positions out of 5000i,500n,10000p,200000b,50000s @@ -1303,9 +1305,9 @@ </usr/local/texlive/2009/texmf-dist/fonts/type1/public/amsfonts/cm/cmss17.pfb>< /usr/local/texlive/2009/texmf-dist/fonts/type1/public/amsfonts/cm/cmssbx10.pfb> </usr/local/texlive/2009/texmf-dist/fonts/type1/public/amsfonts/cm/cmtt10.pfb> -Output written on cactuswriting.pdf (16 pages, 406393 bytes). +Output written on cactuswriting.pdf (17 pages, 414847 bytes). PDF statistics: - 495 PDF objects out of 1000 (max. 8388607) - 33 named destinations out of 1000 (max. 500000) + 521 PDF objects out of 1000 (max. 8388607) + 35 named destinations out of 1000 (max. 500000) 122 words of extra memory for PDF output out of 10000 (max. 10000000) File [modified]: cactuswriting.nav Delta lines: +7 -5 =================================================================== --- cactuswriting/cactuswriting.nav 2012-04-04 12:52:23 UTC (rev 35) +++ cactuswriting/cactuswriting.nav 2012-04-04 13:06:33 UTC (rev 36) @@ -31,8 +31,10 @@ \headcommand {\beamer@framepages {15}{15}} \headcommand {\slideentry {0}{0}{16}{16/16}{}{0}} \headcommand {\beamer@framepages {16}{16}} -\headcommand {\beamer@partpages {1}{16}} -\headcommand {\beamer@subsectionpages {1}{16}} -\headcommand {\beamer@sectionpages {1}{16}} -\headcommand {\beamer@documentpages {16}} -\headcommand {\def \inserttotalframenumber {16}} +\headcommand {\slideentry {0}{0}{17}{17/17}{}{0}} +\headcommand {\beamer@framepages {17}{17}} +\headcommand {\beamer@partpages {1}{17}} +\headcommand {\beamer@subsectionpages {1}{17}} +\headcommand {\beamer@sectionpages {1}{17}} +\headcommand {\beamer@documentpages {17}} +\headcommand {\def \inserttotalframenumber {17}} File [modified]: cactuswriting.pdf Delta lines: +0 -0 =================================================================== (Binary files differ) File [modified]: cactuswriting.tex Delta lines: +30 -1 =================================================================== --- cactuswriting/cactuswriting.tex 2012-04-04 12:52:23 UTC (rev 35) +++ cactuswriting/cactuswriting.tex 2012-04-04 13:06:33 UTC (rev 36) @@ -12,11 +12,32 @@ \frame{\frametitle{Outline} \begin{itemize} + \item Before you begin \item Directory structure \item configuration.ccl +\item interface.ccl +\item param.ccl +\item schedule.ccl +\item Thornlist +\item source code \end{itemize} } + +\frame{\frametitle{Starting out} + +First, ask yourself what you want to do with the code: +\begin{itemize} +\item Modification of existing routine, or new routine? +\item What variables will you need? +\item Which routines can you use as possible templates? +\item What external libraries do you need? This is not a huge problem for the ET! +\item Has this already been done? +\end{itemize} + +Choose a convenient language! +} + \frame{\frametitle{Directory structure} In {\tt arrangements/ThornName} you want the following files @@ -103,7 +124,15 @@ } \frame{\frametitle{Final tips} -\LARGE Plagiarize and copy liberally, especially at first!!! +\begin{itemize} +\item \LARGE Plagiarize and copy liberally, especially at first!!! + +\item Consider compilation time when you write your thorn. You want a thornlist that is as small as possible, but not smaller. + +\item Get version control running first! + +\item Learn how to output. It is a huge timesaver for debugging/running. +\end{itemize} } \end{document}

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/genrel/ (Rev. 35)
by jfaber＠einsteintoolkit.org 04 Apr '12

04 Apr '12

User: jfaber Date: 2012/04/04 07:52 AM Modified: /genrel/ genrel.aux, genrel.log, genrel.nav, genrel.pdf, genrel.tex Log: update for today's talk File Changes: Directory: /genrel/ =================== File [modified]: genrel.aux Delta lines: +7 -5 =================================================================== --- genrel/genrel.aux 2012-04-04 12:15:48 UTC (rev 34) +++ genrel/genrel.aux 2012-04-04 12:52:23 UTC (rev 35) @@ -57,8 +57,10 @@ \@writefile{nav}{\headcommand {\beamer@framepages {21}{21}}} \@writefile{nav}{\headcommand {\slideentry {0}{0}{22}{22/22}{}{0}}} \@writefile{nav}{\headcommand {\beamer@framepages {22}{22}}} -\@writefile{nav}{\headcommand {\beamer@partpages {1}{22}}} -\@writefile{nav}{\headcommand {\beamer@subsectionpages {1}{22}}} -\@writefile{nav}{\headcommand {\beamer@sectionpages {1}{22}}} -\@writefile{nav}{\headcommand {\beamer@documentpages {22}}} -\@writefile{nav}{\headcommand {\def \inserttotalframenumber {22}}} +\@writefile{nav}{\headcommand {\slideentry {0}{0}{23}{23/23}{}{0}}} +\@writefile{nav}{\headcommand {\beamer@framepages {23}{23}}} +\@writefile{nav}{\headcommand {\beamer@partpages {1}{23}}} +\@writefile{nav}{\headcommand {\beamer@subsectionpages {1}{23}}} +\@writefile{nav}{\headcommand {\beamer@sectionpages {1}{23}}} +\@writefile{nav}{\headcommand {\beamer@documentpages {23}}} +\@writefile{nav}{\headcommand {\def \inserttotalframenumber {23}}} File [modified]: genrel.log Delta lines: +15 -13 =================================================================== --- genrel/genrel.log 2012-04-04 12:15:48 UTC (rev 34) +++ genrel/genrel.log 2012-04-04 12:52:23 UTC (rev 35) @@ -1,4 +1,4 @@ -This is pdfTeX, Version 3.1415926-1.40.10 (TeX Live 2009) (format=pdflatex 2009.11.7) 3 APR 2012 21:34 +This is pdfTeX, Version 3.1415926-1.40.10 (TeX Live 2009) (format=pdflatex 2009.11.7) 4 APR 2012 08:50 entering extended mode \write18 enabled. %&-line parsing enabled. @@ -1231,18 +1231,20 @@ File: lapseshift.png Graphic file (type png) <use lapseshift.png> [20 - <./lapseshift.png (PNG copy)>] + <./lapseshift.png (PNG copy)>] [21 -LaTeX Warning: Command \d invalid in math mode on input line 261. +] +LaTeX Warning: Command \d invalid in math mode on input line 282. -LaTeX Warning: Command \d invalid in math mode on input line 261. -[21 +LaTeX Warning: Command \d invalid in math mode on input line 282. -] <gwpol.png, id=594, 803.0pt x 401.5pt> +[22 + +] <gwpol.png, id=620, 803.0pt x 401.5pt> File: gwpol.png Graphic file (type png) - <use gwpol.png> [22 + <use gwpol.png> [23 <./gwpol.png (PNG copy)>] \tf@nav=\write6 @@ -1261,10 +1263,10 @@ ) Here is how much of TeX's memory you used: - 15010 strings out of 493848 - 276828 string characters out of 3149207 + 15014 strings out of 493848 + 276863 string characters out of 3149207 358008 words of memory out of 3000000 - 17806 multiletter control sequences out of 15000+200000 + 17808 multiletter control sequences out of 15000+200000 11753 words of font info for 43 fonts, out of 3000000 for 9000 718 hyphenation exceptions out of 8191 47i,16n,57p,830b,534s stack positions out of 5000i,500n,10000p,200000b,50000s @@ -1282,9 +1284,9 @@ al/texlive/2009/texmf-dist/fonts/type1/public/amsfonts/cm/cmsy6.pfb></usr/local /texlive/2009/texmf-dist/fonts/type1/public/amsfonts/cm/cmsy8.pfb></usr/local/t exlive/2009/texmf-dist/fonts/type1/public/amsfonts/symbols/msam10.pfb> -Output written on genrel.pdf (22 pages, 2415482 bytes). +Output written on genrel.pdf (23 pages, 2422055 bytes). PDF statistics: - 677 PDF objects out of 1000 (max. 8388607) - 45 named destinations out of 1000 (max. 500000) + 703 PDF objects out of 1000 (max. 8388607) + 47 named destinations out of 1000 (max. 500000) 102 words of extra memory for PDF output out of 10000 (max. 10000000) File [modified]: genrel.nav Delta lines: +7 -5 =================================================================== --- genrel/genrel.nav 2012-04-04 12:15:48 UTC (rev 34) +++ genrel/genrel.nav 2012-04-04 12:52:23 UTC (rev 35) @@ -43,8 +43,10 @@ \headcommand {\beamer@framepages {21}{21}} \headcommand {\slideentry {0}{0}{22}{22/22}{}{0}} \headcommand {\beamer@framepages {22}{22}} -\headcommand {\beamer@partpages {1}{22}} -\headcommand {\beamer@subsectionpages {1}{22}} -\headcommand {\beamer@sectionpages {1}{22}} -\headcommand {\beamer@documentpages {22}} -\headcommand {\def \inserttotalframenumber {22}} +\headcommand {\slideentry {0}{0}{23}{23/23}{}{0}} +\headcommand {\beamer@framepages {23}{23}} +\headcommand {\beamer@partpages {1}{23}} +\headcommand {\beamer@subsectionpages {1}{23}} +\headcommand {\beamer@sectionpages {1}{23}} +\headcommand {\beamer@documentpages {23}} +\headcommand {\def \inserttotalframenumber {23}} File [modified]: genrel.pdf Delta lines: +0 -0 =================================================================== (Binary files differ) File [modified]: genrel.tex Delta lines: +22 -1 =================================================================== --- genrel/genrel.tex 2012-04-04 12:15:48 UTC (rev 34) +++ genrel/genrel.tex 2012-04-04 12:52:23 UTC (rev 35) @@ -222,7 +222,10 @@ \end{equation} where $\rho_0$ is the rest-mass density, $\rho = \rho_0[1+\epsilon]$ is the mass-energy density, $\epsilon$ the internal energy density of a fluid, $p$ the pressure, $u^\mu$ the 4-velocity. We are assuming an ideal fluid, measured in its own rest frame, for the matrix expression. - The stress energy tensor encodes not only the density of mass and momentum, but also their flows (the pressure term). + The stress energy tensor encodes not only the density of mass and momentum, but also their flows (the pressure term). The evolution set is: + \begin{equation} + \nabla_\mu (\rho u^\mu) = 0;~~\nabla_\mu T^{\mu\nu} = 0 + \end{equation} } \frame{\frametitle{Einstein's equations} @@ -232,6 +235,11 @@ R_{\mu\nu}-\frac{1}{2}g_{\mu\nu} R = 8\pi G T_{\mu\nu} \end{equation} Put simply, the left hand side involves second derivatives, both in time (accelerations) and space (just like Gauss's law for gravitation). The right hand-side incorporates the gravitational contributions of mass, and momentum, and stresses in general. Because we have a second-order system, we need both the metric and its first derivatives to be given as initial data to evolve the system. + +The left hand side is automatically divergence free, and this fact played a critical role in deriving it in the first place: +\begin{equation} +\nabla_{\mu}(R_{\mu\nu}-\frac{1}{2}g_{\mu\nu} R ) = 0 +\end{equation} } \frame{\frametitle{3+1 formalism} @@ -245,6 +253,19 @@ \includegraphics[width=1.5in]{lapseshift} } +\frame{\frametitle{3+1 formalism} + +Much of numerical relativity comes from determining how to do the following steps: +\begin{itemize} +\item The 3+1 formalism is not numerically stable, but by evolving additional variables whose behavior we should know, it can stabilize the system: +\begin{equation} +\partial_t \gamma_{ij} = ... +\end{equation} +\item We need an evolution for the time derivative of the (spatial) metric. This is captured by the {\em extrinsic curvature} $K_{ij}$, which is closely related to the metric derivatives. Think of it as how the metric evolves from slice to slice. +\item We need gauge evolution equations for $\alpha$ and $\beta^i$. We can choose these freely, but our choice affects the resulting evolution set. This is the choice of coordinates for the spacetime. +\end{itemize} +} + \frame{\frametitle{Black holes} The most familiar black hole solution, for an uncharged BH without any angular momentum, is the Schwarzschild solution:

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/ (Rev. 34)
by knarf＠cct.lsu.edu 04 Apr '12

04 Apr '12

User: knarf Date: 2012/04/04 07:15 AM Added: /cactus_intro/ cactus_intro.pdf /hpc_intro/ hpc_intro.pdf /introduction/ introduction.pdf /pics/ supernova_model.pdf Modified: /cactus_intro/ cactus_intro.tex /hpc_intro/ hpc_intro.tex Log: some changes to introduction, commit pdfs Directory Changes: Directory: /svn:mime-type/ ========================== + application/pdf File Changes: Directory: /cactus_intro/ ========================= File [added]: cactus_intro.pdf Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: cactus_intro/cactus_intro.pdf ___________________________________________________________________ File [modified]: cactus_intro.tex Delta lines: +68 -16 =================================================================== --- cactus_intro/cactus_intro.tex 2012-04-04 07:50:13 UTC (rev 33) +++ cactus_intro/cactus_intro.tex 2012-04-04 12:15:48 UTC (rev 34) @@ -72,24 +72,60 @@ \end{itemize} } -% Copied from the web site. -\begin{frame} \frametitle{History} - \begin{itemize} \justifying - \item Cactus is a direct descendant of many years of code development in Ed -Seidel's group of researchers at the National Center for Supercomputing -Applications (NCSA) as the group wrestled to numerically solve Einstein's -Equations for General Relativity and thus model black holes, neutron stars, and -boson stars. - \item In 1995, Ed Seidel and many of his group moved to the Max Planck -Institute for Gravitational Physics (Albert Einstein Institute or AEI) in -Potsdam, Germany. Frustrated by the difficulties of coordinating the -development and use of several different codes across a large collaborative -group, Paul Walker, Joan MassÃ³, Ed Seidel, and John Shalf, designed and wrote -the original version of Cactus, Cactus 1.0, which provided a collaborative -parallel toolkit for numerical relativity. +\frame{\frametitle{Modularization Motivation} +Simulations are complex, may contain several models at once.\\ +Example: Long Gamma Ray Burst: + \begin{columns}[t] + \column{.4\textwidth} + \begin{center} + \includegraphics[width=0.9\textwidth]{supernova_model} + \end{center} + \column{.5\textwidth} + \begin{itemize} + \item General Relativity (black hole) + \item Relativistic Hydrodynamics (star) + \item Microphysics, equation of state (shock wave) + \item Neutrino radiation (cooling, heating) + \item Magnetic Fields (Jet formation) + \item Photon radiation (afterglow) + \end{itemize} + + \end{columns} +} + +\frame{\frametitle{Modularization Motivation} + Added Problems: + \begin{itemize} + \item Example: Einstein Toolkit (not untypical) + \item Code 12+ years old + \item Grad students leave after 3 productive years + \item Most original authors not available anymore + \item Developers distributed over many places in several continents + \item Most physicists are not good programmers \end{itemize} -\end{frame} +} +\frame{\frametitle{Component Architecture} + \begin{itemize} + \item Split program into independent components (as much as possible) + \item \textit{Framework} provides glue between these + \item Components are developed independently by small groups + \item End user assembles all code: no central control, no authorative version + \item Framework itself does no real work + \item Components don't interact directly with each other + \end{itemize} +} + +\frame { \frametitle{Cactus History} + \begin{itemize} + \item Direct descendant of many years of code development in Ed Seidel's group of + researchers at NCSA + \item 1995, Paul Walker, Joan Masso, Ed Seidel, and John Shalf: Cactus~1.0 + \item Originally for numerical relativity + \item Over the years generalized for use by scientists in other domains + \end{itemize} +} + \frame { \frametitle{Current Users and Developers} \centering \includegraphics[height=7cm]{groups_map}\\ @@ -224,6 +260,22 @@ \end{columns} \end{frame} +\frame { \frametitle{Typical list of component tasks} + \abspic{puzzle}{0.0}{0.6}{0.4} + \begin{itemize} + \item Evolution systems + \item Boundary conditions + \item Initial conditions + \item Time stepping method + \item Finite Difference methods + \item Simulation grid (distributed arrays) + \item I/O (more on next slide) + \item Simulation domain specification + \item Termination condition + \item Twitter client + \end{itemize} +} + \frame { \frametitle{I/O Capabilities} Usual I/O and checkpointing in different formats: \begin{columns}[t] Directory: /hpc_intro/ ====================== File [added]: hpc_intro.pdf Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: hpc_intro/hpc_intro.pdf ___________________________________________________________________ File [modified]: hpc_intro.tex Delta lines: +31 -0 =================================================================== --- hpc_intro/hpc_intro.tex 2012-04-04 07:50:13 UTC (rev 33) +++ hpc_intro/hpc_intro.tex 2012-04-04 12:15:48 UTC (rev 34) @@ -249,4 +249,35 @@ \end{itemize} } +\frame{\frametitle{Allocations} +\begin{itemize} + \item Supercomputer time is expensive. + \item 1 CPU hour costs about 5 US-cents + \item 100k cores: \$5,000/h + \item Allocation: Amount of computer time a group of + users can use within some time (a year typically) + \item How to get allocation: proposal +\end{itemize} +} + +\frame{\frametitle{Batch Processing} +Simulations take too long +\begin{itemize} + \item Cannot be awake at all times + \item Each user error can destroy weeks of data + \item SC Expensive + \begin{itemize} + \item cannot wait for user input + \item should be utilized close to 100\% + \end{itemize} +\end{itemize} +$\rightarrow$ +\begin{itemize} + \item Need to plan simulations carefully ahead of time + \item Then let them run automatically + \item Often don't start right away + \item Might need to be \textit{checkpointed} +\end{itemize} +} + \end{document} Directory: /introduction/ ========================= File [added]: introduction.pdf Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: introduction/introduction.pdf ___________________________________________________________________ Directory: /pics/ ================= File [added]: supernova_model.pdf Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: pics/supernova_model.pdf ___________________________________________________________________

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 04 Apr '12

04 Apr '12

changeset: 3505:ab6a712c7eb6 tag: tip user: Ian Hinder <ian.hinder(a)aei.mpg.de> date: Wed Apr 04 13:14:38 2012 +0200 summary: Make timer tree output precision and threshold configurable with parameters

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/numerical_relativity/ (Rev. 33)
by bcmsma＠astro.rit.edu 04 Apr '12

04 Apr '12

User: bmundim Date: 2012/04/04 02:50 AM Added: /numerical_relativity/figs/ all_waves.eps, all_waves.pdf, finite_volume.eps, finite_volume.pdf, rarefaction.eps, rarefaction.pdf, shock.eps, shock.pdf Modified: /numerical_relativity/ numerical_relativity.tex Log: valencia formulation plus figs for conservation laws... (must sleep nw) File Changes: Directory: /numerical_relativity/figs/ ====================================== File [added]: all_waves.eps Delta lines: <Unavailable> File [added]: all_waves.pdf Delta lines: <Unavailable> File [added]: finite_volume.eps Delta lines: <Unavailable> File [added]: finite_volume.pdf Delta lines: <Unavailable> File [added]: rarefaction.eps Delta lines: <Unavailable> File [added]: rarefaction.pdf Delta lines: <Unavailable> File [added]: shock.eps Delta lines: <Unavailable> File [added]: shock.pdf Delta lines: <Unavailable> Directory: /numerical_relativity/ ================================= File [modified]: numerical_relativity.tex Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/cactuswriting/ (Rev. 32)
by jfaber＠einsteintoolkit.org 03 Apr '12

03 Apr '12

User: jfaber Date: 2012/04/03 10:43 PM Added: /cactuswriting/ cactuswriting.aux, cactuswriting.log, cactuswriting.nav, cactuswriting.out, cactuswriting.pdf, cactuswriting.snm, cactuswriting.tex, cactuswriting.toc, config.png, config.tiff, interface1.png, interface1.tiff, interface2.png, interface2.tiff, interface3.png, interface3.tiff, make.png, make.tiff, param1.png, param1.tiff, param2.png, param2.tiff, schedule1.png, schedule1.tiff, schedule2.png, schedule2.tiff, schedule3.png, schedule3.tiff, srccode.png, srccode.tiff, srccode2.png, srccode2.tiff, thornlist.png, thornlist.tiff Log: Added notes on writing a cactus thorn (low level) File Changes: Directory: /cactuswriting/ ========================== File [added]: cactuswriting.aux Delta lines: <Unavailable> File [added]: cactuswriting.log Delta lines: <Unavailable> File [added]: cactuswriting.nav Delta lines: <Unavailable> File [added]: cactuswriting.out Delta lines: <Unavailable> File [added]: cactuswriting.pdf Delta lines: <Unavailable> File [added]: cactuswriting.snm Delta lines: <Unavailable> File [added]: cactuswriting.tex Delta lines: <Unavailable> File [added]: cactuswriting.toc Delta lines: <Unavailable> File [added]: config.png Delta lines: <Unavailable> File [added]: config.tiff Delta lines: <Unavailable> File [added]: interface1.png Delta lines: <Unavailable> File [added]: interface1.tiff Delta lines: <Unavailable> File [added]: interface2.png Delta lines: <Unavailable> File [added]: interface2.tiff Delta lines: <Unavailable> File [added]: interface3.png Delta lines: <Unavailable> File [added]: interface3.tiff Delta lines: <Unavailable> File [added]: make.png Delta lines: <Unavailable> File [added]: make.tiff Delta lines: <Unavailable> File [added]: param1.png Delta lines: <Unavailable> File [added]: param1.tiff Delta lines: <Unavailable> File [added]: param2.png Delta lines: <Unavailable> File [added]: param2.tiff Delta lines: <Unavailable> File [added]: schedule1.png Delta lines: <Unavailable> File [added]: schedule1.tiff Delta lines: <Unavailable> File [added]: schedule2.png Delta lines: <Unavailable> File [added]: schedule2.tiff Delta lines: <Unavailable> File [added]: schedule3.png Delta lines: <Unavailable> File [added]: schedule3.tiff Delta lines: <Unavailable> File [added]: srccode.png Delta lines: <Unavailable> File [added]: srccode.tiff Delta lines: <Unavailable> File [added]: srccode2.png Delta lines: <Unavailable> File [added]: srccode2.tiff Delta lines: <Unavailable> File [added]: thornlist.png Delta lines: <Unavailable> File [added]: thornlist.tiff Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/numerical_relativity/ (Rev. 31)
by diener＠cct.lsu.edu 03 Apr '12

03 Apr '12

User: diener Date: 2012/04/03 10:41 PM Modified: /numerical_relativity/ numerical_relativity.tex Log: First slide on MoL. File Changes: Directory: /numerical_relativity/ ================================= File [modified]: numerical_relativity.tex Delta lines: +30 -5 =================================================================== --- numerical_relativity/numerical_relativity.tex 2012-04-04 01:40:55 UTC (rev 30) +++ numerical_relativity/numerical_relativity.tex 2012-04-04 03:41:18 UTC (rev 31) @@ -15,8 +15,8 @@ % \item 3+1 spacetime decomposition. \item The BSSN formulation. \item Puncture data. - \item Finite Differencing. -% \item The method of lines. + \item Finite differencing. + \item The method of lines. \end{itemize} } @@ -214,7 +214,7 @@ that $\psi$ has a unique solution. } -\frame{\frametitle{Finite Differencing} +\frame{\frametitle{Finite differencing} With finite differencing we discretize a function by sampling it at a collection of grid points.\pause @@ -231,7 +231,7 @@ \left. \frac{df}{dx} \right |_{x_i}\approx \frac{1}{\Delta x}\sum_{j=-1}^{j=1} a_j f_{i+j}. \end{equation} } -\frame{\frametitle{Finite Differencing (continued)} +\frame{\frametitle{Finite differencing (continued)} The coefficients $a_j$ can be found be expanding $f$ in a Taylor series around $x_i$ for the grid points in the stencil \begin{eqnarray} @@ -253,7 +253,7 @@ with the solution $a_{-1}=-1/2, a_0=0, a_1=1/2$. } -\frame{\frametitle{Finite Differencing (continued)} +\frame{\frametitle{Finite differencing (continued)} Thus we find that \begin{equation} \left. \frac{df}{dx} \right |_{x_i} = \frac{f_{i+1}-f_{i-1}}{2\Delta x} @@ -276,4 +276,29 @@ It is also clear from either approach that the error estimates are only correct if $f$ is smooth enough. } + +\frame{\frametitle{The method of lines} +Consider the set of hyperbolic PDE's +\begin{equation} +\partial_t \mathbf{q}+\mathbf{A}^i(\mathbf{q})\partial_i +\mathbf{B}(\mathbf{q}) = \mathbf{S}(\mathbf{q}). +\end{equation}\pause +The idea then is to discretize in space first, i.e.\ write the equations as +\begin{equation} +\partial_t \mathbf{q} = \mathbf{L}(\mathbf{q}), +\end{equation} +where $\mathbf{L}(\mathbf{q})$ is a discrete approximation to the equations +(e.g.\ using finite differencing).\pause + +This then turns the equations into a set of coupled ODE's with respect to +time.\pause + +If the spatial discretization (including boundary conditions) is stable +we we can then evolve the system of equations using any stable ODE time +integrator.\pause + +Often Runge-Kutta integration schemes are the scheme of choice. + + +} \end{document}

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits April 2012