def METHOD_NAME(array: f64[:, :, :, :], result: f64[:, :, :, :], size1: i32, size2: i32, size3: i32, size4: i32):
    i: i32; size: i32;
    eps: f64 = 1e-12
    shape: i32[1] = empty((1,), dtype=int32)
    size = size1 * size2 * size3 * size4
    shape[0] = size
    array1d: f64[12800] = reshape(array, shape)
    result1d: f64[12800] = reshape(result, shape)
    for i in range(size):
        assert abs((sinh(array1d[i]) + 2.0)/2.0 - result1d[i]) <= eps