def orthoDen(self): """Routine to orthogonalize the AO Density matrix to orthonormal basis""" self.PO = np.dot(self.U,np.dot(self.P,self.U.T))