| Commit message (Collapse) | Author | Age |
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@420 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
|
|
|
|
|
|
|
|
| |
which used to be called PrimitiveBoundaries
From: Roland Haas <roland.haas@physics.gatech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@400 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
| |
MoL_PostStep rather than HydroBase_PostStep, since the latter is scheduled IN MoL_PostStep
-> rely on correct scheduling rather than on luck to get the correct
atmosphere map after recovery.
From: Christian Ott <cott@bethe.tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@399 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
| |
this only matters if wk_atmosphere is used and is based on the argument
that it makes no sense to do this in postrestrictinitial (at least when
vertex centering is used) but makes sense in postregridinitial
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@384 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
| |
this is required since the real valued atmosphere mask is not
checkpointed. It is valid since atmosphere_mask_real is non-zero
wherever atmosphere_mask is non-zero due to GRHydroUpdateAtmosphereMask
and atmosphere_mask is non-zero wherever atmosphere_mask_real is
non-zero because of GRHydroPostSyncAtmosphereMask, ie. the two masks are
non-zero/zero at the same points.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@383 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@381 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@379 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
variables are known on the finer grid when sync_conserved_only = yes. This is fix only affects runs which use a microphysical EOS.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@378 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
| |
This reverts commit f9f8bffea9d960fb9f20cd3fe09ca9460b247adb.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@369 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@362 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
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
| |
* 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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@308 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
| |
* 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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@296 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
bcom
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@295 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
used/allocated already anyway
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@292 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
| |
GRHydro now compiles and runs without MP thorns present and behaves
identical to the old version
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@274 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
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
Schedule all those items which were scheduled "before MoL_PostStep"
also "before MoL_PostStepModify".
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@265 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@264 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
| |
a routine that is scheduled in BASEGRID after ADMBase sets the shift
state.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@261 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
| |
Do not schedule group SetTmunu from GRHydro; group SetTmunu is now
scheduled by thorn TmunuBase.
Correct scheduling of GRHydro_AtmosphereReset, and
GRHydro_AtmosphereResetM, and the associated HydroBase_Boundaries.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@260 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@244 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
|
|
|
|
|
|
|
|
| |
* test suites pass
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@223 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
| |
Lorentz factor variable from GRHydro to HydroBase
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@218 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@206 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
| |
schedule entries for the case with or without evolved tracers -- the
tracers won't get synced if they are not evolved, because they don't
have storage in that case.
* add entropy to the list of variables to be synced
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@203 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Revert previous commit changes to schedule.ccl
Merge MHD and GRHydro routines:
Boundaries, CalcUpdate, PPM, Reconstruct(poly), RegisterGZ,
RegisterVars
Solve a few race conditions arising in the "magnetic" routines.
This is still being tested, since it passes tests with intel
compilers but fails for gnu ones.
Alias function GRHydro_Con2PrimM_pt for C2P2CM test.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@202 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Fix to GRHydro_RegisterVarsM.cc and allows MHD to work without Y_E bing evolved.
Add C2P polytype interface.
Change schedule.ccl to represent a more detailed version of choosing what
gridfucntions get sync'd in Boundaries depending on what is being evolved.
Un-hardwire fixed values of gamma for all MHD Con2Prim routines.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@199 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Several bugfixes, including properly isolating divergence cleaning
code and initialization of the "psidc" gridfunction in
GRHydro_DivergenceClean.F90, as well as properly initializing rhoenth
in GRHydro_HLLEM (thanks to Roland for catching that). Minor
changes to scheduler to properly call MHD and non-MHD routines in
parallel for Riemann solving. Remove GRHydro_ParamCheckM.F90 and
add GRHydro_DivergenceClean.F90.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@196 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@186 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@161 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
Add the magnetized counterparts for several GRHydro routines.
Adjust interface.ccl, param.ccl and schedule.ccl appropriately.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@158 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
| |
to worry about it. This changes the data slitely for EOSG testsuites because this now uses EOSG calls instead of the otherwise fixed calculations of the TOVSolver, thus the change of two testsuites
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@156 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@142 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
-> add comoving b-field: bcom[3] and bcom0
-> add integer variable MHD which is set in ParamCheck to 0 or 1
according to HydroBase::Bvec_evolution_method being GRHydro or something
else (don't want to do string comparisons very often)
-> skeleton routine Conservative2PrimitiveM in separate .F90 file
-> extend keyword Bvec_evolution_method with "GRHydro"
* code compiles, should not hurt anything already there, does not do
anyting (yet)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@115 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@112 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|