aboutsummaryrefslogtreecommitdiff
path: root/src/GRHydro_Con2PrimM.F90
Commit message (Collapse)AuthorAge
* GRHydro: disable warnings in Prim2ConM and only set mask in Con2PrimMrhaas2014-04-15
| | | | | | | We don't abort in ptwise C2P call with MHD so that we can check for the mask and abort if necessary. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@638 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Better handling of excision for MHD.rhaas2013-12-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@593 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
* 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: 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: 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: 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: remove commented out code, update docsrhaas2013-02-11
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@481 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: correct tau when resetting to atmosphererhaas2013-01-14
| | | | | | | | before a term rho*b2 was present when the correct one is only b2 From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@462 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: 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: disable in-AH polytrope call by defaultrhaas2013-01-14
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@452 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: call polytrope EOS in a region dominated by mag. field.rhaas2013-01-14
| | | | | | From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@451 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Call pointwise polytype upon con2primM failure. If it fails just ↵rhaas2013-01-14
| | | | | | | | use previous primitive values and get the cons from them. Basically freezing their values at previous time step. From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@450 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: call pointwise con2primM with polytrope eos handle deep inside AH ↵rhaas2013-01-14
| | | | | | | | upon inversion failure of gamma law eos. From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@449 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Introduce inequalities cons variables needs to satisfy when B/=0.rhaas2012-11-09
| | | | | | | | | | Introduce temporary prim variables as well. * Please refer to Illinois' paper (arXiv:1112.0568) appendices for details. From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@439 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Check if 3-metric is definite positive before con2prim inversion.rhaas2012-11-09
| | | | | | | | | Only do this deep inside the horizon. Reset it to conformally flat if it isn't. From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@438 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: MP pointer stuff for MHDrhaas2012-08-27
| | | | | | | | this time with minimal replacements in SourceM.F90 From: Roland Haas <roland.haas@physics.gatech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@411 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: remove unused variablesrhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@386 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Error checking for C2P_failed > 1, not only equals 1rhaas2012-05-13
| | | | | | | | and check for atmosphere using density after Con2PrimM_pt was performed. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@320 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Fix Con2Prim EOS calls in MHD version.tbode2011-11-04
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@297 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
* 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
* MERGE divergence cleaning feature into trunk.tbode2011-04-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@244 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
* RIT MHD development:bmundim2011-01-06
| | | | | | | | | | Fix race conditions arising in the Con2Prim routines. Now both GNU and Intel compilers agree on the results for the magnetic tests. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@208 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
* * make EOS Omni calls thread safecott2010-12-21
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@194 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* merge branch hot_and_MHD_temp_dev into branch at revision r185knarf2010-11-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@186 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
* * tinycott2010-05-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@116 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
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-14 14:24:39 +0000