@@ -21,17 +21,17 @@ def itebd_heisenberg(chi = 32, J = 1.0, dt = 0.1, CvgCrit = 1.0e-12):
2121
2222 ## Create Si Sj local H with symmetry:
2323 ## SzSz + S+S- + h.c.
24- bdi = cytnx .Bond (BD_IN ,[Qs (1 )>> 1 ,Qs (- 1 )>> 1 ]);
25- bdo = bdi .clone ().set_type (BD_OUT );
26- H = cytnx .UniTensor ([bdi ,bdi ,bdo ,bdo ],labels = ['2' ,'3' ,'1' ,'0' ]);
24+ bdi = cytnx .Bond (BD_IN ,[Qs (1 )>> 1 ,Qs (- 1 )>> 1 ])
25+ bdo = bdi .clone ().set_type (BD_OUT )
26+ H = cytnx .UniTensor ([bdi ,bdi ,bdo ,bdo ],labels = ['2' ,'3' ,'1' ,'0' ])
2727
2828 ## assign:
29- H .at ([0 ,0 ,0 ,0 ]).value = 1 ;
30- H .at ([0 ,1 ,0 ,1 ]).value = - 1 ;
31- H .at ([1 ,0 ,1 ,0 ]).value = - 1 ;
32- H .at ([0 ,1 ,1 ,0 ]).value = 1 ;
33- H .at ([1 ,0 ,0 ,1 ]).value = 1 ;
34- H .at ([1 ,1 ,1 ,1 ]).value = 1 ;
29+ H .at ([0 ,0 ,0 ,0 ]).value = 1
30+ H .at ([0 ,1 ,0 ,1 ]).value = - 1
31+ H .at ([1 ,0 ,1 ,0 ]).value = - 1
32+ H .at ([0 ,1 ,1 ,0 ]).value = 1
33+ H .at ([1 ,0 ,0 ,1 ]).value = 1
34+ H .at ([1 ,1 ,1 ,1 ]).value = 1
3535
3636 ## create gate:
3737 eH = cytnx .linalg .ExpH (H ,- dt )
@@ -42,14 +42,14 @@ def itebd_heisenberg(chi = 32, J = 1.0, dt = 0.1, CvgCrit = 1.0e-12):
4242 # | |
4343 # --A-la-B-lb--
4444 #
45- bd_mid = bdi .combineBond (bdi , True );
46- A = cytnx .UniTensor ([bdi ,bdi ,bd_mid .redirect ()],labels = ['a' ,'0' ,'b' ]);
47- B = cytnx .UniTensor ([bd_mid ,bdi ,bdo ],labels = ['c' ,'1' ,'d' ]);
45+ bd_mid = bdi .combineBond (bdi , True )
46+ A = cytnx .UniTensor ([bdi ,bdi ,bd_mid .redirect ()],labels = ['a' ,'0' ,'b' ])
47+ B = cytnx .UniTensor ([bd_mid ,bdi ,bdo ],labels = ['c' ,'1' ,'d' ])
4848
4949 for b in range (len (B .get_blocks_ ())):
50- cytnx .random .normal_ (B .get_block_ (b ), mean = 0 , std = 0.2 , seed = 0 );
50+ cytnx .random .normal_ (B .get_block_ (b ), mean = 0 , std = 0.2 , seed = 0 )
5151 for a in range (len (A .get_blocks_ ())):
52- cytnx .random .normal_ (A .get_block_ (a ), mean = 0 , std = 0.2 , seed = 0 );
52+ cytnx .random .normal_ (A .get_block_ (a ), mean = 0 , std = 0.2 , seed = 0 )
5353
5454 A .print_diagram ()
5555 B .print_diagram ()
@@ -81,7 +81,7 @@ def itebd_heisenberg(chi = 32, J = 1.0, dt = 0.1, CvgCrit = 1.0e-12):
8181 ## contract all
8282 tmpA = cytnx .Contract (A ,la )
8383 tmpB = cytnx .Contract (B ,lb )
84- X = cytnx .Contract (tmpA ,tmpB ); # << "this line cause problem!\n";
84+ X = cytnx .Contract (tmpA ,tmpB ) # << "this line cause problem!\n"
8585 #X = cytnx.Contract(cytnx.Contract(A,la),cytnx.Contract(B,lb))
8686 #exit(1)
8787 lb .set_label ("e" ,new_label = 'a' )
@@ -138,8 +138,8 @@ def itebd_heisenberg(chi = 32, J = 1.0, dt = 0.1, CvgCrit = 1.0e-12):
138138 # again, but A' and B' are updated
139139 lb_inv = lb .clone ()
140140 for b in range (len (lb_inv .get_blocks_ ())):
141- T = lb_inv .get_block_ (b );
142- lb_inv .put_block_ (1. / T ,b );
141+ T = lb_inv .get_block_ (b )
142+ lb_inv .put_block_ (1. / T ,b )
143143
144144 lb_inv .set_labels (['e' ,'d' ])
145145
0 commit comments