//Newton-Raphson root polishing (we don't care about bounds, assumes very near the desired root)
static Real polish_cubicroot(Real a, Real b, Real c, Real d, Real t, Real *dpdt)
{
const Real cn[4] = {a,b,c,d};
Real p,dp,newt,oldpmag=FLT_MAX;
//Newton-Raphson root polishing (we don't care about bounds, assumes very near the desired root)
static Real polish_cubicroot(Real a, Real b, Real c, Real d, Real t, Real *dpdt)
{
const Real cn[4] = {a,b,c,d};
Real p,dp,newt,oldpmag=FLT_MAX;