aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
* use an early atmosphere reset to ensure past timelevels are initialized ↵knarf2012-10-23
| | | | | | correctly. This does not work for relative settings, but at least for absolute ones. See ticket #1071 for more information git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@423 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: fix uninitialized pointerrhaas2012-10-01
| | | | | | Found by Frank Loeffler. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@422 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Correct format specifiereschnett2012-09-12
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@421 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* put { on new line after elseknarf2012-09-04
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@420 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: clear atmosphere_mask_real in MHD AtmosphereResetrhaas2012-08-27
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@419 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add OpenMP parallelization to IntialAtmosphereResetrhaas2012-08-27
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@418 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove unused filerhaas2012-08-27
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@417 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Option to reconstruct W*vel (W = Lorentz factor) instead of just ↵rhaas2012-08-27
| | | | | | | | | | | | | the primitive velocity 'vel'. This ensures that the Lorentz factor can never become unphysical during reconstruction. This patch is necessary for NS collapse when ePPM+Refluxing is used. The idea is due to Roland. From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@416 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Changes fixing issues with hot EOS treatment:rhaas2012-08-27
| | | | | | | | | | | | | | | 1) Set correct keytemp in C2P 2) Fix OpenMP private variable declaration in HLLEM 3) Fix velocity arguments for Y_e 1D PPM reconstruction 4) If in atmosphere and if evolve a Y_e also reset Y_e with cell averaged value in Reconstruction 5) Fix index permutation bug in ReconstructPoly routines for divergence cleaning field 6) Remove extra comment character in UpdateMaskM 7) Fix a missing velocity permutation when calling 1D reconstruction for Y_e. 8) Insert missing sqrt(det) factor for pressure term in resetting tau when temperature got too cold with a hot EOS. From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@415 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Adding back OpenMP support to HLLEM.F90rhaas2012-08-27
| | | | | | From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@414 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: improve error handling in C2P hot routinerhaas2012-08-27
| | | | | | | | | | | * improve error handling in C2P hot routine * further bugfix to tau reset treatment in hot part * fix an inconsistency in hot-EOS Con2Prim routine -- when resetting utau, also reset tau for consistency From: Christian Ott <cott@bethe.tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@413 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: use inverse Jacobian to transform Bvec back to global basisrhaas2012-08-27
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@412 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: MP pointer stuff for MHDrhaas2012-08-27
| | | | | | | | this time with minimal replacements in SourceM.F90 From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@411 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Changes to make GRHydro MHD work with nuclear/hot equation of state:rhaas2012-08-27
| | | | | | | | | | | | | | | | | | 1) Switch to EOSOmni pointwise C2P routine and modify where necessary. 2) Modify Con2PrimM.F90 to allow for the evolution of temperature and adjust the wrapper routine. 3) Create EigenProblemM_hot pointwise routine and call that from HLLEM.F90 when temperature is evolved. Additionally adjust HLLEM where necessary. 4) Adjust InterfacesM.h to incorporate the newly created functions. 5) Fix a loop problem (not taking into account constraint transport) in PPM reconstruction of Y_e 6) Introduce Prim2ConM_hot and call this pointwise routine from Prim2ConM.F90 when temperature is evolved. Additionally also make this routine available to initial data routine in GRHydro_InitData 7) Adjust loops in GRHydro_PoloidalMagFieldM.F90 to not set boundary points it cannot set but instead call boundary group afterwards! Pay attention as this will not work with boundary conditions set to "none" in MHD case anymore but is the correct thing to do. 8) Allow StarMapper to extend HydroBase::initial_hydro = "starmapper". 9) Smaller fixes. From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@410 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: fix inverted logic in ParamWarn about divergence cleaningrhaas2012-08-14
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@409 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Conditionally register Tmunu with MoL based on Tmunu storage being set.bmundim2012-08-09
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@408 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add forgotten label in error message, whitespace changerhaas2012-08-04
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@407 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add test case for slow MoL sector of RK4-RK2rhaas2012-08-04
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@406 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Option to use the slow sector of multirate RK methods.rhaas2012-08-04
| | | | | | | | NOTE: This patch requires support from MoL. Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@405 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Move MPI support from flesh to a thorneschnett2012-07-21
| | | | | | | | | | Move MPI support from flesh to thorn ExternalLibraries/MPI. This also requires thorns that call MPI directly to declare this in their configuration.ccl. Existing configurations using MPI need to include ExternalLibraries/MPI into their thorn list. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@404 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: only loop over points where rhs is defined for UpdateMaskrhaas2012-07-17
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@403 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove comparison of logical to .true.rhaas2012-07-17
| | | | | | | | | | | | | | | | | I appreciate the effort to make the code compile with gcc by no longer using integers in place of logicals. I really do. However it seems that F0RTRAN in its infinite wisdom requires you to write ... .and. mustbisect .eqv. .true. (note the all important 'v') when testing logicals. I have instead opted to remove the whole .eqv. term in favor of ... .and. mustbisect From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@402 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Con2Prim improvement to bisection algorithm. This fixes some NANs ↵rhaas2012-07-17
| | | | | | | | associated with NS collapse test case when full GR is used. From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@401 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: schedule AtmosphereReset before PrimitiveInitialGuessesBoundariesrhaas2012-07-17
| | | | | | | | which used to be called PrimitiveBoundaries From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@400 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * schedule GRHydroCopyIntegerMask in POST_RECOVER_VARIABLES before ↵rhaas2012-07-17
| | | | | | | | | | | MoL_PostStep rather than HydroBase_PostStep, since the latter is scheduled IN MoL_PostStep -> rely on correct scheduling rather than on luck to get the correct atmosphere map after recovery. From: Christian Ott <cott@bethe.tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@399 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: ePPM: Allow for tolerance in vel^2 limiting.rhaas2012-07-17
| | | | | | From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@398 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: ePPM: Corrected a logic bug in vel2 constraint. Also, use averaged ↵rhaas2012-07-17
| | | | | | | | metric. From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@397 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: ePPM: Constrain reconstructed vel to below speed of light!rhaas2012-07-17
| | | | | | From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@396 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * ensure consistency in calls to prim2con_hot (fixes issue with my previous ↵rhaas2012-07-17
| | | | | | | | patch) From: Christian Ott <cott@bethe.tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@395 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * add parameter used in Prim2Con_hotrhaas2012-07-17
| | | | | | From: Christian Ott <cott@bethe.tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@394 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * improve treatment of hot EOS warningrhaas2012-07-17
| | | | | | From: Christian Ott <cott@bethe.tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@393 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: limit Carpet based test to run with two processes onlyrhaas2012-07-11
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@392 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove headers from tov_carpet_refined_nosync test datarhaas2012-07-09
| | | | | | | | | No actual data was updated. data was generated without sync_only_primities but we compare against sync_only_primitives which was always the intend of the test. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@391 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove OpenMP options from Riemann1D buildrhaas2012-07-08
| | | | | | | | | | to maek it compiel on Kraken without specifying -openmp to the linker and complaining about ifcore vs ifcoremt issues sort object file before library files for static linkage avoid compiler warning about using libm math rather than mkl math git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@390 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: clear atmosphere_mask_real in MHD AtmosphereResetrhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@389 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove SAVE attribtute from _p and _pp varsrhaas2012-07-05
| | | | | | | this was harmless since the pointers to grid functions affected were OMP shared anyway git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@388 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove unused filesrhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@387 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove unused variablesrhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@386 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove superflous code, unused variablesrhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@385 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: fix apparent typos in scheduling of GRHydro_InitAtmosMaskrhaas2012-07-05
| | | | | | | | this only matters if wk_atmosphere is used and is based on the argument that it makes no sense to do this in postrestrictinitial (at least when vertex centering is used) but makes sense in postregridinitial git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@384 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: copy integer atmosphere mask in Post_Recover_Variablesrhaas2012-07-05
| | | | | | | | | | | this is required since the real valued atmosphere mask is not checkpointed. It is valid since atmosphere_mask_real is non-zero wherever atmosphere_mask is non-zero due to GRHydroUpdateAtmosphereMask and atmosphere_mask is non-zero wherever atmosphere_mask_real is non-zero because of GRHydroPostSyncAtmosphereMask, ie. the two masks are non-zero/zero at the same points. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@383 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * improve error message formatting in hot con2primrhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@382 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: text-only typos in schedule descriptionrhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@381 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * change comparison of double with integer to double with doublerhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@380 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Fixed mislabeled scheduling comment (cosmetic change).rhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@379 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Sync ALL prinmitives after regridding. This ensures that all ↵rhaas2012-07-05
| | | | | | variables are known on the finer grid when sync_conserved_only = yes. This is fix only affects runs which use a microphysical EOS. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@378 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Fixes to atmosphere mask and the way we sync.rhaas2012-07-05
| | | | | | | | | | | | | These are the changes as discussed with Roland. In particular, we sync the atmo mask in an extra call before all other syncs and then do atmo reset. Furthermore, primitives must always be synced in last post step in order to get initial guesses in the buffer zones via prolongation. Atmo reset must be done before prolongation to ensure we don't overwrite "good" data from the coarse grid! git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@377 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * fix range error in GRHydro_Y_e_maxrhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@376 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * add limit on values of Y_e also in Prim2Conrhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@375 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove "use" for con2prim interfaces from pointwise routinesrhaas2012-06-26
| | | | | | | | | since the Fortran standard does not allow this and gfortran 4.7 actually chokes on it. Pointed out by Erik Schnetter. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@374 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45