def METHOD_NAME(self):
    s = self.options['s']
    epsilon = self.options['epsilon']
    uh0 = self.uh0
    uh1 = self.uh1
    @barycentric
    def f0(bcs):
        val0 = uh0(bcs)
        val1 = uh1(bcs)
        val = (val0 + val1)**2/2
        val += epsilon*val0**2/2
        val -= s*val0**3/6
        val += val0**4/24
        return val
    H = self.space.integralalg.mesh_integral(f0)
    self.H.append(H)
    @barycentric
    def f1(bcs):
        val = uh0(bcs)
        return s*val**2/2 - val**3/6
    grad = -self.M*uh0 
    grad -= epsilon*self.M*uh0 
    grad += 2*self.A*uh0 
    grad += self.space.source_vector(f1) 
    grad += self.A*uh1
    self.G.append(np.linalg.norm(grad))