Skip to content

Commit

Permalink
Remove debug prints
Browse files Browse the repository at this point in the history
  • Loading branch information
basilib committed Dec 6, 2024
1 parent 6748cea commit 26eefad
Showing 1 changed file with 1 addition and 14 deletions.
15 changes: 1 addition & 14 deletions vayesta/ewf/ufragment.py
Original file line number Diff line number Diff line change
Expand Up @@ -47,7 +47,6 @@ def get_fragment_energy(self, c1, c2, hamil=None, fock=None, axis1="fragment", c
pxa, pxb = self.get_overlap("proj|cluster-occ")
oaf, obf = np.s_[: pxa.shape[0]], np.s_[: pxb.shape[0]]


# Doubles energy
# TODO: loop to reduce memory?
if hamil is None:
Expand All @@ -66,25 +65,12 @@ def get_fragment_energy(self, c1, c2, hamil=None, fock=None, axis1="fragment", c
assert len(c2) == 4
caa, cab, cba, cbb = c2

print('caa: ', caa.shape)
print('cab: ', cab.shape)
print('cba: ', cba.shape)
print('cbb: ', cbb.shape)

# Remove padding
caa = caa[oaf, oa, va, va]
cab = cab[oaf, ob, va, vb]
cba = cba[obf, oa, vb, va]
cbb = cbb[obf, ob, vb, vb]

print("pxa: ", pxa.shape)
print("gaa: ", gaa.shape)
print("caa: ", caa.shape)
print("gab: ", gab.shape)
print("cab: ", cab.shape)
print("gbb: ", gbb.shape)
print("cbb: ", cbb.shape)
print("pxb: ", pxb.shape)
if c2ba_order == "ab":
cba = cba.transpose(1, 0, 3, 2)
e_doubles = (
Expand All @@ -95,6 +81,7 @@ def get_fragment_energy(self, c1, c2, hamil=None, fock=None, axis1="fragment", c
+ einsum("xi,xjab,iajb", pxa, cab, gab) / 2
+ einsum("xi,xjab,jbia", pxb, cba, gab) / 2
)

else:
assert len(c2) == 3
caa, cab, cbb = c2
Expand Down

0 comments on commit 26eefad

Please sign in to comment.