code: public static QRP factorize(Matrix A) {
          return new QRP(A.numRows(), A.numColumns()).factor(A);
      }
docs: Convenience method to compute a QR decomposition
  
  @param A
  the matrix to decompose (not modified)
  @return Newly allocated decomposition