aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
* 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: 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
* 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: 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
* * 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: 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: 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
* * 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: 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
* * 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
* Revert "* improve handling of cs2 < 0 in the "standard" (non-MHD, non-hot ↵rhaas2012-06-22
| | | | | | | | | | EOS) eigenvalue routine." actually this introduce the real valued atmosphere mask. This reverts commit deccf528b463820b615dfa73e46e137cd59c65be. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@372 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Revert "Revert "GRHydro: do not access whichpsidcspeed in ParamCheck unless ↵rhaas2012-06-22
| | | | | | | | | | evolve_MHD"" Ahem. I had reverted the wrong commit. This one was blameless. This reverts commit 0b5a561872841b50ea32da2b665a04b2984a8c16. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@371 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Revert "GRHydro: do not access whichpsidcspeed in ParamCheck unless evolve_MHD"rhaas2012-06-22
| | | | | | This reverts commit c2122aa2b4762388eb90cc159750aa4d418b112d. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@370 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRydro: avoid implict SAVE attribute in Con2PrimPolytyperhaas2012-06-21
| | | | | | | | | initializing a variable makes it SAVE which is incomaptible with OpenMP. Also fix some 1-element-array vs. scalar variable type issues. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@364 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * introduce parameters GRHydro_Y_e_min and GRHydro_Y_e_max to fix issues ↵rhaas2012-06-21
| | | | | | with EOS tables that don't extend sufficiently far. Enforcement happens in Con2Prim. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@363 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* really improve handling of cs2 < 0 in the "standard" (non-MHD, non-hot EOS)rhaas2012-06-20
| | | | | | | | | | | | | | | | | | | | | | | eigenvalue routine. the previous patch with this same commit message (r359) has the wrong commit message. It should instead read: --8<-- Introduce real-valued atmo mask that can be interpolated. Synchronize real-valued atmo mask and update integer-valued atmo mask based on synchronized value. This fix ensures that the results are independent of the number of processors used. It also gets rid of some extra primitive syncs that are otherwise necessary when sync_conserved_only = yes. Does not cause any new test failures (Roland). Patch by Christian Reisswig. --8<-- Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@361 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: do not access whichpsidcspeed in ParamCheck unless evolve_MHDrhaas2012-06-20
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@360 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * improve handling of cs2 < 0 in the "standard" (non-MHD, non-hot EOS) ↵rhaas2012-06-20
| | | | | | | | | eigenvalue routine. Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@359 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Change the characteristic speed for psidc.rhaas2012-06-20
| | | | | | | | | | | | | | It was set to +-c when it should be the maximum coming from the constraint that the normal vector to the characteristic hypersurface of a wave propagating along an arbitrary x direction be spacelike. This results into the following bound for the wavespeed: +- \alpha sqrt(\gamma^{xx}) - \beta^x) Patch by Bruno Mundim. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@358 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: do not use integers as booleansrhaas2012-06-20
| | | | | | | gfortran (and the Fortran standard I guess) does not allow this. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@357 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* C2P_hot: Resort to bisection as well (if user requested so).rhaas2012-06-20
| | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@355 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: C2P: resort to bisection if Newton-Raphson fails. This is onlyrhaas2012-06-20
| | | | | | | | | activated if user requested this. Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@354 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add file missing from r343rhaas2012-06-06
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@352 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Fix slope limiters and get rid of broken cases;rhaas2012-06-05
| | | | | | | | | | | there is a related patch in TRAC at the moment Patch by Josh Faber. Some additional checks by Roland Haas. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@351 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: register TmunuBase variables as MoL constraint.rhaas2012-06-05
| | | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@350 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * allow disabling check of the GRHydro_C2P_failed maskrhaas2012-06-05
| | | | | | | | | | | on coarse reflevels via GRHydro_c2p_warn_from_level parameter. Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@348 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: fix format statements in hot C2P.rhaas2012-06-05
| | | | | | | | | | Note: This should be propagated to ET! Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@347 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Added option to constrain evolution to 1D byrhaas2012-06-05
| | | | | | | | | | projecting out non-radial parts of the conserved velocity RHSs. Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@346 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: do not register any SandR variables if MaxNumSandR ==0rhaas2012-06-05
| | | | | | | | | this is useful when running with a spacetime code that will register the ADMBase variables as constrained anyway, so one can save the memory used for SandR scratch space. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@345 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Verbosity on which boundary selction routine is executed.rhaas2012-06-05
| | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@344 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Schedule AtmosphereReset in LAST Mol_PostSteprhaas2012-06-05
| | | | | | | | | | | | This saves some syncs. This also gets rid of the akward parameter "disable_sync_and_tmunu_after_mol". Original patch by Christian Reisswig. Adapted to not rely on MoL support by Roland Haas. Bug fixed by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@343 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Added verbosity parameter.rhaas2012-06-05
| | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@342 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: sync those primitives which are needed for initial Con2Prim guessesrhaas2012-06-05
| | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@341 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Added parameter to perform sync on conserved quantities only.rhaas2012-06-05
| | | | | | | | | | | | | | Primitives are only sync'ed after regridding to get initial guesses. Schedule Primitive boundary call in PostRegridInitial if sync_only_coserved == yes. Fixed issue with the sync call in CCTK_EVOL after MoL_Evolution when sync_conserved_only==yes: Conserved vars must be synced here as well, otherwise conserved vars are reset to atmospshere whereas primitives are not! This is wrong. Patches by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@340 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: a bit of a clean-up of unused variables to avoid compiler complaintsbmundim2012-05-17
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@336 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45