fix determinant calculation

beast-dev · Dec 17, 2023 · fd2a747 · fd2a747
1 parent db1fa38
commit fd2a747
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/dr/math/distributions/GaussianMarkovRandomField.java b/src/dr/math/distributions/GaussianMarkovRandomField.java
@@ -243,7 +243,7 @@ private double getLogDeterminant() {
 
         if (!isImproper()) {
             double lambda = lambdaParameter.getParameterValue(0);
-            logDet += (1 - dim) * Math.log(1 - lambda * lambda);
+            logDet += Math.log(1 - lambda * lambda);
         }
 
         if (CHECK_DETERMINANT) {
@@ -260,7 +260,7 @@ private double getLogDeterminant() {
             }
 
             if (Math.abs(sum - logDet) > 1E-6) {
-                throw new RuntimeException("Incorrect pseudo-determinant");
+                throw new RuntimeException("Incorrect (pseudo-) determinant");
             }
         }