more conservative scattering controls

lanl · Dec 23, 2024 · 29cf5ef · 29cf5ef
1 parent f2544f7
commit 29cf5ef
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 5 deletions.
diff --git a/core/decs.h b/core/decs.h
@@ -98,7 +98,7 @@
 
 // Chosen to give some wiggle room
 // in scattering stability criterion
-#define SCATT_BIAS_SAFETY (0.75) // (0.9/RAD_SCATT_TYPES)
+#define SCATT_BIAS_SAFETY (0.5) // (0.9/RAD_SCATT_TYPES)
 
 //#define NUMIN     (1.e10)
 //#define NUMAX     (1.e25)

diff --git a/core/interact.c b/core/interact.c
@@ -17,8 +17,6 @@
 
 #if RADIATION == RAD_TYPE_LIGHT
 #define CONSERVATIVE_BOUND (1)
-#elif RAD_NUM_TYPES >= 4
-#define CONSERVATIVE_BOUND (1)
 #else
 #define CONSERVATIVE_BOUND (0)
 #endif

diff --git a/core/rad_utils.c b/core/rad_utils.c
@@ -137,8 +137,8 @@ void update_superphoton_resolution(
     correction = ideal / real;
 
     // Limit strength of correction
-    correction = MY_MIN(correction, 1.5); // MY_MIN(correction, 2.0);
-    correction = MY_MAX(correction, 0.5);
+    correction = MY_MIN(correction, 1.25); // MY_MIN(correction, 2.0);
+    correction = MY_MAX(correction, 0.25);
 
     // If no superphotons are being emitted (yet) don't modify emission strength
     if (real < SMALL)