From 48e1ebc120b5b3c068b7a7f86950a4b106ac89eb Mon Sep 17 00:00:00 2001 From: Michael Orlitzky Date: Mon, 28 Sep 2020 11:08:37 -0400 Subject: [PATCH] mjo/ldlt.py: actually permute with the permutation matrix. --- mjo/ldlt.py | 1 + 1 file changed, 1 insertion(+) 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] -- 2.44.2