aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
* GRHydro: Let the driver functions to take care of retval initializationrhaas2013-01-14
| | | | | | From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@448 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: reset points inside the excision region to atmosphere inrhaas2013-01-14
| | | | | | | | | InitialAtmosphereReset(M) so that the ID thorns do not have to know about he excised region. From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@447 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove old "hopefullyfixed" con2primM_pt version.rhaas2013-01-14
| | | | | | From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@446 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* only define validate_x() if DEBUG_CON2PRIMM is set because it is currently ↵knarf2012-12-17
| | | | | | unused and trigger warnings, but might be useful for debugging git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@441 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* remove unused variable; instead of using aint which truncates but doesn't ↵knarf2012-11-23
| | | | | | convert type, use int which does both -> silences compiler warning git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@440 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Introduce inequalities cons variables needs to satisfy when B/=0.rhaas2012-11-09
| | | | | | | | | | Introduce temporary prim variables as well. * Please refer to Illinois' paper (arXiv:1112.0568) appendices for details. From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@439 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Check if 3-metric is definite positive before con2prim inversion.rhaas2012-11-09
| | | | | | | | | Only do this deep inside the horizon. Reset it to conformally flat if it isn't. From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@438 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Introduce GRHydro_tau_min = 10^-10 max(tau(t=0)) to be used in tau ↵rhaas2012-11-09
| | | | | | | | rescaling in con2prim failures. From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@437 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* This routine hopefully fixes the polytrope calls to magnetized Con2Prim ↵rhaas2012-11-09
| | | | | | | | while leaving the hybrid EOS behavior unchanged From: Josh Faber <jfaberrit@jfaber3.local> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@436 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRhydro: add recovstruct_Wv to magnetized PPM routinerhaas2012-11-09
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@434 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: make sure to pass distinct dummy arguments to Fortran routinesrhaas2012-11-09
| | | | | | | | since Fortran does not allow aliased arguments From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@433 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add extensive information output to check_C2P_failedrhaas2012-11-09
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@432 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove usage of EOS_Module, poly_gamma_inirhaas2012-11-09
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@431 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: undo revision 426 since it causes test failures since poltytropes userhaas2012-10-25
| | | | | | | pressure as an initial guess git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@428 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Unify function interface for GRHydro_Con2PrimM_pt.c and ↵rhaas2012-10-25
| | | | | | | | | | | | | | | GRHydro_Con2PrimM_pt_EOSOmni.c * Interestingly Cactus doesn't complain about the function interface in interface.ccl. It has the old set of arguments for that function and declare it as written in Fortran instead of C. * Cylindrical explosion works again with GRHydro_Con2PrimM_pt.c routine. You can just swap both routines in make.code.defn From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@427 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Don't sync Lorentz factor and pressure since they are always ↵rhaas2012-10-25
| | | | | | | | computed in Con2Prim and are not needed for initial guesses. From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@426 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
* 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
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-10 02:03:21 +0000