Herbie run

Date:Saturday, May 10th, 2025
Commit:363acd86 on unsound-last
Seed:2025130
Parameters:256 points for 4 iterations
Flags:
reduce:regimesreduce:binary-searchreduce:branch-expressionssetup:searchrules:arithmeticrules:polynomialsrules:fractionsrules:exponentsrules:trigonometryrules:hyperbolicrules:numericsrules:specialrules:boolsrules:branchesgenerate:rrgenerate:taylorgenerate:proofs
default
Memory:634 069.5 MB

Time bar (total: 12.4min)

sample3.6min (28.9%)

Memory
4 427.8MiB live, 201 196.2MiB allocated; 1.5min collecting garbage
Samples
1.9min11 400×3exit
51.2s267 836×0valid
5.9s28 674×1valid
2.2s11 288×0invalid
691.0ms2 501×2valid
185.0ms1 645×1exit
150.0ms676×1invalid
49.0ms2×5exit
26.0ms2×4exit
0.0ms1×4valid
0.0ms1×3valid
Precisions
Click to see histograms. Total time spent on operations: 2.6min
ival-mult!: 25.7s (16.5% of total)
ival-exp: 21.1s (13.6% of total)
adjust: 19.6s (12.6% of total)
ival-pow2: 19.3s (12.4% of total)
ival-log: 15.7s (10.1% of total)
ival-div!: 10.6s (6.8% of total)
ival-add!: 7.7s (5% of total)
ival-sqrt: 6.7s (4.3% of total)
ival-sub!: 6.5s (4.2% of total)
ival-neg: 5.9s (3.8% of total)
ival-fmax: 2.7s (1.7% of total)
ival-floor: 2.4s (1.6% of total)
ival-cos: 2.4s (1.5% of total)
ival-fabs: 2.0s (1.3% of total)
ival-sin: 1.4s (0.9% of total)
ival-if: 1.1s (0.7% of total)
ival-fmin: 881.0ms (0.6% of total)
ival-sinu: 724.0ms (0.5% of total)
ival-<=: 513.0ms (0.3% of total)
ival-expm1: 416.0ms (0.3% of total)
ival-cosu: 387.0ms (0.2% of total)
ival-tan: 365.0ms (0.2% of total)
ival->: 326.0ms (0.2% of total)
ival-and: 293.0ms (0.2% of total)
ival-log2: 195.0ms (0.1% of total)
ival-log1p: 170.0ms (0.1% of total)
ival->=: 142.0ms (0.1% of total)
ival-sinh: 115.0ms (0.1% of total)
ival-asin: 102.0ms (0.1% of total)
ival-atan: 80.0ms (0.1% of total)
ival-assert: 29.0ms (0% of total)
ival-<: 23.0ms (0% of total)
const: 3.0ms (0% of total)
ival-pi: 0.0ms (0% of total)
Bogosity

series2.8min (22.7%)

Memory
1 392.3MiB live, 152 550.9MiB allocated; 1.1min collecting garbage
Counts
5 940 → 27 146
Calls

1500 calls:

