c = self.rho * c + (1 - self.rho) * np.power(g, 2)