Separate HDF5 I/O from algorithm.

[danielballan]

This will have been achieved when it is possible to make Python bindings that pass in a numpy array and get back a numpy array.

[iamfaisalkhan]

@danielballan Wouldn't it be simpler first to start using the C++ interface and see how far we can go with it. Later we can invest some effort into a python binding.

[danielballan]

Sure. The goal here is to separate the I/O from the scientific logic internally. Actually writing the Python bindings is not in scope for this issue; laying track from bindings in general is.