RIT GRMHD dev:

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

1 files changed, 17 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index 399ff7d..a2226f1 100644
--- a/param.ccl
+++ b/param.ccl
@@ -516,6 +516,23 @@ CCTK_REAL cp_dc "The c_p parameter for divergence cleaning"
   0:* :: "Any value, but one to 12 is preferred"
 } 1.0
 
+KEYWORD psidcspeed "Which speed to set for psidc"
+{
+  "char speed"  :: "Based on the characteristic speeds"
+  "light speed" :: "Set the characteristic speeds to speed of light"
+  "set speed"   :: "Manually set the characteristic speeds: [setcharmin,setcharmax]"
+} "light speed"
+
+CCTK_REAL setcharmax "Maximum characteristic speed for psidc if psidcspeed is set"
+{
+  0:1 :: "Any value smaller than speed of light"
+} 1.0
+
+CCTK_REAL setcharmin "Minimum characteristic speed for psidc if psidcspeed is set"
+{
+  -1:0 :: "Any value smaller than speed of light - sign should be negative"
+} -1.0
+
 boolean track_divB "Track the magnetic field constraint violations"
 {
 } "no"