aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
* 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
* 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
* GRHydro: add GRHydro_max_temp parameter for max allowed Trhaas2013-07-06
| | | | | | From: Christian David Ott <cott@zwicky-b.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@538 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add EOS_Omni multi-var call for dpde and dprhorhaas2013-07-06
| | | | | | From: Christian David Ott <cott@zwicky-b.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@537 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 check for stencil width for MP5rhaas2013-05-29
| | | | | | | | regenerate test that was using wrong stencil before From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@535 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: avoid implict save for ierr in boundary registrationrhaas2013-05-29
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@534 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add PARAMCHECK for grhydro_stencil for WENO5rhaas2013-05-29
| | | | | | | | WENO5 needs at least a stencil width of 3 From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@533 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: initialy set y_e_plus/minus to cell center value in reconstructrhaas2013-05-29
| | | | | | | | This ensures that all EOS calls have some sensible y_e to work with. From: Christian Ott <cott@bethe.tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@531 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: correct parameter description textrhaas2013-05-29
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@530 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45