interpolate(
            interpolate(
                cube[0, 0, 0], cube[1, 0, 0], r[0]
            ),  # trilinear interpolation => http://en.wikipedia.org/wiki/Trilinear_interpolation
            interpolate(cube[0, 1, 0], cube[1, 1, 0], r[0]),