| Commit message (Collapse) | Author | Age |
|
|
|
|
|
|
|
| |
the C-like coding had some problems, in particular it would compile the same
module more than once and would circumvent Cactus' module dependency tracking.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@520 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
| |
This re-introduces routines that work for hybrid/hot EOS and
corresponding changes in pointwise routines for hot EOS error checking
and temperature treatment by adding old EOSOmni pointwise routine.
From: Philipp Moesta <pmoesta@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@511 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
| |
* Tested successfully pointwisely. It still needs to be tested on evolution.
From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@500 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
| |
this is just a copy of the shock tube test using WENO
From: Roland Haas <rhaas@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@489 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
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
GRHydro_Con2PrimM_pt_EOSOmni.c
* Interestingly Cactus doesn't complain about the function interface in
interface.ccl. It has the old set of arguments for that function and declare it
as written in Fortran instead of C.
* Cylindrical explosion works again with GRHydro_Con2PrimM_pt.c routine. You
can just swap both routines in make.code.defn
From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@427 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
1) Switch to EOSOmni pointwise C2P routine and modify where necessary.
2) Modify Con2PrimM.F90 to allow for the evolution of temperature and adjust the wrapper routine.
3) Create EigenProblemM_hot pointwise routine and call that from HLLEM.F90 when temperature is evolved.
Additionally adjust HLLEM where necessary.
4) Adjust InterfacesM.h to incorporate the newly created functions.
5) Fix a loop problem (not taking into account constraint transport) in PPM reconstruction of Y_e
6) Introduce Prim2ConM_hot and call this pointwise routine from Prim2ConM.F90 when temperature is evolved.
Additionally also make this routine available to initial data routine in GRHydro_InitData
7) Adjust loops in GRHydro_PoloidalMagFieldM.F90 to not set boundary points it cannot set but instead call boundary group afterwards! Pay attention as this will not work with boundary conditions set to "none" in MHD case anymore but is the correct thing to do.
8) Allow StarMapper to extend HydroBase::initial_hydro = "starmapper".
9) Smaller fixes.
From: Philipp Moesta <pmoesta@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@410 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
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Comments:
1) The reconstruct routine is the most expensive routine
in GRHydro and may require some serious thought in the
future regarding optimization.
2) This routine was merged by the end of the last year.
With too many loops, the compiler wasn't able to fully
optimize it (and it would warn about that).
3) With this split, the code became easier to read,
the compiler was able to fully optimize it, but we only
have a marginal gain in the routine performance (3.2%).
Further optimization and clean-up is still desirable.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@252 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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
true content
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@189 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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
-> 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
|
|
|
|
|
|
| |
funtions because apparently unused - could also be implemented in extra thorn if requested
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@34 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|
|
code right now. I have put in the public version of Whisky to start from.
Everybody with commit rights should get commit messages (and the other
way around). It should not be a problem to add people to that list, just
ask. I don't want to get into political problems because someone feels
excluded, but I also don't want to give everyone access per se.
Frank
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@3 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
|