SimpleMatrix P = SimpleMatrix.wrap(alg.getPivotMatrix(null));
SimpleMatrix A_found = Q.mult(R);
SimpleMatrix A = SimpleMatrix.wrap(orig);
return A.mult(P).minus(A_found).normF()/A.normF();
}
private static void runAlgorithms( DenseMatrix64F mat )
{
System.out.println("qr = "+ evaluate(new QRDecompositionHouseholder(),mat));