wire = 2.000;
# number of wires
-kv = 24
-kw = 16
+kv = 7
+kw = 4
# resolution
-nv = 80
-nw = 80
+nv = 7
+nw = 20
ns = 4 # around tube, should be even
nv += -nv % kv
#print('VWS',v,w,s, sigmas, file=sys.stderr)
if not w % each_w:
extents[w % nw].quad([ cnr for cnr in qc(v,w,sigmas,False) ])
- if not v % each_v and w < nw:
+ if False and not v % each_v and w < nw:
travers[v].rquad([ cnr for cnr in qc(v,w,sigmas,True) ])
- if not v % each_v:
+ if False and not v % each_v:
for w in 0, nw:
p, norm, extt, trav = m.details(v, w)
cnrs = points_round_wire(p, norm, extt, calc_sigmas(range(0,ns)))