NAN from normalisation

Lukasz Laniewski-Wollk · Lukasz Laniewski-Wollk · commit 1408a4b4e474 · 2025-05-19T16:54:10.000+02:00
diff --git a/models/PDE/lif/Dynamics.c.Rt b/models/PDE/lif/Dynamics.c.Rt
@@ -132,7 +132,13 @@ CudaDeviceFunction void Collision()
                 }
             }
         }
-
+        real_t fac = m[0]/mat[0];
+        lambda[0] = lambda[0]*log(fac);
+        for (int j=0; j<k; j++) {
+            for (int g=j; g<k; g++) {
+                mat[j+g*k] = mat[j+g*k] * fac;
+            }
+        }        
         for (int j=0; j<k; j++) {
             res[j] = mat[0+j*k] - m[j];
             for (int g=0; g<j; g++) {

Original file line number	Diff line number	Diff line change
`@@ -132,7 +132,13 @@ CudaDeviceFunction void Collision()`
`132`	`132`	`}`
`133`	`133`	`}`
`134`	`134`	`}`
`135`		`-`
	`135`	`+ real_t fac = m[0]/mat[0];`
	`136`	`+ lambda[0] = lambda[0]*log(fac);`
	`137`	`+ for (int j=0; j<k; j++) {`
	`138`	`+ for (int g=j; g<k; g++) {`
	`139`	`+ mat[j+gk] = mat[j+gk] * fac;`
	`140`	`+ }`
	`141`	`+ }`
`136`	`142`	`for (int j=0; j<k; j++) {`
`137`	`143`	`res[j] = mat[0+j*k] - m[j];`
`138`	`144`	`for (int g=0; g<j; g++) {`