From: Michael Orlitzky 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]