Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 82)
by roland.haas＠physics.gatech.edu 08 May '11

08 May '11

User: rhaas Date: 2011/05/08 12:23 AM Modified: / ET.tex Log: work a bit on language in simulation domain section File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +57 -36 =================================================================== --- ET.tex 2011-04-28 19:13:37 UTC (rev 81) +++ ET.tex 2011-05-08 05:23:51 UTC (rev 82) @@ -1597,11 +1597,12 @@ \subsection{Simulation Domain, Symmetries, Boundaries\pages{3 Roland}} \subsubsection{Domains and Coordinates.} -The Einstein Toolkit provides a module \codename{CoordBase} to facilate +The Einstein Toolkit provided thorn \codename{CoordBase} facilates 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 +evolution thorn used. +The simulation domain is specified at run time via a parameter file at the +same time as parameters describing the physical system are specified. +Cactus distinguishes between the \emph{physical} domain, which lives in the continuum, and \emph{discrete} domain, which consists of a discrete set of grid points. The physical domain is defined via its coordinate extent and @@ -1618,33 +1619,39 @@ 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. -\codename{Carpet} queries \codename{CoordBase} for the discrete -grid when creating the hierarchy of grids. Evolution modules use this +\codename{CoordBase} exposes a public runtime interface that allows +other +thorns to query the domain description in a uniform way. This is used +by +\codename{Carpet} to query \codename{CoordBase} for the discrete +grid when creating the hierarchy of grids, automatically ensuring a +consistent grid description between the two thorns. +Evolution thorns such as \codename{McLachlan} use the domain +information 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. \subsubsection{Symmetries, Boundary Conditions.} -The Einstein Toolkit includes a set of modules, \codename{Boundary} -and \codename{SymBase}, to provide a generic interface to boundary and +The Einstein Toolkit includes a set of thorns, \codename{Boundary} +and \codename{SymBase}, to provide a generic interface to specify and +implement boundary and symmetry conditions. -The toolkit includes built-in suppport for a set of reflecting or +For symmetry conditions the toolkit includes built-in suppport for a +set of reflecting or rotating symmetry conditions that can be used to reduce the size of the simulation domain. These symmetries include periodicity in any of the coordinate directions -(via the \codename{Periodic} module), reflections about the +(via the \codename{Periodic} module), reflections across the coordinate planes (via the \codename{Reflection} module), -and $90^{\circ}$ and $180^{\circ}$ rotational +$90^{\circ}$ and $180^{\circ}$ rotational (via the \codename{RotatingSymmetry90} and \codename{RotatingSymmetry180} modules respectively) -symmetries about the $z$ axis. -\codename{Cartoon2D} provides a continuous -rotational symmetry rather than a discrete -symmetry~\cite{Alcubierre1999ab}. \codename{Cartoon2D} allows fully +symmetries about the $z$ axis, and a continuous rotational symmetry (via +the \codename{Cartoon2D} thorn)~\cite{Alcubierre1999ab}. +\codename{Cartoon2D} allows fully three dimensional codes to be used in axissymmtric problems by evolving a slice in the $y=0$ plane and using the rotational symmetry to populate ghost points off the plane (see Figure~\ref{fig:cartoon-plane}). @@ -1658,15 +1665,16 @@ \label{fig:cartoon-plane} \end{figure} -In applying symmetries, the +In applying symmetries to populate ghost zones, 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 +Similarly 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 are handled transparently +tensorial quantities. Thus symmetries are handled transparently from the point of view of user modules (see Figure~\ref{fig:faces} for an illustration). \begin{figure}[htbp] @@ -1682,16 +1690,20 @@ \label{fig:faces} \end{figure} -Thorn \codename{Boundary} provides basic boundary conditions. A boundary +Thorn \codename{Boundary} provides basic physical boundary conditions. +A boundary condition suitable for matter fields which approach a constant -(atmosphere) value is provided via either the ``flat'' or ``scalar'' +(atmosphere) value on the boundary is provided via either the ``flat'' +or ``scalar'' (Dirichlet) boundary conditions. These boundary conditions set the values on the boundary to either those close by in the direction of the -normal to the boundary or to a fixed given value, respectively. +normal to the boundary (``flat'') or to a fixed given value +(``scalar''). The ``copy'' boundary condition uses a secondary grid variable to fill in the boundary values, allowing arbitrary boundary conditions to be imposed. -A fall-off (``Robin'') boundary condition is provided suitable for field +\codename{Boundary} provides a fall-off (``Robin'') boundary condition +suitable for fields whose behaviour near the boundary is described by \begin{equation} f(r) = f_0 + \frac{k}{r^n}\text{,} @@ -1709,25 +1721,31 @@ the paper is almost ready)?} \subsubsection{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 +The \codename{Carpet} driver thorn provides adaptive mesh refinement +optionally using subcycling in time and level-dependend spatial +refinement factors as well as +Multi-Block capabilities to thorns 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. -Carpet supports feature-based mesh refinement, which is based on +\codename{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 on the locations, sizes, and speeds of these objects. This allows tracking objects as they move through the domain. One can also add or -remove stacks if e.g.\ the number of objects changes. If initial +remove stacks if e.g.\ the number of objects changes. Fully AMR based on +a local error estimate is supported by \codename{Carpet} however the +Einstein Toolkit does not provide a suitable redgridding thorn to create +such a grid. If initial conditions are constructed outside of Carpet (which is often the case), then the initial mesh hierarchy has to be defined manually. -The Einstein toolkit currently provides two ``regridding'' modules in +The Einstein toolkit currently provides two regridding modules in the \codename{CarpetRegrid} and \codename{CarpetRegrid2} thorns. Both thorns primarily support box-in-box type refined meshes, which is well suited to current binary black hole simulations in which the -high-resolution region is centered on the individual black holes. +high-resolution regions are centered on the individual black holes. Figure~\ref{fig:bbh-boxes} shows a typical set of nested boxes during the inspiral phase of a binary black hole merger simulation. \begin{figure}[htbp] @@ -1742,18 +1760,20 @@ \codename{RotatingSymmetry180} to reduce the computational domain.} \label{fig:bbh-boxes} \end{figure} + \codename{CarpetRegrid} provides a number of different ways to specify the refined regions, either as a set of boxes centered around the origin or as an explicit list of boxes that make up the grid hierarchy. Traditionally groups using \codename{CarpetRegrid} have employed auxiliary thorns, that are not part of the Einstein Toolkit, to create -the list of boxes based on information obtained e.g. from apparent +this list of boxes based on information obtained e.g.\ from apparent horizon tracking. \codename{CarpetRegrid2} provides a user friendly interface to define sets of nested boxes that follow black holes or other tracked objects. \codename{CarpetRegrid2} supports up to 10 sets -of nested regions which can be either moving or stationary. The number -of refinement levels in each region as well as the radii of each nested -box are allowed to change during runtime, making it possible to e.g. +of nested boxes which can be either moving or stationary. The number +of refinement levels in each set of boxes as well as the radii of each +nested +box are allowed to change during runtime, making it possible to e.g.\ adapt the shape of the refined region to the surface of a star. \codename{CarpetRegrid2} contains code to handle the $\pi$-symmetry provided by \codename{RotatingSymmetry180}, enforcing the symmetry on @@ -1777,12 +1797,13 @@ only as a source for the boundary condition.} \label{fig:rot180-grid} \end{figure} + \codename{CarpetTracker} provides a simple interface to slave \codename{CarpetRegrid2}'s regions to the object trackers \codename{PunctureTracker} and \codename{NSTracker} (see -section~\ref{sec:object-tracking}) by copying the position of tracker +section~\ref{sec:object-tracking}) by copying the position of tracked objects -positions out of a \codename{SphericalSurface} into +from a \codename{SphericalSurface} into \codename{CarpetRegrid2}.

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 587)
by knarf＠cct.lsu.edu 06 May '11

