aboutsummaryrefslogtreecommitdiff
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
* GRHydro: use timelevels parameter in STORAGE statementrhaas2013-01-14
| | | | | | | | rather than if cascade to test for each value From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@445 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove any traces of non-functional Fisheye supportrhaas2013-01-14
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@444 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: clean up schedule - no-op changesrhaas2013-01-14
| | | | | | | | these mostly refer to groups that used to exist but no longer do From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@443 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: make more paramters steerablerhaas2013-01-14
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@442 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
* ng the interface for magentized con2prim for consistencyrhaas2012-11-09
| | | | | | From: Josh Faber <jfaberrit@jfaber3.local> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@435 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: restrict Carpet based tests to 2 MPI processesrhaas2012-10-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@429 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: Get rid of unused ePPM parameter.rhaas2012-10-25
| | | | | | From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@425 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Add test suites for MHD coderhaas2012-10-25
| | | | | | | | these are all based on the shocktube and cylindrical explosion tests of the ET MHD paper of the same name. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@424 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* 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