Break up pybind11 interface into separate compilation units (not as beneficial as hoped); still very slow

And JSON I/O

Also add the accessors as methods of the class in the python interface (overload resolution has large runtime overhead)

Refactor what is implemented

Bug in mcx with negative arguments to pow(x,2.0), square fixes it

Not especially fast

Lots of testing remains for other numerical types