06 May '11

User: knarf Date: 2011/05/06 11:38 PM Modified: /about/members/ members.txt Log: add David Garrison File Changes: Directory: /about/members/ ========================== File [modified]: members.txt Delta lines: +3 -0 =================================================================== --- about/members/members.txt 2011-05-07 04:31:31 UTC (rev 586) +++ about/members/members.txt 2011-05-07 04:38:08 UTC (rev 587) @@ -94,6 +94,9 @@ Universidad Nacional AutÃ³noma de MÃ©xico Jose Manuel Torres +University of Houston Clear Lake + David Garrison + University of Maryland Bruno Giacomazzo

1 0

[svn:einsteintoolkit] www/documentation/ (Rev. 586)
by knarf＠cct.lsu.edu 06 May '11

06 May '11

User: knarf Date: 2011/05/06 11:31 PM Modified: /documentation/ request_email Log: reformat File Changes: Directory: /documentation/ ========================== File [modified]: request_email Delta lines: +14 -13 =================================================================== --- documentation/request_email 2011-05-07 04:20:25 UTC (rev 585) +++ documentation/request_email 2011-05-07 04:31:31 UTC (rev 586) @@ -1,11 +1,13 @@ Thank you for your request for an Einstein Toolkit demonstration account. Please go to <https://allocations.loni.org/login_request.php> -and request a LONI account. Enter "Frank Loeffler" as your LONI -Contact/Collaborator. Enter "Einstein Toolkit" as your research area. -You do not need to enter a Globus DN. +and request a LONI account. Use an email address connected to your +afficilation, independent providers like yahoo or google will not be +accepted. Enter "Frank Loeffler" as your LONI Contact/Collaborator. +Enter "Einstein Toolkit" as your research area. You do not need to +enter a Globus DN. -(If you read XX here, this email was not properly created, it should - be two numbers!) +(If you read XX here, this email was not properly created, it should be two +numbers!) Please use the user name 'ettestXX'. This account is only meant for testing the Einstein Toolkit and will be removed within one to two months, without warning, and possibly including all stored data. @@ -13,9 +15,9 @@ Once you have received your account (hopefully within one or two days), go to <https://allocations.loni.org/balances.php> and log in. Select "Request Allocation" on the right hand side, and then "Join -Allocation". Enter "Frank Loeffler" and then select the -"loni_cactus06" allocation. This will be approved within one or two -days, and you will receive an email confirming the allocation. +Allocation". Enter "Frank Loeffler" and then select the "loni_cactus06" +allocation. This will be approved within one or two days, and you will +receive an email confirming the allocation. Note that this allocation is reserved for running the examples for the Einstein Toolkit, please limit your total usage to 2 kSU. If you need @@ -26,9 +28,8 @@ and files. The user guide for Queen Bee is at <http://www.loni.org/teragrid/users_guide.php>. -If you continue to use the Einstein Toolkit, please join our -consortium by signing up at -<http://einsteintoolkit.org/about/members/join>. Please also write to -<users(a)einsteintoolkit.org> if you have problems or questions about -the examples. +If you continue to use the Einstein Toolkit, please join our consortium +by signing up at <http://einsteintoolkit.org/about/members/join>. +Please also write to <users(a)einsteintoolkit.org> if you have problems +or questions about the examples.

