Esempio n. 1
0
    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])
Esempio n. 2
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])
Esempio n. 3
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])
Esempio n. 4
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)