aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
* 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
* 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: 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
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-06 18:39:36 +0000