aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
...
* GRhydro: Pressure in polytypeM_pt was off by factor inv_sqrt_detgrhaas2013-05-11
| | | | | | From: srichers <srichers@caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@510 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add SAVE attribute to module variablesrhaas2013-04-18
| | | | | | | | | | | | I forgot the SAVE options since my initial google search seemed to indicate that they are redundant. It turns out that is only true for Fortran 2008 and an earlier Fortran compiler would in theory be free do have variables use their values once the last procedure using the module goes out of scope. This seems to be not done by any compiler though. From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@508 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove make.code.deps, coeffs_allocated, unused includes rhaas2013-04-18
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@507 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Set InLastMoLPostStep during initial data setupeschnett2013-04-13
| | | | | | | | The previous code accessed MoL_Intermediate_Step, which is undefined during initial data setup. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@506 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove superfluous "divb(i,j,k) +" from divb computationrhaas2013-04-10
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@505 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Critical bugfix: Con2Prim was probably never executed inrhaas2013-04-10
| | | | | | | | | | | | post_recover_variables or for new points in post_regrid since excision mask does typically not get initialized before MoL_PostStep! Also: Instead of aborting Con2Prim, set points to atmopshere for a point where hydro_excision_mask > 0. From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@504 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: support reconstruct_Wv in TVD reconstructorrhaas2013-03-28
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@501 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Finish the implementation of con2prim using the entropy equation.rhaas2013-03-28
| | | | | | | | * Tested successfully pointwisely. It still needs to be tested on evolution. From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@500 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: move critical section in prim2conM into if statementrhaas2013-03-28
| | | | | | | | | otherwise the code will enter and leave the critical section during each loop iteration From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@499 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: explicitly cast from CCTK_REAL to int to avoid warningrhaas2013-03-28
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@498 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: change CCTK_WARN(0,...) to CCTK_ERRORrhaas2013-03-28
| | | | | | | | | this should tell the compiler that the function never returns and help optimizing the code by indicating eg slow branches of if statements From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@497 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: pass pressure into eigenvalue routine for MHDrhaas2013-03-28
| | | | | | | | bugfix for hot eos case, just more efficient for cold eos From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@496 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove unused variablesrhaas2013-03-28
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@495 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: reduce size of critical section in eigenbalues routinerhaas2013-03-28
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@494 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: convert more tabs to (8) spaces. No other changes.rhaas2013-03-28
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@493 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: implemented H viscosity for HLLE solver to eliminate the carbuncle ↵rhaas2013-03-28
| | | | | | | | instability. Tested with shocktube and TOV. Not yet tested for problem where carbuncles occur. From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@491 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add test case for WENO5rhaas2013-03-06
| | | | | | | | this is just a copy of the shock tube test using WENO From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@489 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Added MP5 reconstruction method. Tested with TOV and shocktube.rhaas2013-03-06
| | | | | | From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@487 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Reconstruct Wv for WENO5.rhaas2013-03-06
| | | | | | From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@486 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* convert tabs to (8) spaces. No other changes.knarf2013-03-04
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@485 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* remove unused variablesknarf2013-03-04
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@484 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove commented out code, update docsrhaas2013-02-11
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@481 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: declare charXXX_dc variables in HLLEM OMP privaterhaas2013-02-11
| | | | | | From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@477 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add fuzzy comparison for evolution_maskrhaas2013-02-11
| | | | | | | | | evolution_mask is every only set to 0d0 or 1d0 (explicitly) right now. Should this change, this tests needs to be amended. From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@474 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * add option to reconstruct temperature to GRHydrorhaas2013-02-11
| | | | | | | | | | | | This does not yet work with the MHD part of the code, but Philipp is going to change this. * Improve TOVSolverHot; in particular, give functions unique names to avoid duplicate symbols with TOVSolverC From: Christian Ott <cott@bethe.tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@473 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: disable in-AH polytrope call by default for Avecrhaas2013-02-11
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@472 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add option to automatically use Evolution_Mask if possiblerhaas2013-01-20
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@469 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* 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
* 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