Commits February 2011

commits@lists.einsteintoolkit.org

8 participants
34 discussions

[svn:einsteintoolkit] IDAnalyticBH/trunk/src/ (Rev. 178)
by bcmsma＠astro.rit.edu 28 Feb '11

28 Feb '11

User: bmundim Date: 2011/02/28 02:09 PM Modified: /trunk/src/ Kerr.c Log: Segfault with kerr initial data using IDAnalyticBH Yosef's patch: If initial_data=kerr and metric_type=physical, IDAnalyticBH writes to psi[x-z] and psi[x-z][x-z] (it tries to set these GFs to zero) even if no storage is allocated. The problem is with the code: IDAnalyticBH_zero_CCTK_REAL_array(npoints, psix); etc.., which is only executed when the metric is physical. I was able to use the Kerr data with this patch, but just deleting the offending lines should be sufficient. File Changes: Directory: /trunk/src/ ====================== File [modified]: Kerr.c Delta lines: +21 -9 =================================================================== --- trunk/src/Kerr.c 2010-11-21 03:29:37 UTC (rev 177) +++ trunk/src/Kerr.c 2011-02-28 20:09:09 UTC (rev 178) @@ -103,8 +103,13 @@ conformal_storage); /*NOTREACHED*/ } } + else + { + *conformal_state = 0; + } + /* printf("npoints: %i\n",npoints); */ for(i = 0; i < npoints; i++) { @@ -311,15 +316,22 @@ } IDAnalyticBH_zero_CCTK_REAL_array(npoints, gxz); IDAnalyticBH_zero_CCTK_REAL_array(npoints, gyz); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psix); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiy); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiz); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixx); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixy); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixz); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiyy); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiyz); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psizz); + + if (psix && *conformal_state < 2) + { + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psix); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiy); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiz); + } + if (psixx && *conformal_state < 3) + { + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixx); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixy); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixz); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiyy); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiyz); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psizz); + } } return;

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 32)
by knarf＠cct.lsu.edu 28 Feb '11

28 Feb '11

User: knarf Date: 2011/02/28 09:42 AM Modified: / ET.tex Log: some equations File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +12 -4 =================================================================== --- ET.tex 2011-02-27 23:06:15 UTC (rev 31) +++ ET.tex 2011-02-28 15:42:55 UTC (rev 32) @@ -422,7 +422,6 @@ \subsection{Base Modules} \label{sec:base_modules} -\todo{1/3 page Frank} Modular designs have proven to be essential for distributed development of complex software systems, and require the use of well defined interfaces. Low-level interoperability within the Einstein Toolkit is provided by @@ -463,11 +462,20 @@ In particular, it defines the rest-mass density $\varrho$, the gas pressure $p$, internal energy density $\epsilon$, the contravariant gas three velocity defined as -$v^i = \frac{u^i}{\alpha u^0} + \frac{\beta^i}{\alpha}$, the electron +\begin{equation} + v^i = \frac{u^i}{\alpha u^0} + \frac{\beta^i}{\alpha} +\end{equation} +the electron fraction $Y_e$, the temperature $T$, the specific entropy per particle $s$ and the contravariant magnetic field vector $B^i$, defined -as $B^i = \frac{1}{\sqrt{4\pi}} n_{\nu} F^{*\nu i}$, in terms of the dual -$F^{*\mu\nu} = \frac{1}{2}\varepsilon^{\mu\nu\alpha\beta}F_{\alpha\beta}$ +as +\begin{equation} + B^i = \frac{1}{\sqrt{4\pi}} n_{\nu} F^{*\nu i} +\end{equation} +in terms of the dual +\begin{equation} + F^{*\mu\nu} = \frac{1}{2}\varepsilon^{\mu\nu\alpha\beta}F_{\alpha\beta} +\end{equation} to the Faraday tensor and the unit normal of the foliation of spacetime $n^\mu$.

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/cosmology/ (Rev. 31)
by bentivegna＠einsteintoolkit.org 27 Feb '11

27 Feb '11

User: bentivegna Date: 2011/02/27 05:06 PM Modified: /examples/cosmology/ gowdy.par Log: Fixed a few details in gowdy.par. File Changes: Directory: /examples/cosmology/ =============================== File [modified]: gowdy.par Delta lines: +3 -3 =================================================================== --- examples/cosmology/gowdy.par 2011-02-15 14:31:53 UTC (rev 30) +++ examples/cosmology/gowdy.par 2011-02-27 23:06:15 UTC (rev 31) @@ -9,8 +9,8 @@ Cactus::cctk_timer_output = "full" Cactus::terminate = "time" -Cactus::cctk_initial_time = 1 -Cactus::cctk_final_time = 10.0 +Cactus::cctk_initial_time = 1.0 +Cactus::cctk_final_time = 2.0 ############################################################################################# ActiveThorns = "Boundary CartGrid3D CoordBase Slab SymBase Periodic" @@ -237,7 +237,7 @@ Activethorns = "CarpetIOHDF5" ################################# -IOHDF5::out_every = 1 +IOHDF5::out_every = 32 IOHDF5::one_file_per_group = yes IOHDF5::compression_level = 0 IOHDF5::out_vars = "

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 219)
by knarf＠cct.lsu.edu 26 Feb '11

26 Feb '11

