docs: public boolean computeDirect( DMatrixRMaj A ) {
  
          initPower(A);
  
          boolean converged = false;
  
          for( int i = 0; i < maxIterations && !converged; i++ ) {
  //            q0.print();
              
              CommonOps_DDRM.mult(A,q0,q1);
              double s = NormOps_DDRM.normPInf(q1);
              CommonOps_DDRM.divide(q1,s,q2);
  
              converged = checkConverged(A);
          }
  
          return converged;
      }