Commits March 2011

commits@lists.einsteintoolkit.org

6 participants
48 discussions

[svn:einsteintoolkit] manifest/trunk/ (Rev. 76)
by knarf＠cct.lsu.edu 29 Mar '11

29 Mar '11

User: knarf Date: 2011/03/29 09:07 PM Modified: /trunk/ einsteintoolkit.th Log: don't use carpetcode.dyndns.org File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.th Delta lines: +5 -5 =================================================================== --- trunk/einsteintoolkit.th 2011-02-14 18:12:08 UTC (rev 75) +++ trunk/einsteintoolkit.th 2011-03-30 02:07:30 UTC (rev 76) @@ -261,8 +261,8 @@ # McLachlan, the spacetime code !TARGET = $ARR !TYPE = git -!AUTH_URL = carpetgit@carpetcode.dyndns.org:McLachlan -!URL = git://carpetcode.dyndns.org/McLachlan +!AUTH_URL = carpetgit@carpetcode.org:McLachlan +!URL = git://carpetcode.org/McLachlan !REPO_PATH= $2 !CHECKOUT = McLachlan/doc McLachlan/m McLachlan/par McLachlan/ML_BSSN @@ -278,12 +278,12 @@ # Carpet, the AMR driver !TARGET = $ARR !TYPE = git -!AUTH_URL = carpetgit@carpetcode.dyndns.org:carpet -!URL = git://carpetcode.dyndns.org/carpet +!AUTH_URL = carpetgit@carpetcode.org:carpet +!URL = git://carpetcode.org/carpet # For experimental version, comment previous three lines and uncomment next two # Note that you will need to get access first (ask Erik Schnetter) #!TYPE = hg -#!URL = ssh://carpetmercurial@carpetcode.dyndns.org/carpet +#!URL = ssh://carpetmercurial@carpetcode.org/carpet !CHECKOUT = Carpet/doc Carpet/Carpet Carpet/CarpetEvolutionMask

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 226)
by cott＠tapir.caltech.edu 28 Mar '11

28 Mar '11

User: cott Date: 2011/03/28 09:59 PM Modified: /trunk/src/ GRHydro_Con2Prim.F90 Log: * fix an issue with hot EOS and inner symmetry boundaries File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: +2 -3 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2011-03-28 15:15:34 UTC (rev 225) +++ trunk/src/GRHydro_Con2Prim.F90 2011-03-29 02:59:13 UTC (rev 226) @@ -541,9 +541,8 @@ (x .lt. -1.0d-12 .or.& y .lt. -1.0d-12 .or.& z .lt. -1.0d-12)) then - return -! failwarnmode = 2 -! failinfomode = 2 + failwarnmode = 2 + failinfomode = 2 endif !!$ Undensitize the variables

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 468)
by knarf＠cct.lsu.edu 28 Mar '11

28 Mar '11

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 467)
by knarf＠cct.lsu.edu 28 Mar '11

28 Mar '11

User: knarf Date: 2011/03/28 03:10 PM Modified: /release-info/ run_tests Log: update runscript File Changes: Directory: /release-info/ ========================= File [modified]: run_tests Delta lines: +17 -8 =================================================================== --- release-info/run_tests 2011-03-28 20:05:57 UTC (rev 466) +++ release-info/run_tests 2011-03-28 20:10:29 UTC (rev 467) @@ -1,11 +1,11 @@ #!/bin/bash # -# philip qsub -I -l nodes=1:ppn=8 -l walltime=01:00:00 -q workq -# queenbee qsub -I -l walltime=01:00:00 -l nodes=1:ppn=8 -A loni_numrel04 -q workq +# philip qsub -I -l nodes=1:ppn=8 -l walltime=02:00:00 -q workq +# queenbee qsub -I -l walltime=02:00:00 -l nodes=1:ppn=8 -A loni_numrel05 -q workq set -e # Abort on errors -cd ET_2010_11 +cd ET_2011_05 export CCTK_TESTSUITE_RUN_COMMAND="mpirun -np \$nprocs \$exe \$parfile" hostname | grep ^philip && @@ -15,9 +15,18 @@ hostname | grep macbook-pro.local && export CCTK_TESTSUITE_RUN_COMMAND="openmpirun -np \$nprocs \$exe \$parfile" -export OMP_NUM_THREADS=1 -for i in 1 2; do - export CCTK_TESTSUITE_RUN_PROCESSORS=$i - make sim-testsuite PROMPT=no - mv TEST/sim/summary.log ../${HOSTNAME}__${i}.log +for i in 1 2 3; do + case $i in + 1 ) + export CCTK_TESTSUITE_RUN_PROCESSORS=1 + export OMP_NUM_THREADS=1;; + 2 ) + export CCTK_TESTSUITE_RUN_PROCESSORS=2 + export OMP_NUM_THREADS=1;; + 3 ) + export CCTK_TESTSUITE_RUN_PROCESSORS=1 + export OMP_NUM_THREADS=2;; + esac + echo make sim-testsuite PROMPT=no + echo mv TEST/sim/summary.log ../${HOSTNAME}__${CCTK_TESTSUITE_RUN_PROCESSORS}_${OMP_NUM_THREADS}.log done

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 466)
by knarf＠cct.lsu.edu 28 Mar '11

