def METHOD_NAME(self, p):
    """ The right hand side of Possion equation
    INPUT:
        p: array object,  
    """
    f = np.array([0.0], np.float64)
    shape = len(p.shape[:-1])*(1, )
    return f.reshape(shape) 