RZ good zeros to t=1M + anti-non BH convergence failsafe

Author: MathisSinet

src/Theories/CTs/RZ.ts

@@ -288,6 +288,8 @@ class rzSim extends theoryClass<theory> {
     bhRewindNorm: number;
     bhRewindDeriv: number;
 
+    bhProcessCounter: number = 0;
+
     getBuyingConditions(): conditionFunction[] {
         const activeStrat = [
             () => {
@@ -437,6 +439,10 @@ class rzSim extends theoryClass<theory> {
 
     bhProcess(zResult: ComplexValue | null = null, tmpZ: ComplexValue | null = null) {
         this.offGrid = true;
+        if (this.strat != "RZdBHRewind") this.bhProcessCounter++;
+
+        // This is a method to prevent the non-convergence of the algorithm
+        const bhLockPreventionCoeff = this.bhProcessCounter < 1000 ? 1 : Math.E ** (-0.002*(this.bhProcessCounter - 1000)); 
 
         if (zResult === null){
             zResult = zeta(this.t_var, this.ticks, this.offGrid, lookups.zetaLookup);
@@ -451,11 +457,11 @@ class rzSim extends theoryClass<theory> {
         if(this.bhSearchingRewind && this.t_var > 14.5 && bhdt > 0)
         {
             let srdt = -Math.min(0.125 / bhdt, 0.125);
-            this.t_var += srdt;
+            this.t_var += srdt * bhLockPreventionCoeff;
         }
         else
         {
-            this.t_var += bhdt;
+            this.t_var += bhdt * bhLockPreventionCoeff;
             this.bhSearchingRewind = false;
             if(Math.abs(bhdt) < 1e-9)
             {

src/Theories/CTs/helpers/RZgoodzeros.json

@@ -51,7 +51,11 @@
         148190.29,
         175030.58,
         210283.59,
-        321842.6
+        321842.6,
+        454350.9,
+        757720.12,
+        833464.81,
+        964594.73
     ],
     "longZeros": [
         7623.42,
@@ -120,7 +124,7 @@
             "toRho": 2300, "from": 50000, "to": 150000
         },
         {
-            "toRho": 2500, "from": 90000, "to": 999999999
+            "toRho": 2500, "from": 90000, "to": 400000
         },
         {
             "toRho": 999999, "from": 200000, "to": 999999999
@@ -155,7 +159,7 @@
             "toRho": 2300, "from": 50000, "to": 150000
         },
         {
-            "toRho": 2500, "from": 90000, "to": 999999999
+            "toRho": 2500, "from": 90000, "to": 400000
         },
         {
             "toRho": 999999, "from": 200000, "to": 999999999