use scipy.special.cbrt in the tests instead of **(1.0/3.0)