Flux point method

Implements #50 .
Why call `np.round` here ... shouldn't be necessary.