def rhochange(self):
        """Re-factorise matrix when rho changes."""

        self.lu, self.piv = sl.cho_factor(self.D, self.rho)
        self.lu = np.asarray(self.lu, dtype=self.dtype)
Re-factorise matrix when rho changes.