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
