def get_S(self): """Calculates the matrix for the coloured noise.""" SST = Constants.kb * (self.C - np.dot(self.T, np.dot(self.C, self.T.T))) # Uses a symetric decomposition rather than Cholesky, since it is more stable return root_herm(SST)
def get_S(self): """Calculates the matrix for the coloured noise.""" SST = Constants.kb*(self.C - np.dot(self.T,np.dot(self.C,self.T.T))) # Uses a symetric decomposition rather than Cholesky, since it is more stable return root_herm(SST)