def METHOD_NAME(self, m):
    self.tmpvec.assign(m)
    # Compute the integral of the control over the domain
    integral = self.smass.inner(self.tmpvec.vector())
    vecmax   = m.vector().max()
    vecmin   = m.vector().min()
    if MPI.rank(MPI.comm_world) == 0:
        print("Current control integral: ", integral)
        print("Maximum of control: ", vecmax)
        print("Minimum of control: ", vecmin)
    return [self.Vol - integral]