aboutsummaryrefslogtreecommitdiff
path: root/src/GRHydro_Con2Prim.F90
Commit message (Collapse)AuthorAge
* GRHydro: use ashX instead of lsh(X) when declaring array sizesrhaas2013-10-29
| | | | | | | | 1.) one should use ash when declaring array sizes 2.) intel11 chokes on lsh(X), this was already in the first Cray pointer patch git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@585 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: replace Fortan pointers with Cray pointersrhaas2013-10-27
| | | | | | | Intel 12 seems to like these better. Otherwise GRHydro_Con2PrimM.F90 takes very long to compile. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@580 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: re-implement prim2con routine in Crhaas2013-07-06
| | | | | | | | | | | | | | | | * new experimental prim2con routine that is about twice as fast as the old one. * Right now tested only for simple EOS. * Right now handles only prim2con call after reconstruction and must be enabled by changing comments in GRHydro_Reconstruct.F90 * remove duplicate code in all reconstruction *_drv.F90 routines and just have the initialization of plus and minus variables in the main GRHydro_Reconstruct.F90 routine From: Christian Ott <cott@bethe.tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@556 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add grid function for sqrt(detg)rhaas2013-07-06
| | | | | | | | | | * add new 1 tl grid function sdetg that stores the sqrt of the determinent of the 3-metric. * replace lots of re-computation of det by use of this grid function From: Christian Ott <cott@bethe.tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@555 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: call dedicated HOT C2P routine in Con2Primrhaas2013-07-06
| | | | | | From: Christian David Ott <cott@zwicky-b.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@539 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Do not declare Fortran pointers explicitlyeschnett2013-06-13
| | | | | | | | | Some compilers do not accept Fortran POINTER declarations if the pointer variable has already been declared to be a particular integer type. Instead, declare them (only) via a POINTER statement. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@536 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add option to control action taken when con2prim failedrhaas2013-05-29
| | | | | | From: Roland Haas <rhaas@caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@529 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: relace C-style include for module with proper modulerhaas2013-05-17
| | | | | | | | the C-like coding had some problems, in particular it would compile the same module more than once and would circumvent Cactus' module dependency tracking. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@520 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: 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
* 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: 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: 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: 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
* Correct format specifiereschnett2012-09-12
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@421 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: 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: 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
* * 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
* GRHydro: remove "use" for con2prim interfaces from pointwise routinesrhaas2012-06-26
| | | | | | | | | since the Fortran standard does not allow this and gfortran 4.7 actually chokes on it. Pointed out by Erik Schnetter. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@374 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRydro: avoid implict SAVE attribute in Con2PrimPolytyperhaas2012-06-21
| | | | | | | | | initializing a variable makes it SAVE which is incomaptible with OpenMP. Also fix some 1-element-array vs. scalar variable type issues. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@364 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * introduce parameters GRHydro_Y_e_min and GRHydro_Y_e_max to fix issues ↵rhaas2012-06-21
| | | | | | with EOS tables that don't extend sufficiently far. Enforcement happens in Con2Prim. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@363 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: do not use integers as booleansrhaas2012-06-20
| | | | | | | gfortran (and the Fortran standard I guess) does not allow this. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@357 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* C2P_hot: Resort to bisection as well (if user requested so).rhaas2012-06-20
| | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@355 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: C2P: resort to bisection if Newton-Raphson fails. This is onlyrhaas2012-06-20
| | | | | | | | | activated if user requested this. Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@354 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * allow disabling check of the GRHydro_C2P_failed maskrhaas2012-06-05
| | | | | | | | | | | on coarse reflevels via GRHydro_c2p_warn_from_level parameter. Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@348 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: fix format statements in hot C2P.rhaas2012-06-05
| | | | | | | | | | Note: This should be propagated to ET! Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@347 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Added parameter to perform sync on conserved quantities only.rhaas2012-06-05
| | | | | | | | | | | | | | Primitives are only sync'ed after regridding to get initial guesses. Schedule Primitive boundary call in PostRegridInitial if sync_only_coserved == yes. Fixed issue with the sync call in CCTK_EVOL after MoL_Evolution when sync_conserved_only==yes: Conserved vars must be synced here as well, otherwise conserved vars are reset to atmospshere whereas primitives are not! This is wrong. Patches by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@340 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * improve error handling/reporting in hot C2Prhaas2012-05-13
| | | | | | | | | * set up directory for performance workshop Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@323 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Error check for C2P_failed mask greater or equal 1, do not abort code if in therhaas2012-05-13
| | | | | | | | | atmosphere though. patch by Philipp Moesta git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@322 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
* dump the atmosphere bitmask in favor of the atmosphere integer mask which is ↵knarf2011-11-01
| | | | | | used/allocated already anyway git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@292 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* move TARGET statements in front of CCTK_DECLARE_ARGUMENTSrhaas2011-09-25
| | | | | | | | | this is need for gfortran 4.1 which otherwise assumes that the arrays have already been used (in the kind(gaa) initializations inside of CCTK_DECLARE_ARGUMENTS). git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@278 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* remove explicit dependence on Multipatch thornsrhaas2011-09-15
| | | | | | | | | GRHydro now compiles and runs without MP thorns present and behaves identical to the old version git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@274 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* add Multipatch support to GRHydrorhaas2011-09-15
| | | | | | | | | | | * not all features of GRHydro are supported yet, in particular only the HLLE solver supports Mulitpatch yet. Original commit by Christian Reisswig and Christian Ott git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@273 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* clean up variable definitions in Eigenproblem_Marquinarhaas2011-09-15
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@272 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * improve behavior of CoreCollapseControl upon recoveryrhaas2011-09-15
| | | | | | | | * improve handling of too-low temperatures in GRHydro git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@270 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* include Y_e_con in atmosphere handlingrhaas2011-09-15
| | | | | | | | original commit by Christian Ott (cott) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@268 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* add more workarounds in Con2Prim for hot EOSrhaas2011-09-15
| | | | | | | | | | | * fix con2prim issues (pertaining to error messages) * fix OMP CRITICAL issues in Prim2Con and Con2Prim original commits by Christian Ott (cott) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@267 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Declare all OpenMP private variables explicitly as privateeschnett2011-06-13
| | | | | | | | | In Fortran OpenMP, loop variables of parallel loops are implicitly private and do not need to be declared as private. To improve code clarity, all such variables are now explicitly declared to be private. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@251 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * improve atmosphere handling so that it is now possible to usecott2011-06-11
| | | | | | | | hot EOS with an atmosphere git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@249 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * fix an issue with hot EOS and inner symmetry boundariescott2011-03-29
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@226 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * remove support for "General" EOS interface (EOSG_*)cott2011-03-28
| | | | | | | | | | * test suites pass with Intel 11 (no optimizaton) on bethe and gcc 4.4.5 git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@225 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * fix OMP bug in _hot part of c2pcott2011-03-12
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@221 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* patch from Roland Haas: in Con2PrimGeneral inside a doknarf2011-02-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@219 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-27 17:58:07 +0000