Optimisations

Forgot to push the actual code I used in the thesis. Here it is
aethrvmn · Sep 18, 2021 · dff5b10 · dff5b10
1 parent b9440c7
commit dff5b10
Show file tree

Hide file tree

Showing 24 changed files with 14 additions and 168 deletions.
diff --git a/.gitignore b/.gitignore
@@ -1,3 +1 @@
-
-__pycache__/chain.cpython-39.pyc
-__pycache__/logchain.cpython-39.pyc
+__pycache__
diff --git a/Figures/Fixed/eta100.png b/Figures/Fixed/eta100.png
diff --git a/Figures/Fixed/eta1000.png → Figures/Fixed/eta_distribution.png b/Figures/Fixed/eta1000.png → Figures/Fixed/eta_distribution.png
diff --git a/Figures/Fixed/near_singlet100.png b/Figures/Fixed/near_singlet100.png
diff --git a/Figures/Fixed/near_singlet1000.png b/Figures/Fixed/near_singlet1000.png
diff --git a/Figures/Matrix/midrg1.png b/Figures/Matrix/midrg1.png
diff --git a/Figures/Matrix/midrg10.png b/Figures/Matrix/midrg10.png
diff --git a/Figures/Matrix/midrg13.png b/Figures/Matrix/midrg13.png
diff --git a/Figures/Matrix/midrg5.png b/Figures/Matrix/midrg5.png
diff --git a/Figures/Matrix/singlet-phase.png b/Figures/Matrix/singlet-phase.png
diff --git a/Figures/distance.png b/Figures/distance.png
diff --git a/Figures/frog.png b/Figures/frog.png
diff --git a/Figures/rg.gif b/Figures/rg.gif
diff --git a/README.md b/README.md
@@ -1,5 +1,3 @@
 # Ongoing Thesis Project
 
 The goal of this code is to model the RG transformation for the 1d Random Antiferromagnetic Heisenberg Quantum Spin Chain.
-
-# WIP
diff --git a/Thesis.pdf b/Thesis.pdf
diff --git a/__pycache__/chain.cpython-37.pyc b/__pycache__/chain.cpython-37.pyc
diff --git a/__pycache__/testchain.cpython-39.pyc b/__pycache__/testchain.cpython-39.pyc
diff --git a/chain.py b/chain.py
@@ -8,9 +8,9 @@ class Random_chain:
 
  version = "v5"
 
- def __init__(self, bond_number, ceiling):
+ def __init__(self, bond_number):
  self.length = int(bond_number)
- self.ceiling = ceiling
+ self.ceiling = 1
  self.matrix_creator()
  self.bonds()
 
@@ -73,9 +73,4 @@ def rescale(self): # Second Step of the RG
  def renormalization(self): # Self-explanatory
  self.decimate()
  self.rescale()
- return self
-
- def parameters(self): # Different physical parameters that can be calculated
- self.distance = int(np.abs(self.right_index[1]-self.left_index[0]))
- self.system_energy = -(1/4)*np.sum(self.bond_matrix[self.bond_matrix > 0])
  return self
diff --git a/test.py → flow.py b/test.py → flow.py
@@ -3,8 +3,10 @@
 from chain import Random_chain
 from tqdm import tqdm
 
+plt.ion
+
 length = 1000
-matrix = Random_chain(length,1)
+matrix = Random_chain(length)
 eta = []
 
 for i in tqdm(range(int(length/2))):
@@ -15,24 +17,15 @@
 range=np.arange(len(eta))
 sum = []
 
-for i in tqdm(range):
+for i in range:
  sum = np.append(sum, np.sum(eta[i:]))
 
 plt.figure(1)
 plt.plot(np.sort(sum/sum.max()), linewidth=1)
 plt.plot(1-np.exp(-range), '--r', linewidth=1)
-plt.ylim(0.9, 1.01)
-plt.xlim(420, 500)
-plt.yscale('log')
-plt.yticks([1,1])
-plt.ylabel('Value of $\eta$ near the random singlet phase (rescaled)')
-plt.xlabel('Iteration step')
-plt.savefig('Figures/Fixed/near_singlet1000.png')
-
-plt.figure(2)
-plt.plot(np.sort(sum/sum.max()), linewidth=1)
-plt.plot(1-np.exp(-range), '--r', linewidth=1)
 plt.ylabel('Value of $\eta$ during the SDRG (rescaled)')
 plt.xlabel('Iteration step')
 plt.yscale('log')
-plt.savefig('Figures/Fixed/eta1000.png')
+plt.savefig('Figures/Fixed/eta_distribution.png')
+
+plt.show()
diff --git a/flow100.py b/flow100.py
diff --git a/matrix.py b/matrix.py
@@ -6,7 +6,7 @@
 from tqdm import tqdm
 
 length = 30
-matrix = Random_chain(length,1)
+matrix = Random_chain(length)
 
 fig = plt.figure()
 
@@ -31,8 +31,8 @@
  img = plt.imshow(matrix.bond_matrix, interpolation='nearest', cmap=cmap, norm=norm, animated=True)
  ims.append([img])
 
+img.figure.savefig('Figures/Matrix/singlet-phase.png', pad_inches=0, transparent=True)
+
 ani = animation.ArtistAnimation(fig, ims, interval=200, blit=True, repeat_delay=10000)
 ani.save('Figures/rg.gif')
-
-img.figure.savefig('Figures/Matrix/singlet-phase.png', pad_inches=0, transparent=True)
 plt.show()
diff --git a/parameters.py b/parameters.py
diff --git a/requirements.txt b/requirements.txt
@@ -1,6 +1,4 @@
-scipy
 numpy
 matplotlib
 tqdm
-seaborn
-statsmodels
+PyQt5
diff --git a/test2.py b/test2.py