28 Mar '11

User: knarf Date: 2011/03/28 03:05 PM Modified: /release-info/ compile_tests Log: update compile script File Changes: Directory: /release-info/ ========================= File [modified]: compile_tests Delta lines: +11 -3 =================================================================== --- release-info/compile_tests 2011-03-09 16:11:15 UTC (rev 465) +++ release-info/compile_tests 2011-03-28 20:05:57 UTC (rev 466) @@ -1,19 +1,27 @@ #! /bin/bash set -e # Abort on errors -CROOT=ET_2010_11 +CROOT=ET_2011_05 +#LOC=branches/$CROOT +LOC=trunk # Ensure that Cactus does not already exist test ! -e $CROOT -wget http://svn.cactuscode.org/Utilities/branches/$CROOT/Scripts/GetComponents -O GetComponents -wget http://svn.einsteintoolkit.org/manifest/branches/$CROOT/einsteintoolkit.th -O einsteintoolkit.th +wget --no-check-certificate https://github.com/gridaphobe/CRL/raw/master/GetComponents -O GetComponents +wget http://svn.einsteintoolkit.org/manifest/$LOC/einsteintoolkit.th -O einsteintoolkit.th chmod a+x GetComponents time ./GetComponents -root $CROOT -a einsteintoolkit.th cd $CROOT cp simfactory/udb.example.pm simfactory/udb.pm + +# 'configure' simfactory hostname | grep ^topf && cp ~/Cactus/simfactory/udb.pm simfactory/udb.pm + +# main build ./simfactory/sim build sim --thornlist=../einsteintoolkit.th + +# additional builds hostname | grep ^numrel && ./simfactory/sim build gcc --thornlist=../einsteintoolkit.th --optionlist=numrel-gcc.cfg

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 225)
by cott＠tapir.caltech.edu 28 Mar '11

28 Mar '11

User: cott Date: 2011/03/28 10:15 AM Removed: /trunk/src/ GRHydro_HLLE_EOSG.F90, GRHydro_HLLE_EOSGM.F90, GRHydro_RegisterGZ.cc, GRHydro_RegisterGZM.cc Modified: /trunk/ interface.ccl, param.ccl, schedule.ccl /trunk/src/ GRHydro_Con2Prim.F90, GRHydro_Con2PrimM.F90, GRHydro_Marquina.F90, GRHydro_Prim2Con.F90, GRHydro_Prim2ConM.F90, GRHydro_Reconstruct.F90, GRHydro_ReconstructM.F90, GRHydro_ReconstructPoly.F90, GRHydro_ReconstructPolyM.F90, GRHydro_RiemannSolve.F90, GRHydro_RiemannSolveM.F90, GRHydro_RoeSolver.F90, GRHydro_Scalars.F90, GRHydro_Startup.F90, make.code.defn Log: * remove support for "General" EOS interface (EOSG_*) * test suites pass with Intel 11 (no optimizaton) on bethe and gcc 4.4.5 File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Con2PrimM.F90 Delta lines: <Unavailable> File [removed]: GRHydro_HLLE_EOSG.F90 Delta lines: <Unavailable> File [removed]: GRHydro_HLLE_EOSGM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Marquina.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Prim2Con.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Prim2ConM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Reconstruct.F90 Delta lines: <Unavailable> File [modified]: GRHydro_ReconstructM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_ReconstructPoly.F90 Delta lines: <Unavailable> File [modified]: GRHydro_ReconstructPolyM.F90 Delta lines: <Unavailable> File [removed]: GRHydro_RegisterGZ.cc Delta lines: <Unavailable> File [removed]: GRHydro_RegisterGZM.cc Delta lines: <Unavailable> File [modified]: GRHydro_RiemannSolve.F90 Delta lines: <Unavailable> File [modified]: GRHydro_RiemannSolveM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_RoeSolver.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Scalars.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Startup.F90 Delta lines: <Unavailable> File [modified]: make.code.defn Delta lines: <Unavailable> Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: <Unavailable> File [modified]: param.ccl Delta lines: <Unavailable> File [modified]: schedule.ccl Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 41)
by roland.haas＠physics.gatech.edu 28 Mar '11

28 Mar '11

User: rhaas Date: 2011/03/28 09:56 AM Modified: / ET.tex Log: remove undefined command call File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +0 -1 =================================================================== --- ET.tex 2011-03-28 14:54:05 UTC (rev 40) +++ ET.tex 2011-03-28 14:56:57 UTC (rev 41) @@ -904,7 +904,6 @@ \subsection{Simulation Domain, Symmetries, Boundaries\pages{3 Roland}} -\unknowasdag \paragraph{Domains and Coordinates.} The Einstein Toolkit provides a module \codename{CoordBase} to facilate the specification of the simulation domain independend from the actual

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 40)
by roland.haas＠physics.gatech.edu 28 Mar '11

28 Mar '11

