mjo/ldlt.py: fix a dimension in block_ldlt_naive().

diff --git a/mjo/ldlt.py b/mjo/ldlt.py
index c8483c707afa6f3b67935e451b7c979401661b09..4955d9fa9248f2bde240f70a526062f927aae9b5 100644 (file)
--- a/mjo/ldlt.py
+++ b/mjo/ldlt.py
@@ -325,7 +325,7 @@ def block_ldlt_naive(A, check_hermitian=False):
 
     # The top-left 2x2 submatrix is now our pivot.
     E = A1[:2,:2]
-    C = A1[2:n,0]
+    C = A1[2:n,0:2]
     B = A1[2:,2:]
 
     if B.nrows() == 0: