aboutsummaryrefslogtreecommitdiff
path: root/src/GRHydro_Reconstruct.F90
Commit message (Collapse)AuthorAge
* GRHydro: remove unused variablesrhaas2013-03-28
| | | | | | From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@495 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Added MP5 reconstruction method. Tested with TOV and shocktube.rhaas2013-03-06
| | | | | | From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@487 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Improved WENO implementation: (i) Pick adaptive epsilon-parameter ↵rhaas2013-01-14
| | | | | | | | 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
* GRHydro: Implemented WENO5. Tested with a shock tube.rhaas2013-01-14
| | | | | | From: Christian Reisswig <reisswig@scriwalker.(none)> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@465 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: Implement entropy evolution (separated from temp evolution)rhaas2013-01-14
| | | | | | 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
* GRHydro: Introduce variable atmopshere level / tolerance as function of ↵rhaas2013-01-14
| | | | | | | | 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
* GRHydro: Option to reconstruct W*vel (W = Lorentz factor) instead of just ↵rhaas2012-08-27
| | | | | | | | | | | | | 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
* GRHydro: Changes fixing issues with hot EOS treatment:rhaas2012-08-27
| | | | | | | | | | | | | | | 1) Set correct keytemp in C2P 2) Fix OpenMP private variable declaration in HLLEM 3) Fix velocity arguments for Y_e 1D PPM reconstruction 4) If in atmosphere and if evolve a Y_e also reset Y_e with cell averaged value in Reconstruction 5) Fix index permutation bug in ReconstructPoly routines for divergence cleaning field 6) Remove extra comment character in UpdateMaskM 7) Fix a missing velocity permutation when calling 1D reconstruction for Y_e. 8) Insert missing sqrt(det) factor for pressure term in resetting tau when temperature got too cold with a hot EOS. From: Philipp Moesta <pmoesta@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@415 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
* move TARGET statement in Reconstruction to compile with gcc 4.1rhaas2011-11-28
| | | | | | | gcc 4.1 requires TARGET statements to be before DECLARE_CCTK_ARGUMENTS git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@303 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Implemented enhanced PPM scheme by Colella & Sekora 2008, McCorquodale & ↵reisswig2011-11-26
| | | | | | | | | | | | 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
* * Optimize: remove support for shift_state = 0 (except for shock tubes andcott2011-08-02
| | | | | | | | | | Cowling calculations of spherically symmetric objects, there is no reason not to have storage for the shift). git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@259 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * modify reconstruction:cott2011-08-02
| | | | | | | | | | | * get rid of nonsensical re-computation of the conserved tracer * remove support for shift_state = 0 and get rid of unnecessary assignments. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@256 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* RIT GRMHD dev:bmundim2011-08-02
| | | | | | | | | | | | | | | | | | | | | | | | | | | 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
* RIT GRMHD dev: split reconstruct routine.bmundim2011-07-20
| | | | | | | | | | | | | | | | | | | | | 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
* Declare all OpenMP private variables explicitly as privateeschnett2011-06-13
| | | | | | | | | In Fortran OpenMP, loop variables of parallel loops are implicitly private and do not need to be declared as private. To improve code clarity, all such variables are now explicitly declared to be private. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@251 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* JOsh's fixes.bmundim2011-04-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@246 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* MERGE divergence cleaning feature into trunk.tbode2011-04-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@244 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * remove support for "General" EOS interface (EOSG_*)cott2011-03-28
| | | | | | | | | | * 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
* more OpenMP WORKSHARE sections in GRHydro Reconstructrhaas2011-03-02
| | | | | | | | multiple WORKSHARE sections are used since the standard does only seem to allow assignments and WHERE statements in WORKSHARE sections git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@220 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * fix a very nasty bug in PPM reconstruction that sneeked in with the very ↵cott2011-01-08
| | | | | | | | | | recent PPM changes. Old version of the code were not affected. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@209 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* RIT MHD development:bmundim2010-12-31
| | | | | | | | | | | | | | | | | | 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
* merge branch hot_and_MHD_temp_dev into branch at revision r185knarf2010-11-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@186 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * remove dependence on StaticConformalcott2010-08-27
| | | | | | | | | * change calculation of the determinant of the 3-metric from a subroutine call to a macro. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@152 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* file/parameter string replacement from whisky to GRHydrobmundim2010-05-02
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@112 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-26 23:19:15 +0000