Bug in mcx with negative arguments to pow(x,2.0), square fixes it
Polishing missing still, otherwise done?
Doesn't work so far
Not especially fast
Lots of testing remains for other numerical types
