aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
...
* 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: register TmunuBase variables as MoL constraint.rhaas2012-06-05
| | | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@350 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * allow disabling check of the GRHydro_C2P_failed maskrhaas2012-06-05
| | | | | | | | | | | on coarse reflevels via GRHydro_c2p_warn_from_level parameter. Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@348 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: fix format statements in hot C2P.rhaas2012-06-05
| | | | | | | | | | Note: This should be propagated to ET! Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@347 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Added option to constrain evolution to 1D byrhaas2012-06-05
| | | | | | | | | | projecting out non-radial parts of the conserved velocity RHSs. Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@346 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: do not register any SandR variables if MaxNumSandR ==0rhaas2012-06-05
| | | | | | | | | this is useful when running with a spacetime code that will register the ADMBase variables as constrained anyway, so one can save the memory used for SandR scratch space. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@345 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Verbosity on which boundary selction routine is executed.rhaas2012-06-05
| | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@344 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: Added verbosity parameter.rhaas2012-06-05
| | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@342 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: sync those primitives which are needed for initial Con2Prim guessesrhaas2012-06-05
| | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@341 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Added parameter to perform sync on conserved quantities only.rhaas2012-06-05
| | | | | | | | | | | | | | Primitives are only sync'ed after regridding to get initial guesses. Schedule Primitive boundary call in PostRegridInitial if sync_only_coserved == yes. Fixed issue with the sync call in CCTK_EVOL after MoL_Evolution when sync_conserved_only==yes: Conserved vars must be synced here as well, otherwise conserved vars are reset to atmospshere whereas primitives are not! This is wrong. Patches by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@340 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: a bit of a clean-up of unused variables to avoid compiler complaintsbmundim2012-05-17
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@336 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: decouple Bvec(flux_direction) from other two for divcleanrhaas2012-05-16
| | | | | | | | since it and psidc form a 2d subsystem of the characteristic system. See eg. Mignonoe et al. arXiv:1001.2832v2 section 3.1. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@335 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Fix missing OpenMP variable private declarations.rhaas2012-05-14
| | | | | | | patch by Philipp Moesta. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@334 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Cleanup of PPM routine. Remove unused PPM parameter. Only apply stricter ↵reisswig2012-05-14
| | | | | | "beyond Colella" limiting scheme to epsilon (to avoid that correction in one step leads to negative values). git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@333 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: correct atmosphere handling after primitve reconstructionrhaas2012-05-14
| | | | | | | | | | | | | | | | | | | | | I have fixed a bug in GRHydro's atmosphere handling after reconstruction. After reconstruction, the reconstructed plus and minus face values are tested for whether they drop below atmosphere level and are then reset. In particular, this means that plus and minus face values for cell i don't generally coincide with the minus/plus face values of the corresponding cells i-1 and i+1. If the reconstruction previously found that the Riemann problem was trivial between, say rhominus(i+1) and rhoplus(i), when changing rhoplus(i), it is not anymore! However, the mask "trivial_rp" testing for a trivial Riemann problem is not changed! This leads to inconsistent behavior at the boundary of the atmosphere. The symptom was a drift in the center of mass of a static and perfectly symmetric TOV star. Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@332 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add file missing from r330rhaas2012-05-14
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@331 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
* GRHydro: Correct syntax errorrhaas2012-05-13
| | | | | | | | | Detected by gcc 4.7. Patch by Erik Schnetter git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@328 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: OpenMP parallelize constraint transport looprhaas2012-05-13
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@327 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: fix OpenMP data race condition in UpdateCalculationrhaas2012-05-13
| | | | | | | for constraint transport git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@326 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Fixed minor bugs in registration in Boundariesrhaas2012-05-13
| | | | | | | | | | Some variables were stil set to "Flat" boundaries even though "None" is supposed to be set. Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@325 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Adding Josh's latest changes to the general EOS Con2PrimM routine.rhaas2012-05-13
| | | | | | | patch by Philipp Moesta git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@324 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * improve error handling/reporting in hot C2Prhaas2012-05-13
| | | | | | | | | * set up directory for performance workshop Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@323 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Error check for C2P_failed mask greater or equal 1, do not abort code if in therhaas2012-05-13
| | | | | | | | | atmosphere though. patch by Philipp Moesta git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@322 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Consistent atmosphere treatment.rhaas2012-05-13
| | | | | | | patch by Philipp Moesta git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@321 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
* Support real*16 (and real*4) in GRHydroeschnett2012-05-11
| | | | | | | | | | | | Support using precisions other than real*8 in GRHydro by removing all explicit references to double precision functions and constants, and using type-generic functions and constants instead. In particular, use "one" and "half" as constants in some places, and use "abs", "max" etc. instead of "dabs", "dmax" etc. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@319 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: correct reconstruction of velocities in ePPMrhaas2012-05-11
| | | | | | | | | | | | This is a bug fix to the ePPM scheme. The velocity components were reconstructed from the plus face values! This is clearly incorrect. Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@318 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Cosmetic change to ePPM3.reisswig2012-02-16
| | | | | | | | Improved handling of cs2 < 0. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@314 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* make EOS change work with MPcott2012-02-09
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@313 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Set cs2 to zero if cs2 < 0 and abs(cs2) < 1e-4reisswig2012-02-03
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@312 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Fixed a problem with gxx vs g11.reisswig2012-01-25
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@311 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
* * introduce check that cs2 > 0.cott2011-12-07
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@309 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* expand region where fluxes are computed for constraint transportrhaas2011-11-29
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@307 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* remove unused variable from GRHydro_Boundariesrhaas2011-11-29
| | | | | | | | remove unused file GRHydro_BoundariesM.F90 remove unused file GRHydro_ReconstructPolyM.F90 remove unused file GRHydro_ReconstructM.F90 git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@306 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* fix error checking in GRHydro_Boundariesrhaas2011-11-29
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@305 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* remove extraneous (and faulty) Bconrhs access introduced be me in r298rhaas2011-11-28
| | | | | | | | could actually cause a segfault since it uses unitiailised values in i,j,k. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@304 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* move TARGET statement in Reconstruction to compile with gcc 4.1rhaas2011-11-28
| | | | | | | gcc 4.1 requires TARGET statements to be before DECLARE_CCTK_ARGUMENTS git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@303 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Implemented enhanced PPM scheme by Colella & Sekora 2008, McCorquodale & ↵reisswig2011-11-26
| | | | | | | | | | | | Colella 2011. Can be actived by setting use_enhanced_ppm = yes. Also fixed a bug in GRHydro_Reconstruct.F90 when using MP. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@302 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* allocate/deallocate ENO scalars in global moderhaas2011-11-21
| | | | | | | | | | really all I need is that this happens only once. Allocation was already protected by a grid scalar, unfortunately deallocation did not check/reset this scalar git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@301 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* protect Bvec in ENO polytype reconstruction by evolve_mhdrhaas2011-11-21
| | | | | | | otherwise unallocated memory is accessed git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@300 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* second iteration of constraint transportrhaas2011-11-04
| | | | | | | | | | * fix some indices * move poison loop to proper spot * compute conserved divergence in track_divB * add Maple worksheet to check constraint transport indices git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@298 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
* First try of a constraint transport scheme for unigridrhaas2011-11-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@296 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* remove (now) unused old Tmunu interface)knarf2011-11-01
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@294 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* remove unused macros SpaceMask_CheckStateBitsF90_1Dknarf2011-11-01
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@293 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
* RIT dev: Minor change in sign convention of an equation.bmundim2011-11-01
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@289 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-12 12:13:55 +0000