MC: Faster User-RPL HILBERT

[Thomas Okken]

To elaborate a bit...

A naive Java implementation would be like this:

Code:

double[][] hilbert(int n) {
    double[][] h = new double[n][n];
    for (int i = 0; i < n; i++)
        for (int j = 0; j < n; j++)
            h[i][j] = 1.0 / (i + j + 1);
    return h;
}

while the division-minimizing approach could also be coded with nested loops, like this:

Code:

double[][] hilbert(int n) {
    double[][] h = new double[n][n];
    for (int j = 0; j < n; j++)
        h[0][j] = 1.0 / (j + 1);
    for (int i = 1; i < n; i++) {
        for (int j = 0; j < n - 1; j++)
            h[i][j] = h[i - 1][j + 1];
        h[i][n - 1] = 1.0 / (i + n);
    }
    return h;
}