aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
...
* GRHydro: Fixes to atmosphere mask and the way we sync.rhaas2012-07-05
| | | | | | | | | | | | | These are the changes as discussed with Roland. In particular, we sync the atmo mask in an extra call before all other syncs and then do atmo reset. Furthermore, primitives must always be synced in last post step in order to get initial guesses in the buffer zones via prolongation. Atmo reset must be done before prolongation to ensure we don't overwrite "good" data from the coarse grid! git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@377 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * add limit on values of Y_e also in Prim2Conrhaas2012-07-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@375 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: remove "use" for con2prim interfaces from pointwise routinesrhaas2012-06-26
| | | | | | | | | since the Fortran standard does not allow this and gfortran 4.7 actually chokes on it. Pointed out by Erik Schnetter. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@374 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Revert "* improve handling of cs2 < 0 in the "standard" (non-MHD, non-hot ↵rhaas2012-06-22
| | | | | | | | | | EOS) eigenvalue routine." actually this introduce the real valued atmosphere mask. This reverts commit deccf528b463820b615dfa73e46e137cd59c65be. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@372 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Revert "Revert "GRHydro: do not access whichpsidcspeed in ParamCheck unless ↵rhaas2012-06-22
| | | | | | | | | | evolve_MHD"" Ahem. I had reverted the wrong commit. This one was blameless. This reverts commit 0b5a561872841b50ea32da2b665a04b2984a8c16. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@371 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Revert "GRHydro: do not access whichpsidcspeed in ParamCheck unless evolve_MHD"rhaas2012-06-22
| | | | | | This reverts commit c2122aa2b4762388eb90cc159750aa4d418b112d. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@370 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRydro: avoid implict SAVE attribute in Con2PrimPolytyperhaas2012-06-21
| | | | | | | | | initializing a variable makes it SAVE which is incomaptible with OpenMP. Also fix some 1-element-array vs. scalar variable type issues. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@364 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * introduce parameters GRHydro_Y_e_min and GRHydro_Y_e_max to fix issues ↵rhaas2012-06-21
| | | | | | with EOS tables that don't extend sufficiently far. Enforcement happens in Con2Prim. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@363 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* really improve handling of cs2 < 0 in the "standard" (non-MHD, non-hot EOS)rhaas2012-06-20
| | | | | | | | | | | | | | | | | | | | | | | eigenvalue routine. the previous patch with this same commit message (r359) has the wrong commit message. It should instead read: --8<-- Introduce real-valued atmo mask that can be interpolated. Synchronize real-valued atmo mask and update integer-valued atmo mask based on synchronized value. This fix ensures that the results are independent of the number of processors used. It also gets rid of some extra primitive syncs that are otherwise necessary when sync_conserved_only = yes. Does not cause any new test failures (Roland). Patch by Christian Reisswig. --8<-- Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@361 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: do not access whichpsidcspeed in ParamCheck unless evolve_MHDrhaas2012-06-20
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@360 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * improve handling of cs2 < 0 in the "standard" (non-MHD, non-hot EOS) ↵rhaas2012-06-20
| | | | | | | | | eigenvalue routine. Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@359 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Change the characteristic speed for psidc.rhaas2012-06-20
| | | | | | | | | | | | | | It was set to +-c when it should be the maximum coming from the constraint that the normal vector to the characteristic hypersurface of a wave propagating along an arbitrary x direction be spacelike. This results into the following bound for the wavespeed: +- \alpha sqrt(\gamma^{xx}) - \beta^x) Patch by Bruno Mundim. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@358 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: do not use integers as booleansrhaas2012-06-20
| | | | | | | gfortran (and the Fortran standard I guess) does not allow this. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@357 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* C2P_hot: Resort to bisection as well (if user requested so).rhaas2012-06-20
| | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@355 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: C2P: resort to bisection if Newton-Raphson fails. This is onlyrhaas2012-06-20
| | | | | | | | | activated if user requested this. Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@354 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: add file missing from r343rhaas2012-06-06
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@352 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: 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
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-25 20:06:04 +0000