+/*
+ s*(2+s) is negative for s = [-2..0]
+ slightly negative values for s due to floating point rounding errors
+ cause nan for sqrt(s*(2+s))
+ We can probably clamp the resulting parable to positive numbers
+*/
+ if( s >= -2.0 && s <= 0.0 )
+ s = 0.0;