Commits June 2012

commits@lists.einsteintoolkit.org

4 participants
73 discussions

New changesets in Carpet Mercurial
by Carpet mercurial repositories 28 Jun '12

28 Jun '12

changeset: 3593:ca7508a01347 tag: tip user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Thu Jun 28 14:31:29 2012 -0700 summary: Carpet: call POSTREGRID on all finer levels of a regridded level

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 374)
by rhaas＠tapir.caltech.edu 26 Jun '12

26 Jun '12

User: rhaas Date: 2012/06/26 06:48 PM Modified: /trunk/src/ GRHydro_Con2Prim.F90 Log: GRHydro: remove "use" for con2prim interfaces from pointwise routines since the Fortran standard does not allow this and gfortran 4.7 actually chokes on it. Pointed out by Erik Schnetter. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: +0 -10 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2012-06-26 15:29:13 UTC (rev 373) +++ trunk/src/GRHydro_Con2Prim.F90 2012-06-26 23:48:01 UTC (rev 374) @@ -349,8 +349,6 @@ uyz, uzz, det, x, y, z, r, epsnegative, GRHydro_rho_min, pmin, epsmin, & GRHydro_reflevel, GRHydro_C2P_failed) - use Con2Prim_fortran_interfaces - implicit none DECLARE_CCTK_PARAMETERS @@ -614,8 +612,6 @@ uyz, uzz, det, x, y, z, r, epsnegative, GRHydro_rho_min, pmin, epsmin, & GRHydro_reflevel, GRHydro_C2P_failed, local_perc_ptol) - use Con2Prim_fortran_interfaces - implicit none DECLARE_CCTK_PARAMETERS @@ -1476,8 +1472,6 @@ velx, vely, velz, epsilon, press, w_lorentz, uxx, uxy, uxz, uyy, & uyz, uzz, det, x, y, z, r, GRHydro_rho_min, GRHydro_reflevel, GRHydro_C2P_failed) - use Con2Prim_fortran_interfaces - implicit none DECLARE_CCTK_PARAMETERS @@ -1980,8 +1974,6 @@ subroutine Con2Prim_ptTracer(cons_tracer, tracer, dens) - use Con2Prim_fortran_interfaces - implicit none DECLARE_CCTK_PARAMETERS @@ -2010,8 +2002,6 @@ subroutine Con2Prim_ptBoundsTracer(cons_tracer, tracer, rho, one_over_w_lorentz, det) - use Con2Prim_fortran_interfaces - implicit none DECLARE_CCTK_PARAMETERS

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 26 Jun '12

26 Jun '12

changeset: 3590:f18d797901f1 user: Roland Haas <rhaas(a)caltech.edu> date: Thu Jun 21 09:53:24 2012 -0700 summary: Carpet: back out of 18d226c24a35 (again) changeset: 3591:a3ad683a0607 user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Mon Jun 18 13:31:11 2012 -0700 summary: CarpetLib: remove comment about not supporting overlaps in vertex changeset: 3592:55759e108807 tag: tip user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Mon Jun 25 21:03:35 2012 -0700 summary: Carpet: ensure scheduling of global-early routines in POSTREGRID

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 373)
by rhaas＠tapir.caltech.edu 26 Jun '12

26 Jun '12

