aboutsummaryrefslogtreecommitdiff
path: root/src/make.code.defn
Commit message (Collapse)AuthorAge
* 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: 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: 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: 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: 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: 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
* 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: add test case for WENO5rhaas2013-03-06
| | | | | | | | this is just a copy of the shock tube test using WENO From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@489 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Added MP5 reconstruction method. Tested with TOV and shocktube.rhaas2013-03-06
| | | | | | From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@487 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Implemented WENO5. Tested with a shock tube.rhaas2013-01-14
| | | | | | From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@465 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Add basic vector potential supportrhaas2013-01-14
| | | | | | | | | | Basic cell-centered, algebraic gauge vector potential method with place-holders for lorenz gauge. Initial Avec constrained to poloidal at the moment. From: Tanja Bode <tanja.bode@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@456 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: Unify function interface for GRHydro_Con2PrimM_pt.c and ↵rhaas2012-10-25
| | | | | | | | | | | | | | | GRHydro_Con2PrimM_pt_EOSOmni.c * Interestingly Cactus doesn't complain about the function interface in interface.ccl. It has the old set of arguments for that function and declare it as written in Fortran instead of C. * Cylindrical explosion works again with GRHydro_Con2PrimM_pt.c routine. You can just swap both routines in make.code.defn From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@427 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Changes to make GRHydro MHD work with nuclear/hot equation of state:rhaas2012-08-27
| | | | | | | | | | | | | | | | | | 1) Switch to EOSOmni pointwise C2P routine and modify where necessary. 2) Modify Con2PrimM.F90 to allow for the evolution of temperature and adjust the wrapper routine. 3) Create EigenProblemM_hot pointwise routine and call that from HLLEM.F90 when temperature is evolved. Additionally adjust HLLEM where necessary. 4) Adjust InterfacesM.h to incorporate the newly created functions. 5) Fix a loop problem (not taking into account constraint transport) in PPM reconstruction of Y_e 6) Introduce Prim2ConM_hot and call this pointwise routine from Prim2ConM.F90 when temperature is evolved. Additionally also make this routine available to initial data routine in GRHydro_InitData 7) Adjust loops in GRHydro_PoloidalMagFieldM.F90 to not set boundary points it cannot set but instead call boundary group afterwards! Pay attention as this will not work with boundary conditions set to "none" in MHD case anymore but is the correct thing to do. 8) Allow StarMapper to extend HydroBase::initial_hydro = "starmapper". 9) Smaller fixes. From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@410 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Fix slope limiters and get rid of broken cases;rhaas2012-06-05
| | | | | | | | | | | there is a related patch in TRAC at the moment Patch by Josh Faber. Some additional checks by Roland Haas. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@351 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Schedule AtmosphereReset in LAST Mol_PostSteprhaas2012-06-05
| | | | | | | | | | | | This saves some syncs. This also gets rid of the akward parameter "disable_sync_and_tmunu_after_mol". Original patch by Christian Reisswig. Adapted to not rely on MoL support by Roland Haas. Bug fixed by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@343 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: fixes for GRMHDrhaas2012-05-14
| | | | | | | | | | | | | | | | | | | | All by Philipp Moesta. 1) Fix parity of psidc and divb 2) Fix a wrong index in the source terms of scon 3) Fix wrong indices of derivatives of space-time metric in the source of the divergence cleaning scalar. 4) Calculate divergence of B in MoL PseudoEvolution and set its Prolongation="Restrict". 5) Correct the source terms and fluxes for the Bfield and the divergence cleaning field when having a non-flat space-time. 6) Make sure alpha factors match between UpdateCalculation and fluxes definition. 7) Include 1/sqrt(detg) factor in calculation of \epsilon^{\muijk} in the cross product to obtain the Bfield form the vector potential. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@330 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Add a routine to calculate mag pressure/comoving fieldbmundim2012-01-03
| | | | | | | | | | * When Tmunu is not calculated, ie when no storage is selected, the comoving magnetic field/magnetic pressure wasn't calculated at all. This commit allows that now. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@310 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
* RIT GRMHD dev:bmundim2011-08-02
| | | | | | | | | | | | | | | | | | | | | | | | | | | 1) Further splitting of PPM Reconstruction routine into magnetic and non-magnetic part. 2) Merge the divergence cleaning loop into the non-divergence one, since profiling indicated no substantial difference between the two cases. Indeed it was a little bit better after the merger (0.14%), but that it is not significant. We decided then to apply Occam's razor and choose the simplest form. Branch prediction seems to work fine in this case. 3) Move reconstruction initialization statements to their repective drivers. 4) Get rid of WORKSHARE in the reconstruction routines. Profiling showed a 4.16% performance improvement for the hydro ppm reconstruction routine when using 1 processor and 2 threads only. Expect a better improvement for a larger number of threads. 5) Introduce a parameter to control characteristic speeds for psidc in HLLE. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@255 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* RIT GRMHD dev: split reconstruct routine.bmundim2011-07-20
| | | | | | | | | | | | | | | | | | | | | Comments: 1) The reconstruct routine is the most expensive routine in GRHydro and may require some serious thought in the future regarding optimization. 2) This routine was merged by the end of the last year. With too many loops, the compiler wasn't able to fully optimize it (and it would warn about that). 3) With this split, the code became easier to read, the compiler was able to fully optimize it, but we only have a marginal gain in the routine performance (3.2%). Further optimization and clean-up is still desirable. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@252 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
* * remove GZ patch system codecott2011-03-12
| | | | | | | | * test suites pass git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@223 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* RIT MHD development:bmundim2010-12-31
| | | | | | | | | | | | | | | | | | Revert previous commit changes to schedule.ccl Merge MHD and GRHydro routines: Boundaries, CalcUpdate, PPM, Reconstruct(poly), RegisterGZ, RegisterVars Solve a few race conditions arising in the "magnetic" routines. This is still being tested, since it passes tests with intel compilers but fails for gnu ones. Alias function GRHydro_Con2PrimM_pt for C2P2CM test. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@202 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* RIT MHD dev:bmundim2010-12-25
| | | | | | | | | | | | | | | | Fix to GRHydro_RegisterVarsM.cc and allows MHD to work without Y_E bing evolved. Add C2P polytype interface. Change schedule.ccl to represent a more detailed version of choosing what gridfucntions get sync'd in Boundaries depending on what is being evolved. Un-hardwire fixed values of gamma for all MHD Con2Prim routines. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@199 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* RIT MHD development: an update.bmundim2010-12-22
| | | | | | | | | | | | | | Several bugfixes, including properly isolating divergence cleaning code and initialization of the "psidc" gridfunction in GRHydro_DivergenceClean.F90, as well as properly initializing rhoenth in GRHydro_HLLEM (thanks to Roland for catching that). Minor changes to scheduler to properly call MHD and non-MHD routines in parallel for Riemann solving. Remove GRHydro_ParamCheckM.F90 and add GRHydro_DivergenceClean.F90. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@196 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * rename HLLE routine source code files to reflectcott2010-12-06
| | | | | | | | true content git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@189 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Current RIT GRMHD code contributions:bmundim2010-09-29
| | | | | | | | | | Add the magnetized counterparts for several GRHydro routines. Adjust interface.ccl, param.ccl and schedule.ccl appropriately. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@158 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * initial structural modifications to start with MHD development (Noble/Ott)cott2010-05-26
| | | | | | | | | | | | | | | | | | | -> add comoving b-field: bcom[3] and bcom0 -> add integer variable MHD which is set in ParamCheck to 0 or 1 according to HydroBase::Bvec_evolution_method being GRHydro or something else (don't want to do string comparisons very often) -> skeleton routine Conservative2PrimitiveM in separate .F90 file -> extend keyword Bvec_evolution_method with "GRHydro" * code compiles, should not hurt anything already there, does not do anyting (yet) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@115 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
* remove analysis routines (moved to Hydro_Analysis), removed SetDriftCorrect ↵knarf2010-01-12
| | | | | | funtions because apparently unused - could also be implemented in extra thorn if requested git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@34 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* This is a _temporary_ repository to be able to start to work on theknarf2009-11-18
code right now. I have put in the public version of Whisky to start from. Everybody with commit rights should get commit messages (and the other way around). It should not be a problem to add people to that list, just ask. I don't want to get into political problems because someone feels excluded, but I also don't want to give everyone access per se. Frank git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@3 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45