aboutsummaryrefslogtreecommitdiff
path: root/src/GRHydro_PPM.F90
Commit message (Collapse)AuthorAge
* GRHydro: write floating point operations in left/right symmetric mannerrhaas2013-07-06
| | | | | | | | | | | | | | | | * write DIFF_X_4 in an explicitly symmetric manner this ensures that the order of evaluation is such that symmetry in i+/-1 is preserved by the numerical derivative * add code to make ePPM and Source more symmetric this is an attempt to better preserve the CoM in long runs by eliminating asymmetries between left and right. This is fragile in that it does not make the code more robust against existing roundoff-level asymmetries, it only prevents them from developping. From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@553 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove unreachable coderhaas2013-07-06
| | | | | | From: Christian Reisswig <reisswig@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@548 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add option to switch to oPPM on a given refinement levelrhaas2013-07-06
| | | | | | | | | | | | When reconstruction method is enhanced PPM, this change allows using old PPM on one refinement level, whereas the rest of the reflevels will use enhanced PPM. It is only for hydro, but could be extended to MHD also. From: abdik <abdik@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@545 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: fixed critical bug in temperature reconstruction in ePPMrhaas2013-05-29
| | | | | | From: Christian David Ott <cott@zwicky.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@532 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: 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
* 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
* * 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: 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: 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
* Cleanup of PPM routine. Remove unused PPM parameter. Only apply stricter ↵reisswig2012-05-14
| | | | | | "beyond Colella" limiting scheme to epsilon (to avoid that correction in one step leads to negative values). git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@333 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Support real*16 (and real*4) in GRHydroeschnett2012-05-11
| | | | | | | | | | | | Support using precisions other than real*8 in GRHydro by removing all explicit references to double precision functions and constants, and using type-generic functions and constants instead. In particular, use "one" and "half" as constants in some places, and use "abs", "max" etc. instead of "dabs", "dmax" etc. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@319 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: correct reconstruction of velocities in ePPMrhaas2012-05-11
| | | | | | | | | | | | This is a bug fix to the ePPM scheme. The velocity components were reconstructed from the plus face values! This is clearly incorrect. Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@318 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Cosmetic change to ePPM3.reisswig2012-02-16
| | | | | | | | Improved handling of cs2 < 0. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@314 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Set cs2 to zero if cs2 < 0 and abs(cs2) < 1e-4reisswig2012-02-03
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@312 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Implemented enhanced PPM scheme by Colella & Sekora 2008, McCorquodale & ↵reisswig2011-11-26
| | | | | | | | | | | | Colella 2011. Can be actived by setting use_enhanced_ppm = yes. Also fixed a bug in GRHydro_Reconstruct.F90 when using MP. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@302 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* remove unused macros SpaceMask_CheckStateBitsF90_1Dknarf2011-11-01
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@293 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * remove psi4 (optimize)cott2011-08-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@257 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* merge branch hot_and_MHD_temp_dev into branch at revision r185knarf2010-11-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@186 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * remove dependence on StaticConformalcott2010-08-27
| | | | | | | | | * change calculation of the determinant of the 3-metric from a subroutine call to a macro. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@152 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * remove unused EOS code from GRHydro_PPM.F90cott2010-06-30
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@129 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Fix error in ppm when using GRHydro_enable_internal_excision == no (which is ↵knarf2010-06-20
| | | | | | _not_ the default setting). Thanks go to Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@123 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* file/parameter string replacement from whisky to GRHydrobmundim2010-05-02
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@112 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45