def __init__(self, inputfile): self.loc_n = pyalps.loadEigenstateMeasurements([inputfile], what='N')[0][0].y self.norb = len(self.loc_n) props = pyalps.getParameters([inputfile]) orbital_order = [] if "orbital_order" in props: orbital_order = list(map(int, props["orbital_order"].split(','))) else: orbital_order = list(map(int, range(1, self.norb + 1))) inv_order = [] for i in range(self.norb): inv_order.append(orbital_order.index(i + 1)) self.orb_order = inv_order empty_diag = np.zeros(self.norb) self.corr_cdag_c = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='dm')[0][0]) self.corr_docc = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='doccdocc')[0][0])
def __init__(self, inputfile): # load data from the HDF5 result file self.nup = assy_vec( pyalps.loadEigenstateMeasurements([inputfile], what='Nup')[0][0]) self.ndown = assy_vec( pyalps.loadEigenstateMeasurements([inputfile], what='Ndown')[0][0]) self.dmup = assy_hc( self.nup, pyalps.loadEigenstateMeasurements([inputfile], what='dm_up')[0][0]) self.dmdown = assy_hc( self.ndown, pyalps.loadEigenstateMeasurements([inputfile], what='dm_down')[0][0])
def __init__(self, inputfile): self.loc_n = pyalps.loadEigenstateMeasurements([inputfile], what='N')[0][0].y[0] self.norb = len(self.loc_n) DMRG_Parms = pyalps.getParameters([inputfile]) orbital_order = map(int, DMRG_Parms[0]['orbital_order'].split(',')) inv_order = [] for i in range(self.norb): inv_order.append(orbital_order.index(i + 1)) self.orb_order = inv_order empty_diag = np.zeros(self.norb) self.corr_cdag_c = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='dm')[0][0]) self.corr_docc = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='doccdocc')[0][0])
def __init__(self, inputfile): self.loc_nup = assy_vec( pyalps.loadEigenstateMeasurements([inputfile], what='Nup')[0][0]) self.loc_ndown = assy_vec( pyalps.loadEigenstateMeasurements([inputfile], what='Ndown')[0][0]) self.loc_docc = assy_vec( pyalps.loadEigenstateMeasurements([inputfile], what='Nupdown')[0][0]) self.norb = len(self.loc_nup) empty_diag = np.zeros(self.norb) self.corr_cdag_up_c_up = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='dm_up')[0][0]) self.corr_cdag_down_c_down = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='dm_down')[0][0]) self.corr_nupnup = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='nupnup')[0][0]) self.corr_nupndown = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='nupndown')[0][0]) self.corr_ndownnup = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='ndownnup')[0][0]) self.corr_ndownndown = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='ndownndown')[0][0]) self.corr_docc = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='doccdocc')[0][0]) self.corr_trans_up = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements( [inputfile], what='transfer_up_while_down')[0][0]) self.corr_trans_down = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements( [inputfile], what='transfer_down_while_up')[0][0]) u1 = pyalps.loadEigenstateMeasurements( [inputfile], what='transfer_up_while_down_at_2')[0][0] u2 = pyalps.loadEigenstateMeasurements( [inputfile], what='transfer_up_while_down_at_1')[0][0] d1 = pyalps.loadEigenstateMeasurements( [inputfile], what='transfer_down_while_up_at_2')[0][0] d2 = pyalps.loadEigenstateMeasurements( [inputfile], what='transfer_down_while_up_at_1')[0][0] self.corr_trans_up_down2 = assy_c(empty_diag, u1, u2) self.corr_trans_up_down1 = assy_c(empty_diag, u2, u1) self.corr_trans_down_up2 = assy_c(empty_diag, d1, d2) self.corr_trans_down_up1 = assy_c(empty_diag, d2, d1) self.corr_trans_pair = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='transfer_pair')[0][0]) self.corr_spinflip = assy_hc( empty_diag, pyalps.loadEigenstateMeasurements([inputfile], what='spinflip')[0][0]) u1 = pyalps.loadEigenstateMeasurements([inputfile], what='nupdocc')[0][0] u2 = pyalps.loadEigenstateMeasurements([inputfile], what='doccnup')[0][0] self.corr_nupdocc = assy_c(empty_diag, u1, u2) self.corr_doccnup = assy_c(empty_diag, u2, u1) u1 = pyalps.loadEigenstateMeasurements([inputfile], what='ndowndocc')[0][0] u2 = pyalps.loadEigenstateMeasurements([inputfile], what='doccndown')[0][0] self.corr_ndowndocc = assy_c(empty_diag, u1, u2) self.corr_doccndown = assy_c(empty_diag, u2, u1)