docs: public static double[][] identity(final int m, final int n) {
      final double[][] A = new double[m][n];
      final int dim = m < n ? m : n;
      for(int i = 0; i < dim; i++) {
        A[i][i] = 1.0;
      }
      return A;
    }