docs: protected double estimateEigenvalue(double[][] mat, double[] in) {
      double de = 0., di = 0.;
      // Matrix multiplication:
      for(int d1 = 0; d1 < in.length; d1++) {
        final double[] row = mat[d1];
        double t = 0.;
        for(int d2 = 0; d2 < in.length; d2++) {
          t += row[d2] * in[d2];
        }
        final double s = in[d1];
        de += t * s;
        di += s * s;
      }
      return de / di;
    }