User: knarf Date: 2011/02/26 02:51 PM Modified: /trunk/src/ GRHydro_Con2Prim.F90 Log: patch from Roland Haas: in Con2PrimGeneral inside a do File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2011-02-21 17:16:59 UTC (rev 218) +++ trunk/src/GRHydro_Con2Prim.F90 2011-02-26 20:50:58 UTC (rev 219) @@ -2086,7 +2086,7 @@ if (SpaceMask_CheckStateBitsF90(space_mask, i, j, k, type_bits, atmosphere) .or. & GRHydro_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle - f = press_new(i,j,k) - press(i,j,k) + f(i,j,k) = press_new(i,j,k) - press(i,j,k) udens = dens(i,j,k) / sqrt(GRHydro_det(i,j,k)) usx = scon(i,j,k,1) / sqrt(GRHydro_det(i,j,k))

1 0

[svn:einsteintoolkit] Multipole/trunk/src/ (Rev. 75)
by ian.hinder＠aei.mpg.de 23 Feb '11

23 Feb '11

User: hinder Date: 2011/02/23 08:23 AM Modified: /trunk/src/ multipole.cc Log: Correct error check in CCTK_TraverseString Previously the code would abort if there were no variables to decompose File Changes: Directory: /trunk/src/ ====================== File [modified]: multipole.cc Delta lines: +1 -1 =================================================================== --- trunk/src/multipole.cc 2011-01-03 18:50:53 UTC (rev 74) +++ trunk/src/multipole.cc 2011-02-23 14:23:41 UTC (rev 75) @@ -89,7 +89,7 @@ vars.vars = v; int ierr = CCTK_TraverseString(var_string.c_str(), fill_variable, &vars, CCTK_GROUP_OR_VAR); - assert(ierr > 0); + assert(ierr >= 0); *n_variables = vars.n_vars; }

1 0

[svn:einsteintoolkit] HydroBase/trunk/ (Rev. 39)
by bcmsma＠astro.rit.edu 21 Feb '11

21 Feb '11

User: bmundim Date: 2011/02/21 01:13 PM Modified: /trunk/ interface.ccl Log: Lorentz factor tensor type alias fixed. File Changes: Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: +1 -1 =================================================================== --- trunk/interface.ccl 2011-02-21 17:14:39 UTC (rev 38) +++ trunk/interface.ccl 2011-02-21 19:13:43 UTC (rev 39) @@ -35,7 +35,7 @@ CCTK_REAL vel[3] type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="U" interpolator="matter"' "velocity v^i" -CCTK_REAL w_lorentz type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="U" interpolator="matter"' "Lorentz Factor" +CCTK_REAL w_lorentz type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" interpolator="matter"' "Lorentz Factor" CCTK_REAL Y_e type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" interpolator="matter"' "Electron Fraction" CCTK_REAL temperature type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" interpolator="matter"' "Temperature [MeV]"

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 218)
by cott＠tapir.caltech.edu 21 Feb '11

21 Feb '11