User: rhaas Date: 2011/03/28 09:54 AM Modified: / ET.tex Log: add some text to Simulation Domain section mostly only name the thorns involved, what services they offer File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +38 -5 =================================================================== --- ET.tex 2011-03-28 13:38:16 UTC (rev 39) +++ ET.tex 2011-03-28 14:54:05 UTC (rev 40) @@ -903,9 +903,12 @@ -\subsection{Simulation Domain, Symmetries, Boundaries\pages{3 Roland?}} - +\subsection{Simulation Domain, Symmetries, Boundaries\pages{3 Roland}} +\unknowasdag \paragraph{Domains and Coordinates.} +The Einstein Toolkit provides a module \codename{CoordBase} to facilate +the specification of the simulation domain independend from the actual +evolution module used. The simulation domain is specified via the parameter file \todo{ES: ensure this has been defined above}. Cactus distinguishes between the \emph{physical} domain, which lives in the continuum, and @@ -924,23 +927,53 @@ the numerical resolution do not affect the extent of the physical domain, i.e.\ that the discrete domains converge to the physical domain in the limint of infinite resolution. +\codename{CoordBase} exposes a public interface that allows other +modules to query the domain and description in a uniform way. This +information is used e.g.\ by modules computing the total barionic mass of +the system to determine the region over which to intergrate the mass +density. Evolution modules use this information to decide which points +are evolved and therefore require the evaluation of the righ-hand-side +expression and which ones are set via boundary or symmetry conditions. \paragraph{Symmetries, Boundary Conditions.} +The Einstein Toolkit includes a set of modules, \codename{Boundary} +and \codename{SymBase}, to provide a generic interface to boundary and +symmetry conditions. The Einstein Toolkit offers a set of reflecting or rotating symmetry conditions that can be used to reduce the size of the simulation -domain. These symmetries include periodicity, reflections about the -coordinate planes, and $90^{\circ}$ and $180^{\circ}$ rotational +domain. These symmetries include periodicity +(via the \codename{Periodic} module), reflections about the +coordinate planes (via the \codename{Reflection} module), +and $90^{\circ}$ and $180^{\circ}$ rotational +(via the \codename{RotatingSymmetry90} and +\codename{RotatingSymmetry180} modules respectively) symmetries about the $z$ axis. In applying these symmetries, the transformation behaviour of tensorial quantities (including tensor densities and non-tensors such as Christoffel symbols) is correctly taken into account. - +The interpolation routines present in \codename{Cactus} automatically +take the existing symmetries into account, transparently mapping points +to the numerical domain and respecting the transformation behaviour of +tensorial quantities. Symmetries therefore are handled transparently +from the point of view of user modules. The Einstein Toolkit also offers a set of boundary conditions, including Dirichlet, von Neumann, and Robin (fall-off) boundary conditions, as well as extrapolation and radiative boundary conditions. +\comment{RH: Christian already describes NewRad's boundary conditions in +his section (these should be used for extrapolation and radiative +boundary conditions I think). Should they eventually be moved here (once +the paper is almost ready)?} \paragraph{Adaptive Mesh Refinement} +The \codename{Carpet} driver module provides adaptive mesh refinement and +Multi-Block capabilities to modules in the Einstein Toolkit. As in most +adaptive mesh refinement frameworks, care has been taken to separate +driver related code from physics codes. A typical evolution or analysis +code like \codename{McLachlan} or \codename{WeylScal4} is agnostic to +the mesh refinement. +\todo{RH: add description and motivation of/for Berger-Rigoustos scheme?} +\todo{RH: add interaction with MoL?} Carpet supports feature-based mesh refinement, which is based on extracting e.g.\ the locations of black holes or neutron stars, and then constructing a mesh hierarchy (stacks of refined regions) based

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 39)
by roland.haas＠physics.gatech.edu 28 Mar '11

28 Mar '11

