def testSolveUHDU(self): x = ls.solveUHDU(self.A, self.b1, self.n) self.assertMatrixAlmosteEqual(self.x1, x, 14) x = ls.solveUHDU(self.complexA, self.complexb, self.n) self.assertMatrixAlmosteEqual(self.complexx, x, 14)
def testSolveUHDUvsCholesky(self): x_uhdu = ls.solveUHDU(self.randA, self.randb, self.L) x_chol = ls.solveCholesky(self.randA, self.randb, self.L) self.assertMatrixAlmosteEqual(x_uhdu, x_chol, 15)