def METHOD_NAME(array: f64[:, :, :], result: f64[:, :, :], size1: i32, size2: i32, size3: i32):
    i: i32
    j: i32
    k: i32
    eps: f64
    eps = 1e-12
    for i in range(size1):
        for j in range(size2):
            for k in range(size3):
                assert abs( (1.0 + arcsinh(array[i, j, k])) - result[i, j, k]) <= eps