User: rhaas Date: 2011/03/28 08:38 AM Added: / cuted.sty, flushend.sty Log: add missing style files required by widetext.sty missing at least on RH5 machines (from /macros/latex/contrib/sttools). File Changes: Directory: / ============ File [added]: cuted.sty Delta lines: +474 -0 =================================================================== --- cuted.sty (rev 0) +++ cuted.sty 2011-03-28 13:38:16 UTC (rev 39) @@ -0,0 +1,474 @@ +%% cuted.sty +%% Copyright 1997 Sigitas Tolu\v sis +%% VTeX Ltd., Akademijos 4, Vilnius, Lithuania +%% e-mail sigitas(a)vtex.lt +%% http://www.vtex.lt/tex/download/macros/ +%% +% This program can redistributed and/or modified under the terms +% of the LaTeX Project Public License Distributed from CTAN +% archives in directory macros/latex/base/lppl.txt; either +% version 1 of the License, or (at your option) any later version. +% +% PURPOSE: Insert some onecolumn material in full-width on double-column page. +% +% SHORT DESCRIPTION: +% +% It is denoted to place some material in full-width at any place +% on double-column page (\twocolumn mode). The part of material +% can be placed on the next page. Existing text is placed +% under inserted material in two columns (it is reformatted). +% +% Provides new environment: +% +% \begin{strip} ... \end{strip} +% ----------------------------- +% +% Puts \vbox{...} in full-width on double-column page. +% Adds \stripsep glue before and after +% (By default \stripsep=15pt plus2pt minus 2pt). +% It can be problems with floats and footnotes. In the case of +% problems with footnotes it is possible to correct the situation +% using commands \footnotemark and \footnotetext. +% +% Provides new command: +% \oldcolsbreak{#1} +% +% Adds #1 in place of original column break if strip is placed +% on the right column. +% +% \changes{1997/09/02}{improved column balancing} +% \changes{1997/09/03}{float control for \@colroom setting} +% \changes{1997/09/04}{improved \@colroom setting} +% \changes{1997/09/05}{checking for \@viper and \hold@viper} +% \changes{1997/09/09}{compatibility with flushend.sty} +% \changes{1997/09/22}{improved \output ir \@specialoutput} +% \changes{1997/10/01}{improved \output ir \@addviper} +% +\NeedsTeXFormat{LaTeX2e} +\ProvidesPackage{cuted}[1997/10/01] +% +\newtoks\At@ViperColsBreak \At@ViperColsBreak{} +\def\oldcolsbreak#1{\global\At@ViperColsBreak{#1}} +\newbox\v@aaa +\newbox\v@ccc +\@ifundefined{stripsep}{\newskip\stripsep\stripsep 15pt plus 2pt minus 2pt}{} +% +\@ifundefined{@viper}{\newbox\@viper}{} +\newif\if@lastviper \@lastviperfalse +\def\@outputdblcol{% + \if@firstcolumn + \global \@firstcolumnfalse + \global \setbox\@leftcolumn \box\@outputbox + \else + \global \@firstcolumntrue + \setbox\@outputbox \vbox {% + \hb@xt@\textwidth {% + \hb@xt@\columnwidth {% + \box\@leftcolumn \hss}% + \hfil + \vrule \@width\columnseprule + \hfil + \hb@xt@\columnwidth {% + \box\@outputbox \hss}% + }% + }% + \ifvoid\hold@viper + \else + \setbox\@outputbox \vbox{\box\hold@viper\box\@outputbox}% + \fi + \@combinedblfloats + \@outputpage + \begingroup + \@dblfloatplacement + \@startdblcolumn + \@whilesw\if@fcolmade \fi + {\@outputpage + \@startdblcolumn}% + \ifvoid\@viper + \else + \global\setbox\@viper\vbox{% + \vskip-\stripsep\unvbox\@viper}\@viperoutput + \fi + \endgroup + \fi +} +\mathchardef\@Mv=10005 +\newcount\viper@penalty +\def\strip{% + \ifinner + \@parmoderr\viper@penalty\z@ + \else + \viper@penalty=-\@Mv\relax + \fi + \global \setbox\@viper\vbox\bgroup% + \hsize\textwidth + \@parboxrestore + \col@number \@ne + \vrule height\topskip width0pt depth0pt +} +\def\endstrip{% + \egroup + \ifnum\viper@penalty <-\@Mii + \penalty -\@Miv + \@tempdima\prevdepth %% saving and restoring \prevdepth added + \vbox{}% %% 26 May 87 to prevent extra vertical + \prevdepth \@tempdima %% space when used in vertical mode + \penalty-\@Mv\relax + \fi +} +\gdef\@viperoutput{% + \global \setbox\@holdpg \vbox{% + \unvbox\@holdpg + \unvbox\@cclv + \setbox\@tempboxa \lastbox + \unskip + }% + \@pagedp \dp\@holdpg + \@pageht \ht\@holdpg + \ifvoid\@viper\else + \advance \@pageht \@pagedp + \ifvoid\footins \else + \advance \@pageht \ht\footins + \advance \@pageht \skip\footins + \advance \@pageht \dp\footins + \fi + \ifvbox \@kludgeins + \ifdim \wd\@kludgeins=\z@ + \advance \@pageht \ht\@kludgeins + \fi + \fi + \@reinserts + \@addviper + \fi +} +\gdef\@specialoutput{% + \ifnum \outputpenalty>-\@Mii + \@doclearpage + \else + \ifnum \outputpenalty<-\@Miii + \ifnum \outputpenalty<-\@MM \deadcycles \z@ \fi + \ifnum \outputpenalty=-10005 + \@viperoutput + \ifnum \outputpenalty<\z@ + \if@nobreak + \else + \addpenalty \interlinepenalty + \fi + \fi + \else + \global \setbox\@holdpg \vbox {\unvbox\@cclv}% + \fi + \else + \global \setbox\@holdpg \vbox{% + \unvbox\@holdpg + \unvbox\@cclv + \setbox\@tempboxa \lastbox + \unskip + }% + \@pagedp \dp\@holdpg + \@pageht \ht\@holdpg + \unvbox \@holdpg + \@next\@currbox\@currlist{% + \ifnum \count\@currbox>\z@ + \advance \@pageht \@pagedp + \ifvoid\footins \else + \advance \@pageht \ht\footins + \advance \@pageht \skip\footins + \advance \@pageht \dp\footins + \fi + \ifvbox \@kludgeins + \ifdim \wd\@kludgeins=\z@ + \advance \@pageht \ht\@kludgeins + \fi + \fi + \@reinserts + \@addtocurcol + \else + \@reinserts + \@addmarginpar + \fi + }\@latexbug + \ifnum \outputpenalty<\z@ + \if@nobreak + \else + \addpenalty \interlinepenalty + \fi + \fi + \fi + \fi +} +\@ifundefined{hold@viper}{\newbox\hold@viper}{} +\newdimen\ht@viper +\def \@addviper {% + \if@firstcolumn + \@insertfalse + \@reqcolroom .5\@pageht + \ht@viper\@reqcolroom + \advance \@reqcolroom \ht\@viper + \advance \@reqcolroom \stripsep + \ifvoid\hold@viper\@tempdima=\z@\else\@tempdima\ht\hold@viper\fi +% + \@tempdimb \ht@viper + \loop + \setbox\v@aaa\copy\@holdpg + \setbox\v@ccc\vbox to\@tempdimb{% + \vsplit\v@aaa to\@tempdimb\vss\vsplit\v@aaa to\@tempdimb}% + \wlog{Extra height:\the\ht\v@aaa\space when \the\@tempdimb}% + \ifvoid\v@aaa \else \advance\@tempdimb 1pt \repeat% + \loop + \setbox\v@aaa\copy\@holdpg + \setbox\v@ccc\vbox to\@tempdimb{% + \vsplit\v@aaa to\@tempdimb\vss}% + \wlog{(2)Left:\the\ht\v@ccc\space Right:\the\ht\v@aaa\space Output:\the\@tempdimb}% + \ifdim \ht\v@ccc<\ht\v@aaa \@tempdimb \the\ht\v@aaa \repeat% + \wlog{- LAST -^^JLeft:\the\ht\v@ccc^^JRight:\the\ht\v@aaa^^JOutput:\the\@tempdimb}% + \setbox\v@ccc\vbox to\@tempdimb{% + \vsplit\@holdpg to\@tempdimb\vss}% +% + \setbox\@tempboxa \vbox {% + \hb@xt@\textwidth {% + \hb@xt@\columnwidth {% +% \vbox to\ht@viper{\vsplit\@holdpg to\ht@viper\vss} \hss}% + \vbox to\ht@viper{\box\v@ccc\vss} \hss}% + \hfil + \vrule \@width\columnseprule + \hfil + \hb@xt@\columnwidth {% + \vbox to\ht@viper{\box\@holdpg\vss} \hss}% + }% + }% + \ifvoid\hold@viper + \global\setbox\hold@viper \vbox {% + \unvbox\@tempboxa + \penalty\interlinepenalty + \vskip\stripsep + \ifnum\outputpenalty <-\@Mii \vskip -\parskip\fi + \outputpenalty \z@ + }% + \else + \global\setbox\hold@viper \vbox {% + \unvbox\hold@viper + \unvbox\@tempboxa + \penalty\interlinepenalty + \vskip\stripsep + \ifnum\outputpenalty <-\@Mii \vskip -\parskip\fi + \outputpenalty \z@ + }% + \fi + \ifdim \@colroom>\@reqcolroom + \global\setbox\hold@viper \vbox {% + \unvbox\hold@viper + \box\@viper + \penalty\interlinepenalty + \vskip\stripsep + \ifnum\outputpenalty <-\@Mii \vskip -\parskip\fi + \outputpenalty \z@ + }% + \@inserttrue + \advance\@tempdima by-\ht\hold@viper + \global \advance \@colht \@tempdima + \ifdim \@colht<2.5\baselineskip + \@latex@warning@no@line {Optional argument of \noexpand\twocolumn + too tall on page \thepage (1 stulpelis)}% + \@emptycol + \if@firstcolumn + \else + \@emptycol + \fi + \global\@lastvipertrue + \else + \global \vsize \@colht + \global \@colroom \@colht + \check@floats + \ifdim\@tempdima>\z@ + \global\advance\@colroom by-\@tempdima + \fi + \ifdim\@colroom<\z@ + \global\@colroom 1pt + \ifdim\@tempdima>\@colht + \@latex@warning@no@line {TOP ir/ar BOT FLOAT's + netelpa lape \thepage (1 stulpelis)}% + \fi + \@emptycol + \if@firstcolumn + \else + \@emptycol + \fi + \global\@lastvipertrue + \fi + \fi + \fi + \else + \@insertfalse + \@reqcolroom\@colht + \advance\@reqcolroom \@pageht + \ht@viper.5\@reqcolroom + \@reqcolroom\ht@viper + \advance \@reqcolroom \ht\@viper + \advance \@reqcolroom \stripsep + \global\setbox\@holdpg\vbox{\unvbox\@leftcolumn + \setbox\@tempboxa\lastbox\unskip\the\At@ViperColsBreak\unvbox\@holdpg}% + \global\At@ViperColsBreak{}% + \ht@viper .5\ht\@holdpg + \ifvoid\hold@viper\@tempdima=\z@\else\@tempdima\ht\hold@viper\fi +% + \@tempdimb \ht@viper + \loop + \setbox\v@aaa\copy\@holdpg + \setbox\v@ccc\vbox to\@tempdimb{% + \vsplit\v@aaa to\@tempdimb\vss\vsplit\v@aaa to\@tempdimb}% + \wlog{Extra height:\the\ht\v@aaa\space when \the\@tempdimb}% + \ifvoid\v@aaa \else \advance\@tempdimb 1pt \repeat% + \loop + \setbox\v@aaa\copy\@holdpg + \setbox\v@ccc\vbox to\@tempdimb{% + \vsplit\v@aaa to\@tempdimb\vss}% + \wlog{(2)Left:\the\ht\v@ccc\space Right:\the\ht\v@aaa\space Output:\the\@tempdimb}% + \ifdim \ht\v@ccc<\ht\v@aaa \@tempdimb \the\ht\v@aaa \repeat% + \wlog{- LAST -^^JLeft:\the\ht\v@ccc^^JRight:\the\ht\v@aaa^^JOutput:\the\@tempdimb}% + \setbox\v@ccc\vbox to\@tempdimb{% + \vsplit\@holdpg to\@tempdimb\vss}% +% + \setbox\@tempboxa\vbox{% + \hb@xt@\textwidth {% + \hb@xt@\columnwidth {% +% \vbox to\ht@viper{\vsplit\@holdpg to \ht@viper\vss} \hss}% + \vbox to\ht@viper{\box\v@ccc\vss} \hss}% + \hfil + \vrule \@width\columnseprule + \hfil + \hb@xt@\columnwidth {% + \vbox to\ht@viper{\box\@holdpg\vss} \hss}% + }% + }% + \ifvoid\hold@viper + \global\setbox\hold@viper \vbox {% + \unvbox\@tempboxa + \penalty\interlinepenalty + \vskip\stripsep + \ifnum\outputpenalty <-\@Mii \vskip -\parskip\fi + \outputpenalty \z@ + }% + \else + \global\setbox\hold@viper \vbox {% + \unvbox\hold@viper + \unvbox\@tempboxa + \penalty\interlinepenalty + \vskip\stripsep + \ifnum\outputpenalty <-\@Mii \vskip -\parskip\fi + \outputpenalty \z@ + }% + \fi + \ifdim \@colht>\@reqcolroom + \global\setbox\hold@viper \vbox {% + \unvbox\hold@viper + \box\@viper + \penalty\interlinepenalty + \vskip\stripsep + \ifnum\outputpenalty <-\@Mii \vskip -\parskip\fi + \outputpenalty \z@ + }% + \@inserttrue + \global\@firstcolumntrue + \advance\@tempdima by-\ht\hold@viper + \global \advance \@colht \@tempdima + \ifdim \@colht<2.5\baselineskip + \@latex@warning@no@line {Optional argument of \noexpand\twocolumn + too tall on page \thepage (2 stulpelis)}% + \@emptycol + \if@firstcolumn + \else + \@emptycol + \fi + \global\@lastvipertrue + \else + \global \vsize \@colht + \global \@colroom \@colht + \check@floats + \ifdim\@tempdima>\z@ + \global\advance\@colroom by-\@tempdima + \fi + \ifdim\@colroom<\z@ + \global\@colroom 1pt + \ifdim\@tempdima>\@colht + \@latex@warning@no@line {TOP ir/ar BOT FLOAT's + netelpa lape \thepage (2 stulpelis)}% + \fi + \@emptycol + \if@firstcolumn + \else + \@emptycol + \fi + \global\@lastvipertrue + \fi + \fi + \fi + \fi + \if@insert + \else + \advance\@tempdima by-\ht\hold@viper + \global \advance \@colht \@tempdima + \global\setbox\hold@viper \vbox {% + \unvbox\hold@viper + \vsplit\@viper to \@colht + \vss + }% + \vbox{}% + \global\@firstcolumnfalse + \newpage + \fi +} + +\def \check@floats {% + \global\@tempdima\z@ + \let \@elt \ht@comflelt + \ifx \@toplist\@empty \else\@toplist \fi + \ifx \@botlist\@empty \else\@botlist \fi +} +\def\ht@comflelt#1{\setbox\@tempboxa\vbox{\copy #1\vskip\floatsep}% +\global\advance\@tempdima\ht\@tempboxa} + +\output {% + \let \par \@@par + \ifnum \outputpenalty<-\@M + \@specialoutput + \else + \@makecol + \@opcol + \@startcolumn + \@whilesw \if@fcolmade \fi + {% + \@opcol\@startcolumn}% + \fi + \ifnum \outputpenalty>-\@Miv + \ifdim \@colroom<1.5\baselineskip + \ifdim \@colroom<\textheight + \@latex@warning@no@line {Text page \thepage\space + contains only floats}% + \if@lastviper\else\@emptycol\fi + \if@twocolumn + \if@firstcolumn + \else + \@emptycol + \fi + \fi + \global\@lastviperfalse + \else + \global \vsize \@colroom + \fi + \else + \global \vsize \@colroom + \fi + \else + \ifnum\outputpenalty<-100005 + \global \vsize \maxdimen + \else + \global \vsize \@colroom + \fi + \fi +} +\endinput + + + + File [added]: flushend.sty Delta lines: +141 -0 =================================================================== --- flushend.sty (rev 0) +++ flushend.sty 2011-03-28 13:38:16 UTC (rev 39) @@ -0,0 +1,141 @@ +%% flushend.sty +%% Copyright 1997 Sigitas Tolu\v sis +%% VTeX Ltd., Akademijos 4, Vilnius, Lithuania +%% e-mail sigitas(a)vtex.lt +%% http://www.vtex.lt/tex/download/macros/ +%% +% This program can redistributed and/or modified under the terms +% of the LaTeX Project Public License Distributed from CTAN +% archives in directory macros/latex/base/lppl.txt; either +% version 1 of the License, or (at your option) any later version. +% +% PURPOSE: Balanced columns on last page in twocolumn mode. +% +% SHORT DESCRIPTION: +% +% \flushend (loaded by default) +% --------- +% Switches on column balancing at last page +% +% \raggedend +% ---------- +% Switches off column balancing at last page +% +% \atColsBreak={#1} +% ------------------ +% Adds #1 in place of original column break (without balancing) +% Example: \atColsBreak{\vskip-2pt} +% +% \showcolsendrule +% ---------------- +% Adds rule to the bottom of columns (just for debugging) +% +% P.S. To stretch right column by #1 add command \vskip-#1 just before +% command \end{document}. +% TO shrink right column by #1 add command \vskip#1 just before +% command \end{document}. +% Example: \vskip-10pt +% \end{document} +% +% \changes{1997/05/16}{first version} +% \changes{1997/09/09}{support for compatibility with cuted.sty} +% \changes{1997/10/01}{\vipersep changed to \stripsep for compatibility with cuted.sty} +% +\NeedsTeXFormat{LaTeX2e} +\ProvidesPackage{flushend}[1997/10/01] +% +\newbox\@aaa +\newbox\@ccc +\@ifundefined{@viper}{\newbox\@viper}{} +\@ifundefined{hold@viper}{\newbox\hold@viper}{} +\newtoks\atColsBreak \atColsBreak={} +\newdimen\@extra@skip \@extra@skip\z@ +\newdimen\@nd@page@rule \@nd@page@rule\z@ +\def\last@outputdblcol{% + \if@firstcolumn + \global \@firstcolumnfalse + \global \setbox\@leftcolumn \box\@outputbox + \else + \global \@firstcolumntrue + \if@lastpage + \@tempdima\ht\@leftcolumn + \splittopskip\topskip\splitmaxdepth\maxdepth + \setbox\@tempboxa\vbox{% + \unvbox\@leftcolumn\setbox0\lastbox\unskip% + \the\atColsBreak% + \unvbox\@outputbox\setbox0\lastbox\unskip}% + \@tempdimb .5\ht\@tempboxa% + \loop + \setbox\@aaa\copy\@tempboxa% + \setbox\@ccc\vbox to\@tempdimb{% + \vsplit\@aaa to\@tempdimb\vss\vsplit\@aaa to\@tempdimb}% + \wlog{Extra height:\the\ht\@aaa\space when \the\@tempdimb}% + \ifvoid\@aaa \else \advance\@tempdimb 1pt \repeat% + \loop + \setbox\@aaa\copy\@tempboxa% + \setbox\@ccc\vbox to\@tempdimb{% + \vsplit\@aaa to\@tempdimb\vss}% + \wlog{(2)Left:\the\ht\@ccc\space Right:\the\ht\@aaa\space Output:\the\@tempdimb}% + \ifdim \ht\@ccc<\ht\@aaa \@tempdimb \the\ht\@aaa \repeat% + \wlog{- LAST -^^JExtra skip:\the\@extra@skip^^JLeft:\the\ht\@ccc^^JRight:\the\ht\@aaa^^JOutput:\the\@tempdimb}% + \setbox\@ccc\vbox to\@tempdimb{% + \vsplit\@tempboxa to\@tempdimb\vss}% + \setbox\@leftcolumn\vbox to\@tempdima{% + \vbox to\@tempdimb{\unvbox\@ccc}% + \hrule\@height\@nd@page@rule% + \vss}% + \setbox\@outputbox\vbox to\@tempdima{% + \vbox to\@tempdimb{\unvbox\@tempboxa\vfilneg\vskip\@extra@skip}% + \hrule\@height\@nd@page@rule% + \vss}% + \setbox\@outputbox \vbox {% + \hb@xt@\textwidth {% + \hb@xt@\columnwidth {% + \box\@leftcolumn \hss}% + \hfil + \vrule \@width\columnseprule + \hfil + \hb@xt@\columnwidth {% + \box\@outputbox \hss}% + }% + }% + \else + \setbox\@outputbox \vbox {% + \hb@xt@\textwidth {% + \hb@xt@\columnwidth {% + \box\@leftcolumn \hss}% + \hfil + \vrule \@width\columnseprule + \hfil + \hb@xt@\columnwidth {% + \box\@outputbox \hss}% + }% + }% + \fi + \ifvoid\hold@viper + \else + \setbox\@outputbox \vbox{\box\hold@viper\box\@outputbox}% + \fi + \@combinedblfloats + \@outputpage + \begingroup + \@dblfloatplacement + \@startdblcolumn + \@whilesw\if@fcolmade \fi + {\@outputpage + \@startdblcolumn}% + \ifvoid\@viper + \else + \global\setbox\@viper\vbox{% + \vskip-\stripsep\unvbox\@viper}\@viperoutput + \fi + \endgroup + \fi +} +\let\prev@enddocument\enddocument +\newif\if@lastpage \@lastpagefalse +\def\enddocument{\global\@lastpagetrue\let\@outputdblcol\last@outputdblcol\prev@enddocument} +\def\raggedend{\global\let\enddocument\prev@enddocument} +\def\flushend{\gdef\enddocument{\global\@lastpagetrue\let\@outputdblcol\last@outputdblcol\prev@enddocument}} +\def\showcolsendrule{\global\@nd@page@rule=.4pt} +\endinput

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 38)
by cott＠tapir.caltech.edu 26 Mar '11

