docs: public static double computeTauAndDivide(final int j, final int numRows ,
                                               final double[] u , final double max) {
          double tau = 0;
  //        double div_max = 1.0/max;
  //        if( Double.isInfinite(div_max)) {
              for( int i = j; i < numRows; i++ ) {
                  double d = u[i] /= max;
                  tau += d*d;
              }
  //        } else {
  //            for( int i = j; i < numRows; i++ ) {
  //                double d = u[i] *= div_max;
  //                tau += d*d;
  //            }
  //        }
          tau = Math.sqrt(tau);
  
          if( u[j] < 0 )
              tau = -tau;
  
          return tau;
      }