def METHOD_NAME(self, p):
    x = p[..., 0, None]
    y = p[..., 1, None]
    Fx = self.Fx(x, y)
    Fy = self.Fy(x, y)
    if type(Fx) is not np.ndarray:
        Fx = np.ones(x.shape, dtype=np.float_)*Fx
    if type(Fy) is not np.ndarray:
        Fy = np.ones(x.shape, dtype=np.float_)*Fy
    f = np.c_[Fx, Fy] 
    return f 