Only change the permutation using rcm if the matrix isn't block

nlmixr2 · Aug 31, 2024 · 221c538 · 221c538
1 parent c7af1c2
commit 221c538
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/R/lotri.R b/R/lotri.R
@@ -1058,8 +1058,8 @@ NULL
   dimnames(.ret) <- list(env$names, env$names)
   dimnames(.retF) <- list(env$names, env$names)
   dimnames(.retU) <- list(env$names, env$names)
-  if (is.logical(env$rcm) && env$rcm && .n >= 1) {
-    .tmp <- lotriMatInv(.ret)
+  if (is.logical(env$rcm) && env$rcm && .n >= 1 &&
+        !lotriIsBlockMat(.ret)) {
     .ret <- rcm(.ret)
     env$names <- dimnames(.ret)[[1]]
     .retF <- .retF[env$names, env$names]