Flux point method

Implements #50 .