User: cott Date: 2011/02/21 11:16 AM Modified: /trunk/ interface.ccl, schedule.ccl /trunk/src/ GRHydro_Boundaries.F90, GRHydro_BoundariesM.F90, GRHydro_Con2Prim.F90, GRHydro_Prim2Con.F90, GRHydro_RegisterGZ.cc, GRHydro_RegisterGZM.cc, GRHydro_RegisterVars.cc, GRHydro_RegisterVarsM.cc Log: * changes to go along with the move of the Lorentz factor variable from GRHydro to HydroBase File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Boundaries.F90 Delta lines: +3 -3 =================================================================== --- trunk/src/GRHydro_Boundaries.F90 2011-02-10 15:23:57 UTC (rev 217) +++ trunk/src/GRHydro_Boundaries.F90 2011-02-21 17:16:59 UTC (rev 218) @@ -64,7 +64,7 @@ call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::press") call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::dens") call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::tau") - call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::w_lorentz") + call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::w_lorentz") call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::eps") call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::GRHydro_C2P_failed") if(evolve_mhd.ne.0.and.clean_divergence.ne.0) then @@ -157,7 +157,7 @@ ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::scon", "Flat") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & - "GRHydro::w_lorentz", "Flat") + "HydroBase::w_lorentz", "Flat") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::rho", "Flat") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & @@ -207,7 +207,7 @@ ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::scon", "None") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & - "GRHydro::w_lorentz", "None") + "HydroBase::w_lorentz", "None") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::rho", "None") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & File [modified]: GRHydro_BoundariesM.F90 Delta lines: +3 -3 =================================================================== --- trunk/src/GRHydro_BoundariesM.F90 2011-02-10 15:23:57 UTC (rev 217) +++ trunk/src/GRHydro_BoundariesM.F90 2011-02-21 17:16:59 UTC (rev 218) @@ -64,7 +64,7 @@ call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::press") call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::dens") call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::tau") - call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::w_lorentz") + call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::w_lorentz") call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::eps") call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::GRHydro_C2P_failed") if(clean_divergence.ne.0) then @@ -152,7 +152,7 @@ ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::scon", "Flat") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & - "GRHydro::w_lorentz", "Flat") + "HydroBase::w_lorentz", "Flat") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::rho", "Flat") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & @@ -192,7 +192,7 @@ ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::scon", "None") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & - "GRHydro::w_lorentz", "None") + "HydroBase::w_lorentz", "None") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::rho", "None") ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & File [modified]: GRHydro_Con2Prim.F90 Delta lines: +27 -0 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2011-02-10 15:23:57 UTC (rev 217) +++ trunk/src/GRHydro_Con2Prim.F90 2011-02-21 17:16:59 UTC (rev 218) @@ -561,6 +561,7 @@ ! error handling if(anyerr.ne.0) then if (GRHydro_reflevel.ge.GRHydro_c2p_warn_from_reflevel) then + !OMP CRITICAL call CCTK_WARN(failinfomode,"EOS error in c2p 0") write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(failinfomode,warnline) @@ -569,6 +570,8 @@ write(warnline,"(A7,i8)") "code: ",keyerr(1) call CCTK_WARN(failinfomode,warnline) write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel + write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel + !OMP END CRITICAL if(keyerr(1).eq.667.and.GRHydro_eos_hot_eps_fix.ne.0) then ! Handling of the case when no new temperature can be ! found for a given epsilon. The amount of times @@ -582,11 +585,14 @@ rho,epsilon,temp,ye,xpress,keyerr,anyerr) nf = nf + 1 enddo + !OMP CRITICAL write(warnline,"(A30,i5)") "Iterations of heat injection: ",nf call CCTK_WARN(failinfomode,warnline) write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(failinfomode,warnline) + !OMP END CRITICAL if(nf.gt.nfudgemax) then + !OMP CRITICAL call CCTK_WARN(failinfomode,"EOS error in c2p 1: injected heat too many times") write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(failinfomode,warnline) @@ -597,10 +603,13 @@ write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel call CCTK_WARN(failinfomode,warnline) call CCTK_WARN(failwarnmode,"Aborting!!!") + !OMP END CRITICAL endif else + !OMP CRITICAL call CCTK_WARN(failinfomode,warnline) call CCTK_WARN(failwarnmode,"Aborting!!!") + !OMP END CRITICAL endif endif endif @@ -636,6 +645,7 @@ ! error handling if(anyerr.ne.0) then if (GRHydro_reflevel.ge.GRHydro_c2p_warn_from_reflevel) then + !OMP CRITICAL call CCTK_WARN(failinfomode,"EOS error in c2p 1") write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(failinfomode,warnline) @@ -644,6 +654,8 @@ write(warnline,"(A7,i8)") "code: ",keyerr(1) call CCTK_WARN(failinfomode,warnline) write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel + call CCTK_WARN(failinfomode,warnline) + !OMP END CRITICAL if(keyerr(1).eq.667.and.GRHydro_eos_hot_eps_fix.ne.0) then ! Handling of the case when no new temperature can be ! found for a given epsilon. The amount of times @@ -657,11 +669,14 @@ rho,epsilon,temp,ye,xpress,keyerr,anyerr) nf = nf + 1 enddo + !OMP CRITICAL write(warnline,"(A30,i5)") "Iterations of heat injection: ",nf call CCTK_WARN(failinfomode,warnline) write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(failinfomode,warnline) + !OMP END CRITICAL if(nf.gt.nfudgemax) then + !OMP CRITICAL call CCTK_WARN(failinfomode,"EOS error in c2p 1: injected heat too many times") write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(failinfomode,warnline) @@ -672,10 +687,13 @@ write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel call CCTK_WARN(failinfomode,warnline) call CCTK_WARN(failwarnmode,"Aborting!!!") + !OMP END CRITICAL endif else + !OMP CRITICAL call CCTK_WARN(failinfomode,warnline) call CCTK_WARN(failwarnmode,"Aborting!!!") + !OMP END CRITICAL endif endif endif @@ -729,6 +747,7 @@ ! error handling if(anyerr.ne.0) then if (GRHydro_reflevel.ge.GRHydro_c2p_warn_from_reflevel) then + !OMP CRITICAL call CCTK_WARN(failinfomode,"EOS error in c2p 2") write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(failinfomode,warnline) @@ -737,6 +756,8 @@ write(warnline,"(A7,i8)") "code: ",keyerr(1) call CCTK_WARN(failinfomode,warnline) write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel + call CCTK_WARN(failinfomode,warnline) + !OMP END CRITICAL if(keyerr(1).eq.667.and.GRHydro_eos_hot_eps_fix.ne.0) then ! Handling of the case when no new temperature can be ! found for a given epsilon. The amount of times @@ -750,11 +771,14 @@ rho,epsilon,temp,ye,xpress,keyerr,anyerr) nf = nf + 1 enddo + !OMP CRITICAL write(warnline,"(A30,i5)") "Iterations of heat injection: ",nf call CCTK_WARN(failinfomode,warnline) write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(failinfomode,warnline) + !OMP END CRITICAL if(nf.gt.nfudgemax) then + !OMP CRITICAL call CCTK_WARN(failinfomode,"EOS error in c2p 1: injected heat too many times") write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(failinfomode,warnline) @@ -765,10 +789,13 @@ write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel call CCTK_WARN(failinfomode,warnline) call CCTK_WARN(failwarnmode,"Aborting!!!") + !OMP END CRITICAL endif else + !OMP CRITICAL call CCTK_WARN(failinfomode,warnline) call CCTK_WARN(failwarnmode,"Aborting!!!") + !OMP END CRITICAL endif endif endif File [modified]: GRHydro_Prim2Con.F90 Delta lines: +4 -2 =================================================================== --- trunk/src/GRHydro_Prim2Con.F90 2011-02-10 15:23:57 UTC (rev 217) +++ trunk/src/GRHydro_Prim2Con.F90 2011-02-21 17:16:59 UTC (rev 218) @@ -269,7 +269,8 @@ ! This is a way of recovering even on finer refinement levels: ! Use the average temperature at the interface instead of the ! reconstructed specific internal energy. - !OMP CRITICAL +!! !OMP CRITICAL +#if 0 call CCTK_WARN(1,"EOS error in prim2con_hot: NOW using averaged temp!") write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(1,warnline) @@ -279,7 +280,8 @@ call CCTK_WARN(1,warnline) write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel call CCTK_WARN(1,warnline) - !OMP END CRITICAL +#endif +!!! !OMP END CRITICAL keytemp=1 temp = temp0 call EOS_Omni_press(handle,keytemp,GRHydro_eos_rf_prec,n,& File [modified]: GRHydro_RegisterGZ.cc Delta lines: +6 -6 =================================================================== --- trunk/src/GRHydro_RegisterGZ.cc 2011-02-10 15:23:57 UTC (rev 217) +++ trunk/src/GRHydro_RegisterGZ.cc 2011-02-21 17:16:59 UTC (rev 218) @@ -39,7 +39,7 @@ if(clean_divergence) { string var[10] = {"HydroBase::rho", "HydroBase::press", "HydroBase::eps", "HydroBase::vel", - "GRHydro::dens", "GRHydro::tau", "GRHydro::w_lorentz", + "GRHydro::dens", "GRHydro::tau", "HydroBase::w_lorentz", "GRHydro::scon", "HydroBase::Bvec", "GRHydro::psidc"}; for (int i = 0; i < 10; i++) { @@ -58,7 +58,7 @@ } else { string var[9] = {"HydroBase::rho", "HydroBase::press", "HydroBase::eps", "HydroBase::vel", - "GRHydro::dens", "GRHydro::tau", "GRHydro::w_lorentz", + "GRHydro::dens", "GRHydro::tau", "HydroBase::w_lorentz", "GRHydro::scon", "HydroBase::Bvec"}; for (int i = 0; i < 9; i++) { @@ -78,7 +78,7 @@ } else { string var[8] = {"HydroBase::rho", "HydroBase::press", "HydroBase::eps", "HydroBase::vel", - "GRHydro::dens", "GRHydro::tau", "GRHydro::w_lorentz", + "GRHydro::dens", "GRHydro::tau", "HydroBase::w_lorentz", "GRHydro::scon"}; for (int i = 0; i < 8; i++) { @@ -110,7 +110,7 @@ if(clean_divergence) { string var[13] = {"HydroBase::rho", "HydroBase::press", "HydroBase::eps", "HydroBase::vel", - "GRHydro::dens", "GRHydro::tau", "GRHydro::w_lorentz", + "GRHydro::dens", "GRHydro::tau", "HydroBase::w_lorentz", "GRHydro::scon", "HydroBase::Bvec", "GRHydro::psidc", "ADMBase::metric", "ADMBase::curv", "ADMBase::shift"}; for (int i = 0; i < 13; i++) @@ -131,7 +131,7 @@ } else { string var[12] = {"HydroBase::rho", "HydroBase::press", "HydroBase::eps", "HydroBase::vel", - "GRHydro::dens", "GRHydro::tau", "GRHydro::w_lorentz", + "GRHydro::dens", "GRHydro::tau", "HydroBase::w_lorentz", "GRHydro::scon", "HydroBase::Bvec", "ADMBase::metric", "ADMBase::curv", "ADMBase::shift"}; for (int i = 0; i < 12; i++) @@ -154,7 +154,7 @@ } else { string var[11] = {"HydroBase::rho", "HydroBase::press", "HydroBase::eps", "HydroBase::vel", - "GRHydro::dens", "GRHydro::tau", "GRHydro::w_lorentz", + "GRHydro::dens", "GRHydro::tau", "HydroBase::w_lorentz", "GRHydro::scon", "ADMBase::metric", "ADMBase::curv", "ADMBase::shift"}; for (int i = 0; i < 11; i++) File [modified]: GRHydro_RegisterGZM.cc Delta lines: +4 -4 =================================================================== --- trunk/src/GRHydro_RegisterGZM.cc 2011-02-10 15:23:57 UTC (rev 217) +++ trunk/src/GRHydro_RegisterGZM.cc 2011-02-21 17:16:59 UTC (rev 218) @@ -38,7 +38,7 @@ if(clean_divergence) { string var[10] = {"HydroBase::rho", "HydroBase::press", "HydroBase::eps", "HydroBase::vel", - "GRHydro::dens", "GRHydro::tau", "GRHydro::w_lorentz", + "GRHydro::dens", "GRHydro::tau", "HydroBase::w_lorentz", "GRHydro::scon", "HydroBase::Bvec", "GRHydro::psidc"}; for (int i = 0; i < 10; i++) { @@ -57,7 +57,7 @@ } else { string var[9] = {"HydroBase::rho", "HydroBase::press", "HydroBase::eps", "HydroBase::vel", - "GRHydro::dens", "GRHydro::tau", "GRHydro::w_lorentz", + "GRHydro::dens", "GRHydro::tau", "HydroBase::w_lorentz", "GRHydro::scon", "HydroBase::Bvec"}; for (int i = 0; i < 9; i++) { @@ -89,7 +89,7 @@ if(clean_divergence) { string var[13] = {"HydroBase::rho", "HydroBase::press", "HydroBase::eps", "HydroBase::vel", - "GRHydro::dens", "GRHydro::tau", "GRHydro::w_lorentz", + "GRHydro::dens", "GRHydro::tau", "HydroBase::w_lorentz", "GRHydro::scon", "HydroBase::Bvec", "GRHydro::psidc", "ADMBase::metric", "ADMBase::curv", "ADMBase::shift"}; for (int i = 0; i < 13; i++) @@ -110,7 +110,7 @@ } else { string var[12] = {"HydroBase::rho", "HydroBase::press", "HydroBase::eps", "HydroBase::vel", - "GRHydro::dens", "GRHydro::tau", "GRHydro::w_lorentz", + "GRHydro::dens", "GRHydro::tau", "HydroBase::w_lorentz", "GRHydro::scon", "HydroBase::Bvec", "ADMBase::metric", "ADMBase::curv", "ADMBase::shift"}; for (int i = 0; i < 12; i++) File [modified]: GRHydro_RegisterVars.cc Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_RegisterVars.cc 2011-02-10 15:23:57 UTC (rev 217) +++ trunk/src/GRHydro_RegisterVars.cc 2011-02-21 17:16:59 UTC (rev 218) @@ -54,7 +54,7 @@ register_constrained("HydroBase::press"); register_constrained("HydroBase::eps"); register_constrained("HydroBase::vel"); - register_constrained("GRHydro::w_lorentz"); + register_constrained("HydroBase::w_lorentz"); if (CCTK_EQUALS(evolution_method, "GRHydro")) { File [modified]: GRHydro_RegisterVarsM.cc Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_RegisterVarsM.cc 2011-02-10 15:23:57 UTC (rev 217) +++ trunk/src/GRHydro_RegisterVarsM.cc 2011-02-21 17:16:59 UTC (rev 218) @@ -54,7 +54,7 @@ register_constrained("HydroBase::press"); register_constrained("HydroBase::eps"); register_constrained("HydroBase::vel"); - register_constrained("GRHydro::w_lorentz"); + register_constrained("HydroBase::w_lorentz"); if (CCTK_EQUALS(evolution_method, "GRHydro")) { Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: +1 -1 =================================================================== --- trunk/interface.ccl 2011-02-10 15:23:57 UTC (rev 217) +++ trunk/interface.ccl 2011-02-21 17:16:59 UTC (rev 218) @@ -345,7 +345,7 @@ tracer } "Tracers" -real w_lorentz type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" tensorweight=+1.0 interpolator="matter"' "Lorentz factor" +#real w_lorentz type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" tensorweight=+1.0 interpolator="matter"' "Lorentz factor" real psidc type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" tensorweight=+1.0 interpolator="matter"' "Psi parameter for divergence cleaning" File [modified]: schedule.ccl Delta lines: +2 -4 =================================================================== --- trunk/schedule.ccl 2011-02-10 15:23:57 UTC (rev 217) +++ trunk/schedule.ccl 2011-02-21 17:16:59 UTC (rev 218) @@ -16,7 +16,6 @@ STORAGE:dens[3] STORAGE:tau[3] STORAGE:scon[3] - STORAGE:w_lorentz[3] if (number_of_particles) { STORAGE:particles[3] @@ -42,7 +41,6 @@ STORAGE:dens[2] STORAGE:tau[2] STORAGE:scon[2] - STORAGE:w_lorentz[2] if (number_of_particles) { STORAGE:particles[2] @@ -375,7 +373,7 @@ SYNC: dens SYNC: tau SYNC: scon - SYNC: w_lorentz + SYNC: hydrobase::w_lorentz SYNC: hydrobase::rho SYNC: hydrobase::press SYNC: hydrobase::eps @@ -1024,7 +1022,7 @@ SYNC: dens SYNC: tau SYNC: scon - SYNC: w_lorentz + SYNC: HydroBase::w_lorentz SYNC: HydroBase::rho SYNC: HydroBase::press SYNC: HydroBase::eps