1 0

[svn:einsteintoolkit] www/documentation/ (Rev. 585)
by knarf＠cct.lsu.edu 06 May '11

06 May '11

User: knarf Date: 2011/05/06 11:20 PM Added: /documentation/ request_email Log: template for request answer for loni account File Changes: Directory: /documentation/ ========================== File [added]: request_email Delta lines: +34 -0 =================================================================== --- documentation/request_email (rev 0) +++ documentation/request_email 2011-05-07 04:20:25 UTC (rev 585) @@ -0,0 +1,34 @@ +Thank you for your request for an Einstein Toolkit demonstration +account. Please go to <https://allocations.loni.org/login_request.php> +and request a LONI account. Enter "Frank Loeffler" as your LONI +Contact/Collaborator. Enter "Einstein Toolkit" as your research area. +You do not need to enter a Globus DN. + +(If you read XX here, this email was not properly created, it should + be two numbers!) +Please use the user name 'ettestXX'. This account is only meant for +testing the Einstein Toolkit and will be removed within one to two +months, without warning, and possibly including all stored data. + +Once you have received your account (hopefully within one or two +days), go to <https://allocations.loni.org/balances.php> and log in. +Select "Request Allocation" on the right hand side, and then "Join +Allocation". Enter "Frank Loeffler" and then select the +"loni_cactus06" allocation. This will be approved within one or two +days, and you will receive an email confirming the allocation. + +Note that this allocation is reserved for running the examples for the +Einstein Toolkit, please limit your total usage to 2 kSU. If you need +more resources once you reach 2 kSU, feel free to contact the toolkit +maintainers for suggestions on where you might find a real allocation +for your work. We also may remove you from the allocation if you have +used over 2 kSU, but you will still be able to access all your data +and files. The user guide for Queen Bee is at +<http://www.loni.org/teragrid/users_guide.php>. + +If you continue to use the Einstein Toolkit, please join our +consortium by signing up at +<http://einsteintoolkit.org/about/members/join>. Please also write to +<users(a)einsteintoolkit.org> if you have problems or questions about +the examples. +

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 584)
by knarf＠cct.lsu.edu 06 May '11

06 May '11

User: knarf Date: 2011/05/06 01:47 PM Modified: /about/members/ members.txt Log: 'move' Erik File Changes: Directory: /about/members/ ========================== File [modified]: members.txt Delta lines: +1 -1 =================================================================== --- about/members/members.txt 2011-04-21 21:00:46 UTC (rev 583) +++ about/members/members.txt 2011-05-06 18:47:25 UTC (rev 584) @@ -48,7 +48,6 @@ Oleg Korobkin Frank LÃ¶ffler Megan Miller - http://www.cct.lsu.edu/~eschnett/ Erik Schnetter Edward Seidel Jian Tao Michael Thomas @@ -67,6 +66,7 @@ http://www.perimeterinstitute.ca/ Perimeter Institute David Rideout + http://www.cct.lsu.edu/~eschnett/ Erik Schnetter http://ccrg.rit.edu/ Rochester Institute of Technology Manuela Campanelli

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 128)
by roland.haas＠physics.gatech.edu 29 Apr '11

