for b in bests:
along = unit_v(np.cross(b[1], unit_z))
upwards = np.cross(b[1], along)
- matrix = [ b[1],
- upwards,
- along ]
+ matrix = [ along,
+ -upwards,
+ -b[1] ]
print('initial', matrix, file=sys.stderr)
matrix = np.concatenate((matrix, [[0,0,0]]))
print('concatd', matrix, file=sys.stderr)