User: rhaas Date: 2012/06/26 10:29 AM Modified: /trunk/ param.ccl Log: GRHydro: increase maximum and default number of constraint variables to include Tmunu in the count. Extra constraints are cheap since MoL does not allocate any scratch or save storage for them. Bruno Mundim pointed this one out and provided the accurate count. File Changes: Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +15 -2 =================================================================== --- trunk/param.ccl 2012-06-22 20:37:46 UTC (rev 372) +++ trunk/param.ccl 2012-06-26 15:29:13 UTC (rev 373) @@ -79,10 +79,23 @@ 5:10 :: "dens scon[3] tau Bvec[3] psidc ye" } 5 +# 7 primitives (rho,press,eps,w_lorentz,vel) +# 10 Tmunu +# 3 Bvec +# 1 tau (polytype EOS) +# 1 Y_e +# 1 temperature +# 1 entropy +# 3 comoving shift +# 3 lvel +# 3 local_shift +# 6 local_metric +# 6 local_extrinsic_curvature +# 3 lBvec CCTK_INT GRHydro_MaxNumConstrainedVars "The maximum number of constrained variables used by GRHydro" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Constrained_Vars { - 7:33 :: "A small range, depending on testing or not" -} 27 + 7:48 :: "A small range, depending on testing or not" +} 37 CCTK_INT GRHydro_MaxNumSandRVars "The maximum number of save and restore variables used by GRHydro" ACCUMULATOR-BASE=MethodofLines::MoL_Num_SaveAndRestore_Vars {

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 25 Jun '12

25 Jun '12

changeset: 3589:435d3519cb69 tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Mon Jun 25 15:08:09 2012 -0400 summary: Carpet: Output informative message if HAVE_SCHED_GETAFFINITY is missing

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 372)
by rhaas＠tapir.caltech.edu 22 Jun '12

22 Jun '12

User: rhaas Date: 2012/06/22 03:37 PM Modified: /trunk/ interface.ccl, schedule.ccl /trunk/src/ GRHydro_Boundaries.F90, GRHydro_UpdateMask.F90 Log: Revert "* improve handling of cs2 < 0 in the "standard" (non-MHD, non-hot EOS) eigenvalue routine." actually this introduce the real valued atmosphere mask. This reverts commit deccf528b463820b615dfa73e46e137cd59c65be. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Boundaries.F90 Delta lines: +1 -6 =================================================================== --- trunk/src/GRHydro_Boundaries.F90 2012-06-22 20:37:44 UTC (rev 371) +++ trunk/src/GRHydro_Boundaries.F90 2012-06-22 20:37:46 UTC (rev 372) @@ -99,8 +99,6 @@ !endif endif - call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::atmosphere_mask_real") - sym(1) = -1 sym(2) = 1 sym(3) = 1 @@ -212,8 +210,6 @@ "GRHydro::tau", "Flat") ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::scon", "Flat") - ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & - "GRHydro::GRHydro_atmosphere_mask_real", "Flat") if (sync_conserved_only .eq. 0) then ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::w_lorentz", "Flat") @@ -286,8 +282,7 @@ "GRHydro::tau", "None") ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::scon", "None") - ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & - "GRHydro::GRHydro_atmosphere_mask_real", "None") + if (sync_conserved_only .eq. 0) then ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::w_lorentz", "None") File [modified]: GRHydro_UpdateMask.F90 Delta lines: +2 -35 =================================================================== --- trunk/src/GRHydro_UpdateMask.F90 2012-06-22 20:37:44 UTC (rev 371) +++ trunk/src/GRHydro_UpdateMask.F90 2012-06-22 20:37:46 UTC (rev 372) @@ -48,9 +48,7 @@ densrhs(i,j,k) = 0.0d0 srhs(i,j,k,:) = 0.0d0 taurhs(i,j,k) = 0.0d0 - ! Set real-valued mask! This will be sync'ed and right after syncing translated to - ! our standard integer based mask (so that atmosphere_mask is still valid!). - atmosphere_mask_real(i,j,k) = 1 + atmosphere_mask(i,j,k) = 1 end if end do end do @@ -69,9 +67,7 @@ densrhs(i,j,k) = 0.0d0 srhs(i,j,k,:) = 0.0d0 taurhs(i,j,k) = 0.0d0 - ! Set real-valued mask! This will be sync'ed and right after syncing translated to - ! our standard integer based mask (so that atmosphere_mask is still valid!). - atmosphere_mask_real(i,j,k) = 1 + atmosphere_mask(i,j,k) = 1 end if end do end do @@ -81,33 +77,6 @@ end subroutine GRHydroUpdateAtmosphereMask - -subroutine GRHydroPostSyncAtmosphereMask(CCTK_ARGUMENTS) - - implicit none - - DECLARE_CCTK_ARGUMENTS - DECLARE_CCTK_PARAMETERS - - CCTK_INT :: i,j,k - -!! This sets the integer atmo mask based on the real-valued (and sync'ed) atmo mask - - !$OMP PARALLEL DO PRIVATE(k,j,i) - do k = 1, cctk_lsh(3) - do j = 1, cctk_lsh(2) - do i = 1, cctk_lsh(1) - if ( atmosphere_mask_real(i,j,k) .ne. 0) then - atmosphere_mask(i,j,k) = 1 - end if - end do - end do - end do - !$OMP END PARALLEL DO - -end subroutine GRHydroPostSyncAtmosphereMask - - /*@@ @routine GRHydro_SetupMask @date Thu Jun 20 13:27:28 2002 @@ -147,7 +116,6 @@ endif atmosphere_mask = 0 - atmosphere_mask_real = 0 call CCTK_INFO("Setting up the atmosphere mask: all points are not_atmosphere") @@ -295,7 +263,6 @@ velz(i,j,k), eps(i,j,k), press(i,j,k), w_lorentz(i,j,k)) if (wk_atmosphere .eq. 0) then atmosphere_mask(i, j, k) = 0 - atmosphere_mask_real(i, j, k) = 0 end if endif Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: +0 -9 =================================================================== --- trunk/interface.ccl 2012-06-22 20:37:44 UTC (rev 371) +++ trunk/interface.ccl 2012-06-22 20:37:46 UTC (rev 372) @@ -424,15 +424,6 @@ atmosphere_mask } "Flags to say whether a point needs to be reset to the atmosphere" -# This real mask is set during UpdateAtmosphereMask and sync'ed afterwards (including possible interpatch interpolation) -# After syn'ing and before any con2prim fun, we set the integer mask above based on the real-valued mask. -# This ensures that any routine using the int mask is still correctly working. -real GRHydro_atmosphere_mask_real type = GF Timelevels = 1 tags='Prolongation="copy"' -{ - atmosphere_mask_real -} "Flags to say whether a point needs to be reset to the atmosphere. This is sync'ed (and possibly interpolated)!" - - int GRHydro_atmosphere_descriptors type=SCALAR { atmosphere_field_descriptor File [modified]: schedule.ccl Delta lines: +7 -11 =================================================================== --- trunk/schedule.ccl 2012-06-22 20:37:44 UTC (rev 371) +++ trunk/schedule.ccl 2012-06-22 20:37:46 UTC (rev 372) @@ -477,7 +477,6 @@ if (CCTK_EQUALS(evolution_method, "GRHydro")) { STORAGE:GRHydro_atmosphere_mask -STORAGE:GRHydro_atmosphere_mask_real STORAGE:GRHydro_atmosphere_descriptors schedule GRHydro_SetupMask AT CCTK_Initial BEFORE HydroBase_Initial @@ -1025,7 +1024,6 @@ SYNC: Y_e_con SYNC: lvel SYNC: lBvec - SYNC: GRHydro_atmosphere_mask_real # Must be synced otherwise we get inconsistent results when using multiple procs } "Select GRHydro boundary conditions" } else @@ -1041,7 +1039,6 @@ SYNC: psidc SYNC: GRHydro_cons_tracers SYNC: Y_e_con - SYNC: GRHydro_atmosphere_mask_real # Must be synced otherwise we get inconsistent results when using multiple procs } "Select GRHydro boundary conditions" # after a regrid Cactus relies on the boundary thorns to fill in the outer @@ -1053,6 +1050,11 @@ { } "Apply boundary conditions to primitives" + # This is necessary to provide initial guesses for Con2Prim after atmo reset (which happens only at last MoL post-step)! + schedule group GRHydro_PrimitiveBoundaries IN HydroBase_PostStep BEFORE HydroBase_Boundaries IF GRHydro::InLastMoLPostStep + { + } "Apply boundary conditions to primitives" + if(CCTK_IsImplementationActive("Coordinates")) { @@ -1158,20 +1160,14 @@ ### This is executed in the last MoL_PosStep ### ################################################################ -# First we need to convert real (and synchronized mask) to an integer mask! -schedule GRHydroPostSyncAtmosphereMask IN HydroBase_PostStep AFTER HydroBase_Boundaries -{ - LANG: Fortran -} "Set integer atmosphere mask from synchronized real atmosphere mask" - if (CCTK_Equals(Bvec_evolution_method,"GRHydro")) { - schedule GRHydro_AtmosphereResetM IN HydroBase_PostStep BEFORE HydroBase_Con2Prim IF GRHydro::InLastMoLPostStep + schedule GRHydro_AtmosphereResetM IN MoL_PostStep AFTER ADMBase_SetADMVars AFTER GRHydro_SetLastMoLPostStep BEFORE HydroBase_PostStep IF GRHydro::InLastMoLPostStep { LANG: Fortran } "Reset the atmosphere - MHD version" } else { - schedule GRHydro_AtmosphereReset IN HydroBase_PostStep BEFORE HydroBase_Con2Prim IF GRHydro::InLastMoLPostStep + schedule GRHydro_AtmosphereReset IN MoL_PostStep AFTER ADMBase_SetADMVars BEFORE HydroBase_PostStep IF GRHydro::InLastMoLPostStep { LANG: Fortran } "Reset the atmosphere"

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 371)
by rhaas＠tapir.caltech.edu 22 Jun '12

22 Jun '12

User: rhaas Date: 2012/06/22 03:37 PM Modified: /trunk/src/ GRHydro_ParamCheck.F90 Log: Revert "Revert "GRHydro: do not access whichpsidcspeed in ParamCheck unless evolve_MHD"" Ahem. I had reverted the wrong commit. This one was blameless. This reverts commit 0b5a561872841b50ea32da2b665a04b2984a8c16. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_ParamCheck.F90 Delta lines: +5 -1 =================================================================== --- trunk/src/GRHydro_ParamCheck.F90 2012-06-22 02:29:21 UTC (rev 370) +++ trunk/src/GRHydro_ParamCheck.F90 2012-06-22 20:37:44 UTC (rev 371) @@ -126,7 +126,7 @@ end if - if(clean_divergence.ne.0) then + if(evolve_MHD.ne.0 .and. clean_divergence.ne.0) then if (CCTK_EQUALS(psidcspeed,"char speed")) then whichpsidcspeed = 0 else if (CCTK_EQUALS(psidcspeed,"light speed")) then @@ -138,6 +138,10 @@ end if end if + if(clean_divergence.eq.0 .and. evolve_MHD.ne.0) then + call CCTK_WARN(1, "You activated divergence cleaning but do not evolve a magnetic field. At best, nothing should happen, be preapered for a segfault otherwise.") + end if + if (Tmunu_damping_radius_min .gt. Tmunu_damping_radius_max) then call CCTK_WARN(0, "Minimum damping radius is greater than maximum damping radius!") end if

1 0

[svn:einsteintoolkit] EinsteinExamples/trunk/par/ (Rev. 33)
by rhaas＠tapir.caltech.edu 22 Jun '12

22 Jun '12

User: rhaas Date: 2012/06/22 12:22 PM Modified: /trunk/par/ qc0-mclachlan.par Log: add puncture tracking and wave extraction to qc0 example file File Changes: Directory: /trunk/par/ ====================== File [modified]: qc0-mclachlan.par Delta lines: +76 -1 =================================================================== --- trunk/par/qc0-mclachlan.par 2012-05-28 22:55:03 UTC (rev 32) +++ trunk/par/qc0-mclachlan.par 2012-06-22 17:22:37 UTC (rev 33) @@ -4,7 +4,7 @@ Cactus::highlight_warning_messages = no Cactus::terminate = "time" -Cactus::cctk_final_time = 100.0 +Cactus::cctk_final_time = 200.0 @@ -32,8 +32,10 @@ ActiveThorns = "Slab" +ActiveThorns = "TGRtensor" + ActiveThorns = "SummationByParts" SummationByParts::order = 4 @@ -331,6 +333,26 @@ +ActiveThorns = "WeylScal4 Multipole" +WeylScal4::offset = 1e-8 +WeylScal4::fd_order = "4th" +WeylScal4::verbose = 0 + +Multipole::nradii = 8 +Multipole::out_every = 128 +Multipole::radius[0] = 15 +Multipole::radius[1] = 30 +Multipole::radius[2] = 40 +Multipole::radius[3] = 50 +Multipole::radius[4] = 60 +Multipole::radius[5] = 70 +Multipole::radius[6] = 80 +Multipole::radius[7] = 90 +Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='Psi4'}" +Multipole::l_mode = 4 +Multipole::m_mode = 4 + + ActiveThorns = "AHFinderDirect" AHFinderDirect::find_every = 128 @@ -390,6 +412,23 @@ +ActiveThorns = "QuasiLocalMeasures" + +QuasiLocalMeasures::verbose = yes +QuasiLocalMeasures::interpolator = "Lagrange polynomial interpolation" +QuasiLocalMeasures::interpolator_options = "order=4" +QuasiLocalMeasures::spatial_order = 4 + +QuasiLocalMeasures::num_surfaces = 6 +QuasiLocalMeasures::surface_index[0] = 0 +QuasiLocalMeasures::surface_index[1] = 1 +QuasiLocalMeasures::surface_index[2] = 2 +QuasiLocalMeasures::surface_index[3] = 3 +QuasiLocalMeasures::surface_index[4] = 4 +QuasiLocalMeasures::surface_index[5] = 5 + + + ActiveThorns = "CarpetIOBasic" IOBasic::outInfo_every = 128 @@ -397,6 +436,8 @@ IOBasic::outInfo_vars = " Carpet::physical_time_per_hour ML_ADMConstraints::H + SphericalSurface::sf_radius + QuasiLocalMeasures::qlm_spin[0] " @@ -414,9 +455,19 @@ ADMBase::shift ADMBase::dtlapse ADMBase::dtshift + WEYLSCAL4::Psi4r + WEYLSCAL4::Psi4i ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_mom SphericalSurface::sf_radius + QuasiLocalMeasures::qlm_newman_penrose + QuasiLocalMeasures::qlm_weyl_scalars + QuasiLocalMeasures::qlm_ricci_scalars + QuasiLocalMeasures::qlm_twometric + QuasiLocalMeasures::qlm_killing_vector + QuasiLocalMeasures::qlm_killed_twometric + QuasiLocalMeasures::qlm_invariant_coordinates + QuasiLocalMeasures::qlm_3determinant " @@ -438,6 +489,8 @@ ADMBase::shift ADMBase::dtlapse ADMBase::dtshift + WEYLSCAL4::Psi4r + WEYLSCAL4::Psi4i ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_mom SphericalSurface::sf_active @@ -446,6 +499,11 @@ SphericalSurface::sf_radius SphericalSurface::sf_origin SphericalSurface::sf_coordinate_descriptors + QuasiLocalMeasures::qlm_state + QuasiLocalMeasures::qlm_grid_int + QuasiLocalMeasures::qlm_grid_real + QuasiLocalMeasures::qlm_scalars + QuasiLocalMeasures::qlm_multipole_moments " IOASCII::out1D_every = 128 @@ -457,9 +515,24 @@ ADMBase::shift ADMBase::dtlapse ADMBase::dtshift + WEYLSCAL4::Psi4r + WEYLSCAL4::Psi4i ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_mom SphericalSurface::sf_radius + QuasiLocalMeasures::qlm_shapes + QuasiLocalMeasures::qlm_coordinates + QuasiLocalMeasures::qlm_tetrad_l + QuasiLocalMeasures::qlm_tetrad_n + QuasiLocalMeasures::qlm_tetrad_m + QuasiLocalMeasures::qlm_newman_penrose + QuasiLocalMeasures::qlm_weyl_scalars + QuasiLocalMeasures::qlm_ricci_scalars + QuasiLocalMeasures::qlm_twometric + QuasiLocalMeasures::qlm_killing_vector + QuasiLocalMeasures::qlm_killed_twometric + QuasiLocalMeasures::qlm_invariant_coordinates + QuasiLocalMeasures::qlm_3determinant " IOASCII::out2D_every = 128 @@ -485,6 +558,8 @@ ADMBase::shift ADMBase::dtlapse ADMBase::dtshift + WEYLSCAL4::Psi4r + WEYLSCAL4::Psi4i ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_mom "

1 0

[svn:einsteintoolkit] EOS_Omni/trunk/ (Rev. 63)
by rhaas＠tapir.caltech.edu 22 Jun '12

22 Jun '12

User: rhaas Date: 2012/06/22 12:18 PM Added: /trunk/src/nuc_eos/ dumpASCIItable.F90 /trunk/test/ TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.h5, eos_omni.par, test.ccl /trunk/test/eos_omni/ TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.asc Modified: /trunk/ param.ccl, schedule.ccl /trunk/src/nuc_eos/ make.code.defn Log: EOS_Omni: add table reader test This adds a test (plus sample hdf5 table) using a small routine EOS_OMNI_dumptable to output the read in data as ASCII into a user selected file. Code and table kindly provided by Evan O'Connor. Small changes afterwards by Roland Haas. File Changes: Directory: /trunk/test/eos_omni/ ================================ File [added]: TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.asc Delta lines: <Unavailable> Directory: /trunk/src/nuc_eos/ ============================== File [added]: dumpASCIItable.F90 Delta lines: <Unavailable> File [modified]: make.code.defn Delta lines: <Unavailable> Directory: /trunk/test/ ======================= File [added]: TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.h5 Delta lines: <Unavailable> File [added]: eos_omni.par Delta lines: <Unavailable> File [added]: test.ccl Delta lines: <Unavailable> Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: <Unavailable> File [modified]: schedule.ccl Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 370)
by rhaas＠tapir.caltech.edu 21 Jun '12

21 Jun '12

User: rhaas Date: 2012/06/21 09:29 PM Modified: /trunk/src/ GRHydro_ParamCheck.F90 Log: Revert "GRHydro: do not access whichpsidcspeed in ParamCheck unless evolve_MHD" This reverts commit c2122aa2b4762388eb90cc159750aa4d418b112d. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_ParamCheck.F90 Delta lines: +1 -5 =================================================================== --- trunk/src/GRHydro_ParamCheck.F90 2012-06-22 02:29:18 UTC (rev 369) +++ trunk/src/GRHydro_ParamCheck.F90 2012-06-22 02:29:21 UTC (rev 370) @@ -126,7 +126,7 @@ end if - if(evolve_MHD.ne.0 .and. clean_divergence.ne.0) then + if(clean_divergence.ne.0) then if (CCTK_EQUALS(psidcspeed,"char speed")) then whichpsidcspeed = 0 else if (CCTK_EQUALS(psidcspeed,"light speed")) then @@ -138,10 +138,6 @@ end if end if - if(clean_divergence.eq.0 .and. evolve_MHD.ne.0) then - call CCTK_WARN(1, "You activated divergence cleaning but do not evolve a magnetic field. At best, nothing should happen, be preapered for a segfault otherwise.") - end if - if (Tmunu_damping_radius_min .gt. Tmunu_damping_radius_max) then call CCTK_WARN(0, "Minimum damping radius is greater than maximum damping radius!") end if

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits June 2012