29 Apr '11

User: rhaas Date: 2011/04/29 07:22 PM Modified: /trunk/src/ GRHydro_ShockTubeM.F90 Log: GRHydro_InitialData: apply boundary conditions to Bcons for diagonal shocks File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_ShockTubeM.F90 Delta lines: +9 -9 =================================================================== --- trunk/src/GRHydro_ShockTubeM.F90 2011-04-28 01:14:57 UTC (rev 127) +++ trunk/src/GRHydro_ShockTubeM.F90 2011-04-30 00:22:45 UTC (rev 128) @@ -708,9 +708,9 @@ sy(i,j,k) = sy(inew,jnew,knew) sz(i,j,k) = sz(inew,jnew,knew) tau(i,j,k) = tau(inew,jnew,knew) - Bvecx(i,j,k)=Bvecx(inew,jnew,knew) - Bvecy(i,j,k)=Bvecy(inew,jnew,knew) - Bvecz(i,j,k)=Bvecz(inew,jnew,knew) + Bconsx(i,j,k)=Bconsx(inew,jnew,knew) + Bconsy(i,j,k)=Bconsy(inew,jnew,knew) + Bconsz(i,j,k)=Bconsz(inew,jnew,knew) if(clean_divergence.ne.0) then psidc(i,j,k)=psidc(inew,jnew,knew) endif @@ -789,9 +789,9 @@ sy(i,j,k) = sy(inew,jnew,k) sz(i,j,k) = sz(inew,jnew,k) tau(i,j,k) = tau(inew,jnew,k) - Bvecx(i,j,k)=Bvecx(inew,jnew,k) - Bvecy(i,j,k)=Bvecy(inew,jnew,k) - Bvecz(i,j,k)=Bvecz(inew,jnew,k) + Bconsx(i,j,k)=Bconsx(inew,jnew,k) + Bconsy(i,j,k)=Bconsy(inew,jnew,k) + Bconsz(i,j,k)=Bconsz(inew,jnew,k) if(clean_divergence.ne.0) then psidc(i,j,k)=psidc(inew,jnew,k) endif @@ -805,9 +805,9 @@ sy(i,j,k) = sy(i,j,kc) sz(i,j,k) = sz(i,j,kc) tau(i,j,k) = tau(i,j,kc) - Bvecx(i,j,k)=Bvecx(i,j,kc) - Bvecy(i,j,k)=Bvecy(i,j,kc) - Bvecz(i,j,k)=Bvecz(i,j,kc) + Bconsx(i,j,k)=Bconsx(i,j,kc) + Bconsy(i,j,k)=Bconsy(i,j,kc) + Bconsz(i,j,k)=Bconsz(i,j,kc) if(clean_divergence.ne.0) then psidc(i,j,k)=psidc(i,j,kc) endif

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 246)
by bcmsma＠astro.rit.edu 28 Apr '11

28 Apr '11

