@@ -23,7 +23,7 @@ REFERENCES:
23
23
# it under the terms of the GNU General Public License as published by
24
24
# the Free Software Foundation, either version 2 of the License, or
25
25
# (at your option) any later version.
26
- # http ://www.gnu.org/licenses/
26
+ # https ://www.gnu.org/licenses/
27
27
# *****************************************************************************
28
28
29
29
from libc.math cimport log, ceil
@@ -134,7 +134,8 @@ def OP_represent(int n, merges, perm):
134
134
if not OP_find(OP, i) == i:
135
135
print (" Failed at i = %d !" % i)
136
136
good = False
137
- if good: print (" Each element reports itself as its root." )
137
+ if good:
138
+ print (" Each element reports itself as its root." )
138
139
print (" Merging:" )
139
140
for i,j in merges:
140
141
OP_join(OP, i, j)
@@ -296,14 +297,16 @@ cdef int PS_first_smallest(PartitionStack *PS, bitset_t b, int *second_pos=NULL,
296
297
n = i - j + 1
297
298
location = j
298
299
j = i + 1
299
- if PS.levels[i] == - 1 : break
300
+ if PS.levels[i] == - 1 :
301
+ break
300
302
i += 1
301
303
# location now points to the beginning of the first, smallest,
302
304
# nontrivial cell
303
305
i = location
304
306
while 1 :
305
307
bitset_flip(b, PS.entries[i])
306
- if PS.levels[i] <= PS.depth: break
308
+ if PS.levels[i] <= PS.depth:
309
+ break
307
310
i += 1
308
311
309
312
if second_pos != NULL :
@@ -480,7 +483,8 @@ def PS_represent(partition, splits):
480
483
if not PS.degree == n or not PS.depth == 0 :
481
484
print (" Incorrect degree or depth!" )
482
485
good = False
483
- if good: print (" Everything seems in order, deallocating." )
486
+ if good:
487
+ print (" Everything seems in order, deallocating." )
484
488
PS_dealloc(PS)
485
489
print (" Deallocated." )
486
490
print (" Creating PartitionStack from partition %s ." % partition)
@@ -633,11 +637,13 @@ cdef inline int SC_realloc_gens(StabilizerChain *SC, int level, int size):
633
637
cdef int n = SC.degree
634
638
635
639
temp = < int * > sig_realloc( SC.generators[level], n * size * sizeof(int ) )
636
- if temp is NULL : return 1
640
+ if temp is NULL :
641
+ return 1
637
642
SC.generators[level] = temp
638
643
639
644
temp = < int * > sig_realloc( SC.gen_inverses[level], n * size * sizeof(int ) )
640
- if temp is NULL : return 1
645
+ if temp is NULL :
646
+ return 1
641
647
SC.gen_inverses[level] = temp
642
648
643
649
SC.array_size[level] = size
@@ -751,7 +757,8 @@ cdef int SC_realloc_bitsets(StabilizerChain *SC, unsigned long size):
751
757
Returns 1 in case of an allocation failure.
752
758
"""
753
759
cdef unsigned long size_old = SC.gen_used.size
754
- if size <= size_old: return 0
760
+ if size <= size_old:
761
+ return 0
755
762
cdef unsigned long new_size = size_old
756
763
while new_size < size:
757
764
new_size *= 2
@@ -1041,8 +1048,10 @@ cdef int SC_insert_and_sift(StabilizerChain *SC, int level, int *pi, int num_per
1041
1048
break
1042
1049
else :
1043
1050
bitset_set(& SC.gen_is_id, perm_gen_index)
1044
- if b != - 1 : break
1045
- if b == - 1 : return 0
1051
+ if b != - 1 :
1052
+ break
1053
+ if b == - 1 :
1054
+ return 0
1046
1055
if sift and level == SC.base_size:
1047
1056
SC_add_base_point(SC, b)
1048
1057
else :
@@ -1064,18 +1073,22 @@ cdef int SC_insert_and_sift(StabilizerChain *SC, int level, int *pi, int num_per
1064
1073
for i from 0 <= i < SC.orbit_sizes[level]:
1065
1074
x = SC.base_orbits[level][i]
1066
1075
for perm_gen_index from 0 <= perm_gen_index < num_perms:
1067
- if sift and bitset_check(& SC.gen_is_id, perm_gen_index): continue
1076
+ if sift and bitset_check(& SC.gen_is_id, perm_gen_index):
1077
+ continue
1068
1078
perm = pi + n* perm_gen_index
1069
1079
if SC.parents[level][perm[x]] == - 1 :
1070
1080
# now we have an x which maps to a new point under perm,
1071
1081
re_treed = 1
1072
- if sift: bitset_set(& SC.gen_used, perm_gen_index)
1082
+ if sift:
1083
+ bitset_set(& SC.gen_used, perm_gen_index)
1073
1084
if SC_re_tree(SC, level, perm, x):
1074
1085
return 1
1075
1086
start_over = 1 # we must look anew
1076
1087
break
1077
- if start_over: break
1078
- if not re_treed: continue
1088
+ if start_over:
1089
+ break
1090
+ if not re_treed:
1091
+ continue
1079
1092
for perm_gen_index from 0 <= perm_gen_index < old_num_gens:
1080
1093
perm = SC.generators[level] + n* perm_gen_index
1081
1094
if SC.parents[level][perm[x]] == - 1 :
@@ -1084,7 +1097,8 @@ cdef int SC_insert_and_sift(StabilizerChain *SC, int level, int *pi, int num_per
1084
1097
return 1
1085
1098
start_over = 1 # we must look anew
1086
1099
break
1087
- if start_over: break
1100
+ if start_over:
1101
+ break
1088
1102
for j from level < j < SC.base_size:
1089
1103
for perm_gen_index from 0 <= perm_gen_index < SC.num_gens[j]:
1090
1104
perm = SC.generators[j] + n* perm_gen_index
@@ -1361,7 +1375,8 @@ def SC_test_list_perms(list L, int n, int limit, bint gap, bint limit_complain,
1361
1375
if gap:
1362
1376
G = PermutationGroup([[i+ 1 for i in p] for p in L])
1363
1377
if G.order() > limit:
1364
- if limit_complain: print (' TOO BIG' )
1378
+ if limit_complain:
1379
+ print (' TOO BIG' )
1365
1380
return
1366
1381
SC = SC_new(n)
1367
1382
cdef int * perm = < int * > sig_malloc(n * (len (L)+ 3 ) * sizeof(int ))
0 commit comments