Oops C implementation of complex solutions had a bug in the complex sqrt.
(the Haskell implementation was fine)