X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=mjo%2Fldlt.py;h=8f39198d19a2007a91667baa34f89fc72b7c73db;hb=48e1ebc120b5b3c068b7a7f86950a4b106ac89eb;hp=8e5f6930da7630d11704458ca14613d9aab64718;hpb=fdafcf90579fd42c869c669f0575dc5812485bc6;p=sage.d.git

diff --git a/mjo/ldlt.py b/mjo/ldlt.py
index 8e5f693..8f39198 100644
--- a/mjo/ldlt.py
+++ b/mjo/ldlt.py
@@ -97,6 +97,7 @@ def ldlt_naive(A):
     diags = A1.diagonal()
     s = diags.index(max(diags))
     P1 = copy(A1.matrix_space().identity_matrix())
+    P1.swap_rows(0,s)
     A1 = P1.T * A1 * P1
     alpha1 = A1[0,0]