6465666768697071727374
} public void testDenseLU() { int n = A.numRows(); DenseLU lu = new DenseLU(n, n); lu.factor(A.copy()); lu.solve(I); Matrix J = I.mult(A, new DenseMatrix(n, n)); for (int i = 0; i < n; ++i)
8081828384858687888990
} public void testDenseLUtranspose() { int n = A.numRows(); DenseLU lu = new DenseLU(n, n); lu.factor(A.copy()); lu.transSolve(I); Matrix J = I.transAmult(A, new DenseMatrix(n, n)); for (int i = 0; i < n; ++i)
96979899100101102103104105106
} public void testDenseLUrcond() { int n = A.numRows(); DenseLU lu = new DenseLU(n, n); lu.factor(A.copy()); lu.rcond(A, Matrix.Norm.One); lu.rcond(A, Matrix.Norm.Infinity); }