Metrics for bug366, discussion (missed optimization)

start0.0ms (0%)

Memory: 0.1MiB live, 0.1MiB allocated; 0ms collecting garbage

analyze55.0ms (2.7%)

Memory

-2.6MiB live, 41.7MiB allocated; 13ms collecting garbage

Algorithm

1×

Search

Probability	Valid	Unknown	Precondition	Infinite	Domain	Can't	Iter
0%	0%	99.9%	0.1%	0%	0%	0%	0
0%	0%	99.9%	0.1%	0%	0%	0%	1
0%	0%	99.9%	0.1%	0%	0%	0%	2
0%	0%	99.9%	0.1%	0%	0%	0%	3
0%	0%	99.9%	0.1%	0%	0%	0%	4
30.8%	25%	56.2%	0.1%	0%	18.7%	0%	5
32%	25%	53.1%	0.1%	0%	21.9%	0%	6
55.8%	37.5%	29.7%	0.1%	0%	32.8%	0%	7
57.8%	37.5%	27.3%	0.1%	0%	35.1%	0%	8
74.2%	43.7%	15.2%	0.1%	0%	41%	0%	9
75.9%	43.7%	13.9%	0.1%	0%	42.3%	0%	10
85.9%	46.8%	7.7%	0.1%	0%	45.4%	0%	11
87%	46.8%	7%	0.1%	0%	46.1%	0%	12

Compiler

Compiled 10 to 8 computations (20% saved)

sample989.0ms (48.2%)

Memory

23.0MiB live, 949.1MiB allocated; 415ms collecting garbage

Samples

546.0ms	8 256×	0	valid
64.0ms	619×	0	invalid

Precisions

Click to see histograms. Total time spent on operations: 269.0ms

ival-mult!: 150.0ms (55.8% of total)

ival-sqrt: 96.0ms (35.7% of total)

ival-sub!: 21.0ms (7.8% of total)

adjust: 2.0ms (0.7% of total)

Bogosity

preprocess24.0ms (1.2%)

Memory

-12.8MiB live, 25.1MiB allocated; 19ms collecting garbage

Iterations

Useful iterations: 0 (0.0ms)

Iter	Nodes	Cost
0	18	50
1	75	50

Stop Event

1×	node-limit

Alt Table

Click to see full alt table

Status	Accuracy	Program
▶	53.2%	(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))

Symmetry

(abs a)

(abs b)

Compiler

Compiled 16 to 12 computations (25% saved)

series137.0ms (6.7%)

Memory

11.4MiB live, 65.7MiB allocated; 61ms collecting garbage

Counts

6 → 32

Calls

Call 1

Inputs
(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))
(-.f64 (.f64 a a) (.f64 b b))
(*.f64 a a)
a
(*.f64 b b)
b

Outputs
#s(approx (sqrt (- (* a a) (* b b))) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) (*.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) (.f64 (pow.f64 a #s(literal 2 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) #s(literal 3 binary64))) (.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64))))))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) (.f64 (pow.f64 a #s(literal 2 binary64)) (fma.f64 (pow.f64 a #s(literal 2 binary64)) (-.f64 (.f64 #s(literal 1/16 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) #s(literal 5 binary64)))) (.f64 #s(literal 1/8 binary64) (/.f64 #s(literal 1 binary64) (pow.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) #s(literal 3 binary64))))) (.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64))))))))))
#s(approx (- (* a a) (* b b)) (*.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))))
#s(approx (- (* a a) (* b b)) (-.f64 (pow.f64 a #s(literal 2 binary64)) (pow.f64 b #s(literal 2 binary64))))
#s(approx (* a a) (pow.f64 a #s(literal 2 binary64)))
#s(approx a a)
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64))) (.f64 #s(literal -1/16 binary64) (/.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 a #s(literal 6 binary64)))))))))
#s(approx (- (* a a) (* b b)) (.f64 (pow.f64 a #s(literal 2 binary64)) (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))))
#s(approx (sqrt (- (* a a) (* b b))) (*.f64 #s(literal -1 binary64) a))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (*.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64))) (*.f64 #s(literal -1/16 binary64) (/.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 a #s(literal 6 binary64))))))))))
#s(approx (sqrt (- (* a a) (* b b))) (sqrt.f64 (pow.f64 a #s(literal 2 binary64))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 3 binary64)))) (*.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (.f64 #s(literal -1/16 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 5 binary64)))) (.f64 #s(literal 1/8 binary64) (/.f64 #s(literal 1 binary64) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 3 binary64)))))) (*.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))))))))
#s(approx (- (* a a) (* b b)) (fma.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))))
#s(approx (* b b) (pow.f64 b #s(literal 2 binary64)))
#s(approx b b)
#s(approx (sqrt (- (* a a) (* b b))) (*.f64 b (sqrt.f64 #s(literal -1 binary64))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (fma.f64 #s(literal 1/16 binary64) (/.f64 (pow.f64 a #s(literal 6 binary64)) (.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 5 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64))))))))))
#s(approx (- (* a a) (* b b)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (/.f64 (pow.f64 a #s(literal 2 binary64)) (pow.f64 b #s(literal 2 binary64))) #s(literal 1 binary64))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (sqrt.f64 #s(literal -1 binary64)))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64))))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (fma.f64 #s(literal 1/16 binary64) (/.f64 (pow.f64 a #s(literal 6 binary64)) (.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 5 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64)))))))))))

Calls

6 calls:

Time	Variable	Point
15.0ms	a	0
2.0ms	a	inf
2.0ms	b	inf
2.0ms	a	-inf
1.0ms	b	0

rewrite148.0ms (7.2%)

Memory

-11.8MiB live, 81.8MiB allocated; 22ms collecting garbage

Counts

6 → 24

Calls

Call 1

Inputs
(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))
(-.f64 (.f64 a a) (.f64 b b))
(*.f64 a a)
a
(*.f64 b b)
b

