From: Michael Orlitzky <michael@orlitzky.com>
Date: Mon, 28 Sep 2020 15:08:37 +0000 (-0400)
Subject: mjo/ldlt.py: actually permute with the permutation matrix.
X-Git-Url: http://gitweb.michael.orlitzky.com/?a=commitdiff_plain;ds=sidebyside;h=48e1ebc120b5b3c068b7a7f86950a4b106ac89eb;hp=fdafcf90579fd42c869c669f0575dc5812485bc6;p=sage.d.git

mjo/ldlt.py: actually permute with the permutation matrix.
---

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]