aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
...
* 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
* * make parameter c2p_resort_to_bisection steerablerhaas2012-06-20
| | | | | | | Patch by Christian Ott. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@356 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 files missing from r340rhaas2012-06-06
| | | | | | | Thanks to Ian Hinder for pointing these (and the last one) out. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@353 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
* GRHydro: Make enahnced params steerable.rhaas2012-06-05
| | | | | | | | Patch by Christian Reisswig. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@349 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
* Correct S_i conserved variable index.bmundim2012-05-18
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@338 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Change default value of kap_dc from 1.0 to 10.0 as well as the recommendations.bmundim2012-05-17
| | | | | | | | | | | | | This is based on the magnetic monopole tests. It indicates that kap_dc = 10.0 reaches the best balance between damping and advecting. kap_dc = 1.0 advects well but damps much slower. kap_dc = 100.0 seems to turn the problem into a stiff one and it fails to either damp or advect the divergencies away when compared to the other two cases. The recommendation then goes to use values between 1.0 and 10.0. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@337 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 vector type of Bcons to match Bvecrhaas2012-05-13
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@329 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
* remove tests using BSSN_MoL for which ML_BSSN alternatives existrhaas2012-05-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@317 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: switch testsuite files to fdOrder=2 from ML_BSSN_O2rhaas2012-05-01
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@316 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* epsilon: specific internal energy (ie energy/mass)bmundim2012-03-13
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@315 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