User: bmundim Date: 2011/04/28 04:33 PM Modified: /trunk/src/ GRHydro_CalcUpdate.F90, GRHydro_HLLEM.F90, GRHydro_Reconstruct.F90, GRHydro_SourceM.F90 Log: JOsh's fixes. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_CalcUpdate.F90 Delta lines: +2 -8 =================================================================== --- trunk/src/GRHydro_CalcUpdate.F90 2011-04-28 19:15:41 UTC (rev 245) +++ trunk/src/GRHydro_CalcUpdate.F90 2011-04-28 21:33:20 UTC (rev 246) @@ -59,7 +59,7 @@ do k = GRHydro_stencil + 1, cctk_lsh(3) - GRHydro_stencil do j = GRHydro_stencil + 1, cctk_lsh(2) - GRHydro_stencil do i = GRHydro_stencil + 1, cctk_lsh(1) - GRHydro_stencil - + alp_l = 0.5d0 * (alp(i,j,k) + & alp(i-xoffset,j-yoffset,k-zoffset)) alp_r = 0.5d0 * (alp(i,j,k) + & @@ -103,12 +103,6 @@ divB(i,j,k) = divB(i,j,k) + ( alp_l * Bvec_l - alp_r * Bvec_r ) * idx endif -!!$ if(i.eq.75.and.j.eq.5.and.k.eq.5)write(6,*)'update:', & -!!$ Bconsrhs(i,j,k,1),Bconsxflux(i,j,k), & -!!$ Bconsxflux(i-xoffset,j-yoffset,k-zoffset), & -!!$ Bconsrhs(i,j,k,3),Bconszflux(i,j,k), & -!!$ Bconszflux(i-xoffset,j-yoffset,k-zoffset) - endif if (evolve_tracer .ne. 0) then @@ -270,7 +264,7 @@ divB(i,j,k) = divB(i,j,k) + ( Bvec_l - Bvec_r ) * idx endif endif - + enddo enddo enddo File [modified]: GRHydro_HLLEM.F90 Delta lines: +2 -9 =================================================================== --- trunk/src/GRHydro_HLLEM.F90 2011-04-28 19:15:41 UTC (rev 245) +++ trunk/src/GRHydro_HLLEM.F90 2011-04-28 21:33:20 UTC (rev 246) @@ -111,8 +111,6 @@ cons_m(7) = Bconsyminus(i+xoffset,j+yoffset,k+zoffset) cons_m(8) = Bconszminus(i+xoffset,j+yoffset,k+zoffset) -!!$ if(i.eq.75.and.j.eq.5.and.k.eq.5)write(6,*)'HLLEM0:',cons_m(6),cons_p(6),cons_m(7),cons_p(7),cons_m(8),cons_p(8) - prim_p(1) = rhoplus(i,j,k) prim_p(2) = velxplus(i,j,k) prim_p(3) = velyplus(i,j,k) @@ -353,6 +351,7 @@ fminus(7),fminus(8),fminus(6),& vytm,vztm,vxtm,pressstarm,bylowm,bzlowm,bxlowm,ab0m,wm, & avg_det,avg_alp,avg_beta) + if(clean_divergence.ne.0) then fminus(6)=fminus(6) + avg_alp*sdet*uxyh*psidcm - cons_m(7)*avg_betax fminus(7)=fminus(7) + avg_alp*sdet*uyyh*psidcm - cons_m(7)*avg_betay @@ -386,6 +385,7 @@ fminus(8),fminus(6),fminus(7), & vztm,vxtm,vytm,pressstarm,bzlowm,bxlowm,bylowm,ab0m,wm, & avg_det,avg_alp,avg_beta) + if(clean_divergence.ne.0) then fminus(6)=fminus(6) + avg_alp*sdet*uxzh*psidcm - cons_m(8)*avg_betax fminus(7)=fminus(7) + avg_alp*sdet*uyzh*psidcm - cons_m(8)*avg_betay @@ -414,8 +414,6 @@ charpm = charmax - charmin -!!$ if(i.eq.75.and.j.eq.5.and.k.eq.5)write(6,*)'HLLEM:',fplus(6),fminus(6),fplus(8),fminus(8),charmax,charmin,cons_p(6),cons_m(6) - !!$ Calculate flux by standard formula do m = 1,8 @@ -427,9 +425,6 @@ end do -!!$ if(i.eq.75.and.j.eq.5.and.k.eq.5)write(6,*)'HLLEM2:',fplus(6),fminus(6),fplus(8),fminus(8),charmax,charmin,qdiff(6),qdiff(8),f1(6),f1(8) - - if(clean_divergence.ne.0) then psidcdiff = psidcm - psidcp @@ -456,8 +451,6 @@ Bconsyflux(i, j, k) = f1(7) Bconszflux(i, j, k) = f1(8) -!!$ if(i.eq.75.and.j.eq.5.and.k.eq.5)write(6,*)'HLLEM2:',fplus(6),fminus(6),fplus(8),fminus(8),f1(6),f1(8) - if(clean_divergence.ne.0) then psidcflux(i,j,k) = psidcf endif File [modified]: GRHydro_Reconstruct.F90 Delta lines: +0 -4 =================================================================== --- trunk/src/GRHydro_Reconstruct.F90 2011-04-28 19:15:41 UTC (rev 245) +++ trunk/src/GRHydro_Reconstruct.F90 2011-04-28 21:33:20 UTC (rev 246) @@ -220,8 +220,6 @@ Bvec(:,:,:,3), Bveczplus, Bveczminus, trivial_rp, hydro_excision_mask) endif -!!$ write(6,*)'recon0:',Bvecxplus(75,5,5),Bvecxminus(75,5,5) - else if (CCTK_EQUALS(recon_vars,"conservative")) then call tvdreconstruct(nx, ny, nz, xoffset, yoffset, zoffset, & dens, densplus, densminus, trivial_rp, hydro_excision_mask) @@ -924,8 +922,6 @@ if(evolve_mhd.ne.0) then call primitive2conservativeM(CCTK_PASS_FTOF) -!!$ write(6,*)'recon1:',Bvecxplus(75,5,5),Bvecxminus(75,5,5),Bconsxplus(75,5,5),Bconsxminus(75,5,5) - else call primitive2conservative(CCTK_PASS_FTOF) endif File [modified]: GRHydro_SourceM.F90 Delta lines: +0 -6 =================================================================== --- trunk/src/GRHydro_SourceM.F90 2011-04-28 19:15:41 UTC (rev 245) +++ trunk/src/GRHydro_SourceM.F90 2011-04-28 21:33:20 UTC (rev 246) @@ -443,9 +443,6 @@ srhs(i,j,k,3) = alp(i,j,k)*sqrtdet*sz_source taurhs(i,j,k) = alp(i,j,k)*sqrtdet*tau_source -!!$ if(i.eq.75.and.j.eq.5.and.k.eq.5)write(6,*)'source:',i,j,k, & -!!$ Bconsrhs(i,j,k,1),Bconsrhs(i,j,k,3) - if(clean_divergence.ne.0) then psidcrhs(i,j,k) = -1.d0 * ( ch_dc*ch_dc/cp_dc/cp_dc*alp(i,j,k) + & dx_betax + dy_betay + dz_betaz ) * psidc(i,j,k) + & @@ -504,9 +501,6 @@ psidc(i,j,k)*alp(i,j,k)*sqrtdet*( uxz*gdg_x + uyz*gdg_y + & uzz*gdg_z ) -!!$ if(i.eq.5.and.j.eq.5.and.k.eq.5)write(6,*)'source:',i,j,k, & -!!$ bvcx_source,bvcz_source - Bconsrhsx(i,j,k) = bvcx_source Bconsrhsy(i,j,k) = bvcy_source Bconsrhsz(i,j,k) = bvcz_source

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 245)
by bcmsma＠astro.rit.edu 28 Apr '11