26 Mar '11

User: cott Date: 2011/03/26 09:28 PM Modified: / ET.tex Log: * more evolution text File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +19 -11 =================================================================== --- ET.tex 2011-03-26 15:14:54 UTC (rev 37) +++ ET.tex 2011-03-27 02:28:19 UTC (rev 38) @@ -553,6 +553,9 @@ Alcubierre~\cite{alcubierre:08} and Baumgarte \& Shapiro~\cite{baumgarte:10}. GR hydrodynamics has been reviewed by Font~\cite{font:08}. +\todo{state Einstein equations, mention conventions, give stress +energy tensor} + \subsubsection{Spacetime Curvature Evolution} The Einstein Toolkit curvature evolution code \codename{McLachlan}~\cite{brown:09} is auto-generated from tensor equations via \codename{Kranc} @@ -561,26 +564,31 @@ \cite{york:79}. For this, the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) conformal-tracefree reformulation \cite{shibata:95,baumgarte:95,alcubierre:00} of the original -Arnowitt-Deser-Misner formalism~\cite{adm:62} is employed. +Arnowitt-Deser-Misner (ADM) formalism~\cite{adm:62} is employed. -This leads to the following set of evolved variables: -\todo{continue editing here} +The evolved variables are the conformal factor $\Phi$, the conformal +3-metric $\tilde{gamma}_{ij}$, the trace $K$ of the extrinsic curvature, +the trace free extrinsic curvature $A_{ij}$ and the conformal connection +functions $\tilde{\Gamma}^i$. These are defined in terms of the +standard ADM 4-metric $g_{ij}$, 3-metric $\gamma{ij}$, and extrinsic +curvature $K_{ij}$ by + \begin{eqnarray} - \phi & := & \log \left[ \frac{1}{12} \det \gamma_{ij} \right] + \phi & := & \log \left[ \frac{1}{12} \det \gamma_{ij} \right]\,, \\ - \tilde\gamma_{ij} & := & e^{-4\phi}\; \gamma_{ij} + \tilde\gamma_{ij} & := & e^{-4\phi}\; \gamma_{ij}\,, \\ - K & := & g^{ij} K_{ij} + K & := & g^{ij} K_{ij}\,, \\ - \tilde A_{ij} & := & e^{-4\phi} \left[ K_{ij} - \frac{1}{3} g_{ij} K + \tilde A_{ij} & := & e^{-4\phi} \left[ K_{ij} - \frac{1}{3} g_{ij} K\,, \right] \\ \tilde\Gamma^i & := & \tilde\gamma^{jk} \tilde\Gamma^i_{jk} . \end{eqnarray} -Our exact evolution equations are as described by Eqs.~(3) to (10) of -\cite{ES-Brown2007b}, which we list here for completeness: + +The evolution equations are then: \begin{widetext} \begin{eqnarray} \partial_0 \alpha & = & -\alpha^2 f(\alpha, \phi, x^\mu) (K - @@ -661,13 +669,13 @@ variables are lapse $\alpha$, shift $\beta^i$, and a quantity $B^i$ related to the time derivative of the shift. The gauge parameters $f$, $G$, $H$, and $\eta$ are determined by our choice of a $1+\log$ -slicing: +\cite{alcubierre:03a} slicing: \begin{eqnarray} f(\alpha,\phi,x^\mu) & := & 2/\alpha \\ K_0(x^\mu) & := & 0 \end{eqnarray} -and $\Gamma$-driver shift condition: +and $\Gamma$-driver shift condition \cite{alcubierre:03a}: \begin{eqnarray} G(\alpha,\phi,x^\mu) & := & (3/4)\, \alpha^{-2} \\

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits March 2011