SimpleMatrix R = SimpleMatrix.wrap(alg.getR(null,true));
SimpleMatrix A_found = Q.mult(R);
SimpleMatrix A = SimpleMatrix.wrap(orig);
return A.minus(A_found).normF()/A.normF();
}
public static double evaluate( QRPDecomposition<DenseMatrix64F> alg , DenseMatrix64F orig ) {
double maxValue = CommonOps.elementMaxAbs(orig);