Also not clear if this is the best way to include it.  Probably need the std::variant approach here too.

I think g++ getting confused by copy elision on return of Eigen temporary?

Python compilation working again

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