else
b <- d
c <- t / b
- d <- sqrt(0.75) * (b - c)
+ d <- sqrt 0.75 * (b - c)
r.[2, 2] <- d
b <- b + c
c <- -0.5 * b - s
for k = 1 to 3 do
r.[2, k] <- 0.0
-let private biquadroots (p : float[]) (r : float[,]) =
+let inline private biquadroots (p : float[]) (r : float[,]) =
if p.[0] <> 1.0 then
for k = 1 to 4 do
p.[k] <- p.[k] / p.[0]
let x2 = -x1
#if DEBUG_LOG
- printf "\n\tx1=%f, y1=%f, A=%f. B=%f ---> ellipse2tr(x1)= %f\n" x1 ychk.[i] a1 b1 (ellipse2tr x1 ychk.[i] aa bb cc dd ee ff)
- printf "\tx2=%f, y1=%f, A=%f. B=%f ---> ellipse2tr(x2)= %f\n" x2 ychk.[i] a1 b1 (ellipse2tr x2 ychk.[i] aa bb cc dd ee ff)
+ printf "\n\tx1=%f, y1=%f, A=%f. B=%f ---> ellipse2tr (x1)= %f\n" x1 ychk.[i] a1 b1 (ellipse2tr x1 ychk.[i] aa bb cc dd ee ff)
+ printf "\tx2=%f, y1=%f, A=%f. B=%f ---> ellipse2tr (x2)= %f\n" x2 ychk.[i] a1 b1 (ellipse2tr x2 ychk.[i] aa bb cc dd ee ff)
#endif
if abs (ellipse2tr x1 ychk.[i] aa bb cc dd ee ff) < EPS then