NaN error.

[pwp71]

I've tried to use:

u = alpha / 0.00001
or
return false

both work.

I've preferred to use "return false" looking at other code, like:

Code:

        if (D < 0) return false;
        float Q = (-0.5) * (B + ((B < (0.0) ? (-1.0) : (1.0))
          * std::sqrt(D)));

in the same function.

You know this code much better than me, feel free to modify my solution.

pwp71