Herbie run

Date:Monday, May 12th, 2025
Commit:2ad78556 on artem-rules-updates
Seed:2025132
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:467 417.6 MB

Time bar (total: 6.2min)

sample3.5min (55.9%)

Memory
5 873.0MiB live, 263 335.2MiB allocated; 1.6min collecting garbage
Samples
2.1min14 043×3exit
40.7s268 255×0valid
6.6s28 686×1valid
1.4s11 164×0invalid
736.0ms2 490×2valid
277.0ms1 578×1exit
209.0ms694×1invalid
32.0ms5exit
25.0ms4exit
0.0ms3valid
Precisions
Click to see histograms. Total time spent on operations: 2.6min
ival-mult!: 25.4s (16.2% of total)
adjust: 23.5s (15.1% of total)
ival-exp: 19.9s (12.7% of total)
ival-pow2: 18.3s (11.7% of total)
ival-log: 16.7s (10.7% of total)
ival-div!: 9.8s (6.2% of total)
ival-neg: 9.3s (5.9% of total)
ival-add!: 9.0s (5.7% of total)
ival-sub!: 5.9s (3.8% of total)
ival-sqrt: 5.9s (3.8% of total)
ival-fmax: 2.9s (1.8% of total)
ival-sin: 1.5s (0.9% of total)
ival-fabs: 1.4s (0.9% of total)
ival-fmin: 1.1s (0.7% of total)
ival-<=: 983.0ms (0.6% of total)
ival-cos: 931.0ms (0.6% of total)
ival-floor: 756.0ms (0.5% of total)
ival-sinu: 519.0ms (0.3% of total)
ival-cosu: 410.0ms (0.3% of total)
ival-expm1: 390.0ms (0.2% of total)
ival-tan: 355.0ms (0.2% of total)
ival-and: 297.0ms (0.2% of total)
ival-if: 247.0ms (0.2% of total)
ival-log2: 200.0ms (0.1% of total)
ival-asin: 198.0ms (0.1% of total)
ival-log1p: 142.0ms (0.1% of total)
ival-atan: 132.0ms (0.1% of total)
ival-sinh: 90.0ms (0.1% of total)
ival->=: 71.0ms (0% of total)
ival->: 44.0ms (0% of total)
ival-<: 42.0ms (0% of total)
ival-assert: 31.0ms (0% of total)
const: 2.0ms (0% of total)
ival-pi: 0.0ms (0% of total)
Bogosity

rewrite39.9s (10.8%)

Memory
883.3MiB live, 39 367.4MiB allocated; 12.2s collecting garbage
Stop Event
275×iter-limit
118×node-limit
36×saturated
26×unsound
Counts
27 752 → 64 773

analyze28.7s (7.7%)

Memory
486.0MiB live, 39 615.3MiB allocated; 9.6s 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)

series23.6s (6.4%)

Memory
1 052.9MiB live, 33 336.0MiB allocated; 7.8s collecting garbage
Counts
6 343 → 21 409
Calls

1515 calls:

