Update SkatO.cpp

possible fixe for zhanxw#133
lindenb · Jan 26, 2022 · 22f1fa0 · 22f1fa0
1 parent 8defd6f
commit 22f1fa0
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/regression/SkatO.cpp b/regression/SkatO.cpp
@@ -312,7 +312,8 @@ class SkatO::SkatOImpl {
       }
     }
     double temp;
-    if (kappa > lambda.sum() * 10000) {
+    /** VarQ <= 0.0 : fixes: https://github.com/zhanxw/rvtests/issues/133 */
+    if (kappa > lambda.sum() * 10000 || VarQ <= 0.0) {
       temp = 0.0;
     } else {
       double Q = (kappa - MuQ) * sqrt(VarQ - VarZeta) / sqrt(VarQ) + MuQ;