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

        self.lu, self.piv = sl.lu_factor(self.Z, self.rho)
        self.lu = np.asarray(self.lu, dtype=self.dtype)
def rhochange ( self ) : self . lu , self . piv = sl . lu_factor ( self . Z , self . rho ) self . lu = np . asarray ( self . lu , dtype = self . dtype )
Re-factorise matrix when rho changes