entanglement entropy added to ES and tested

SA_v2/main.py

@@ -571,7 +571,7 @@ def run_ES(parameters, model:MODEL, utils):
 
     n_step = parameters['n_step']
     n_protocol = 2**n_step
-    exact_data = np.zeros((n_protocol,2), dtype=np.float64) # 15 digits precision
+    exact_data = np.zeros((n_protocol,3), dtype=np.float64) # 15 digits precision
 
     b2_array = lambda n10 : np.array(list(np.binary_repr(n10, width=n_step)), dtype=np.int)
     st=time.time()
@@ -581,7 +581,7 @@ def run_ES(parameters, model:MODEL, utils):
     model.compute_fidelity(psi_evolve=psi)
     model.compute_Sent(psi_evolve=psi)
     model.compute_energy(psi_evolve=psi)
-    print("Est. run time : \t %.3f s"%(0.5*n_protocol*(time.time()-st)))
+    print("Est. run time : \t %.3f s"%(0.5* n_protocol*(time.time()-st)))
     # ---> Starting real calculation <---
 
     st=time.time()
@@ -593,6 +593,8 @@ def run_ES(parameters, model:MODEL, utils):
     outfile = utils.make_file_name(parameters, root=parameters['root'])
     with open(outfile,'wb') as f:
         pickle.dump(exact_data, f, protocol=4)
+
+    print("Saved results in %s"%outfile)
     print("Total run time : \t %.3f s"%(time.time()-st))
     print("\n Thank you and goodbye !")
     f.close()

SA_v2/para.dat

@@ -1,5 +1,5 @@
-task	SD
-L	2
+task	ES
+L	8
 J	1.0
 hz	1.0
 hx_i	-2.
@@ -8,7 +8,7 @@ hx_max	4.
 hx_min	-4.
 dh	8.
 dt	-1
-n_step	100
+n_step	12
 n_sample	50
 n_quench	1000
 T	2.0