Flux point method

Implements #50 .
Can you get this test to work or do you want to do this in a future pull request?