d = self.rho * d + (1 - self.rho) * np.power(update, 2)