28 Apr '11

User: bmundim Date: 2011/04/28 02:15 PM Modified: /trunk/src/ GRHydro_Prim2ConM.F90 Log: Josh's fix for prim2conM File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Prim2ConM.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_Prim2ConM.F90 2011-04-28 17:55:57 UTC (rev 244) +++ trunk/src/GRHydro_Prim2ConM.F90 2011-04-28 19:15:41 UTC (rev 245) @@ -91,7 +91,7 @@ Bconsxplus(i,j,k),Bconsyplus(i,j,k),Bconszplus(i,j,k), & rhoplus(i,j,k),velxplus(i,j,k),velyplus(i,j,k),& velzplus(i,j,k),epsplus(i,j,k),pressplus(i,j,k),& - Bvecxminus(i,j,k), Bvecyminus(i,j,k), Bveczminus(i,j,k), & + Bvecxplus(i,j,k), Bvecyplus(i,j,k), Bveczplus(i,j,k), & w_lorentzplus(i,j,k)) end do

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 81)
by tanja.bode＠physics.gatech.edu 28 Apr '11

28 Apr '11

User: tbode Date: 2011/04/28 02:13 PM Modified: / ET.tex Log: Analysis section: some rewriting. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +20 -20 =================================================================== --- ET.tex 2011-04-28 18:45:59 UTC (rev 80) +++ ET.tex 2011-04-28 19:13:37 UTC (rev 81) @@ -1233,28 +1233,30 @@ \subsection{Analysis\pages{4 Tanja}} -It is beneficial to evaluate common analysis quantities online -rather than offline, in time-consuming post-processing procedures. -Beyond extracting physics, these quantities are often used as measures -of how accurately the simulation is progressing. Below we detail the -quantities available within Einstein Toolkit modules and the -assumptions and equations used by each. The analysis capabilities of -the Einstein Toolkit broadly fall into three categories: horizons, -masses and momenta, and gravitational waves. Many modules bridge these -these categories and some fall outside them. The latter are described -in the last subsection, including constraint monitoring and tools for -commonly required derived spacetime quantities. The following discussion -is meant as an overview of the most common tools rather than an -exhaustive list of functionality. In most cases, the analysis modules +It is often beneficial and sometimes necessary to evaluate analysis quantities +during the simulation rather than post-processing variable output. Beyond +extracting physics, these quantities are often used as measures of how +accurately the simulation is progressing. In the following, we describe the +common quantities available through Einstein Toolkit modules, and how different +modules approach these quantities with differing assumptions and algorithms. +The most common analysis quantities provided broadly fall into three +categories: horizons, masses and momenta, and gravitational waves. Several +modules bridge these categories and some fall outside them. The latter, are +described in the last subsection, including constraint monitoring and commonly +desired derived spacetime quantities. The following discussion is meant as an +overview of the most common tools rather than an exhaustive list of the +functionality provided by the Einstein Toolkit. In most cases, the analysis modules work on the variables stored in the base modules discussed in Sec.~\ref{sec:base_modules} (\codename{ADMBase}, \codename{TmunuBase}, and \codename{HydroBase}) to create as portable a tool as possible. \subsubsection{Horizons} -For spacetimes which contain a \bh{,} the Einstein Toolkit provides -one module (\codename{EHFinder}) for finding event horizons and -two modules for finding \ahz{s} (\codename{AHFinder} -and \codename{AHFinderDirect}). +When spacetimes contain a \bh{,} localizing a \bh{'s} horizon are necessary +for describing time-dependent quasi-local measures of the \bh{} such as +mass and spin. The Einstein Toolkit provides two modules (\codename{AHFinder} +and \codename{AHFinderDirect}) for locating the \ahz{s}, defined locally on a +hypersurface. The module \codename{EHFinder} is also available to search an +evolved spacetime for the globally defined event horizons. % Event horizon The event horizon module \codename{EHFinder}~\cite{Diener:2003jc} @@ -1530,9 +1532,7 @@ where $S_i=-\frac{1}{\alpha} \left( T_{i0} - \beta^j T_{ij} \right)$. The difference between these modules lies in how they access the stress energy tensor $T_{\mu\nu}$, as the module \codename{ADMConstraints} -uses a \todo{deprecated?} -\todo{ES: yes, deprecated} -functionality which does not require storage +uses a deprecated functionality which does not require storage for $T_{\mu\nu}$. Finally, \codename{ADMAnalysis} calculates a variety of derived spacetime

1 0

[svn:einsteintoolkit] manifest/trunk/ (Rev. 96)
by schnetter＠cct.lsu.edu 28 Apr '11

28 Apr '11

User: eschnett Date: 2011/04/28 01:48 PM Modified: /trunk/ einsteintoolkit.bib Log: Beautify "A" entries File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.bib Delta lines: +66 -54 =================================================================== --- trunk/einsteintoolkit.bib 2011-04-27 22:33:23 UTC (rev 95) +++ trunk/einsteintoolkit.bib 2011-04-28 18:48:33 UTC (rev 96) @@ -3,8 +3,13 @@ not always correct, many have been manually post-processed, e.g. to complete or correct author lists. -We also add a url entry, pointing back to the SPIRES entry. This url +We also add either a DOI entry, if available. + +We also add a URL entry, pointing back to the SPIRES entry. This URl entry uses the DOI (if available), otherwise the arXiv eprint number. +This URL looks like +{http://www.slac.stanford.edu/spires/find/hep/www?rawcmd=find+doi+...} or +{http://www.slac.stanford.edu/spires/find/hep/www?rawcmd=find+eprint+...}. In addition to the standard bibtex keys there are two new ones: requested-for and suggested-for. Both contain space-separated lists of @@ -45,83 +50,90 @@ volume = 77, pages = 104017, year = 2008, - eprint = {arXiv:0710.2335}, - SLACcitation = "%%CITATION = ARXIV:0710.2335;%%" + eprint = {arXiv:0710.2335 [gr-qc]}, + doi = {10.1103/PhysRevD.77.104017}, + url = + {http://www.slac.stanford.edu/spires/find/hep/www?rawcmd=find+doi+10.1103/physrevd.77.104017}, + SLACcitation = "%%CITATION = ARXIV:0710.2335w;%%" } @Article{Alcubierre:1999ab, - author = "Miguel Alcubierre and Steven Brandt and Bernd + author = {Miguel Alcubierre and Steven Brandt and Bernd Br{\"u}gmann and Daniel Holz and Edward Seidel and - Ryoji Takahashi and Jonathan Thornburg", + Ryoji Takahashi and Jonathan Thornburg}, title = "{Symmetry without Symmetry: Numerical Simulation of Axisymmetric Systems using Cartesian Grids}", - journal = "Int. J. Mod. Phys. D", + journal = {Int. J. Mod. Phys. D}, volume = 10, year = 2001, - pages = "273-290", - eprint = "gr-qc/9908012", - archivePrefix ="arXiv", - doi = "10.1142/S0218271801000834", + pages = {273-290}, + eprint = {arXiv:gr-qc/9908012}, + doi = {10.1142/S0218271801000834}, + url = + {http://www.slac.stanford.edu/spires/find/hep/www?rawcmd=find+doi+10.1142/S0218271801000834}, SLACcitation = "%%CITATION = GR-QC/9908012;%%" } @Article{Alcubierre:2004bm, - author = "Miguel Alcubierre and Bernd Br{\"u}gmann and Peter - Diener and Frank Herrmann and Denis Pollney and - Edward Seidel and Ryoji Takahashi", - title = "{Testing excision techniques for dynamical 3D black - hole evolutions}", - year = "2004", - eprint = "gr-qc/0411137", - archivePrefix ="arXiv", - url = {http://www.slac.stanford.edu/spires/find/hep/www?rawcmd=find+eprint+gr-qc/0411137}, - SLACcitation = "%%CITATION = GR-QC/0411137;%%" + author = {Miguel Alcubierre and Bernd Br{\"u}gmann and Peter + Diener and Frank Herrmann and Denis Pollney and + Edward Seidel and Ryoji Takahashi}, + title = {Testing excision techniques for dynamical 3-D black + hole evolutions}, + year = {2004}, + eprint = {arXiv:gr-qc/0411137}, + url = + {http://www.slac.stanford.edu/spires/find/hep/www?rawcmd=find+eprint+gr-qc/0411137}, + SLACcitation = "%%CITATION = GR-QC/0411137;%%" } @Article{Anninos:1996ez, - author = "Peter Anninos and Karen Camarda and Joseph Libson - and Joan Mass{'o} and Edward Seidel and Wai-Mo Suen", - title = "{Finding apparent horizons in dynamic 3-D numerical - space-times}", - journal = "Phys. Rev. D", + author = {Peter Anninos and Karen Camarda and Joseph Libson + and Joan Mass{'o} and Edward Seidel and Wai-Mo Suen}, + title = {Finding apparent horizons in dynamic 3-D numerical + space-times}, + journal = {Phys. Rev. D}, volume = 58, pages = 024003, - doi = "10.1103/PhysRevD.58.024003", year = 1998, - eprint = "gr-qc/9609059", - archivePrefix ="arXiv", - primaryClass = "gr-qc", + eprint = {arXiv:gr-qc/9609059}, + doi = {10.1103/PhysRevD.58.024003}, + url = + {http://www.slac.stanford.edu/spires/find/hep/www?rawcmd=find+doi+10.1103/PhysRevD.58.024003}, + SLACcitation = "%%CITATION = GR-QC/9609059;%%", } @Article{Ansorg:2004ds, - requested-for ="EinsteinInitialData/TwoPunctures", - author = "Ansorg, Marcus and Br{\"u}gmann, Bernd and Tichy, - Wolfgang", - title = "{A single-domain spectral method for black hole - puncture data}", - journal = "Phys. Rev. D", - volume = "70", - year = 2004, - pages = 064011, - eprint = "gr-qc/0404056", - archivePrefix ="arXiv", - url = {http://www.slac.stanford.edu/spires/find/hep/www?rawcmd=find+doi+10.1103/PhysRevD.70.064011}, - doi = "10.1103/PhysRevD.70.064011", - SLACcitation = "%%CITATION = GR-QC/0404056;%%" + requested-for ={EinsteinInitialData/TwoPunctures}, + author = {Ansorg, Marcus and Br{\"u}gmann, Bernd and Tichy, + Wolfgang}, + title = {A single-domain spectral method for black hole + puncture data}, + journal = {Phys. Rev. D}, + volume = 70, + year = 2004, + pages = 064011, + eprint = {arXiv:gr-qc/0404056}, + doi = {10.1103/PhysRevD.70.064011}, + url = + {http://www.slac.stanford.edu/spires/find/hep/www?rawcmd=find+doi+10.1103/PhysRevD.70.064011}, + SLACcitation = "%%CITATION = GR-QC/0404056;%%" } @ARTICLE{Arnowitt:1962hi, - author = "Arnowitt, Richard L. and Deser, Stanley and Misner, Charles W.", - title = "{The dynamics of general relativity}", - journal = {General Relativity and Gravitation}, - year = 2008, - month = sep, - volume = 40, - pages = {1997-2027}, - doi = {10.1007/s10714-008-0661-1}, - eprint = "gr-qc/0405109", - archivePrefix = "arXiv", - SLACcitation = "%%CITATION = GR-QC/0405109;%%" + author = {Arnowitt, Richard L. and Deser, Stanley and Misner, + Charles W.}, + title = {The dynamics of general relativity}, + journal = {General Relativity and Gravitation}, + year = 2008, + month = sep, + volume = 40, + pages = {1997-2027}, + doi = {10.1007/s10714-008-0661-1}, + eprint = {arXiv:gr-qc/0405109}, + url = + {http://www.slac.stanford.edu/spires/find/hep/www?rawcmd=find+doi+10.1007/s10714-008-0661-1}, + SLACcitation = "%%CITATION = GR-QC/0405109;%%" } ##B ############################################################################

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits