diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:33 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:33 +0000 |
| commit | 48b053d41b2aa9ce720d443d288eef535fb9651d (patch) | |
| tree | bb5274de4dfb94de758bafc0ee77e6b5f499fb06 /src/GRHydro_RegisterVars.cc | |
| parent | d6e45e79d92f3d2adf3eb4ebcc7641db07119f30 (diff) | |
GRHydro: Implement entropy evolution (separated from temp evolution)
From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@455 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_RegisterVars.cc')
| -rw-r--r-- | src/GRHydro_RegisterVars.cc | 16 |
1 files changed, 10 insertions, 6 deletions
diff --git a/src/GRHydro_RegisterVars.cc b/src/GRHydro_RegisterVars.cc index 48a8651..95a7053 100644 --- a/src/GRHydro_RegisterVars.cc +++ b/src/GRHydro_RegisterVars.cc @@ -92,9 +92,13 @@ extern "C" void GRHydro_Register(CCTK_ARGUMENTS) } register_evolved("GRHydro::Bcons", "GRHydro::Bconsrhs"); if(clean_divergence) { - register_evolved("GRHydro::psidc" , "GRHydro::psidcrhs"); + register_evolved("GRHydro::psidc" , "GRHydro::psidcrhs"); } } + // entropycons + if(CCTK_EQUALS(entropy_evolution_method,"GRHydro")){ + register_evolved("GRHydro::entropycons" , "GRHydro::entropyrhs"); + } // tau if (CCTK_EQUALS(GRHydro_eos_type, "General")) @@ -118,7 +122,7 @@ extern "C" void GRHydro_Register(CCTK_ARGUMENTS) { register_constrained("admbase::shift"); register_evolved("GRHydro::GRHydro_coords", - "GRHydro::GRHydro_coords_rhs"); + "GRHydro::GRHydro_coords_rhs"); } else if(GRHydro_MaxNumSandRVars != 0) // hack to save some memory since we "know" that someone else will register these as constrained { @@ -129,15 +133,15 @@ extern "C" void GRHydro_Register(CCTK_ARGUMENTS) // tracer if (evolve_tracer != 0) register_evolved("GRHydro::GRHydro_cons_tracers", - "GRHydro::GRHydro_tracer_rhs"); - + "GRHydro::GRHydro_tracer_rhs"); + // note that we set in pararamcheck // evolve_y_e and evolve_temper, but MoLRegister // happens before paramcheck... if (CCTK_EQUALS(Y_e_evolution_method,"GRHydro")) { register_constrained("HydroBase::Y_e"); register_evolved("GRHydro::Y_e_con", - "GRHydro::Y_e_con_rhs"); + "GRHydro::Y_e_con_rhs"); } if (CCTK_EQUALS(temperature_evolution_method,"GRHydro")) { @@ -157,7 +161,7 @@ extern "C" void GRHydro_Register(CCTK_ARGUMENTS) if (CCTK_EQUALS(Bvec_evolution_method, "GRHydro")) { register_constrained("HydroBase::Bcons"); if(clean_divergence) { - register_constrained("GRHydro::psidc"); + register_constrained("GRHydro::psidc"); } } } |