Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

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

04 Apr '12

User: rhaas Date: 2012/04/04 10:32 AM Modified: /example_wavetoy/Cactus/misc/ WaveToy.par Log: keep jpeg files around in wavetoy example par file File Changes: Directory: /example_wavetoy/Cactus/misc/ ======================================== File [modified]: WaveToy.par Delta lines: +1 -1 =================================================================== --- example_wavetoy/Cactus/misc/WaveToy.par 2012-04-04 15:24:28 UTC (rev 42) +++ example_wavetoy/Cactus/misc/WaveToy.par 2012-04-04 15:32:35 UTC (rev 43) @@ -26,7 +26,7 @@ IOJpeg::gridpoints = "interpolate" IOJpeg::out_every = 5 IOJpeg::refinement_factor = 10 -IOJpeg::mode = "remove" +IOJpeg::mode = "standard" IOJpeg::colormap = "custom" IOJpeg::colormap_min = -0.25 IOJpeg::colormap_max = +0.25

1 0

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

04 Apr '12

User: rhaas Date: 2012/04/04 10:24 AM Modified: /example_wavetoy/Cactus/misc/ WaveToy.th Log: update WaveToy thorn list File Changes: Directory: /example_wavetoy/Cactus/misc/ ======================================== File [modified]: WaveToy.th Delta lines: +5 -0 =================================================================== --- example_wavetoy/Cactus/misc/WaveToy.th 2012-04-04 14:46:16 UTC (rev 41) +++ example_wavetoy/Cactus/misc/WaveToy.th 2012-04-04 15:24:28 UTC (rev 42) @@ -117,6 +117,11 @@ !CHECKOUT = Introduction/WaveToy +!TARGET = $ROOT +!TYPE = svn +!URL = https://svn.einsteintoolkit.org/documents/Workshop_Spring_2012/example_wave… +!CHECKOUT = wavetoy-misc + # Private thorns !TARGET = $ARR !TYPE = ignore

1 0

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

04 Apr '12

