aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
...
* GRHydro: Improved WENO implementation: (i) Pick adaptive epsilon-parameter ↵rhaas2013-01-14
| | | | | | | | according to WHAM code paper, (ii) Catch epsilon<0 for non-microphysical EOS. Tested with shock tube and TOV. From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@466 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Implemented WENO5. Tested with a shock tube.rhaas2013-01-14
| | | | | | From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@465 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: output v2 when warning about NaN in w_lorentzrhaas2013-01-14
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@464 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: use distinct dummys for calc_vlow_blowrhaas2013-01-14
| | | | | | | | also check actual timelevels information for storage of bcom From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@463 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: correct tau when resetting to atmosphererhaas2013-01-14
| | | | | | | | before a term rho*b2 was present when the correct one is only b2 From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@462 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: clear divB everywhere rather than just in interiorrhaas2013-01-14
| | | | | | From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@461 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: CCTK_Equals -> CCTK_EQUALSrhaas2013-01-14
| | | | | | | | | since only the later (a macro) contains the required ".ne. 0" to turn the return value into a logical type From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@460 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Fixing comment to reflect what is used in the code (sqrt(\gamma) - ↵rhaas2013-01-14
| | | | | | | | | | > \gamma). not code change From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@459 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Starting to push vector potential changesrhaas2013-01-14
| | | | | | | | | | * incl boundary calls for Avec * disable call to BvecfromAvec from within loop, this one is a schedule function call now. From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@458 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Improved Tmunu damping.rhaas2013-01-14
| | | | | | | | | * faster falloff to make sure to damp to atmosphere at Tmunu_damping_radius_max From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@457 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Add basic vector potential supportrhaas2013-01-14
| | | | | | | | | | Basic cell-centered, algebraic gauge vector potential method with place-holders for lorenz gauge. Initial Avec constrained to poloidal at the moment. From: Tanja Bode <tanja.bode@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@456 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Implement entropy evolution (separated from temp evolution)rhaas2013-01-14
| | | | | | From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@455 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: optionally use CarpetEvolutionMask::evolution_mask to skip validity ↵rhaas2013-01-14
| | | | | | | | tests in regions that are not needed for evolution From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@454 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Introduce variable atmopshere level / tolerance as function of ↵rhaas2013-01-14
| | | | | | | | radius. Off by default. From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@453 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: disable in-AH polytrope call by defaultrhaas2013-01-14
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@452 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: call polytrope EOS in a region dominated by mag. field.rhaas2013-01-14
| | | | | | From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@451 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Call pointwise polytype upon con2primM failure. If it fails just ↵rhaas2013-01-14
| | | | | | | | use previous primitive values and get the cons from them. Basically freezing their values at previous time step. From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@450 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: call pointwise con2primM with polytrope eos handle deep inside AH ↵rhaas2013-01-14
| | | | | | | | upon inversion failure of gamma law eos. From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@449 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* 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