Outputs
(fabs.f64 (sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b))))
(pow.f64 (*.f64 (+.f64 a b) (-.f64 a b)) #s(literal 1/2 binary64))
(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
(*.f64 (+.f64 a b) (-.f64 a b))
(+.f64 (.f64 a a) (neg.f64 (.f64 b b)))
(+.f64 (.f64 a a) (.f64 (neg.f64 b) b))
(-.f64 (.f64 a a) (.f64 b b))
(fma.f64 a a (neg.f64 (*.f64 b b)))
(fma.f64 (neg.f64 a) (neg.f64 a) (neg.f64 (*.f64 b b)))
(fma.f64 (fabs.f64 a) (fabs.f64 a) (neg.f64 (*.f64 b b)))
(fma.f64 a a (*.f64 (neg.f64 b) b))
(fma.f64 (neg.f64 a) (neg.f64 a) (*.f64 (neg.f64 b) b))
(fma.f64 (fabs.f64 a) (fabs.f64 a) (*.f64 (neg.f64 b) b))
(neg.f64 (-.f64 (.f64 b b) (.f64 a a)))
(*.f64 a a)
(*.f64 (neg.f64 a) (neg.f64 a))
(*.f64 (fabs.f64 a) (fabs.f64 a))
(pow.f64 a #s(literal 2 binary64))
a
(*.f64 b b)
(*.f64 (neg.f64 b) (neg.f64 b))
(*.f64 (fabs.f64 b) (fabs.f64 b))
(pow.f64 b #s(literal 2 binary64))
b

eval7.0ms (0.4%)

Memory: 8.0MiB live, 8.0MiB allocated; 0ms collecting garbage
Compiler: Compiled 1 060 to 216 computations (79.6% saved)

prune3.0ms (0.1%)

Memory

3.2MiB live, 3.2MiB allocated; 0ms collecting garbage

Pruning

3 alts after pruning (3 fresh and 0 done)

	Pruned	Kept	Total
New	48	3	51
Fresh	0	0	0
Picked	1	0	1
Done	0	0	0
Total	49	3	52

Accuracy

95.5%

Counts

52 → 3

Alt Table

Click to see full alt table

Status	Accuracy	Program
▶	53.5%	(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
▶	90.6%	#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))))
▶	1.0%	#s(approx (sqrt (- (* a a) (* b b))) (*.f64 #s(literal -1 binary64) a))

Status

Accuracy

Program

▶

53.5%

(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))

▶

90.6%

#s(approx (sqrt (- (* a a) (* b b))) (*.f64 a (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))))

▶

1.0%

#s(approx (sqrt (- (* a a) (* b b))) (*.f64 #s(literal -1 binary64) a))

Compiler

Compiled 42 to 31 computations (26.2% saved)

series142.0ms (6.9%)

Memory

-13.1MiB live, 82.7MiB allocated; 23ms collecting garbage

Counts

19 → 81

Calls

Call 1

Inputs
#s(approx (sqrt (- (* a a) (* b b))) (*.f64 #s(literal -1 binary64) a))
(*.f64 #s(literal -1 binary64) a)
#s(literal -1 binary64)
a
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))))
(.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))))))
(+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))
#s(literal 1 binary64)
(*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))))
#s(literal -1/2 binary64)
(/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))
(pow.f64 b #s(literal 2 binary64))
b
#s(literal 2 binary64)
(pow.f64 a #s(literal 2 binary64))
(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
(*.f64 (+.f64 a b) (-.f64 a b))
(+.f64 a b)
(-.f64 a b)

Outputs
#s(approx (sqrt (- (* a a) (* b b))) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) (*.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) (.f64 (pow.f64 a #s(literal 2 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) #s(literal 3 binary64))) (.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64))))))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) (.f64 (pow.f64 a #s(literal 2 binary64)) (fma.f64 (pow.f64 a #s(literal 2 binary64)) (-.f64 (.f64 #s(literal 1/16 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) #s(literal 5 binary64)))) (.f64 #s(literal 1/8 binary64) (/.f64 #s(literal 1 binary64) (pow.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) #s(literal 3 binary64))))) (.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64))))))))))
#s(approx (* -1 a) (*.f64 #s(literal -1 binary64) a))
#s(approx a a)
#s(approx (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) a)))
#s(approx (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))) (/.f64 (fma.f64 #s(literal -1/2 binary64) (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) a))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) (/.f64 (fma.f64 #s(literal -1/2 binary64) (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (pow.f64 a #s(literal 2 binary64))))
#s(approx (/ (pow b 2) (pow a 2)) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))))
#s(approx (pow a 2) (pow.f64 a #s(literal 2 binary64)))
#s(approx (sqrt (* (+ a b) (- a b))) (sqrt.f64 (*.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64)))))
#s(approx (sqrt (* (+ a b) (- a b))) (+.f64 (sqrt.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (.f64 a (+.f64 b (.f64 #s(literal -1 binary64) b))) (sqrt.f64 (*.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))))))))
#s(approx (sqrt (* (+ a b) (- a b))) (+.f64 (sqrt.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64)))) (.f64 a (fma.f64 #s(literal 1/2 binary64) (/.f64 (.f64 a (-.f64 #s(literal 1 binary64) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64)))) #s(literal 2 binary64)))))) (sqrt.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))))) (.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) (sqrt.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))))))))))
#s(approx (sqrt (* (+ a b) (- a b))) (+.f64 (sqrt.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64)))) (.f64 a (fma.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) (sqrt.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))))) (.f64 a (fma.f64 #s(literal -1/4 binary64) (/.f64 (.f64 a (.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) (-.f64 #s(literal 1 binary64) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64)))) #s(literal 2 binary64))))))) (pow.f64 (sqrt.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64)))) #s(literal 3 binary64))) (.f64 #s(literal 1/2 binary64) (/.f64 (-.f64 #s(literal 1 binary64) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64)))) #s(literal 2 binary64))))) (sqrt.f64 (*.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))))))))))))
#s(approx (* (+ a b) (- a b)) (*.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))))
#s(approx (* (+ a b) (- a b)) (fma.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64)) (.f64 a (+.f64 b (.f64 #s(literal -1 binary64) b)))))
#s(approx (* (+ a b) (- a b)) (fma.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64)) (.f64 a (+.f64 a (+.f64 b (.f64 #s(literal -1 binary64) b))))))
#s(approx (+ a b) b)
#s(approx (+ a b) (+.f64 a b))
#s(approx (- a b) (*.f64 #s(literal -1 binary64) b))
#s(approx (- a b) (-.f64 a b))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64))) (.f64 #s(literal -1/16 binary64) (/.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 a #s(literal 6 binary64)))))))))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 b (*.f64 #s(literal -1 binary64) b)) a)))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) a) (.f64 #s(literal 1/2 binary64) (/.f64 (-.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))) (.f64 #s(literal 1/4 binary64) (pow.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) #s(literal 2 binary64)))) (pow.f64 a #s(literal 2 binary64))))))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/4 binary64) (/.f64 (.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) (-.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))) (.f64 #s(literal 1/4 binary64) (pow.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) #s(literal 2 binary64))))) (pow.f64 a #s(literal 3 binary64))) (fma.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) a) (.f64 #s(literal 1/2 binary64) (/.f64 (-.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))) (.f64 #s(literal 1/4 binary64) (pow.f64 (+.f64 b (*.f64 #s(literal -1 binary64) b)) #s(literal 2 binary64)))) (pow.f64 a #s(literal 2 binary64)))))))))
#s(approx (* (+ a b) (- a b)) (*.f64 (pow.f64 a #s(literal 2 binary64)) (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1 binary64) (/.f64 b a) (/.f64 b a)))))
#s(approx (* (+ a b) (- a b)) (*.f64 (pow.f64 a #s(literal 2 binary64)) (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1 binary64) (/.f64 b a) (fma.f64 #s(literal -1 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (/.f64 b a))))))
#s(approx (+ a b) (*.f64 a (+.f64 #s(literal 1 binary64) (/.f64 b a))))
#s(approx (- a b) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1 binary64) (/.f64 b a)))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (*.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64))) (*.f64 #s(literal -1/16 binary64) (/.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 a #s(literal 6 binary64))))))))))
#s(approx (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))) (.f64 #s(literal -1 binary64) (.f64 a (-.f64 (*.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))) #s(literal 1 binary64)))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) a))))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1 binary64) (/.f64 (fma.f64 #s(literal -1/2 binary64) (+.f64 b (.f64 #s(literal -1 binary64) b)) (.f64 #s(literal -1/2 binary64) (/.f64 (-.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))) (.f64 #s(literal 1/4 binary64) (pow.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) #s(literal 2 binary64)))) a))) a))))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1 binary64) (/.f64 (fma.f64 #s(literal -1 binary64) (/.f64 (fma.f64 #s(literal -1/4 binary64) (/.f64 (.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) (-.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))) (.f64 #s(literal 1/4 binary64) (pow.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) #s(literal 2 binary64))))) a) (.f64 #s(literal 1/2 binary64) (-.f64 (.f64 #s(literal -1 binary64) (pow.f64 b #s(literal 2 binary64))) (.f64 #s(literal 1/4 binary64) (pow.f64 (+.f64 b (.f64 #s(literal -1 binary64) b)) #s(literal 2 binary64)))))) a) (.f64 #s(literal -1/2 binary64) (+.f64 b (.f64 #s(literal -1 binary64) b)))) a))))))
#s(approx (* (+ a b) (- a b)) (.f64 (pow.f64 a #s(literal 2 binary64)) (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1 binary64) (/.f64 (fma.f64 #s(literal -1 binary64) (+.f64 b (*.f64 #s(literal -1 binary64) b)) (/.f64 (pow.f64 b #s(literal 2 binary64)) a)) a)))))
#s(approx (+ a b) (.f64 #s(literal -1 binary64) (.f64 a (-.f64 (*.f64 #s(literal -1 binary64) (/.f64 b a)) #s(literal 1 binary64)))))
#s(approx (- a b) (.f64 #s(literal -1 binary64) (.f64 a (-.f64 (/.f64 b a) #s(literal 1 binary64)))))
#s(approx (sqrt (- (* a a) (* b b))) (sqrt.f64 (pow.f64 a #s(literal 2 binary64))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 3 binary64)))) (*.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (.f64 #s(literal -1/16 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 5 binary64)))) (.f64 #s(literal 1/8 binary64) (/.f64 #s(literal 1 binary64) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 3 binary64)))))) (*.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))))))))
#s(approx (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))) (+.f64 a (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) a))))
#s(approx (pow b 2) (pow.f64 b #s(literal 2 binary64)))
#s(approx (sqrt (* (+ a b) (- a b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (.f64 #s(literal 1/2 binary64) (/.f64 (.f64 b (+.f64 a (*.f64 #s(literal -1 binary64) a))) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))))))
#s(approx (sqrt (* (+ a b) (- a b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (.f64 b (fma.f64 #s(literal -1/2 binary64) (/.f64 (.f64 b (+.f64 #s(literal 1 binary64) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 2 binary64)))))) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))))))))
#s(approx (sqrt (* (+ a b) (- a b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (.f64 b (fma.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))) (.f64 b (fma.f64 #s(literal -1/2 binary64) (/.f64 (+.f64 #s(literal 1 binary64) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 2 binary64))))) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))) (.f64 #s(literal 1/4 binary64) (/.f64 (.f64 b (.f64 (+.f64 #s(literal 1 binary64) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 2 binary64))))) (+.f64 a (*.f64 #s(literal -1 binary64) a)))) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 3 binary64))))))))))
#s(approx (* (+ a b) (- a b)) (fma.f64 b (+.f64 a (*.f64 #s(literal -1 binary64) a)) (pow.f64 a #s(literal 2 binary64))))
#s(approx (* (+ a b) (- a b)) (fma.f64 b (+.f64 a (fma.f64 #s(literal -1 binary64) a (*.f64 #s(literal -1 binary64) b))) (pow.f64 a #s(literal 2 binary64))))
#s(approx (- a b) (+.f64 a (*.f64 #s(literal -1 binary64) b)))
#s(approx (sqrt (- (* a a) (* b b))) (*.f64 b (sqrt.f64 #s(literal -1 binary64))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (fma.f64 #s(literal 1/16 binary64) (/.f64 (pow.f64 a #s(literal 6 binary64)) (.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 5 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64))))))))))
#s(approx (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (/.f64 a (pow.f64 b #s(literal 2 binary64))) (.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) a)))))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (/.f64 #s(literal 1 binary64) (pow.f64 b #s(literal 2 binary64))) (.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (pow.f64 a #s(literal 2 binary64)))))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) (.f64 b (sqrt.f64 #s(literal -1 binary64))))))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) (.f64 b (sqrt.f64 #s(literal -1 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (-.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 2 binary64))))) (*.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64)))))))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/4 binary64) (/.f64 (.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) (-.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 2 binary64)))))) (.f64 (pow.f64 b #s(literal 3 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (fma.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) (.f64 b (sqrt.f64 #s(literal -1 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (-.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 2 binary64))))) (.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64))))))))))
#s(approx (* (+ a b) (- a b)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (fma.f64 #s(literal -1 binary64) (/.f64 a b) (/.f64 a b)) #s(literal 1 binary64))))
#s(approx (* (+ a b) (- a b)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (fma.f64 #s(literal -1 binary64) (/.f64 a b) (+.f64 (/.f64 a b) (/.f64 (pow.f64 a #s(literal 2 binary64)) (pow.f64 b #s(literal 2 binary64))))) #s(literal 1 binary64))))
#s(approx (+ a b) (*.f64 b (+.f64 #s(literal 1 binary64) (/.f64 a b))))
#s(approx (- a b) (*.f64 b (-.f64 (/.f64 a b) #s(literal 1 binary64))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (sqrt.f64 #s(literal -1 binary64)))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64))))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (fma.f64 #s(literal 1/16 binary64) (/.f64 (pow.f64 a #s(literal 6 binary64)) (.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 5 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64)))))))))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (.f64 #s(literal -1/2 binary64) (/.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) (*.f64 b (sqrt.f64 #s(literal -1 binary64)))))))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (.f64 #s(literal -1 binary64) (/.f64 (fma.f64 #s(literal -1/2 binary64) (/.f64 (-.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 2 binary64))))) (.f64 b (sqrt.f64 #s(literal -1 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) (sqrt.f64 #s(literal -1 binary64))))) b))))))
#s(approx (sqrt (* (+ a b) (- a b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (.f64 #s(literal -1 binary64) (/.f64 (fma.f64 #s(literal -1 binary64) (/.f64 (fma.f64 #s(literal 1/4 binary64) (/.f64 (.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) (-.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 2 binary64)))))) (.f64 b (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (-.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 #s(literal 1/4 binary64) (/.f64 (pow.f64 (+.f64 a (.f64 #s(literal -1 binary64) a)) #s(literal 2 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 2 binary64))))) (sqrt.f64 #s(literal -1 binary64))))) b) (.f64 #s(literal 1/2 binary64) (/.f64 (+.f64 a (*.f64 #s(literal -1 binary64) a)) (sqrt.f64 #s(literal -1 binary64))))) b))))))
#s(approx (* (+ a b) (- a b)) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (.f64 #s(literal -1 binary64) (/.f64 (+.f64 a (*.f64 #s(literal -1 binary64) a)) b)) #s(literal 1 binary64))))
#s(approx (* (+ a b) (- a b)) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (.f64 #s(literal -1 binary64) (/.f64 (+.f64 a (fma.f64 #s(literal -1 binary64) a (*.f64 #s(literal -1 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) b)))) b)) #s(literal 1 binary64))))
#s(approx (+ a b) (.f64 #s(literal -1 binary64) (.f64 b (-.f64 (*.f64 #s(literal -1 binary64) (/.f64 a b)) #s(literal 1 binary64)))))
#s(approx (- a b) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1 binary64) (/.f64 a b))))))

Calls

6 calls:

Time	Variable	Point
6.0ms	b	0
3.0ms	a	0
3.0ms	b	inf
2.0ms	a	-inf
2.0ms	a	inf

rewrite143.0ms (7%)

Memory

-2.6MiB live, 134.8MiB allocated; 33ms collecting garbage

Counts

19 → 61

Calls

Call 1

Inputs
#s(approx (sqrt (- (* a a) (* b b))) (*.f64 #s(literal -1 binary64) a))
(*.f64 #s(literal -1 binary64) a)
#s(literal -1 binary64)
a
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))))
(.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))))))
(+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))
#s(literal 1 binary64)
(*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))))
#s(literal -1/2 binary64)
(/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))
(pow.f64 b #s(literal 2 binary64))
b
#s(literal 2 binary64)
(pow.f64 a #s(literal 2 binary64))
(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
(*.f64 (+.f64 a b) (-.f64 a b))
(+.f64 a b)
(-.f64 a b)

Outputs
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))
(*.f64 #s(literal -1 binary64) a)
(*.f64 a #s(literal -1 binary64))
(neg.f64 a)
#s(literal -1 binary64)
a
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (fma.f64 (/.f64 (.f64 b b) (*.f64 a a)) #s(literal -1/2 binary64) #s(literal 1 binary64))))
(.f64 a (fma.f64 (/.f64 (.f64 b b) (*.f64 a a)) #s(literal -1/2 binary64) #s(literal 1 binary64)))
(.f64 (fma.f64 (/.f64 (.f64 b b) (*.f64 a a)) #s(literal -1/2 binary64) #s(literal 1 binary64)) a)
(+.f64 (.f64 #s(literal 1 binary64) a) (.f64 (/.f64 (.f64 #s(literal -1/2 binary64) (.f64 b b)) (*.f64 a a)) a))
(+.f64 (.f64 a #s(literal 1 binary64)) (.f64 a (/.f64 (.f64 #s(literal -1/2 binary64) (.f64 b b)) (*.f64 a a))))
(fma.f64 #s(literal 1 binary64) a (.f64 (/.f64 (.f64 #s(literal -1/2 binary64) (.f64 b b)) (.f64 a a)) a))
(fma.f64 a #s(literal 1 binary64) (.f64 a (/.f64 (.f64 #s(literal -1/2 binary64) (.f64 b b)) (.f64 a a))))
(+.f64 #s(literal 1 binary64) (/.f64 (.f64 #s(literal -1/2 binary64) (.f64 b b)) (*.f64 a a)))
(+.f64 (/.f64 (.f64 #s(literal -1/2 binary64) (.f64 b b)) (*.f64 a a)) #s(literal 1 binary64))
(-.f64 #s(literal 1 binary64) (.f64 #s(literal 1/2 binary64) (/.f64 (.f64 b b) (*.f64 a a))))
(-.f64 #s(literal 1 binary64) (neg.f64 (/.f64 (.f64 #s(literal -1/2 binary64) (.f64 b b)) (*.f64 a a))))
(fma.f64 (/.f64 (.f64 b b) (.f64 a a)) #s(literal -1/2 binary64) #s(literal 1 binary64))
(fma.f64 #s(literal -1/2 binary64) (/.f64 (.f64 b b) (.f64 a a)) #s(literal 1 binary64))
#s(literal 1 binary64)
(cosh.f64 #s(literal 0 binary64))
(exp.f64 #s(literal 0 binary64))
(.f64 #s(literal -1/2 binary64) (/.f64 (.f64 b b) (*.f64 a a)))
(.f64 (/.f64 (.f64 b b) (*.f64 a a)) #s(literal -1/2 binary64))
(/.f64 (.f64 #s(literal -1/2 binary64) (.f64 b b)) (*.f64 a a))
#s(literal -1/2 binary64)
(.f64 (.f64 b b) (/.f64 #s(literal 1 binary64) (*.f64 a a)))
(/.f64 (.f64 b b) (.f64 a a))
(/.f64 (neg.f64 (.f64 b b)) (neg.f64 (.f64 a a)))
(*.f64 b b)
(*.f64 (fabs.f64 b) (fabs.f64 b))
(*.f64 (neg.f64 b) (neg.f64 b))
(pow.f64 b #s(literal 2 binary64))
b
#s(literal 2 binary64)
(+.f64 #s(literal 1 binary64) #s(literal 1 binary64))
(*.f64 a a)
(*.f64 (neg.f64 a) (neg.f64 a))
(*.f64 (fabs.f64 a) (fabs.f64 a))
(pow.f64 a #s(literal 2 binary64))
(fabs.f64 (sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b))))
(pow.f64 (*.f64 (+.f64 a b) (-.f64 a b)) #s(literal 1/2 binary64))
(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
(*.f64 (+.f64 a b) (-.f64 a b))
(*.f64 (-.f64 a b) (+.f64 a b))
(+.f64 (.f64 a a) (.f64 (neg.f64 b) b))
(+.f64 (.f64 a a) (neg.f64 (.f64 b b)))
(-.f64 (.f64 a a) (.f64 b b))
(fma.f64 a a (*.f64 (neg.f64 b) b))
(fma.f64 (neg.f64 a) (neg.f64 a) (*.f64 (neg.f64 b) b))
(fma.f64 (fabs.f64 a) (fabs.f64 a) (*.f64 (neg.f64 b) b))
(fma.f64 a a (neg.f64 (*.f64 b b)))
(fma.f64 (neg.f64 a) (neg.f64 a) (neg.f64 (*.f64 b b)))
(fma.f64 (fabs.f64 a) (fabs.f64 a) (neg.f64 (*.f64 b b)))
(neg.f64 (-.f64 (.f64 b b) (.f64 a a)))
(+.f64 a b)
(+.f64 b a)
(-.f64 a (neg.f64 b))
(+.f64 a (neg.f64 b))
(-.f64 a b)
(neg.f64 (-.f64 b a))

eval64.0ms (3.1%)

Memory: 31.4MiB live, 31.5MiB allocated; 0ms collecting garbage
Compiler: Compiled 3 519 to 674 computations (80.8% saved)

prune15.0ms (0.7%)

Memory

-43.4MiB live, 6.2MiB allocated; 18ms collecting garbage

Pruning

3 alts after pruning (2 fresh and 1 done)

	Pruned	Kept	Total
New	128	2	130
Fresh	0	0	0
Picked	2	1	3
Done	0	0	0
Total	130	3	133

Accuracy

99.9%

Counts

133 → 3

Alt Table

Click to see full alt table

Status	Accuracy	Program
✓	53.5%	(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
▶	99.4%	#s(approx (sqrt (- (* a a) (* b b))) (.f64 a #s(approx (+ 1 ( -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))))
▶	1.0%	#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))

Status

Accuracy

Program

✓

53.5%

(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))

▶

99.4%

#s(approx (sqrt (- (* a a) (* b b))) (*.f64 a #s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))))

▶

1.0%

#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))

Compiler

Compiled 43 to 32 computations (25.6% saved)

series138.0ms (6.7%)

Memory

7.4MiB live, 54.0MiB allocated; 6ms collecting garbage

Counts

7 → 34

Calls

Call 1

Inputs
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))
(neg.f64 a)
a
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a #s(approx (+ 1 ( -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))))
(.f64 a #s(approx (+ 1 ( -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64)))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))
#s(literal 1 binary64)

Outputs
#s(approx (sqrt (- (* a a) (* b b))) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) (*.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) (.f64 (pow.f64 a #s(literal 2 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) #s(literal 3 binary64))) (.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64))))))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) (.f64 (pow.f64 a #s(literal 2 binary64)) (fma.f64 (pow.f64 a #s(literal 2 binary64)) (-.f64 (.f64 #s(literal 1/16 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) #s(literal 5 binary64)))) (.f64 #s(literal 1/8 binary64) (/.f64 #s(literal 1 binary64) (pow.f64 (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64)))) #s(literal 3 binary64))))) (.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (neg.f64 (pow.f64 b #s(literal 2 binary64))))))))))
#s(approx (neg a) (*.f64 #s(literal -1 binary64) a))
#s(approx a a)
#s(approx (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) a)))
#s(approx (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))) (/.f64 (fma.f64 #s(literal -1/2 binary64) (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) a))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) (/.f64 (fma.f64 #s(literal -1/2 binary64) (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (pow.f64 a #s(literal 2 binary64))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64))) (.f64 #s(literal -1/16 binary64) (/.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 a #s(literal 6 binary64)))))))))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (*.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 a (+.f64 #s(literal 1 binary64) (fma.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 a #s(literal 4 binary64))) (*.f64 #s(literal -1/16 binary64) (/.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 a #s(literal 6 binary64))))))))))
#s(approx (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))) (.f64 #s(literal -1 binary64) (.f64 a (-.f64 (*.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))) #s(literal 1 binary64)))))
#s(approx (sqrt (- (* a a) (* b b))) (sqrt.f64 (pow.f64 a #s(literal 2 binary64))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 3 binary64)))) (*.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (+.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (.f64 #s(literal -1/16 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 5 binary64)))) (.f64 #s(literal 1/8 binary64) (/.f64 #s(literal 1 binary64) (pow.f64 (sqrt.f64 (pow.f64 a #s(literal 2 binary64))) #s(literal 3 binary64)))))) (*.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (sqrt.f64 (pow.f64 a #s(literal 2 binary64)))))))))
#s(approx (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))) (+.f64 a (*.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) a))))
#s(approx (sqrt (- (* a a) (* b b))) (*.f64 b (sqrt.f64 #s(literal -1 binary64))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (fma.f64 #s(literal 1/16 binary64) (/.f64 (pow.f64 a #s(literal 6 binary64)) (.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 5 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64))))))))))
#s(approx (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (/.f64 a (pow.f64 b #s(literal 2 binary64))) (.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) a)))))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) (.f64 (pow.f64 b #s(literal 2 binary64)) (-.f64 (/.f64 #s(literal 1 binary64) (pow.f64 b #s(literal 2 binary64))) (.f64 #s(literal 1/2 binary64) (/.f64 #s(literal 1 binary64) (pow.f64 a #s(literal 2 binary64)))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (sqrt.f64 #s(literal -1 binary64)))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64)))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (*.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64))))))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 #s(literal -1 binary64) (.f64 b (+.f64 (sqrt.f64 #s(literal -1 binary64)) (fma.f64 #s(literal -1/8 binary64) (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 3 binary64)))) (fma.f64 #s(literal 1/16 binary64) (/.f64 (pow.f64 a #s(literal 6 binary64)) (.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 (sqrt.f64 #s(literal -1 binary64)) #s(literal 5 binary64)))) (.f64 #s(literal 1/2 binary64) (/.f64 (pow.f64 a #s(literal 2 binary64)) (.f64 (pow.f64 b #s(literal 2 binary64)) (sqrt.f64 #s(literal -1 binary64)))))))))))

Calls

6 calls:

Time	Variable	Point
1.0ms	a	-inf
1.0ms	a	0
1.0ms	b	inf
1.0ms	a	inf
1.0ms	b	-inf

rewrite123.0ms (6%)

Memory

11.2MiB live, 56.9MiB allocated; 8ms collecting garbage

Counts

7 → 9

Calls

Call 1

Inputs
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))
(neg.f64 a)
a
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a #s(approx (+ 1 ( -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))))
(.f64 a #s(approx (+ 1 ( -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64)))
#s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))
#s(literal 1 binary64)

Outputs
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))
(neg.f64 a)
a
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a #s(approx (+ 1 (/ ( -1/2 (* b b)) (* a a))) #s(literal 1 binary64))))
(.f64 a #s(approx (+ 1 (/ ( -1/2 (* b b)) (* a a))) #s(literal 1 binary64)))
#s(approx (+ 1 (/ (* -1/2 (* b b)) (* a a))) #s(literal 1 binary64))
#s(literal 1 binary64)
(cosh.f64 #s(literal 0 binary64))
(exp.f64 #s(literal 0 binary64))

eval4.0ms (0.2%)

Memory: 4.3MiB live, 4.3MiB allocated; 0ms collecting garbage
Compiler: Compiled 1 222 to 208 computations (83% saved)

prune2.0ms (0.1%)

Memory

2.1MiB live, 2.1MiB allocated; 0ms collecting garbage

Pruning

3 alts after pruning (0 fresh and 3 done)

	Pruned	Kept	Total
New	39	0	39
Fresh	0	0	0
Picked	0	2	2
Done	0	1	1
Total	39	3	42

Accuracy

99.9%

Counts

42 → 3

Alt Table

Click to see full alt table

Status	Accuracy	Program
✓	53.5%	(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
✓	99.4%	#s(approx (sqrt (- (* a a) (* b b))) (.f64 a #s(approx (+ 1 ( -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))))
✓	1.0%	#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))

Status

Accuracy

Program

✓

53.5%

(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))

✓

99.4%

#s(approx (sqrt (- (* a a) (* b b))) (*.f64 a #s(approx (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))))

✓

1.0%

#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))

Compiler

Compiled 43 to 32 computations (25.6% saved)

regimes28.0ms (1.4%)

Memory

-20.2MiB live, 26.9MiB allocated; 2ms collecting garbage

Accuracy

Total 0.0b remaining (0%)

Threshold costs 0b (0%)

Counts

2 → 1

6 → 1

Calls

Call 1

Inputs
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))
#s(approx (sqrt (- (* a a) (* b b))) (*.f64 #s(literal -1 binary64) a))

Outputs
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))

Call 2

Inputs
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))
#s(approx (sqrt (- (* a a) (* b b))) (*.f64 #s(literal -1 binary64) a))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a #s(approx (+ 1 ( -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))))
(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a (+.f64 #s(literal 1 binary64) (.f64 #s(literal -1/2 binary64) (/.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64)))))))

Outputs
#s(approx (sqrt (- (* a a) (* b b))) (.f64 a #s(approx (+ 1 ( -1/2 (/ (pow b 2) (pow a 2)))) #s(literal 1 binary64))))

Calls

12 calls:

3.0ms	(-.f64 (.f64 a a) (.f64 b b))
3.0ms	b
3.0ms	a
2.0ms	(-.f64 (.f64 a a) (.f64 b b))
2.0ms	(*.f64 a a)

Results

Accuracy	Segments	Branch
1.0%	1	`(*.f64 b b)`
1.0%	1	`(*.f64 a a)`
1.0%	1	`(-.f64 (.f64 a a) (.f64 b b))`
1.0%	1	`(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))`
1.0%	1	`a`
1.0%	1	`b`
99.4%	1	`(*.f64 b b)`
99.4%	1	`(*.f64 a a)`
99.4%	1	`(-.f64 (.f64 a a) (.f64 b b))`
99.4%	1	`(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))`
99.4%	1	`a`
99.4%	1	`b`

Compiler

Compiled 131 to 82 computations (37.4% saved)

derivations1.0ms (0%)

Memory

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

Stop Event

1×

done

Compiler

Compiled 35 to 20 computations (42.9% saved)

preprocess27.0ms (1.3%)

Memory

23.2MiB live, 23.2MiB allocated; 0ms collecting garbage

Remove

(abs b)

(abs a)

Compiler

Compiled 231 to 166 computations (28.1% saved)

end0.0ms (0%)

Memory: 0.0MiB live, 0.0MiB allocated; 0ms collecting garbage

start0.0ms (0%)

analyze55.0ms (2.7%)

sample989.0ms (48.2%)

preprocess24.0ms (1.2%)

series137.0ms (6.7%)

rewrite148.0ms (7.2%)

eval7.0ms (0.4%)

prune3.0ms (0.1%)

series142.0ms (6.9%)

rewrite143.0ms (7%)

eval64.0ms (3.1%)

prune15.0ms (0.7%)

series138.0ms (6.7%)

rewrite123.0ms (6%)

eval4.0ms (0.2%)

prune2.0ms (0.1%)

regimes28.0ms (1.4%)

derivations1.0ms (0%)

preprocess27.0ms (1.3%)

end0.0ms (0%)

Profiling