TimeVariablePointExpression
1.0s
cos2phi
@0
((/ (neg (log (+ 1 (neg u0)))) (+ (/ cos2phi (* alphax alphax)) (/ sin2phi (* alphay alphay)))) (neg (log (+ 1 (neg u0)))) (log (+ 1 (neg u0))) (neg u0) u0 (+ (/ cos2phi (* alphax alphax)) (/ sin2phi (* alphay alphay))) (/ (+ (* (* alphay alphay) (/ cos2phi (* alphax alphax))) sin2phi) (* alphay alphay)) (+ (* (* alphay alphay) (/ cos2phi (* alphax alphax))) sin2phi) (* alphay alphay) alphay (/ cos2phi (* alphax alphax)) cos2phi (* alphax alphax) alphax sin2phi (/ (neg (log (- 1 u0))) (+ (/ cos2phi (* alphax alphax)) (/ sin2phi (* alphay alphay)))) (* (* alphax alphax) (/ (neg (log (- 1 u0))) cos2phi)) (/ (* alphax (* alphax u0)) cos2phi) (* alphax (* alphax u0)) (* alphax u0) (/ (neg (log (- 1 u0))) (+ (/ cos2phi (* alphax alphax)) (/ sin2phi (* alphay alphay)))) (neg (log (- 1 u0))) (log (- 1 u0)) (- 1 u0) 1 (+ (/ cos2phi (* alphax alphax)) (/ sin2phi (* alphay alphay))) (/ (/ sin2phi alphay) alphay) (/ sin2phi alphay) (/ (neg (log (- 1 u0))) (+ (/ cos2phi (* alphax alphax)) (/ sin2phi (* alphay alphay)))) (/ (* (* (* alphax alphax) alphay) u0) (+ (* alphay cos2phi) (/ (* (* alphax alphax) sin2phi) alphay))) (* (* (* alphax alphax) alphay) u0) (* (* alphax alphax) alphay) (+ (* alphay cos2phi) (/ (* (* alphax alphax) sin2phi) alphay)) (/ (* (* alphax alphax) sin2phi) alphay) (* (* alphax alphax) sin2phi) (/ (neg (log (- 1 u0))) (+ (/ cos2phi (* alphax alphax)) (/ sin2phi (* alphay alphay)))) (* (* alphax alphax) (/ (neg (log (- 1 u0))) cos2phi)) (* (+ (* (/ (+ (* 1/2 (* alphax alphax)) (* 1/3 (* (* alphax alphax) u0))) cos2phi) u0) (/ (* alphax alphax) cos2phi)) u0) (+ (* (/ (+ (* 1/2 (* alphax alphax)) (* 1/3 (* (* alphax alphax) u0))) cos2phi) u0) (/ (* alphax alphax) cos2phi)) (/ (+ (* 1/2 (* alphax alphax)) (* 1/3 (* (* alphax alphax) u0))) cos2phi) (+ (* 1/2 (* alphax alphax)) (* 1/3 (* (* alphax alphax) u0))) 1/2 (* 1/3 (* (* alphax alphax) u0)) 1/3 (* (* alphax alphax) u0) (/ (* alphax alphax) cos2phi))
956.0ms
w
@-inf
((if (< (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) 1) (fmax 1 (* (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))))) (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))))) (< (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) 1) (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (* (* (floor w) (floor w)) dX.u) (* (floor w) (floor w)) (floor w) w dX.u (* (* (floor h) (floor h)) (* dX.v dX.v)) (* (floor h) (floor h)) (floor h) h (* dX.v dX.v) dX.v (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))) (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (floor w) dY.u) dY.u (* (* (floor h) dY.v) (* (floor h) dY.v)) (* (floor h) dY.v) dY.v (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))) (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (floor w) dX.u) (* (* (floor h) dX.v) (* (floor w) dY.u)) (* (floor h) dX.v) (floor maxAniso) maxAniso (/ (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))) 1 (fmax 1 (* (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))))) (* (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))))) (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (< (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) 1) (fmax 1 (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))))) (< (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) 1) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))) (fmax 1 (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))))) (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (< (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) 1) (fmax 1 (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))))) (< (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) 1) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v)) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))) (* (* (floor w) dY.v) (* (floor h) dX.u)) (* (floor w) dY.v) (* (floor h) dX.u) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))) (fmax 1 (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))))) (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (< (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) 1) (fmax 1 (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))))) (< (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) 1) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))) (* (* (* (floor h) (floor w)) (- (/ (* dY.v dX.u) dY.u) dX.v)) dY.u) (* (* (floor h) (floor w)) (- (/ (* dY.v dX.u) dY.u) dX.v)) (* (floor h) (floor w)) (- (/ (* dY.v dX.u) dY.u) dX.v) (/ (* dY.v dX.u) dY.u) (* dY.v dX.u) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))) (fmax 1 (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))))) (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u))))) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (fabs (- (* (* (floor w) dX.u) (* (floor h) dY.v)) (* (* (floor h) dX.v) (* (floor w) dY.u)))))) (if (< (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) 1) (fmax 1 (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)))) (< (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) 1) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1) (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) (- (* dY.u dX.v) (* dY.v dX.u)) (* dY.u dX.v) (/ (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))) (fmax 1 (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))))) (* (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (/ (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1)) (floor maxAniso)) (floor maxAniso) (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (pow (* (- (* dY.u dX.v) (* dY.v dX.u)) (* (floor h) (floor w))) 1))))
692.0ms
sinTheta_i
@-inf
((/ (* (exp (neg (/ (* sinTheta_i sinTheta_O) v))) (/ (* cosTheta_i cosTheta_O) v)) (* (* (sinh (/ 1 v)) 2) v)) (* (exp (neg (/ (* sinTheta_i sinTheta_O) v))) (/ (* cosTheta_i cosTheta_O) v)) (exp (neg (/ (* sinTheta_i sinTheta_O) v))) (neg (/ (* sinTheta_i sinTheta_O) v)) (/ (* sinTheta_i sinTheta_O) v) (* sinTheta_i sinTheta_O) sinTheta_i sinTheta_O v (/ (* cosTheta_i cosTheta_O) v) (* cosTheta_i cosTheta_O) cosTheta_i cosTheta_O (* (* (sinh (/ 1 v)) 2) v) (* (sinh (/ 1 v)) 2) (sinh (/ 1 v)) (/ 1 v) 1 2)
412.0ms
u1
@0
((sqrt (/ 1 (+ (/ u0 (* (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (- 1 u0))) 1))) (/ 1 (+ (/ u0 (* (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (- 1 u0))) 1)) 1 (+ (/ u0 (* (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (- 1 u0))) 1) (/ u0 (* (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (- 1 u0))) u0 (* (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (- 1 u0)) (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2))) (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2) (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) (tan (* (PI) (+ (* 2 u1) 1/2))) (* (PI) (+ (* 2 u1) 1/2)) (PI) (+ (* 2 u1) 1/2) 2 u1 1/2 (/ alphay alphax) alphay alphax (* alphay alphay) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax))) (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)) (* alphax alphax) (- 1 u0) (/ 1 (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))))) (/ 1 (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))))) (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0)))) (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))) (+ (* (* alphax alphax) (/ (* u0 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2))) (- 1 u0))) 1) (/ (* u0 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2))) (- 1 u0)) (* u0 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2))) (/ 1 (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))))) (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0)))) (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))) (+ (* (* alphay alphay) (/ u0 (* (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (- 1 u0)))) 1) (/ u0 (* (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (- 1 u0))) (* (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (- 1 u0)) (/ 1 (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))))) (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0)))) (+ (* (/ u0 (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax))))) 1/2) 1) (/ u0 (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax))))))
370.0ms
dY.w
@inf
((log2 (exp (* (log (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (+ (* (* (floor h) dX.v) (* (floor h) dX.v)) (* (* (floor d) dX.w) (* (floor d) dX.w)))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w))))) 1/2))) (exp (* (log (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (+ (* (* (floor h) dX.v) (* (floor h) dX.v)) (* (* (floor d) dX.w) (* (floor d) dX.w)))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w))))) 1/2)) (* (log (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (+ (* (* (floor h) dX.v) (* (floor h) dX.v)) (* (* (floor d) dX.w) (* (floor d) dX.w)))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w))))) 1/2) (log (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (+ (* (* (floor h) dX.v) (* (floor h) dX.v)) (* (* (floor d) dX.w) (* (floor d) dX.w)))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w))))) (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (+ (* (* (floor h) dX.v) (* (floor h) dX.v)) (* (* (floor d) dX.w) (* (floor d) dX.w)))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w)))) (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (+ (* (* (floor h) dX.v) (* (floor h) dX.v)) (* (* (floor d) dX.w) (* (floor d) dX.w)))) (* (floor w) dX.u) (floor w) w dX.u (+ (* (* (floor h) dX.v) (* (floor h) dX.v)) (* (* (floor d) dX.w) (* (floor d) dX.w))) (* (floor h) dX.v) (floor h) h dX.v (* (* (floor d) dX.w) (* (floor d) dX.w)) (* (floor d) dX.w) (floor d) d dX.w (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w))) (* (* dY.w (floor d)) (* dY.w (floor d))) (* dY.w (floor d)) dY.w 1/2 (log2 (sqrt (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))) (* (* (floor d) dY.w) (* (floor d) dY.w)))))) (sqrt (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))) (* (* (floor d) dY.w) (* (floor d) dY.w))))) (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))) (* (* (floor d) dY.w) (* (floor d) dY.w)))) (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))) (* (* (floor d) dY.w) (* (floor d) dY.w))) (* (* dY.u dY.u) (* (floor w) (floor w))) (* dY.u dY.u) dY.u (* (floor w) (floor w)) (log2 (sqrt (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))) (* (* (floor d) dY.w) (* (floor d) dY.w)))))) (sqrt (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))) (* (* (floor d) dY.w) (* (floor d) dY.w))))) (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))) (* (* (floor d) dY.w) (* (floor d) dY.w)))) (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))) (* (* (floor d) dY.w) (* (floor d) dY.w))) (* (* dY.w dY.w) (* (floor d) (floor d))) (* dY.w dY.w) (* (floor d) (floor d)) (log2 (exp (* (log (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w))))) 1/2))) (exp (* (log (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w))))) 1/2)) (* (log (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w))))) 1/2) (log (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w))))) (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w)))) (+ (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (exp (* (log (* (floor h) dY.v)) 2))) (* (* (floor d) dY.w) (* (floor d) dY.w))) (+ (* (* dY.w (floor d)) (* dY.w (floor d))) (* (* dY.u (floor w)) (* dY.u (floor w)))) (* (* dY.u (floor w)) (* dY.u (floor w))) (* dY.u (floor w)) (log2 (sqrt (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (* (* (* (floor w) (floor w)) dY.u) dY.u) (+ (* 2 (* (/ (- (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2) (/ (+ (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2))) (+ (* (* dY.w (floor d)) (* dY.w (floor d))) (cosh (* (log (* dY.v (floor h))) 2)))))))) (sqrt (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (* (* (* (floor w) (floor w)) dY.u) dY.u) (+ (* 2 (* (/ (- (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2) (/ (+ (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2))) (+ (* (* dY.w (floor d)) (* dY.w (floor d))) (cosh (* (log (* dY.v (floor h))) 2))))))) (fmax (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (+ (* (* (* (floor w) (floor w)) dY.u) dY.u) (+ (* 2 (* (/ (- (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2) (/ (+ (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2))) (+ (* (* dY.w (floor d)) (* dY.w (floor d))) (cosh (* (log (* dY.v (floor h))) 2)))))) (+ (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (* (* (floor d) dX.w) (* (floor d) dX.w))) (* (* (floor w) (floor w)) (* dX.u dX.u)) (* dX.u dX.u) (+ (* (* (* (floor w) (floor w)) dY.u) dY.u) (+ (* 2 (* (/ (- (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2) (/ (+ (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2))) (+ (* (* dY.w (floor d)) (* dY.w (floor d))) (cosh (* (log (* dY.v (floor h))) 2))))) (* (* (floor w) (floor w)) dY.u) (+ (* 2 (* (/ (- (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2) (/ (+ (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2))) (+ (* (* dY.w (floor d)) (* dY.w (floor d))) (cosh (* (log (* dY.v (floor h))) 2)))) 2 (* (/ (- (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2) (/ (+ (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2)) (/ (- (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2) (- (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) (* dY.v (floor h)) dY.v (/ 1 (* dY.v (floor h))) 1 (/ (+ (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) 2) (+ (* dY.v (floor h)) (/ 1 (* dY.v (floor h)))) (+ (* (* dY.w (floor d)) (* dY.w (floor d))) (cosh (* (log (* dY.v (floor h))) 2))) (cosh (* (log (* dY.v (floor h))) 2)) (* (log (* dY.v (floor h))) 2) (log (* dY.v (floor h))))

rewrite1.8min (14.4%)

Memory
-194.7MiB live, 82 255.3MiB allocated; 29.3s collecting garbage
Stop Event
458×iter-limit
136×node-limit
4×unsound
Counts
33 086 → 194 348

derivations1.5min (11.8%)

Memory
-85.6MiB live, 33 116.0MiB allocated; 12.9s collecting garbage
Stop Event
33×fuel
2×done
Compiler

Compiled 13 207 to 4 484 computations (66% saved)

regimes49.1s (6.6%)

Memory
196.7MiB live, 46 502.4MiB allocated; 17.9s collecting garbage
Counts
14 159 → 606
Calls

158 calls:

4.3s
(*.f32 uy #s(literal 2 binary32))
2.3s
dY.u
1.4s
cosTheta_i
1.3s
uy
1.3s
s
Compiler

Compiled 20 641 to 19 739 computations (4.4% saved)

eval47.1s (6.3%)

Memory
1 090.7MiB live, 49 369.5MiB allocated; 16.5s collecting garbage
Compiler

Compiled 7 843 329 to 1 228 607 computations (84.3% saved)

analyze37.9s (5.1%)

Memory
543.6MiB live, 37 247.6MiB allocated; 11.1s collecting garbage
Algorithm
37×search
Search
ProbabilityValidUnknownPreconditionInfiniteDomainCan'tIter
0%0%21%79%0%0%0%0
25.7%5.4%15.6%79%0%0%0%1
26.5%5.6%15.5%79%0%0%0%2
26.5%5.6%15.4%79%0%0%0%3
28.6%6%15%79%0%0%0%4
40.9%8.6%12.4%79%0%0%0%5
46.2%9.7%11.3%79%0%0%0%6
53.2%11%9.7%79%0%0.4%0%7
55.6%11.5%9.2%79%0%0.4%0%8
60.6%12.5%8.1%79%0%0.4%0%9
63.8%13.1%7.4%79%0%0.5%0%10
66.1%13.6%7%79%0%0.5%0%11
67.5%13.9%6.7%79%0%0.5%0%12
Compiler

Compiled 20 293 to 1 873 computations (90.8% saved)

prune18.0s (2.4%)

Memory
-1 966.2MiB live, 19 161.1MiB allocated; 5.6s collecting garbage
Counts
193 481 → 4 494
Compiler

Compiled 404 855 to 225 898 computations (44.2% saved)

preprocess11.4s (1.5%)

Memory
70.0MiB live, 11 190.1MiB allocated; 2.5s collecting garbage
Stop Event
36×node-limit
Compiler

Compiled 107 482 to 55 412 computations (48.4% saved)

bsearch1.4s (0.2%)

Memory
-77.0MiB live, 1 386.7MiB allocated; 224ms collecting garbage
Algorithm
127×binary-search
71×left-value
Stop Event
122×narrow-enough
5×predicate-same
Samples
535.0ms2 514×0valid
193.0ms989×1valid
32.0ms120×1invalid
18.0ms81×2valid
15.0ms160×0invalid
Compiler

Compiled 34 253 to 26 326 computations (23.1% saved)

Precisions
Click to see histograms. Total time spent on operations: 665.0ms
ival-fmax: 157.0ms (23.6% of total)
ival-exp: 119.0ms (17.9% of total)
ival-sin: 72.0ms (10.8% of total)
ival-pow2: 70.0ms (10.5% of total)
ival-mult!: 50.0ms (7.5% of total)
ival-sub!: 49.0ms (7.4% of total)
ival-cos: 36.0ms (5.4% of total)
adjust: 25.0ms (3.8% of total)
ival-sqrt: 23.0ms (3.5% of total)
ival-log1p: 18.0ms (2.7% of total)
ival-add!: 14.0ms (2.1% of total)
ival-fmin: 14.0ms (2.1% of total)
ival-fabs: 10.0ms (1.5% of total)
ival-neg: 8.0ms (1.2% of total)

start34.0ms (0%)

Memory
-4.8MiB live, 93.0MiB allocated; 19ms collecting garbage

end0.0ms (0%)

Memory
0.7MiB live, 0.7MiB allocated; 0ms collecting garbage
Calls

1 calls:

TimeVariablePointExpression
2.3min
u1
@0
((/ 1 (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))))) 1 (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0)))) (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))) (+ 1 (/ u0 (* (+ (/ (- 1 (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2)) (* alphax alphax)) (/ (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay))) (- 1 u0)))) (/ u0 (* (+ (/ (- 1 (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2)) (* alphax alphax)) (/ (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay))) (- 1 u0))) u0 (* (+ (/ (- 1 (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2)) (* alphax alphax)) (/ (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay))) (- 1 u0)) (+ (/ (- 1 (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2)) (* alphax alphax)) (/ (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay))) (/ (- 1 (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2)) (* alphax alphax)) (- 1 (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2)) (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (/ alphay alphax) alphay alphax (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))) (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))) 1/2 (PI) (* 2 (* u1 (PI))) 2 (* u1 (PI)) u1 (* alphax alphax) (/ (pow (tanh (asinh (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay)) (* alphay alphay) (- 1 u0) (/ 1 (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))))) (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0)))) (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))) (+ (* (* alphay alphay) (/ u0 (* (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (- 1 u0)))) 1) (* (* alphay alphay) (/ u0 (- 1 u0))) (/ u0 (- 1 u0)) (sqrt (/ 1 (+ (/ u0 (* (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (- 1 u0))) 1))) (/ 1 (+ (/ u0 (* (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (- 1 u0))) 1)) (+ (/ u0 (* (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (- 1 u0))) 1) (/ u0 (* (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (- 1 u0))) (* (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (- 1 u0)) (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2))) (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax))) (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)) (/ 1 (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))))) (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0)))) (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))) (+ (* (* alphay alphay) (/ u0 (* (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (- 1 u0)))) 1) (/ u0 (* (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (- 1 u0))) (* (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (- 1 u0)) (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2))) (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2) (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) (tan (* (PI) (+ (* 2 u1) 1/2))) (* (PI) (+ (* 2 u1) 1/2)) (+ (* 2 u1) 1/2) (- 1 u0) (/ 1 (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0))))) (sqrt (+ 1 (/ (* (/ 1 (+ (/ (* (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (cos (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphax alphax)) (/ (* (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI))))))) (sin (atan (* (/ alphay alphax) (tan (+ (* (* 2 (PI)) u1) (* 1/2 (PI)))))))) (* alphay alphay)))) u0) (- 1 u0)))) (+ (* (/ u0 (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax))))) 1/2) 1) (/ u0 (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax))))) (+ (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))) (/ (- 1 (/ 1 (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)))) (* alphay alphay)) (/ 1 (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax))) (* (+ 1 (pow (* (tan (* (PI) (+ (* 2 u1) 1/2))) (/ alphay alphax)) 2)) (* alphax alphax)))

Profiling

Loading profile data...