def METHOD_NAME(self, rho, mesh, smooth=True):
    self.rho = rho
    self.mesh = mesh
    self.area = mesh.entity_measure('cell')
    if smooth is True:
        self.smooth_rho()
        #self.loop_smooth_rho()
    self.eta = self.compute_eta()