aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
* GRHydro: increase test tolerance for balsara4_1d testrhaas2013-11-22
| | | | | | | | this is a very highly relativistic shock (v=0.999) and apparently different compilers produce very different results. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@589 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove misleading commentrhaas2013-11-18
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@588 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: work around Intel 11 failures when using Fortran pointersrhaas2013-11-15
| | | | | | | | | | | | this rewrites the Reconstruction driver routines since Intel 11 miscompiles the code using Fortran pointers and generates segfault errors at runtime. Nested subroutines (another solution to make intel 11 work) do not work with intel 13. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@587 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Declare loop variables OpenMP privateeschnett2013-10-30
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@586 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* 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
* Add parentheses around macro definitionseschnett2013-10-29
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@584 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: correct usage of sqrt(detg) in prim2con wrappersrhaas2013-10-29
| | | | | | found by Ian Hawke git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@583 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Make some macros safereschnett2013-10-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@582 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Correct array indices in etaXYZeschnett2013-10-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@581 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
* Adjust loops for Prim2Con call (needs to run in the ghost zones as well ↵HEADmasterrhaas2013-08-13
| | | | | | | | after being called from ZelmaniLeak). From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@579 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Get rid of old heat injection in C2P MHDrhaas2013-08-13
| | | | | | From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@578 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Notes on how to work out the stsggered A-field evolutionrhaas2013-08-13
| | | | | | From: Josh Faber <jfaberrit@jfaber3.local> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@577 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: use C prim2con routinerhaas2013-08-13
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@576 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: finish up call for templated C++ PPM routinerhaas2013-08-13
| | | | | | | | activated by paramater use_optimized_PPM From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@575 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: First draft of templated version of PPM Reconstruction call,rhaas2013-08-13
| | | | | | | | | Mostly involves slicing and unslicing 3-d grids into 1-d segments and calling the new templated routines in GRHydro_PPM.cc From: Josh Faber <jfaberrit@jfaber3.local> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@574 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add templated (yay\!) PPM routine; can do old ppm, but for MHD, ↵rhaas2013-08-13
| | | | | | | | hot, cold whatever From: Christian Ott <cott@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@573 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: in c PPM routine, remove mistakerhaas2013-08-13
| | | | | | From: Christian D. Ott <cott@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@572 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Move Evec calculation in seperate loopsrhaas2013-08-13
| | | | | | | | | this avoids accessing unitinialized values in the ghost zones when computing the hydro varaible updates. This also saves time in running those loops since they are smaller. From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@571 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Adding parfile for vector potential TOV.rhaas2013-08-13
| | | | | | From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@570 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* add interface use OMP version of ESO_Omni_press functionrhaas2013-08-13
| | | | | | From: Christian Ott <cott@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@569 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: avoid allocating large arrays on the stackrhaas2013-08-13
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@568 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* expand warnings in GRHydro_Con2PrimHotrhaas2013-08-13
| | | | | | | | | * compute reflevel in GRHydro_Con2PrimHot to be sure it is current * only output warnings if we want them From: Christian Ott <cott@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@567 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* don't specify boundaries twicerhaas2013-08-13
| | | | | | From: Frank Loeffler <knarf@cct.lsu.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@566 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: explicitly recompute SqrtDetg after RECOVERYrhaas2013-08-13
| | | | | | | | | this is usually not required, unless MoL_run_MoL_PostStep_in_PostRecoverVariables is set to false. From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@565 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: fix metric GF used for sqrtdetg in multipatch runsrhaas2013-08-13
| | | | | | | | fix some scheduling From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@564 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Implemented WENO-z. Ability to choose between standard WENO5 and ↵rhaas2013-08-13
| | | | | | | | "adaptive epsilon"-variant from WHAM. From: Christian Reisswig <reisswig@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@563 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: always set up SqrtDetg GF in INITIALrhaas2013-07-15
| | | | | | | | | rather than only when evolution_method == "GRHydro" so that ID only runs (or static) runs can use InitialAtmosphereReset From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@562 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * update par files using TOVSolver to not use [0] after TOV_K and cott2013-07-07
| | | | | | | | TOV_Gamma git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@561 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: reimplement PPM in Crhaas2013-07-06
| | | | | | | | | | | | | | | | * add prototype of a re-implementation of PPM in C. So far we can do old PPM without and with temperature/Y_e reconstruction. * also add a version of the reconstruction driver that is optimized for the use of the new PPM routine. * no change to standard code version; one needs to manually replace source files to get the new stuff. * new version is factor 3 faster than old version (based on preliminary measurements). From: Christian Ott <cott@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@560 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: make atmosphere threshold steerable on recoveryrhaas2013-07-06
| | | | | | From: Christian David Ott <cott@zwicky.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@559 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: activate storage for Y_e, temperature, entropy when we evolve itrhaas2013-07-06
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@558 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: provide wrapper routines that take detg as argumentsrhaas2013-07-06
| | | | | | | | | | then call internal routines taking sqrt(detg) as argument clean up and correct function prototypes in interface.ccl From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@557 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: pass correct temperature o prim2con call on coarse levelsrhaas2013-07-06
| | | | | | | | | Fix a bug in prim2con_hot that was causing temp to go negative near the shock in 3D CCSN simulations with WENO reconstruction. From: abdik <abdik@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@554 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* 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: fix compiler warningsrhaas2013-07-06
| | | | | | | | | | | | | * remove unused variable * pacify compiler warning about possibly uninitialized value not sure if one could somehow get the compiler to realize that the condition on which access is based is false in both case, ie that evolve_MHD is an effective constant. From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@552 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Con2PrimM: Set points to atmosphere if in excised region.rhaas2013-07-06
| | | | | | From: Christian Reisswig <reisswig@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@551 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: whitespace change, tab->8 spacesrhaas2013-07-06
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@550 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: small cosmetic change to Con2PrimHotrhaas2013-07-06
| | | | | | From: Christian David Ott <cott@zwicky.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@549 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: Don't compute Bcom twice.rhaas2013-07-06
| | | | | | From: Christian Reisswig <reisswig@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@547 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Make "CalcBcom" compatible with multipatch by enforcing scheduling ↵rhaas2013-07-06
| | | | | | | | after local->global trafo of primitives. From: Christian Reisswig <reisswig@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@546 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: Make MHD version of Prim2Con compatible with multipatch.rhaas2013-07-06
| | | | | | From: Christian Reisswig <reisswig@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@544 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Fixed a problem with excision mask and atmopshere treatment in hot ↵rhaas2013-07-06
| | | | | | | | Con2Prim. From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@543 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Critical bugfix: Con2Prim was probably never executed in ↵rhaas2013-07-06
| | | | | | | | | | | | 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@542 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: limit temperature to valid range in prim2conrhaas2013-07-06
| | | | | | | | add some more error checking From: Christian David Ott <cott@zwicky-b.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@541 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: change C2P interfacerhaas2013-07-06
| | | | | | From: Christian David Ott <cott@zwicky-b.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@540 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45