| Commit message (Collapse) | Author | Age |
... | |
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@360 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
| |
eigenvalue routine.
Patch by Christian Ott.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@359 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
| |
Patch by Christian Ott.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@356 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
| |
Patch by Christian Reisswig.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@355 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
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@352 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@350 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
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
| |
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 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
|
|
|
|
|
|
|
| |
Patch by Christian Reisswig.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@344 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
|
|
|
|
|
|
|
| |
Patch by Christian Reisswig.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@341 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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@338 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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@336 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
|
|
|
|
|
|
|
| |
patch by Philipp Moesta.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@334 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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@331 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@329 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@327 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
| |
for constraint transport
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@326 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
| |
patch by Philipp Moesta
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@324 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
| |
* 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
|
|
|
|
|
|
|
|
|
| |
atmosphere though.
patch by Philipp Moesta
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@322 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
| |
patch by Philipp Moesta
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@321 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
| |
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 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
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@317 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@316 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@315 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
| |
Improved handling of cs2 < 0.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@314 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@313 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@312 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@311 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
* 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
|