def get_di_coef(i, n, i0):
    """ 
    Provide the index/indices (di) and coef/coefs (cx)
    for the interpolation at point 'i' in the 'i' direction
    """
    if (i == 0) or (i == n-1):
        cx = [1.]
        di = [0]
    else:
        direc = -1 if (i-i0) % 2 == 0 else +1
        cx = [0.25, 0.75]
        di = [direc, 0]

    # if (i-i0) % 2 == 0:
    #     if (i == 0) or (i==n-1):
    #         cx = [1.]
    #         di = [0]
    #     else:
    #         cx = [0.25, 0.75]
    #         di = [-1, 0]
    # else:
    #     if (i==0) or (i == n-1):
    #         cx = [1.]
    #         di = [0]
    #     else:
    #         cx = [0.25, 0.75]
    #         di = [+1, 0]
    return di, cx