User: diener Date: 2012/04/04 09:46 AM Modified: /numerical_relativity/ numerical_relativity.tex Log: Add command /mypause to easily turn on or off all pauses. Change the date. Remove the Dr's. File Changes: Directory: /numerical_relativity/ ================================= File [modified]: numerical_relativity.tex Delta lines: +51 -49 =================================================================== --- numerical_relativity/numerical_relativity.tex 2012-04-04 14:42:26 UTC (rev 40) +++ numerical_relativity/numerical_relativity.tex 2012-04-04 14:46:16 UTC (rev 41) @@ -3,12 +3,14 @@ \usecolortheme[RGB={200,200,200}]{structure} \subtitle[Introduction to numerical methods \& numerical Relativity]{{\large Einstein Toolkit Workshop}\\*[0.3em] Numerical methods for spacetime and matter evolution} -\author[\mbox{Peter Diener \& Bruno Mundim}]{Dr Peter Diener \& Dr Bruno Mundim} -\date{Apr 05 2012} +\author[\mbox{Peter Diener \& Bruno Mundim}]{Peter Diener \& Bruno Mundim} +\date{Apr 4 2012} \begin{document} \newcommand{\code}[1]{\texttt{#1}} +%\newcommand{\mypause}{\pause} +\newcommand{\mypause}{} \frame{\titlepage} @@ -34,11 +36,11 @@ \begin{align} {\cal H} & \equiv R + K^2 - K_{ij} K^{ij} = 0, \\ {\cal M}^i & \equiv D_j (K^{ij} - \gamma^{ij} K) = 0. -\end{align}\pause -These equations were used a lot in the early years of numerical relativity.\pause +\end{align}\mypause +These equations were used a lot in the early years of numerical relativity.\mypause However, it was later discovered that the ADM evolution equations are only -weakly hyperbolic.\pause +weakly hyperbolic.\mypause In the mid to late 90's a new formulation was introduced that proved to be much more robust and stable: The BSSN formulation. @@ -49,18 +51,18 @@ Introduce a conformal rescaling of the three metric\vspace{\vs} \begin{equation} \gamma_{ij} = \psi^4 \tg_{ij}.\vspace{\vs} -\end{equation}\pause +\end{equation}\mypause We choose $\psi = \gamma^{1/12}$ such that the determinant of $\tg_{ij}$ is 1%\vspace{\vs} \begin{equation} \mathrm{det}(\tg_{ij}) = \mathrm{det}(\psi^{-4}\gamma_{ij}) = \mathrm{det}(\gamma^{-1/3}\gamma_{ij}) = \gamma^{-1} \mathrm{det}(\gamma_{ij}) = 1.\vspace{\vs} -\end{equation}\pause +\end{equation}\mypause In addition we introduce a trace decomposition of the extrinsic curvature.\vspace{\vs} \begin{align} K & = \gamma^{ij} K_{ij}, \\ A_{ij} & = K_{ij} - \frac{1}{3} \gamma_{ij} K.\vspace{\vs} -\end{align}\pause +\end{align}\mypause We then promote the following variables to evolution variables\vspace{\vs} \begin{align} \phi &= \ln \psi = \frac{1}{12} \ln \gamma, @@ -78,17 +80,17 @@ \begin{equation} \tG^i = \tg^{jk} \tG^i{}_{jk} = - \partial_j \tg^{ij}, \end{equation} -to evolved variables as well.\pause +to evolved variables as well.\mypause The final set of evolution variables are $\phi$, $K$, $\tg_{ij}$, $\tA_{ij}$ -and $\tG^i$.\pause +and $\tG^i$.\mypause The BSSN evolution equations can be derived from the ADM equations. As an example take the equation for $\phi$. \begin{equation} \dt \phi = \partial_0\phi = \partial_0 \left (\frac{1}{12}\ln\gamma\right) = \frac{1}{12}\frac{1}{\gamma}\partial_0\gamma. -\end{equation}\pause +\end{equation}\mypause Using the expression for the derivative of the determinant of the metric in terms of the derivatives of the metric ($\partial_0\gamma = \gamma \gamma^{ij} \partial_0\gamma_{ij}$) we find @@ -144,18 +146,18 @@ \tilde{{\cal G}} & \equiv \tg-1 = 0, \\ \tilde{{\cal A}} & \equiv \tg^{ij}\tA_{ij} = 0, \\ \tilde{{\cal L}}^{i} & \equiv \tG^{i}+\partial_j\tg^{ij} = 0. -\end{align}\pause +\end{align}\mypause The constraints $\tilde{{\cal G}}$ and $\tilde{{\cal A}}$ are enforced actively -at each timestep.\pause +at each timestep.\mypause The other constraints ($\tilde{{\cal H}}$, $\tilde{{\cal M}}^i$ and -$\tilde{{\cal L}}^{i}$) are not enforced.\pause +$\tilde{{\cal L}}^{i}$) are not enforced.\mypause To improve stability and to help maintain $\tilde{{\cal L}}^{i}$ at a low level the following rule is employed in an implementation \begin{itemize} \item Where derivatives of $\tG^{i}$ are needed the evolved $\tG^{i}$ are used - directly.\pause + directly.\mypause \item Where $\tG^{i}$ are needed without taking derivatives $\tg^{jk}\tG^{i}{}_{jk}$ are used instead. \end{itemize} @@ -163,7 +165,7 @@ \frame{\frametitle{The BSSN formulation (continued)} In order to evolve a spacetime with the BSSN equations, you have to specify -the gauges $\alpha$ and $\beta^{i}$.\pause +the gauges $\alpha$ and $\beta^{i}$.\mypause Most codes use the ``moving puncture'' gauges \begin{align} @@ -181,12 +183,12 @@ \gamma^{\mathrm{ph}}_{ij} = \psi^4\gamma_{ab}, \mbox{\hspace{2em}} K^{\mathrm{ph}}_{ij} = \psi^{-2} K_{ij}, \end{equation} where $\gamma_{ij}$ is chosen to be the flat metric and $K_{ij}$ is assumed -to be tracefree.\pause \\ +to be tracefree.\mypause \\ The constraint equations become\vspace{\vs} \begin{align} 0 & = \Delta\psi +\frac{1}{8}K^{ij}K_{ij}\psi^{-7} \\ 0 & = D_j K^{ij}. -\end{align}\pause +\end{align}\mypause The momentum constraint has an analytic solution\vspace{\vs} \begin{align} K^{ij}_{\mathrm{BY}} = {} & \frac{3}{2r^2}\left ( P^i n^j + P^j n^i - @@ -203,29 +205,29 @@ black holes. \begin{equation} \psi = 1+\sum_{i=1}^N\frac{m_i}{2|\vec{r}-\vec{r}_i|}. -\end{equation}\pause +\end{equation}\mypause Inspired by this, for $K_{ij}\ne 0$ we make the ansatz \begin{equation} \psi=\frac{1}{\alpha}+u, \mbox{\hspace{2em}} \frac{1}{\alpha} = \sum_{i=1}^N\frac{m_i}{2|\vec{r}-\vec{r}_i|}, -\end{equation}\pause +\end{equation}\mypause In which case the Hamiltonian constraint becomes an equation for $u$ \begin{equation} \Delta u+\frac{1}{8}\alpha^7 K^{ij}K_{ij}(1+\alpha u)^{-7} = 0. -\end{equation}\pause +\end{equation}\mypause It can be shown that $u$ will be $C^2$ at the location of the `punctures' and that $\psi$ has a unique solution. } \frame{\frametitle{Finite differencing} With finite differencing we discretize a function by sampling it at a -collection of grid points.\pause +collection of grid points.\mypause -The grid points are usually (but not necessarily) equally spaced. \pause +The grid points are usually (but not necessarily) equally spaced. \mypause We can then approximate derivatives of a function at a grid point by a weigthed sum of function values at grid points in the neighbourhood (the -stencil) of the grid point. \pause +stencil) of the grid point. \mypause As an example consider a stencil containing the grid point ($f_i$) and it's two nearest neighbors ($f_{i-1}$ and $f_{i+1}$) with @@ -245,7 +247,7 @@ f(x_{i+1}) & = & f(x_i) + \left. \frac{df}{dx}\right |_{x_i} \Delta x + \frac{1}{2} \left. \frac{d^2 f}{dx^2}\right |_{x_i} \Delta x^2 + O((\Delta x)^3) \nonumber -\end{eqnarray} \pause +\end{eqnarray} \mypause Requiring that the weighted sum approximates the derivative yields the following equations for $a_{-1}$, $a_0$ and $a_1$ \begin{eqnarray} @@ -261,20 +263,20 @@ \begin{equation} \left. \frac{df}{dx} \right |_{x_i} = \frac{f_{i+1}-f_{i-1}}{2\Delta x} + O((\Delta x)^2). -\end{equation}\pause +\end{equation}\mypause Similarly we find for the second derivative that \begin{equation} \left. \frac{d^2f}{dx^2} \right |_{x_i} = \frac{f_{i+1}-2 f_i+f_{i-1}}{(\Delta x)^2} + O((\Delta x)^2). -\end{equation}\pause -These finite difference operators are second order accurate.\pause +\end{equation}\mypause +These finite difference operators are second order accurate.\mypause -Higher order accuracy or higher order derivatives require larger stencils.\pause +Higher order accuracy or higher order derivatives require larger stencils.\mypause Another way of looking at finite differencing operators is through -interpolating polynomials.\pause +interpolating polynomials.\mypause -Either approach gives the same coefficients for the same stencil.\pause +Either approach gives the same coefficients for the same stencil.\mypause It is also clear from either approach that the error estimates are only correct if $f$ is smooth enough. @@ -285,20 +287,20 @@ \begin{equation} \partial_t \mathbf{q}+\mathbf{A}^i(\mathbf{q})\partial_i \mathbf{B}(\mathbf{q}) = \mathbf{S}(\mathbf{q}). -\end{equation}\pause +\end{equation}\mypause 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 +(e.g.\ using finite differencing).\mypause This then turns the equations into a set of coupled ODE's with respect to -time.\pause +time.\mypause 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 +integrator.\mypause Often Runge-Kutta integration schemes are the scheme of choice. } @@ -307,37 +309,37 @@ 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 + order).\mypause \item It is easy to couple different evolution equations maintaining - high order coupling.\pause + high order coupling.\mypause \end{itemize} The method of lines is implemented in Cactus in the thorn -\code{CactusNumerical/MoL}.\pause +\code{CactusNumerical/MoL}.\mypause 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 + variables with \code{MoL}.\mypause + \item Schedule routines to set the RHS variables.\mypause \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 the Einstein Toolkit implementation of BSSN.\mypause \code{McLachlan} is named in honor of the Canadian Singer/Song writer Sarah -McLachlan.\pause +McLachlan.\mypause -It is based on finite differencing and the method of lines.\pause +It is based on finite differencing and the method of lines.\mypause -It supports high order finite differencing (8th order).\pause +It supports high order finite differencing (8th order).\mypause Since it is generated from the tensor equations by Kranc it is easy to -maintain and modify.\pause +maintain and modify.\mypause 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 +ongoing to make it be able to run on GPU's as well.\mypause In the future it may be necessary to generate different versions optimized for specific computer architectures. @@ -348,7 +350,7 @@ Once we have established the spacetime evolution and initial data equations, we need to obtain the evolution equations for the matter fields and the magnetic field evolution equation in -ideal MHD case. \pause +ideal MHD case. \mypause These equations can be expressed as the local conservation laws of baryon number and energy momentum. For baryon number we have: @@ -357,7 +359,7 @@ \nabla_\nu J^\nu=0, \end{equation} where $J^\mu=\rho u^\mu$ is the rest-mass current, $\rho$ the rest-mass -density and $u^\mu$ is the four-velocity of a fluid comoving observer.\pause +density and $u^\mu$ is the four-velocity of a fluid comoving observer.\mypause The conservation of energy-momentum is given by: % @@ -385,7 +387,7 @@ \end{equation} where we define $p^*=p+b^2/2$, $h^*=h+b^2/\rho$, $\varepsilon^*=\varepsilon+b^2/(2\rho)$, that results into -$h^*=1+\varepsilon^*+p^*/\rho$. \pause +$h^*=1+\varepsilon^*+p^*/\rho$. \mypause Note that in the expression above we have used the magnetic field $b^{\mu}$ as measured by the comoving observer. It can expressed

1 0

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

04 Apr '12

User: knarf Date: 2012/04/04 09:42 AM Modified: /cactuswriting/ cactuswriting.pdf /genrel/ genrel.pdf Log: change mime type of pdfs Directory Changes: Directory: /svn:mime-type/ ========================== - application/octet-stream + application/pdf File Changes: Directory: /cactuswriting/ ========================== File [modified]: cactuswriting.pdf Delta lines: +0 -0 =================================================================== Property changes on: cactuswriting/cactuswriting.pdf ___________________________________________________________________ Directory: /genrel/ =================== File [modified]: genrel.pdf Delta lines: +0 -0 =================================================================== Property changes on: genrel/genrel.pdf ___________________________________________________________________

1 0

[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

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits