Switch back to choleski inverse

haziqj · May 29, 2024 · 5e2a9ea · 5e2a9ea
1 parent 15a288b
commit 5e2a9ea
Showing 1 changed file with 4 additions and 5 deletions.
diff --git a/R/20-rgeneric.R b/R/20-rgeneric.R
@@ -130,11 +130,10 @@ inla_sem <- function(
       front <- Lambda %*% solve(IminB)
     }
     Sigma <- front %*% tcrossprod(Psi, front) + Theta
-    # Sigma <- Sigma + diag(1e-10, nrow(Sigma))  # for stability
-    # chol2inv(chol(Sigma))
-    # MASS::ginv(Sigma)
-    Sigma <- Matrix::forceSymmetric(Matrix::Matrix(Sigma))
-    Matrix::solve(Sigma)
+    Sigma <- Sigma + diag(1e-5, nrow(Sigma))  # for stability
+    chol2inv(chol(Sigma))
+    # Sigma <- Matrix::forceSymmetric(Matrix::Matrix(Sigma))
+    # Matrix::solve(Sigma)
   }
 
   mu <- function() { numeric(0) }