aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
* 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
* 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: MP5: Make epsilon parameter scale dependent on reconstructed field.rhaas2013-05-29
| | | | | | From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@528 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add error checking after call for cs2 in Eigenproblemrhaas2013-05-29
| | | | | | From: Christian David Ott <cott@zwicky-b.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@527 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Fixing index offset bug in tempplus(i,j,k) for x direction flux.rhaas2013-05-29
| | | | | | From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@526 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Roland Haas: replace the diagonal shocktube in the WENO test by a shock ↵rhaas2013-05-22
| | | | | | | | | | along the y direction. The diagonal shock causes issues near the boundary where the shock fails due to differences between gfortran and intel compilers. One issue is obviously that the boundary condition is incorrect and one really would want to use a boundary condition similar to the ones implemented for the MHD shocks by Josh and Bruno. For now I simplified the test which also makes to work more reliably on the different compilers. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@524 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add missing filesrhaas2013-05-17
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@522 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: remove duplicate (and incorrect) GRHydro_InterfacesAM.hrhaas2013-05-13
| | | | | | | sync up GRHydro_Con2PrimAM with GRHydro_Con2PrimM git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@517 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: do not declare grid scalar as OpenMP PRIVATErhaas2013-05-13
| | | | | | | instead leave it SHARED since it is a global read-only variable. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@516 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: fix use of uninitialized variable in eigenproblem_leftrightrhaas2013-05-13
| | | | | | | this is only used when doing flux split FD which is not the default. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@515 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Fix a problem with ghost zones in Prim2ConCells calls in MHD case.rhaas2013-05-11
| | | | | | From: Christian Reisswig <reisswig@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@512 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Fixing Con2PrimM to call the proper point-wise routinesrhaas2013-05-11
| | | | | | | | | | This re-introduces routines that work for hybrid/hot EOS and corresponding changes in pointwise routines for hot EOS error checking and temperature treatment by adding old EOSOmni pointwise routine. From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@511 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* 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
* Although probably often not wanted, there is no reason to have ↵knarf2013-03-28
| | | | | | riemann_solver not steerable=always. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@503 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Mark GRHydro_NaN_verbose to be steerable. There is no reason this shouldn't ↵knarf2013-03-28
| | | | | | be and odds are that you would like to enable this at some point for debugging... git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@502 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
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-05 07:48:25 +0000