def METHOD_NAME(self, p): # 这里在边界上，始终知道边界外面的浓度
    x = p[..., 0]
    y = p[..., 1]
    val = np.zeros(p.shape[:-1], dtype=np.float64)
    flag0 = (np.abs(x) < 1e-13) & (y < 1/16)
    flag1 = (np.abs(y) < 1e-13) & (x < 1/16)
    val[flag0 | flag1] = 1
    return val