TimeVariablePointExpression
875.0ms
dX.u
@0
((log2 (if (> (/ (fmax (+ (* (* (* (floor h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (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 h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (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 h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))))) (if (> (/ (fmax (+ (* (* (* (floor h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (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 h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (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 h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))))) (> (/ (fmax (+ (* (* (* (floor h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (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 h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (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 h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))) (+ (* (* (* (floor h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (* (* (floor h) (floor h)) dX.v) (* (floor h) (floor h)) (floor h) h dX.v (* (* (floor w) (floor w)) (* dX.u dX.u)) (* (floor w) (floor w)) (floor w) w (* dX.u dX.u) dX.u (+ (* (* (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 h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v))))) (floor maxAniso)) (sqrt (fmax (+ (* (* (* (floor h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (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 h) (floor h)) dX.v) dX.v) (* (* (floor w) (floor w)) (* dX.u dX.u))) (+ (* (* (floor w) dY.u) (* (floor w) dY.u)) (* (* (floor h) dY.v) (* (floor h) dY.v)))))) (log2 (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)) (/ (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))) (* (* (neg dX.v) dY.u) (* (floor h) (floor w))) (* (neg dX.v) dY.u) (neg dX.v) (* (floor h) (floor w)) (/ (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)))))) (log2 (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w))))) (fabs (* (* (floor h) (floor w)) (- (* dY.u dX.v) (* dY.v dX.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w)))))) (floor maxAniso)) (/ (fabs (* (* (floor h) (floor w)) (- (* dY.u dX.v) (* dY.v dX.u)))) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w))))))))) (if (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w))))) (fabs (* (* (floor h) (floor w)) (- (* dY.u dX.v) (* dY.v dX.u))))) (floor maxAniso)) (/ (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w)))))) (floor maxAniso)) (/ (fabs (* (* (floor h) (floor w)) (- (* dY.u dX.v) (* dY.v dX.u)))) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w)))))))) (> (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w))))) (fabs (* (* (floor h) (floor w)) (- (* dY.u dX.v) (* dY.v dX.u))))) (floor maxAniso)) (/ (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w))))) (fabs (* (* (floor h) (floor w)) (- (* dY.u dX.v) (* dY.v dX.u))))) (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w))))) (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (* (* (floor w) (floor w)) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v)) (* dX.v dX.v) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w)))) (* (* (floor h) (floor h)) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w))) (* dY.u dY.u) (fabs (* (* (floor h) (floor w)) (- (* dY.u dX.v) (* dY.v dX.u)))) (* (* (floor h) (floor w)) (- (* dY.u dX.v) (* dY.v dX.u))) (- (* dY.u dX.v) (* dY.v dX.u)) (* dY.u dX.v) (* dY.v dX.u) (/ (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w)))))) (floor maxAniso)) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w)))))) (/ (fabs (* (* (floor h) (floor w)) (- (* dY.u dX.v) (* dY.v dX.u)))) (sqrt (fmax (+ (* (* (* (floor w) (floor w)) dX.u) dX.u) (* (* (floor h) (floor h)) (* dX.v dX.v))) (+ (* (* (* (floor h) (floor h)) dY.v) dY.v) (* (* dY.u dY.u) (* (floor w) (floor w))))))) (log2 (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h)))))) (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))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h))))))) (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))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h)))))))))) (if (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h)))))) (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))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h))))))) (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))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h))))))))) (> (/ (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h)))))) (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))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h)))))) (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))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h)))))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h))))) (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (pow (* (floor w) dY.u) 4) 4 (pow (* (floor h) dY.v) 4) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h)))) (* (* dY.v dY.v) (* (floor h) (floor h))) (* dY.v dY.v) (/ (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h))))))) (floor maxAniso)) (sqrt (fmax (+ (* (* (floor w) dX.u) (* (floor w) dX.u)) (* (* (floor h) dX.v) (* (floor h) dX.v))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h))))))) (/ (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))) (/ (- (pow (* (floor w) dY.u) 4) (pow (* (floor h) dY.v) 4)) (- (* (* dY.u dY.u) (* (floor w) (floor w))) (* (* dY.v dY.v) (* (floor h) (floor h)))))))) (log2 (if (> (/ (fmax (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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 (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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 (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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 (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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 (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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 (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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 (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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 (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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 (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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)))) (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (floor w) (floor w)) (* dX.u dX.u)) (* (* (floor h) (floor h)) (* dX.v dX.v)))) (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4) (- (* (* (floor w) (floor w)) (* dX.u dX.u)) (* (* (floor h) (floor h)) (* dX.v dX.v))) (/ (sqrt (fmax (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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 (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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 (/ (- (pow (* (floor w) dX.u) 4) (pow (* (floor h) dX.v) 4)) (- (* (* (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)))))))
301.0ms
alphay
@inf
((/ 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 (/ u0 (* (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay))) (- 1 u0)))))) (/ 1 (+ 1 (/ u0 (* (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay))) (- 1 u0))))) 1 (+ 1 (/ u0 (* (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay))) (- 1 u0)))) (/ u0 (* (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay))) (- 1 u0))) u0 (* (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay))) (- 1 u0)) (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay))) (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ 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 (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (* alphay alphay)) (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))) (* 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))))) (- 1 (/ (* 1/2 (* (* alphax alphax) u0)) (* (pow (cos (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (- 1 u0)))) (/ (* 1/2 (* (* alphax alphax) u0)) (* (pow (cos (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 2) (- 1 u0))) (* 1/2 (* (* alphax alphax) u0)) (* (* alphax alphax) u0) (* (pow (cos (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) 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)))) (+ 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 (+ (/ (* (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 (+ (/ (* (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 (+ (/ (* (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)))) (* (* alphay alphay) (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) -2)) (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) -2) (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))) (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))) (* 1/2 (PI)) -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 (* (* (* alphay alphay) (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) -2)) (/ u0 (- 1 u0)))))) (/ 1 (+ 1 (* (* (* alphay alphay) (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) -2)) (/ u0 (- 1 u0))))) (+ 1 (* (* (* alphay alphay) (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) -2)) (/ u0 (- 1 u0)))) (* (* (* alphay alphay) (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) -2)) (/ u0 (- 1 u0))) (* (* alphay alphay) (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) -2)) (pow (sin (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))) -2) (/ 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 (/ u0 (* (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (- 1/2 (* 1/2 (cos (* 2 (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))))))) (* alphay alphay))) (- 1 u0)))))) (/ 1 (+ 1 (/ u0 (* (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (- 1/2 (* 1/2 (cos (* 2 (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))))))) (* alphay alphay))) (- 1 u0))))) (+ 1 (/ u0 (* (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (- 1/2 (* 1/2 (cos (* 2 (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))))))) (* alphay alphay))) (- 1 u0)))) (/ u0 (* (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (- 1/2 (* 1/2 (cos (* 2 (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))))))) (* alphay alphay))) (- 1 u0))) (* (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (- 1/2 (* 1/2 (cos (* 2 (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))))))) (* alphay alphay))) (- 1 u0)) (+ (/ (pow (/ 1 (sqrt (+ (* 1 1) (* (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))) (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) 2) (* alphax alphax)) (/ (- 1/2 (* 1/2 (cos (* 2 (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))))))) (* alphay alphay))) (/ (- 1/2 (* 1/2 (cos (* 2 (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))))))) (* alphay alphay)) (- 1/2 (* 1/2 (cos (* 2 (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))))))) (* 1/2 (cos (* 2 (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))) (cos (* 2 (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI))))))))) (* 2 (atan (* (/ alphay alphax) (tan (+ (* 1/2 (PI)) (* 2 (* u1 (PI)))))))))
283.0ms
u
@-inf
((* (* s 3) (neg (log (+ 1 (neg (+ (* u 4/3) -1/3)))))) (* s 3) s 3 (neg (log (+ 1 (neg (+ (* u 4/3) -1/3))))) (log (+ 1 (neg (+ (* u 4/3) -1/3)))) (neg (+ (* u 4/3) -1/3)) (+ (* -4/3 u) 1/3) -4/3 u 1/3 (* (* 3 s) (log (/ 1 (- 1 (/ (- u 1/4) 3/4))))) (* 3 s) (log (/ 1 (- 1 (/ (- u 1/4) 3/4)))) (/ 1 (- 1 (/ (- u 1/4) 3/4))) 3/4 (* s (* 3 (neg (log (- 1 (+ (* u 4/3) -1/3)))))) (* 3 (neg (log (- 1 (+ (* u 4/3) -1/3))))) (neg (log (- 1 (+ (* u 4/3) -1/3)))) (log (- 1 (+ (* u 4/3) -1/3))) (- 1 (+ (* u 4/3) -1/3)) 1 (+ (* u 4/3) -1/3) 4/3 -1/3 (* (* s 3) (log (/ 1 (- 1 (+ (* u 4/3) -1/3))))) (log (/ 1 (- 1 (+ (* u 4/3) -1/3)))) (/ 1 (- 1 (+ (* u 4/3) -1/3))) (* (* s 3) (neg (log (+ 1 (neg (/ (- (* (* 4/3 u) (* 4/3 u)) 1/9) (+ (* 4/3 u) 1/3))))))) (neg (log (+ 1 (neg (/ (- (* (* 4/3 u) (* 4/3 u)) 1/9) (+ (* 4/3 u) 1/3)))))) (log (+ 1 (neg (/ (- (* (* 4/3 u) (* 4/3 u)) 1/9) (+ (* 4/3 u) 1/3))))) (neg (/ (- (* (* 4/3 u) (* 4/3 u)) 1/9) (+ (* 4/3 u) 1/3))) (/ (- (* (* 4/3 u) (* 4/3 u)) 1/9) (+ (* 4/3 u) 1/3)) (- (* (* 4/3 u) (* 4/3 u)) 1/9) (* (* 4/3 u) (* 4/3 u)) (* 4/3 u) 1/9 (+ (* 4/3 u) 1/3))
258.0ms
s
@-inf
((* (neg s) (log (- (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1))) (neg s) s (log (- (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1)) (- (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1) (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s))))) (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (+ 1 (exp (/ (neg (PI)) s))) (exp (/ (neg (PI)) s)) (/ (neg (PI)) s) (neg (PI)) (PI) (/ 1 (+ 1 (exp (/ (PI) s)))) (+ 1 (exp (/ (PI) s))) (exp (/ (PI) s)) (/ (PI) s) u (* (neg s) (log (/ (- (* (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s))))))) 1) (+ (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1)))) (log (/ (- (* (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s))))))) 1) (+ (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1))) 0 (* (neg s) (log (- (/ 1 (+ (* u (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s)))))) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1))) (log (- (/ 1 (+ (* u (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s)))))) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1)) (* (/ (+ (* (* (PI) 1/2) u) (* -1/4 (PI))) s) -4) (/ (+ (* (* (PI) 1/2) u) (* -1/4 (PI))) s) (+ (* (* (PI) 1/2) u) (* -1/4 (PI))) (* u (+ (* -1/4 (/ (PI) u)) (* 1/2 (PI)))) (+ (* -1/4 (/ (PI) u)) (* 1/2 (PI))) -1/4 (/ (PI) u) (* 1/2 (PI)) 1/2 -4 (* (neg s) (log (- (/ 1 (+ (* u (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s)))))) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1))) (+ (* (neg s) 0) (/ s (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u))) (/ s (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u)) (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (* (neg s) (log (- (/ (pow (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 2) (+ (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1)) (/ 1 (+ (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1))))) (log (- (/ (pow (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 2) (+ (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1)) (/ 1 (+ (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1)))) (- (/ (pow (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 2) (+ (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1)) (/ 1 (+ (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1))) (/ (pow (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 2) (+ (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1)) (pow (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 2) 2 (+ (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1) (/ 1 (+ (/ 1 (+ (* (- (/ 1 (+ 1 (exp (/ (neg (PI)) s)))) (/ 1 (+ 1 (exp (/ (PI) s))))) u) (/ 1 (+ 1 (exp (/ (PI) s)))))) 1)))
237.0ms
s
@0
((* (* s 3) (neg (log (+ 1 (neg (+ (* u 4/3) -1/3)))))) (* s 3) s 3 (neg (log (+ 1 (neg (+ (* u 4/3) -1/3))))) (log (+ 1 (neg (+ (* u 4/3) -1/3)))) (neg (+ (* u 4/3) -1/3)) (+ (* u 4/3) -1/3) u 4/3 -1/3 (* (* 3 s) (log (/ 1 (- 1 (/ (- u 1/4) 3/4))))) (* (* (log 3/4) s) 3) (* (log 3/4) s) (log 3/4) 3/4 (* (* s 3) (neg (log (- 1 (+ (* 4/3 u) -1/3))))) (neg (log (- 1 (+ (* 4/3 u) -1/3)))) (log (- 1 (+ (* 4/3 u) -1/3))) (- 1 (+ (* 4/3 u) -1/3)) (+ (* -4/3 u) 4/3) -4/3 (* (* s 3) (neg (log (- 1 (- (* u 4/3) 1/3))))) (neg (log (- 1 (- (* u 4/3) 1/3)))) (log (- 1 (- (* u 4/3) 1/3))) (- 1 (- (* u 4/3) 1/3)) 1 (- (* u 4/3) 1/3) (* u 4/3) 1/3 (* (* 3 s) (log (/ 1 (- 1 (/ (- u 1/4) 3/4))))) (* 3 s) (log (/ 1 (- 1 (/ (- u 1/4) 3/4)))) (/ 1 (- 1 (/ (- u 1/4) 3/4))) (- 1 (/ (- u 1/4) 3/4)) (* (+ (* (/ 1 u) 4/3) -4/3) u) (+ (* (/ 1 u) 4/3) -4/3) (/ 1 u))

regimes21.9s (5.9%)

Memory
432.1MiB live, 34 068.4MiB allocated; 7.6s collecting garbage
Counts
9 025 → 481
Calls

167 calls:

1.6s
(+.f32 (/.f32 (*.f32 #s(literal 1/4 binary32) (exp.f32 (/.f32 (neg.f32 r) s))) (*.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) s) r)) (/.f32 (*.f32 #s(literal 3/4 binary32) (exp.f32 (/.f32 (neg.f32 r) (*.f32 #s(literal 3 binary32) s)))) (*.f32 (*.f32 (*.f32 #s(literal 6 binary32) (PI.f32)) s) r)))
739.0ms
(-.f32 #s(literal 1 binary32) u0)
655.0ms
x
570.0ms
s
554.0ms
(*.f32 (/.f32 (sin.f32 (*.f32 (*.f32 x (PI.f32)) tau)) (*.f32 (*.f32 x (PI.f32)) tau)) (/.f32 (sin.f32 (*.f32 x (PI.f32))) (*.f32 x (PI.f32))))
Compiler

Compiled 20 381 to 18 660 computations (8.4% saved)

derivations16.1s (4.3%)

Memory
-247.0MiB live, 13 039.7MiB allocated; 3.7s collecting garbage
Stop Event
33×fuel
done
Compiler

Compiled 10 422 to 4 079 computations (60.9% saved)

eval14.6s (3.9%)

Memory
241.6MiB live, 20 052.4MiB allocated; 8.2s collecting garbage
Compiler

Compiled 2 989 902 to 472 361 computations (84.2% saved)

preprocess8.8s (2.4%)

Memory
-435.7MiB live, 10 213.5MiB allocated; 5.4s collecting garbage
Stop Event
36×node-limit
Compiler

Compiled 72 893 to 43 355 computations (40.5% saved)

prune8.6s (2.3%)

Memory
-432.8MiB live, 12 736.3MiB allocated; 4.6s collecting garbage
Counts
66 936 → 3 422
Compiler

Compiled 309 530 to 190 585 computations (38.4% saved)

bsearch1.1s (0.3%)

Memory
191.0MiB live, 1 601.2MiB allocated; 135ms collecting garbage
Algorithm
93×binary-search
47×left-value
Stop Event
92×narrow-enough
predicate-same
Samples
446.0ms2 546×0valid
226.0ms803×1valid
19.0ms193×0invalid
17.0ms95×1invalid
15.0ms59×2valid
Compiler

Compiled 32 423 to 26 403 computations (18.6% saved)

Precisions
Click to see histograms. Total time spent on operations: 603.0ms
ival-log1p: 90.0ms (14.9% of total)
ival-pow2: 71.0ms (11.8% of total)
ival-exp: 70.0ms (11.6% of total)
ival-fmax: 63.0ms (10.4% of total)
ival-sin: 51.0ms (8.5% of total)
ival-sub!: 47.0ms (7.8% of total)
ival-cos: 46.0ms (7.6% of total)
ival-mult!: 43.0ms (7.1% of total)
adjust: 29.0ms (4.8% of total)
ival-add!: 27.0ms (4.5% of total)
ival-sqrt: 24.0ms (4% of total)
ival-fmin: 22.0ms (3.6% of total)
ival-fabs: 15.0ms (2.5% of total)
ival-neg: 8.0ms (1.3% of total)

start23.0ms (0%)

Memory
6.2MiB live, 51.3MiB allocated; 5ms collecting garbage

end0.0ms (0%)

Memory
0.6MiB live, 0.6MiB allocated; 0ms collecting garbage

Profiling

Loading profile data...