def METHOD_NAME(ts, coors, mode="qp", **kwargs):
    # t = ts.dt * ts.step
    eps = diffcoef
    sin = nm.sin
    cos = nm.cos
    exp = nm.exp
    sqrt = nm.sqrt
    pi = nm.pi
    if mode == "qp":
        x_1 = coors[..., 0]
        x_2 = coors[..., 1]
        res = -2*pi*(x_2**2 - x_2)*cos(2*pi*x_1)\
              - 2*(2*pi**2*(x_2**2 - x_2)*sin(2*pi*x_1) - sin(2*pi*x_1))*eps\
              - (2*x_2 - 1)*sin(2*pi*x_1)
        return {"val": res[..., None, None]}