1 0

[svn:einsteintoolkit] HydroBase/trunk/ (Rev. 38)
by cott＠tapir.caltech.edu 21 Feb '11

21 Feb '11

User: cott Date: 2011/02/21 11:14 AM Modified: /trunk/ interface.ccl, schedule.ccl Log: * add variable w_lorentz for the Lorentz factor File Changes: Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: +2 -0 =================================================================== --- trunk/interface.ccl 2011-02-21 15:49:06 UTC (rev 37) +++ trunk/interface.ccl 2011-02-21 17:14:39 UTC (rev 38) @@ -35,6 +35,8 @@ CCTK_REAL vel[3] type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="U" interpolator="matter"' "velocity v^i" +CCTK_REAL w_lorentz type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="U" interpolator="matter"' "Lorentz Factor" + CCTK_REAL Y_e type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" interpolator="matter"' "Electron Fraction" CCTK_REAL temperature type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" interpolator="matter"' "Temperature [MeV]" CCTK_REAL entropy type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" interpolator="matter"' "Specific Entropy [k_b/baryon]" File [modified]: schedule.ccl Delta lines: +3 -0 =================================================================== --- trunk/schedule.ccl 2011-02-21 15:49:06 UTC (rev 37) +++ trunk/schedule.ccl 2011-02-21 17:14:39 UTC (rev 38) @@ -6,6 +6,7 @@ STORAGE: press[1] STORAGE: eps[1] STORAGE: vel[1] + STORAGE: w_lorentz[1] if (!CCTK_EQUALS(initial_Y_e, "none")) { STORAGE: Y_e[1] @@ -29,6 +30,7 @@ STORAGE: press[2] STORAGE: eps[2] STORAGE: vel[2] + STORAGE: w_lorentz[2] if (!CCTK_EQUALS(initial_Y_e, "none")) { STORAGE: Y_e[2] @@ -52,6 +54,7 @@ STORAGE: press[3] STORAGE: eps[3] STORAGE: vel[3] + STORAGE: w_lorentz[3] if (!CCTK_EQUALS(initial_Y_e, "none")) { STORAGE: Y_e[3]

1 0

[svn:einsteintoolkit] HydroBase/trunk/src/ (Rev. 37)
by roland.haas＠physics.gatech.edu 21 Feb '11

21 Feb '11

User: rhaas Date: 2011/02/21 09:49 AM Modified: /trunk/src/ StartUp.c Log: declare argument list of HydroBase_StartUp as (void) otherwise C (both C89 and C99) assume that it is as being and old K&R style function definition (C differs from C++ here) File Changes: Directory: /trunk/src/ ====================== File [modified]: StartUp.c Delta lines: +1 -1 =================================================================== --- trunk/src/StartUp.c 2011-02-04 14:23:52 UTC (rev 36) +++ trunk/src/StartUp.c 2011-02-21 15:49:06 UTC (rev 37) @@ -1,6 +1,6 @@ #include <cctk.h> -int HydroBase_StartUp() +int HydroBase_StartUp(void) { CCTK_RegisterBanner("HydroBase: Let it flow.");

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/cosmology/ (Rev. 30)
by bentivegna＠einsteintoolkit.org 15 Feb '11

15 Feb '11

User: bentivegna Date: 2011/02/15 08:31 AM Added: /examples/cosmology/ gowdy.par, kasner.par Log: Added two cosmological example solutions. File Changes: Directory: /examples/cosmology/ =============================== File [added]: gowdy.par Delta lines: +257 -0 =================================================================== --- examples/cosmology/gowdy.par (rev 0) +++ examples/cosmology/gowdy.par 2011-02-15 14:31:53 UTC (rev 30) @@ -0,0 +1,257 @@ +###################################### +# Flesh parameters +###################################### + +Cactus::cctk_run_title = "Gowdy spacetime" + +Cactus::cctk_full_warnings = yes +Cactus::highlight_warning_messages = no +Cactus::cctk_timer_output = "full" + +Cactus::terminate = "time" +Cactus::cctk_initial_time = 1 +Cactus::cctk_final_time = 10.0 + +############################################################################################# +ActiveThorns = "Boundary CartGrid3D CoordBase Slab SymBase Periodic" +############################################################################################# + +CoordBase::domainsize = "minmax" + +CoordBase::xmin = -0.51 +CoordBase::ymin = -0.51 +CoordBase::zmin = -0.51 +CoordBase::xmax = +0.51 +CoordBase::ymax = +0.51 +CoordBase::zmax = +0.51 +CoordBase::dx = 0.02 +CoordBase::dy = 0.02 +CoordBase::dz = 0.02 + +CoordBase::boundary_size_x_lower = 3 +CoordBase::boundary_size_y_lower = 3 +CoordBase::boundary_size_z_lower = 3 +CoordBase::boundary_size_x_upper = 3 +CoordBase::boundary_size_y_upper = 3 +CoordBase::boundary_size_z_upper = 3 + +CoordBase::boundary_shiftout_x_lower = 0 +CoordBase::boundary_shiftout_y_lower = 0 +CoordBase::boundary_shiftout_z_lower = 0 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 0 +CoordBase::boundary_shiftout_z_upper = 0 + +CartGrid3D::type = "coordbase" + +Periodic::periodic = "yes" +Periodic::periodic_x = "yes" +Periodic::periodic_y = "yes" +Periodic::periodic_z = "yes" + +################################################################################## +ActiveThorns = "IOUtil Carpet CarpetLib CarpetInterp LoopControl CarpetReduce CarpetSlab" +################################################################################## + +Carpet::schedule_barriers = no +CarpetLib::output_bboxes = no +Carpet::init_fill_timelevels = "yes" + +Carpet::domain_from_coordbase = yes +Carpet::max_refinement_levels = 1 + +driver::ghost_size = 3 +Carpet::use_buffer_zones = yes + +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 2 + +Carpet::convergence_level = 0 + +Carpet::poison_new_timelevels = yes +CarpetLib::poison_new_memory = yes + +Carpet::output_timers_every = 512 +CarpetLib::print_timestats_every = 512 +CarpetLib::print_memstats_every = 512 + +CarpetLib::max_allowed_memory_MB = 2048 + +################################################################## +ActiveThorns = "SphericalSurface ADMBase" +################################################################## + +Carpet::grid_structure_filename = "grid_structure.asc" +Carpet::grid_coordinates_filename = "grid_coordinates.asc" + + +################################# +ActiveThorns = "AEILocalInterp" +ActiveThorns = "BLAS LAPACK" +ActiveThorns = "Fortran" +ActiveThorns = "GenericFD" +ActiveThorns = "HDF5" +ActiveThorns = "InitBase" +ActiveThorns = "LocalInterp" +ActiveThorns = "TGRtensor" +################################# + +InitBase::initial_data_setup_method = "init_some_levels" + +################################# +ActiveThorns = "MoL Time" +################################# + +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +Time::dtfac = 0.05 + +################################################################################################### +ActiveThorns = "ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase" +################################################################################################### + +ADMMacros::spatial_order = 4 + +################################# +ActiveThorns = "Exact" +################################# + +ADMBase::metric_type = "physical" + +Exact::exact_model = "Gowdy-wave" +Exact::Gowdy_wave__amplitude = 1.0 + +################################# +ActiveThorns = "ML_BSSN ML_BSSN_Helper ML_ADMConstraints" +################################# + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" +ADMBase::lapse_timelevels = 1 +ADMBase::shift_timelevels = 1 + +ADMBase::evolution_method = "ML_BSSN" +ADMBase::lapse_evolution_method = "exact" +ADMBase::shift_evolution_method = "exact" +ADMBase::dtlapse_evolution_method = "static" +ADMBase::dtshift_evolution_method = "static" + +ML_BSSN::ML_log_confac_bound = "none" +ML_BSSN::ML_metric_bound = "none" +ML_BSSN::ML_Gamma_bound = "none" +ML_BSSN::ML_trace_curv_bound = "none" +ML_BSSN::ML_curv_bound = "none" +ML_BSSN::ML_lapse_bound = "none" +ML_BSSN::ML_dtlapse_bound = "none" +ML_BSSN::ML_shift_bound = "none" +ML_BSSN::ML_dtshift_bound = "none" + + +################################# +ActiveThorns = "NaNChecker" +################################# + +NaNChecker::check_every = 1 +NaNChecker::report_max = 1 + +NaNChecker::action_if_found = "terminate" + +NaNChecker::verbose = "all" +NaNChecker::check_vars = " + ML_BSSN::ML_curv + ML_BSSN::ML_Gamma + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_trace_curv +" + +################################# +ActiveThorns = "CarpetIOBasic" +################################# + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_reductions = "norm2" +IOBasic::outInfo_vars = " + ADMBase::lapse + ADMBase::gxx + ADMBase::Kxx + ML_ADMConstraints::ML_ham + Carpet::physical_time_per_hour +" + + +################################# +ActiveThorns = "CarpetIOScalar" +################################# + +IOScalar::one_file_per_group = yes + +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ML_ADMConstraints::ML_ham + ML_ADMConstraints::ML_mom + ADMBase::lapse + ADMBase::shift + ADMBase::metric + ADMBase::curv +" + +################################# +ActiveThorns = "CarpetIOASCII" +################################# + +IOASCII::one_file_per_group = yes + +IOASCII::output_symmetry_points = no +IOASCII::out3D_ghosts = no + +IOASCII::out0D_every = 0 +IOASCII::out0D_vars = " + Carpet::timing +" + +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ML_ADMConstraints::ML_ham + ML_ADMConstraints::ML_mom + ADMBase::lapse + ADMBase::shift + ADMBase::metric + ADMBase::curv +" + +IOASCII::out2D_every = 0 +IOASCII::out2D_vars = " +" + +################################# +Activethorns = "CarpetIOHDF5" +################################# + +IOHDF5::out_every = 1 +IOHDF5::one_file_per_group = yes +IOHDF5::compression_level = 0 +IOHDF5::out_vars = " + ML_ADMConstraints::ML_ham + ML_ADMConstraints::ML_mom + ADMBase::lapse + ADMBase::shift + ADMBase::metric + ADMBase::curv +" + +################################# +ActiveThorns = "TimerReport" +################################# + +TimerReport::out_every = 512 +TimerReport::out_filename = "TimerReport" File [added]: kasner.par Delta lines: +259 -0 =================================================================== --- examples/cosmology/kasner.par (rev 0) +++ examples/cosmology/kasner.par 2011-02-15 14:31:53 UTC (rev 30) @@ -0,0 +1,259 @@ +###################################### +# Flesh parameters +###################################### + +Cactus::cctk_run_title = "Kasner spacetime" + +Cactus::cctk_full_warnings = yes +Cactus::highlight_warning_messages = no +Cactus::cctk_timer_output = "full" + +Cactus::terminate = "time" +Cactus::cctk_initial_time = 1 +Cactus::cctk_final_time = 10.0 + + +############################################################################################# +ActiveThorns = "Boundary CartGrid3D CoordBase Slab SymBase Periodic" +############################################################################################# + +CoordBase::domainsize = "minmax" + +CoordBase::xmin = -5.00 +CoordBase::ymin = -5.00 +CoordBase::zmin = -5.00 +CoordBase::xmax = +5.00 +CoordBase::ymax = +5.00 +CoordBase::zmax = +5.00 +CoordBase::dx = 0.20 +CoordBase::dy = 0.20 +CoordBase::dz = 0.20 + +CoordBase::boundary_size_x_lower = 3 +CoordBase::boundary_size_y_lower = 3 +CoordBase::boundary_size_z_lower = 3 +CoordBase::boundary_size_x_upper = 3 +CoordBase::boundary_size_y_upper = 3 +CoordBase::boundary_size_z_upper = 3 + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 +CoordBase::boundary_shiftout_x_upper = 1 +CoordBase::boundary_shiftout_y_upper = 1 +CoordBase::boundary_shiftout_z_upper = 1 + +CartGrid3D::type = "coordbase" + +Periodic::periodic = "yes" +Periodic::periodic_x = "yes" +Periodic::periodic_y = "yes" +Periodic::periodic_z = "yes" + +################################################################################## +ActiveThorns = "IOUtil Carpet CarpetLib CarpetInterp LoopControl CarpetReduce CarpetSlab" +################################################################################## + +Carpet::schedule_barriers = no +CarpetLib::output_bboxes = no +Carpet::init_fill_timelevels = "yes" + +Carpet::domain_from_coordbase = yes +Carpet::max_refinement_levels = 1 + +driver::ghost_size = 3 +Carpet::use_buffer_zones = yes + +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 2 + +Carpet::convergence_level = 0 + +Carpet::poison_new_timelevels = yes +CarpetLib::poison_new_memory = yes + +Carpet::output_timers_every = 512 +CarpetLib::print_timestats_every = 512 +CarpetLib::print_memstats_every = 512 + +CarpetLib::max_allowed_memory_MB = 2048 + +################################################################## +ActiveThorns = "SphericalSurface ADMBase" +################################################################## + +Carpet::grid_structure_filename = "grid_structure.asc" +Carpet::grid_coordinates_filename = "grid_coordinates.asc" + + +################################# +ActiveThorns = "AEILocalInterp" +ActiveThorns = "BLAS LAPACK" +ActiveThorns = "Fortran" +ActiveThorns = "GenericFD" +ActiveThorns = "HDF5" +ActiveThorns = "InitBase" +ActiveThorns = "LocalInterp" +ActiveThorns = "TGRtensor" +################################# + +InitBase::initial_data_setup_method = "init_some_levels" + +################################# +ActiveThorns = "MoL Time" +################################# + +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +Time::dtfac = 0.05 + +################################################################################################### +ActiveThorns = "ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase" +################################################################################################### + +ADMMacros::spatial_order = 4 + +################################# +ActiveThorns = "Exact" +################################# + +ADMBase::metric_type = "physical" + +Exact::exact_model = "Kasner-like" +Exact::Kasner_generalized__p1 = 0.666666666666666666666666666666666666666666666666666666666666 +Exact::Kasner_generalized__p2 = 0.666666666666666666666666666666666666666666666666666666666666 + +################################# +ActiveThorns = "ML_BSSN ML_BSSN_Helper ML_ADMConstraints" +################################# + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "one" +ADMBase::initial_shift = "zero" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" +ADMBase::lapse_timelevels = 1 +ADMBase::shift_timelevels = 1 + +ADMBase::evolution_method = "ML_BSSN" +ADMBase::lapse_evolution_method = "static" +ADMBase::shift_evolution_method = "static" +ADMBase::dtlapse_evolution_method = "static" +ADMBase::dtshift_evolution_method = "static" + +ML_BSSN::ML_log_confac_bound = "none" +ML_BSSN::ML_metric_bound = "none" +ML_BSSN::ML_Gamma_bound = "none" +ML_BSSN::ML_trace_curv_bound = "none" +ML_BSSN::ML_curv_bound = "none" +ML_BSSN::ML_lapse_bound = "none" +ML_BSSN::ML_dtlapse_bound = "none" +ML_BSSN::ML_shift_bound = "none" +ML_BSSN::ML_dtshift_bound = "none" + + +################################# +ActiveThorns = "NaNChecker" +################################# + +NaNChecker::check_every = 1 +NaNChecker::report_max = 1 + +NaNChecker::action_if_found = "terminate" + +NaNChecker::verbose = "all" +NaNChecker::check_vars = " + ML_BSSN::ML_curv + ML_BSSN::ML_Gamma + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_trace_curv +" + +################################# +ActiveThorns = "CarpetIOBasic" +################################# + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_reductions = "norm2" +IOBasic::outInfo_vars = " + ADMBase::lapse + ADMBase::gxx + ADMBase::Kxx + ML_ADMConstraints::ML_ham + Carpet::physical_time_per_hour +" + + +################################# +ActiveThorns = "CarpetIOScalar" +################################# + +IOScalar::one_file_per_group = yes + +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ML_ADMConstraints::ML_ham + ML_ADMConstraints::ML_mom + ADMBase::lapse + ADMBase::shift + ADMBase::metric + ADMBase::curv +" + +################################# +ActiveThorns = "CarpetIOASCII" +################################# + +IOASCII::one_file_per_group = yes + +IOASCII::output_symmetry_points = no +IOASCII::out3D_ghosts = no + +IOASCII::out0D_every = 0 +IOASCII::out0D_vars = " + Carpet::timing +" + +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ML_ADMConstraints::ML_ham + ML_ADMConstraints::ML_mom + ADMBase::lapse + ADMBase::shift + ADMBase::metric + ADMBase::curv +" + +IOASCII::out2D_every = 0 +IOASCII::out2D_vars = " +" + +################################# +Activethorns = "CarpetIOHDF5" +################################# + +IOHDF5::out_every = 1 +IOHDF5::one_file_per_group = yes +IOHDF5::compression_level = 0 +IOHDF5::out_vars = " + ML_ADMConstraints::ML_ham + ML_ADMConstraints::ML_mom + ADMBase::lapse + ADMBase::shift + ADMBase::metric + ADMBase::curv +" + +################################# +ActiveThorns = "TimerReport" +################################# + +TimerReport::out_every = 512 +TimerReport::out_filename = "TimerReport"

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits February 2011