| Commit message (Collapse) | Author | Age |
|
|
|
|
|
| |
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@492 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
| |
instability. Tested with shocktube and TOV. Not yet tested for problem where carbuncles occur.
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@491 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@487 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
| |
This does not yet work with the MHD part of the code,
but Philipp is going to change this.
* Improve TOVSolverHot; in particular, give functions unique
names to avoid duplicate symbols with TOVSolverC
From: Christian Ott <cott@bethe.tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@473 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
From: Roland Haas <rhaas@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@469 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
| |
according to WHAM code paper, (ii) Catch epsilon<0 for non-microphysical EOS. Tested with shock tube and TOV.
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@466 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@465 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
radius. Off by default.
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@453 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
From: Roland Haas <rhaas@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@452 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
| |
From: Roland Haas <roland.haas@physics.gatech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@442 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
rescaling in con2prim failures.
From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@437 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
From: Roland Haas <roland.haas@physics.gatech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@431 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@425 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
the primitive velocity 'vel'.
This ensures that the Lorentz factor can never become unphysical during reconstruction.
This patch is necessary for NS collapse when ePPM+Refluxing is used.
The idea is due to Roland.
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@416 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
| |
NOTE: This patch requires support from MoL.
Patch by Christian Reisswig.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@405 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@398 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
From: Christian Ott <cott@bethe.tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@394 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@376 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
to include Tmunu in the count. Extra constraints are cheap since MoL
does not allocate any scratch or save storage for them.
Bruno Mundim pointed this one out and provided the accurate count.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@373 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
| |
Patch by Christian Ott.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@356 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
Patch by Christian Reisswig.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@349 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
| |
Patch by Christian Reisswig.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@342 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
| |
"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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@315 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@296 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@291 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@281 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
| |
* 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
|
|
|
|
|
|
|
|
| |
* improve handling of too-low temperatures in GRHydro
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@270 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
additional sync/settmunu after MoL is done. This will
work only in runs that do not use an atmosphere.
original commit by Christian Ott
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@269 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@266 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
Change legacy divergence cleaning parameters ch_dc and cp_dc,
introduced for the special relativistic case, to kap_dc, the
one introduced in the covariant formulation (Liebling et al.).
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@263 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@250 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
1) Fix a bug in the divergence cleaning implementation:
a psidc term in the induction equation was implemented
as a source term where it was supposed to be coded as
part of the flux calculation.
2) Introduce divB as a diagnostic grid function.
3) Remove an old file: GRHydro_CalcUpdateM.F90
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@228 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
* 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
|