for v in range(0, nv):
for w in range(0, nw+1):
for s in range(0, ns):
- sigmas = [ (s + sx + 0.)/ns * tau for sx in 0,1 ]
+ sigmas = [ (s + sx + 0.5)/ns * tau for sx in 0,1 ]
#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) ])