aboutsummaryrefslogtreecommitdiff
path: root/src/GRHydro_ENOReconstruct_drv.F90
Commit message (Collapse)AuthorAge
* GRHydro: use ash instead of lsh when allocating arraysrhaas2013-12-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@594 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: work around Intel 11 failures when using Fortran pointersrhaas2013-11-15
| | | | | | | | | | | | this rewrites the Reconstruction driver routines since Intel 11 miscompiles the code using Fortran pointers and generates segfault errors at runtime. Nested subroutines (another solution to make intel 11 work) do not work with intel 13. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@587 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* Declare loop variables OpenMP privateeschnett2013-10-30
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@586 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: use ashX instead of lsh(X) when declaring array sizesrhaas2013-10-29
| | | | | | | | 1.) one should use ash when declaring array sizes 2.) intel11 chokes on lsh(X), this was already in the first Cray pointer patch git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@585 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: replace Fortan pointers with Cray pointersrhaas2013-10-27
| | | | | | | Intel 12 seems to like these better. Otherwise GRHydro_Con2PrimM.F90 takes very long to compile. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@580 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: re-implement prim2con routine in Crhaas2013-07-06
| | | | | | | | | | | | | | | | * new experimental prim2con routine that is about twice as fast as the old one. * Right now tested only for simple EOS. * Right now handles only prim2con call after reconstruction and must be enabled by changing comments in GRHydro_Reconstruct.F90 * remove duplicate code in all reconstruction *_drv.F90 routines and just have the initialization of plus and minus variables in the main GRHydro_Reconstruct.F90 routine From: Christian Ott <cott@bethe.tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@556 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* GRHydro: initialy set y_e_plus/minus to cell center value in reconstructrhaas2013-05-29
| | | | | | | | This ensures that all EOS calls have some sensible y_e to work with. From: Christian Ott <cott@bethe.tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@531 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* * add option to reconstruct temperature to GRHydrorhaas2013-02-11
| | | | | | | | | | | | 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
* 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
* expand region where fluxes are computed for constraint transportrhaas2011-11-29
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@307 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* move TARGET statements in front of CCTK_DECLARE_ARGUMENTSrhaas2011-09-25
| | | | | | | | | this is need for gfortran 4.1 which otherwise assumes that the arrays have already been used (in the kind(gaa) initializations inside of CCTK_DECLARE_ARGUMENTS). git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@278 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* fix non-multipatch support for MHDrhaas2011-09-19
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@276 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* fix non-multipatch velocities in ENO reconstructionrhaas2011-09-19
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@275 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
* add Multipatch support to GRHydrorhaas2011-09-15
| | | | | | | | | | | * 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
* * 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