- Real l = a*tl*tl*tl + b*tl*tl + c*tl + d;
- Real u = a*tu*tu*tu + b*tu*tu + c*tu + d;
- if ((l < 0 && u < 0) || (l > 0 && u > 0))
- Debug("Discarding segment (no roots) l = %f (%f), u = %f (%f)", tl, l, tu, u);
+ BReal l = a*tl*tl*tl + b*tl*tl + c*tl + d;
+ BReal u = a*tu*tu*tu + b*tu*tu + c*tu + d;
+ if ((l < BReal(0) && u < BReal(0)) || (l > BReal(0) && u > BReal(0)))
+ {
+ //Debug("Discarding segment (no roots) l = %f (%f), u = %f (%f)", Double(tl), Double(l), Double(tu), Double(u));