Metrics for bug366, discussion (missed optimization)

start0.0ms (0%)

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

analyze34.0ms (2%)

Memory

14.5MiB live, 60.4MiB allocated; 5ms 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)

sample798.0ms (46%)

Memory

10.2MiB live, 950.5MiB allocated; 350ms collecting garbage

Samples

479.0ms	8 256×	0	valid
43.0ms	583×	0	invalid

Precisions

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

ival-mult: 168.0ms (59.6% of total)

ival-sqrt: 58.0ms (20.6% of total)

ival-sub: 51.0ms (18.1% of total)

ival-assert: 3.0ms (1.1% of total)

adjust: 2.0ms (0.7% of total)

Bogosity

explain106.0ms (6.1%)

Memory

-0.5MiB live, 89.4MiB allocated; 20ms collecting garbage

FPErrors

Click to see full error table

Ground Truth	Example	Example	Subexpression
132	-	-	`(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))`
12	-	-	`(-.f64 (.f64 a a) (.f64 b b))`
0	-	-	`(*.f64 a a)`
0	-	-	`a`
0	-	-	`b`
0	-	-	`(*.f64 b b)`

Explanations

Click to see full explanations table

Operator	Subexpression	Explanation	Count
`sqrt.f64`	`(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))`	oflow-rescue	120	0
↳	`(*.f64 b b)`	overflow	12
↳	`(-.f64 (.f64 a a) (.f64 b b))`	overflow	120
↳	`(*.f64 a a)`	overflow	120
`-.f64`	`(-.f64 (.f64 a a) (.f64 b b))`	nan-rescue	12	0
↳	`(*.f64 a a)`	overflow	120
↳	`(*.f64 b b)`	overflow	12
`sqrt.f64`	`(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))`	uflow-rescue	12	0
↳	`(*.f64 b b)`	underflow	97
↳	`(-.f64 (.f64 a a) (.f64 b b))`	underflow	12
↳	`(*.f64 a a)`	underflow	12

Confusion

	Predicted +	Predicted -
+	132	0
-	0	124

Precision

1.0

Recall

1.0

Confusion?

	Predicted +	Predicted Maybe	Predicted -
+	132	0	0
-	0	0	124

Precision?

1.0

Recall?

1.0

Freqs

test

number	freq
0	124
1	120
2	12

Total Confusion?

	Predicted +	Predicted Maybe	Predicted -
+	1	0	0
-	0	0	0

Precision?

1.0

Recall?

1.0

Samples

33.0ms

512×

valid

Compiler

Compiled 58 to 22 computations (62.1% saved)

Precisions

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

ival-mult: 7.0ms (52.8% of total)

ival-sqrt: 3.0ms (22.6% of total)

ival-sub: 2.0ms (15.1% of total)

ival-true: 0.0ms (0% of total)

adjust: 0.0ms (0% of total)

ival-assert: 0.0ms (0% of total)

preprocess119.0ms (6.9%)

Memory

18.0MiB live, 122.3MiB allocated; 11ms collecting garbage

Iterations

Useful iterations: 0 (0.0ms)

Iter	Nodes	Cost
0	18	50
1	55	50
2	214	50
3	344	50
4	875	50
5	1839	50
6	3017	50
7	3359	50
8	3437	50
9	3437	50
10	3633	50
0	6	8
0	10	8
1	23	8
2	83	8
3	225	8
4	405	8
5	836	8
6	1387	8
7	2071	8
8	2896	8
9	3206	8
10	3268	8
11	3268	8
12	3346	8
0	3346	8

Stop Event

1×	iter limit
1×	saturated
1×	iter limit
1×	saturated

Calls

Call 1

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

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

Symmetry

(abs a)

(abs b)

Compiler

Compiled 8 to 6 computations (25% saved)

eval0.0ms (0%)

Memory: 0.0MiB live, 0.0MiB allocated; 0ms collecting garbage
Compiler: Compiled 0 to 2 computations (-∞% saved)

prune0.0ms (0%)

Memory

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

Alt Table

Click to see full alt table

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

Compiler

Compiled 8 to 6 computations (25% saved)

series7.0ms (0.4%)

Memory

9.9MiB live, 9.9MiB allocated; 0ms 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))) #s(hole binary64 (* b (sqrt -1))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ (* 1/2 (/ (pow a 2) (* b (sqrt -1)))) (* b (sqrt -1)))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ (* b (sqrt -1)) (* (pow a 2) (+ (* -1/8 (/ (pow a 2) (* (pow b 3) (pow (sqrt -1) 3)))) (* 1/2 (/ 1 (* b (sqrt -1)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ (* b (sqrt -1)) (* (pow a 2) (+ (* (pow a 2) (- (* 1/16 (/ (pow a 2) (* (pow b 5) (pow (sqrt -1) 5)))) (* 1/8 (/ 1 (* (pow b 3) (pow (sqrt -1) 3)))))) (* 1/2 (/ 1 (* b (sqrt -1)))))))))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (* -1 (pow b 2))))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (- (pow a 2) (pow b 2))))
#s(approx (* a a) #s(hole binary64 (pow a 2)))
#s(approx a #s(hole binary64 a))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 a))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (* -1/8 (/ (pow b 4) (pow a 4))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (+ (* -1/8 (/ (pow b 4) (pow a 4))) (* -1/16 (/ (pow b 6) (pow a 6)))))))))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (pow a 2)))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (* (pow a 2) (+ 1 (* -1 (/ (pow b 2) (pow a 2)))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 a)))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (* -1/8 (/ (pow b 4) (pow a 4)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (+ (* -1/8 (/ (pow b 4) (pow a 4))) (* -1/16 (/ (pow b 6) (pow a 6))))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ a (* -1/2 (/ (pow b 2) a)))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ a (* (pow b 2) (- (* -1/8 (/ (pow b 2) (pow a 3))) (* 1/2 (/ 1 a)))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ a (* (pow b 2) (- (* (pow b 2) (- (* -1/16 (/ (pow b 2) (pow a 5))) (* 1/8 (/ 1 (pow a 3))))) (* 1/2 (/ 1 a)))))))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (+ (* -1 (pow b 2)) (pow a 2))))
#s(approx (* b b) #s(hole binary64 (pow b 2)))
#s(approx b #s(hole binary64 b))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* b (+ (sqrt -1) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (+ (* 1/16 (/ (pow a 6) (* (pow b 6) (pow (sqrt -1) 5)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1))))))))))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (* (pow b 2) (- (/ (pow a 2) (pow b 2)) 1))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (sqrt -1)))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1))))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (+ (* 1/16 (/ (pow a 6) (* (pow b 6) (pow (sqrt -1) 5)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1)))))))))))

Calls

6 calls:

Time	Variable		Point	Expression
1.0ms	a	@	0	((sqrt (- (* a a) (* b b))) (- (* a a) (* b b)) (* a a) a (* b b) b)
1.0ms	b	@	inf	((sqrt (- (* a a) (* b b))) (- (* a a) (* b b)) (* a a) a (* b b) b)
1.0ms	a	@	inf	((sqrt (- (* a a) (* b b))) (- (* a a) (* b b)) (* a a) a (* b b) b)
1.0ms	a	@	-inf	((sqrt (- (* a a) (* b b))) (- (* a a) (* b b)) (* a a) a (* b b) b)
1.0ms	b	@	0	((sqrt (- (* a a) (* b b))) (- (* a a) (* b b)) (* a a) a (* b b) b)

rewrite50.0ms (2.9%)

Memory

-6.6MiB live, 39.4MiB allocated; 4ms collecting garbage

Iterations

Useful iterations: 0 (0.0ms)

Iter	Nodes	Cost
0	159	856
0	164	821
0	818	795

Stop Event

1×	iter limit
1×	iter limit
1×	node limit
1×	iter limit

Counts

38 → 57

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
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* b (sqrt -1))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ (* 1/2 (/ (pow a 2) (* b (sqrt -1)))) (* b (sqrt -1)))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ (* b (sqrt -1)) (* (pow a 2) (+ (* -1/8 (/ (pow a 2) (* (pow b 3) (pow (sqrt -1) 3)))) (* 1/2 (/ 1 (* b (sqrt -1)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ (* b (sqrt -1)) (* (pow a 2) (+ (* (pow a 2) (- (* 1/16 (/ (pow a 2) (* (pow b 5) (pow (sqrt -1) 5)))) (* 1/8 (/ 1 (* (pow b 3) (pow (sqrt -1) 3)))))) (* 1/2 (/ 1 (* b (sqrt -1)))))))))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (* -1 (pow b 2))))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (- (pow a 2) (pow b 2))))
#s(approx (* a a) #s(hole binary64 (pow a 2)))
#s(approx a #s(hole binary64 a))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 a))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (* -1/8 (/ (pow b 4) (pow a 4))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (+ (* -1/8 (/ (pow b 4) (pow a 4))) (* -1/16 (/ (pow b 6) (pow a 6)))))))))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (pow a 2)))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (* (pow a 2) (+ 1 (* -1 (/ (pow b 2) (pow a 2)))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 a)))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (* -1/8 (/ (pow b 4) (pow a 4)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (+ (* -1/8 (/ (pow b 4) (pow a 4))) (* -1/16 (/ (pow b 6) (pow a 6))))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ a (* -1/2 (/ (pow b 2) a)))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ a (* (pow b 2) (- (* -1/8 (/ (pow b 2) (pow a 3))) (* 1/2 (/ 1 a)))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ a (* (pow b 2) (- (* (pow b 2) (- (* -1/16 (/ (pow b 2) (pow a 5))) (* 1/8 (/ 1 (pow a 3))))) (* 1/2 (/ 1 a)))))))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (+ (* -1 (pow b 2)) (pow a 2))))
#s(approx (* b b) #s(hole binary64 (pow b 2)))
#s(approx b #s(hole binary64 b))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* b (+ (sqrt -1) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (+ (* 1/16 (/ (pow a 6) (* (pow b 6) (pow (sqrt -1) 5)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1))))))))))
#s(approx (- (* a a) (* b b)) #s(hole binary64 (* (pow b 2) (- (/ (pow a 2) (pow b 2)) 1))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (sqrt -1)))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1))))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (+ (* 1/16 (/ (pow a 6) (* (pow b 6) (pow (sqrt -1) 5)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1)))))))))))

Outputs
(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 (pow.f64 a #s(literal 4 binary64)) (pow.f64 b #s(literal 4 binary64))) (fma.f64 a a (*.f64 b b)))
(/.f64 (-.f64 (pow.f64 (.f64 a a) #s(literal 3 binary64)) (pow.f64 (.f64 b b) #s(literal 3 binary64))) (+.f64 (pow.f64 a #s(literal 4 binary64)) (+.f64 (pow.f64 b #s(literal 4 binary64)) (.f64 (.f64 a a) (*.f64 b b)))))
(fma.f64 (pow.f64 a #s(literal 1 binary64)) (pow.f64 a #s(literal 1 binary64)) (*.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 (*.f64 (neg.f64 b) b))
(-.f64 (.f64 a a) (.f64 b b))
(+.f64 (.f64 a a) (.f64 (neg.f64 b) b))
(*.f64 (pow.f64 a #s(literal 1 binary64)) (pow.f64 a #s(literal 1 binary64)))
(*.f64 (neg.f64 a) (neg.f64 a))
(*.f64 (fabs.f64 a) (fabs.f64 a))
(*.f64 a a)
(pow.f64 a #s(literal 2 binary64))
(exp.f64 (*.f64 (log.f64 a) #s(literal 2 binary64)))
a
(*.f64 (pow.f64 b #s(literal 1 binary64)) (pow.f64 b #s(literal 1 binary64)))
(*.f64 (neg.f64 b) (neg.f64 b))
(*.f64 (fabs.f64 b) (fabs.f64 b))
(*.f64 b b)
(pow.f64 b #s(literal 2 binary64))
(exp.f64 (*.f64 (log.f64 b) #s(literal 2 binary64)))
b
#s(approx (sqrt (- (* a a) (* b b))) (*.f64 (sqrt.f64 #s(literal -1 binary64)) b))
#s(approx (sqrt (- (* a a) (* b b))) (fma.f64 (/.f64 (/.f64 (.f64 a a) b) (sqrt.f64 #s(literal -1 binary64))) #s(literal 1/2 binary64) (.f64 (sqrt.f64 #s(literal -1 binary64)) b)))
#s(approx (sqrt (- (* a a) (* b b))) (fma.f64 (fma.f64 (/.f64 (.f64 a a) (pow.f64 (.f64 (sqrt.f64 #s(literal -1 binary64)) b) #s(literal 3 binary64))) #s(literal -1/8 binary64) (/.f64 #s(literal 1/2 binary64) (.f64 (sqrt.f64 #s(literal -1 binary64)) b))) (.f64 a a) (*.f64 (sqrt.f64 #s(literal -1 binary64)) b)))
#s(approx (sqrt (- (* a a) (* b b))) (fma.f64 (fma.f64 (-.f64 (/.f64 (.f64 #s(literal 1/16 binary64) (.f64 a a)) (pow.f64 (.f64 (sqrt.f64 #s(literal -1 binary64)) b) #s(literal 5 binary64))) (/.f64 #s(literal 1/8 binary64) (pow.f64 (.f64 (sqrt.f64 #s(literal -1 binary64)) b) #s(literal 3 binary64)))) (.f64 a a) (/.f64 #s(literal 1/2 binary64) (.f64 (sqrt.f64 #s(literal -1 binary64)) b))) (.f64 a a) (.f64 (sqrt.f64 #s(literal -1 binary64)) b)))
#s(approx (- (* a a) (* b b)) (neg.f64 (*.f64 b b)))
#s(approx (- (* a a) (* b b)) (*.f64 (+.f64 a b) (-.f64 a b)))
#s(approx (* a a) (*.f64 a a))
#s(approx a a)
#s(approx (sqrt (- (* a a) (* b b))) a)
#s(approx (sqrt (- (* a a) (* b b))) (.f64 (fma.f64 (/.f64 (.f64 b b) (*.f64 a a)) #s(literal -1/2 binary64) #s(literal 1 binary64)) a))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 (+.f64 (fma.f64 (/.f64 (.f64 b b) (.f64 a a)) #s(literal -1/2 binary64) #s(literal 1 binary64)) (/.f64 (.f64 #s(literal -1/8 binary64) (pow.f64 b #s(literal 4 binary64))) (pow.f64 a #s(literal 4 binary64)))) a))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 (+.f64 (fma.f64 (/.f64 (.f64 b b) (.f64 a a)) #s(literal -1/2 binary64) (fma.f64 (/.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 a #s(literal 6 binary64))) #s(literal -1/16 binary64) (/.f64 (.f64 #s(literal -1/8 binary64) (pow.f64 b #s(literal 4 binary64))) (pow.f64 a #s(literal 4 binary64))))) #s(literal 1 binary64)) a))
#s(approx (- (* a a) (* b b)) (*.f64 a a))
#s(approx (- (* a a) (* b b)) (.f64 (fma.f64 (/.f64 (.f64 b b) (.f64 a a)) #s(literal -1 binary64) #s(literal 1 binary64)) (.f64 a a)))
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 a))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 (neg.f64 a) (fma.f64 (/.f64 (.f64 b b) (*.f64 a a)) #s(literal -1/2 binary64) #s(literal 1 binary64))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 (neg.f64 a) (+.f64 (fma.f64 (/.f64 (.f64 b b) (.f64 a a)) #s(literal -1/2 binary64) #s(literal 1 binary64)) (/.f64 (.f64 #s(literal -1/8 binary64) (pow.f64 b #s(literal 4 binary64))) (pow.f64 a #s(literal 4 binary64))))))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 (neg.f64 a) (+.f64 (fma.f64 (/.f64 (.f64 b b) (.f64 a a)) #s(literal -1/2 binary64) (fma.f64 (/.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 a #s(literal 6 binary64))) #s(literal -1/16 binary64) (/.f64 (.f64 #s(literal -1/8 binary64) (pow.f64 b #s(literal 4 binary64))) (pow.f64 a #s(literal 4 binary64))))) #s(literal 1 binary64))))
#s(approx (sqrt (- (* a a) (* b b))) (fma.f64 (/.f64 (*.f64 b b) a) #s(literal -1/2 binary64) a))
#s(approx (sqrt (- (* a a) (* b b))) (fma.f64 (-.f64 (/.f64 (.f64 #s(literal -1/8 binary64) (.f64 b b)) (pow.f64 a #s(literal 3 binary64))) (/.f64 #s(literal 1/2 binary64) a)) (*.f64 b b) a))
#s(approx (sqrt (- (* a a) (* b b))) (fma.f64 (-.f64 (.f64 (fma.f64 (/.f64 (.f64 b b) (pow.f64 a #s(literal 5 binary64))) #s(literal -1/16 binary64) (.f64 #s(literal -1/8 binary64) (pow.f64 a #s(literal -3 binary64)))) (.f64 b b)) (/.f64 #s(literal 1/2 binary64) a)) (*.f64 b b) a))
#s(approx (- (* a a) (* b b)) (fma.f64 a a (neg.f64 (*.f64 b b))))
#s(approx (* b b) (*.f64 b b))
#s(approx b b)
#s(approx (sqrt (- (* a a) (* b b))) (.f64 (fma.f64 (/.f64 (.f64 a a) (.f64 (.f64 b b) (sqrt.f64 #s(literal -1 binary64)))) #s(literal 1/2 binary64) (sqrt.f64 #s(literal -1 binary64))) b))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 (+.f64 (fma.f64 (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 #s(literal -1 binary64) #s(literal 3/2 binary64)))) #s(literal -1/8 binary64) (/.f64 (.f64 #s(literal 1/2 binary64) (.f64 a a)) (.f64 (.f64 b b) (sqrt.f64 #s(literal -1 binary64))))) (sqrt.f64 #s(literal -1 binary64))) b))
#s(approx (sqrt (- (* a a) (* b b))) (.f64 (+.f64 (fma.f64 (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 #s(literal -1 binary64) #s(literal 3/2 binary64)))) #s(literal -1/8 binary64) (fma.f64 (/.f64 (pow.f64 a #s(literal 6 binary64)) (.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 #s(literal -1 binary64) #s(literal 5/2 binary64)))) #s(literal 1/16 binary64) (/.f64 (.f64 #s(literal 1/2 binary64) (.f64 a a)) (.f64 (*.f64 b b) (sqrt.f64 #s(literal -1 binary64)))))) (sqrt.f64 #s(literal -1 binary64))) b))
#s(approx (- (* a a) (* b b)) (.f64 (-.f64 (/.f64 (.f64 a a) (.f64 b b)) #s(literal 1 binary64)) (.f64 b b)))
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 (*.f64 (sqrt.f64 #s(literal -1 binary64)) b)))
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 (.f64 (fma.f64 (/.f64 (.f64 a a) (.f64 (.f64 b b) (sqrt.f64 #s(literal -1 binary64)))) #s(literal 1/2 binary64) (sqrt.f64 #s(literal -1 binary64))) b)))
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 (.f64 (+.f64 (fma.f64 (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 #s(literal -1 binary64) #s(literal 3/2 binary64)))) #s(literal -1/8 binary64) (/.f64 (.f64 #s(literal 1/2 binary64) (.f64 a a)) (.f64 (.f64 b b) (sqrt.f64 #s(literal -1 binary64))))) (sqrt.f64 #s(literal -1 binary64))) b)))
#s(approx (sqrt (- (* a a) (* b b))) (neg.f64 (.f64 (+.f64 (fma.f64 (/.f64 (pow.f64 a #s(literal 4 binary64)) (.f64 (pow.f64 b #s(literal 4 binary64)) (pow.f64 #s(literal -1 binary64) #s(literal 3/2 binary64)))) #s(literal -1/8 binary64) (fma.f64 (/.f64 (pow.f64 a #s(literal 6 binary64)) (.f64 (pow.f64 b #s(literal 6 binary64)) (pow.f64 #s(literal -1 binary64) #s(literal 5/2 binary64)))) #s(literal 1/16 binary64) (/.f64 (.f64 #s(literal 1/2 binary64) (.f64 a a)) (.f64 (*.f64 b b) (sqrt.f64 #s(literal -1 binary64)))))) (sqrt.f64 #s(literal -1 binary64))) b)))

eval6.0ms (0.4%)

Memory: 4.4MiB live, 4.4MiB allocated; 0ms collecting garbage
Compiler: Compiled 1 095 to 230 computations (79% saved)

prune2.0ms (0.1%)

Memory

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

Pruning

3 alts after pruning (2 fresh and 1 done)

	Pruned	Kept	Total
New	59	2	61
Fresh	0	0	0
Picked	0	1	1
Done	0	0	0
Total	59	3	62

Accuracy

99.5%

Counts

62 → 3

Alt Table

Click to see full alt table

Status	Accuracy	Program
✓	51.5%	(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))
▶	51.9%	(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
▶	98.7%	#s(approx (sqrt (- (* a a) (* b b))) a)

Compiler

Compiled 26 to 19 computations (26.9% saved)

series9.0ms (0.5%)

Memory

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

Counts

7 → 74

Calls

Call 1

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

Outputs
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* b (sqrt -1))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ (* 1/2 (/ (pow a 2) (* b (sqrt -1)))) (* b (sqrt -1)))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ (* b (sqrt -1)) (* (pow a 2) (+ (* -1/8 (/ (pow a 2) (* (pow b 3) (pow (sqrt -1) 3)))) (* 1/2 (/ 1 (* b (sqrt -1)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ (* b (sqrt -1)) (* (pow a 2) (+ (* (pow a 2) (- (* 1/16 (/ (pow a 2) (* (pow b 5) (pow (sqrt -1) 5)))) (* 1/8 (/ 1 (* (pow b 3) (pow (sqrt -1) 3)))))) (* 1/2 (/ 1 (* b (sqrt -1)))))))))
#s(approx a #s(hole binary64 a))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* b (sqrt -1))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (+ (* 1/2 (/ (* a (+ b (* -1 b))) (* b (sqrt -1)))) (* b (sqrt -1)))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (+ (* a (+ (* 1/2 (/ (* a (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2)))))) (* b (sqrt -1)))) (* 1/2 (/ (+ b (* -1 b)) (* b (sqrt -1)))))) (* b (sqrt -1)))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (+ (* a (+ (* 1/2 (/ (+ b (* -1 b)) (* b (sqrt -1)))) (* a (+ (* -1/4 (/ (* a (* (+ b (* -1 b)) (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2))))))) (* (pow b 3) (pow (sqrt -1) 3)))) (* 1/2 (/ (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2))))) (* b (sqrt -1)))))))) (* b (sqrt -1)))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (* -1 (pow b 2))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (+ (* -1 (pow b 2)) (* a (+ b (* -1 b))))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (+ (* -1 (pow b 2)) (* a (+ a (+ b (* -1 b)))))))
#s(approx (+ a b) #s(hole binary64 b))
#s(approx (+ a b) #s(hole binary64 (+ a b)))
#s(approx (- a b) #s(hole binary64 (* -1 b)))
#s(approx (- a b) #s(hole binary64 (- a b)))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 a))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2)))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (* -1/8 (/ (pow b 4) (pow a 4))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (+ (* -1/8 (/ (pow b 4) (pow a 4))) (* -1/16 (/ (pow b 6) (pow a 6)))))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 a))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* a (+ 1 (* 1/2 (/ (+ b (* -1 b)) a))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* a (+ 1 (+ (* 1/2 (/ (+ b (* -1 b)) a)) (* 1/2 (/ (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2))) (pow a 2))))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* a (+ 1 (+ (* -1/4 (/ (* (+ b (* -1 b)) (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2)))) (pow a 3))) (+ (* 1/2 (/ (+ b (* -1 b)) a)) (* 1/2 (/ (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2))) (pow a 2)))))))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (pow a 2)))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (* (pow a 2) (+ 1 (+ (* -1 (/ b a)) (/ b a))))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (* (pow a 2) (+ 1 (+ (* -1 (/ b a)) (+ (* -1 (/ (pow b 2) (pow a 2))) (/ b a)))))))
#s(approx (+ a b) #s(hole binary64 a))
#s(approx (+ a b) #s(hole binary64 (* a (+ 1 (/ b a)))))
#s(approx (- a b) #s(hole binary64 a))
#s(approx (- a b) #s(hole binary64 (* a (+ 1 (* -1 (/ b a))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 a)))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* a (+ 1 (* -1/2 (/ (pow b 2) (pow a 2))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (* -1/8 (/ (pow b 4) (pow a 4)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* a (+ 1 (+ (* -1/2 (/ (pow b 2) (pow a 2))) (+ (* -1/8 (/ (pow b 4) (pow a 4))) (* -1/16 (/ (pow b 6) (pow a 6))))))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* -1 (* a (pow (sqrt -1) 2)))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1/2 (/ (+ b (* -1 b)) a)) (pow (sqrt -1) 2))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1 (/ (+ (* -1/2 (/ (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2)))) a)) (* 1/2 (+ b (* -1 b)))) a)) (pow (sqrt -1) 2))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1 (/ (+ (* -1 (/ (+ (* 1/4 (/ (* (+ b (* -1 b)) (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2))))) (* a (pow (sqrt -1) 2)))) (* 1/2 (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2)))))) a)) (* 1/2 (+ b (* -1 b)))) a)) (pow (sqrt -1) 2))))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (* (pow a 2) (+ 1 (* -1 (/ (+ (* -1 (+ b (* -1 b))) (/ (pow b 2) a)) a))))))
#s(approx (+ a b) #s(hole binary64 (* -1 (* a (- (* -1 (/ b a)) 1)))))
#s(approx (- a b) #s(hole binary64 (* -1 (* a (- (/ b a) 1)))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ a (* -1/2 (/ (pow b 2) a)))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ a (* (pow b 2) (- (* -1/8 (/ (pow b 2) (pow a 3))) (* 1/2 (/ 1 a)))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (+ a (* (pow b 2) (- (* (pow b 2) (- (* -1/16 (/ (pow b 2) (pow a 5))) (* 1/8 (/ 1 (pow a 3))))) (* 1/2 (/ 1 a)))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (+ a (* 1/2 (/ (* b (+ a (* -1 a))) a)))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (+ a (* b (+ (* -1/2 (/ (* b (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2))))) a)) (* 1/2 (/ (+ a (* -1 a)) a)))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (+ a (* b (+ (* 1/2 (/ (+ a (* -1 a)) a)) (* b (+ (* -1/2 (/ (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2)))) a)) (* 1/4 (/ (* b (* (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2)))) (+ a (* -1 a)))) (pow a 3))))))))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (+ (* b (+ a (* -1 a))) (pow a 2))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (+ (* b (+ a (+ (* -1 a) (* -1 b)))) (pow a 2))))
#s(approx b #s(hole binary64 b))
#s(approx (- a b) #s(hole binary64 (+ a (* -1 b))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* b (+ (sqrt -1) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (+ (* 1/16 (/ (pow a 6) (* (pow b 6) (pow (sqrt -1) 5)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1))))))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1))))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* (pow b 2) (sqrt -1)))))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* -1/4 (/ (* (+ a (* -1 a)) (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2))))) (* (pow b 3) (pow (sqrt -1) 3)))) (+ (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* (pow b 2) (sqrt -1))))))))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (* (pow b 2) (- (+ (* -1 (/ a b)) (/ a b)) 1))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (* (pow b 2) (- (+ (* -1 (/ a b)) (+ (/ a b) (/ (pow a 2) (pow b 2)))) 1))))
#s(approx (+ a b) #s(hole binary64 (* b (+ 1 (/ a b)))))
#s(approx (- a b) #s(hole binary64 (* b (- (/ a b) 1))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (sqrt -1)))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1)))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1))))))))))
#s(approx (sqrt (- (* a a) (* b b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (+ (* -1/8 (/ (pow a 4) (* (pow b 4) (pow (sqrt -1) 3)))) (+ (* 1/16 (/ (pow a 6) (* (pow b 6) (pow (sqrt -1) 5)))) (* 1/2 (/ (pow a 2) (* (pow b 2) (sqrt -1)))))))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* -1 (* b (sqrt -1)))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1 (/ (+ (* -1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* b (sqrt -1)))) (* 1/2 (/ (+ a (* -1 a)) (sqrt -1)))) b)))))))
#s(approx (sqrt (* (+ a b) (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1 (/ (+ (* -1 (/ (+ (* 1/4 (/ (* (+ a (* -1 a)) (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2))))) (* b (pow (sqrt -1) 3)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (sqrt -1)))) b)) (* 1/2 (/ (+ a (* -1 a)) (sqrt -1)))) b)))))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (* (pow b 2) (- (* -1 (/ (+ a (* -1 a)) b)) 1))))
#s(approx (* (+ a b) (- a b)) #s(hole binary64 (* (pow b 2) (- (* -1 (/ (+ a (+ (* -1 a) (* -1 (/ (pow a 2) b)))) b)) 1))))
#s(approx (+ a b) #s(hole binary64 (* -1 (* b (- (* -1 (/ a b)) 1)))))
#s(approx (- a b) #s(hole binary64 (* -1 (* b (+ 1 (* -1 (/ a b)))))))

Calls

6 calls:

Time	Variable		Point	Expression
2.0ms	b	@	inf	((sqrt (- (* a a) (* b b))) a (sqrt (* (+ a b) (- a b))) (* (+ a b) (- a b)) (+ a b) b (- a b))
2.0ms	a	@	0	((sqrt (- (* a a) (* b b))) a (sqrt (* (+ a b) (- a b))) (* (+ a b) (- a b)) (+ a b) b (- a b))
2.0ms	a	@	-inf	((sqrt (- (* a a) (* b b))) a (sqrt (* (+ a b) (- a b))) (* (+ a b) (- a b)) (+ a b) b (- a b))
1.0ms	a	@	inf	((sqrt (- (* a a) (* b b))) a (sqrt (* (+ a b) (- a b))) (* (+ a b) (- a b)) (+ a b) b (- a b))
1.0ms	b	@	-inf	((sqrt (- (* a a) (* b b))) a (sqrt (* (+ a b) (- a b))) (* (+ a b) (- a b)) (+ a b) b (- a b))

rewrite131.0ms (7.6%)

Memory

-33.4MiB live, 69.1MiB allocated; 31ms collecting garbage

Iterations

Useful iterations: 0 (0.0ms)

Iter	Nodes	Cost
0	420	2236
0	423	1985
0	2598	1911

Stop Event

1×	iter limit
1×	iter limit
1×	node limit
1×	iter limit

Counts

81 → 97

Calls

Call 1

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

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

eval14.0ms (0.8%)

Memory: -12.2MiB live, 32.0MiB allocated; 6ms collecting garbage
Compiler: Compiled 2 260 to 463 computations (79.5% saved)

prune2.0ms (0.1%)

Memory

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

Pruning

4 alts after pruning (1 fresh and 3 done)

	Pruned	Kept	Total
New	98	1	99
Fresh	0	0	0
Picked	0	2	2
Done	0	1	1
Total	98	4	102

Accuracy

100.0%

Counts

102 → 4

Alt Table

Click to see full alt table

Status	Accuracy	Program
▶	99.3%	(*.f64 (sqrt.f64 (+.f64 b a)) (sqrt.f64 (-.f64 a b)))
✓	51.5%	(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))
✓	51.9%	(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
✓	98.7%	#s(approx (sqrt (- (* a a) (* b b))) a)

Compiler

Compiled 35 to 26 computations (25.7% saved)

series11.0ms (0.7%)

Memory

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

Counts

7 → 83

Calls

Call 1

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

Outputs
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (sqrt -1))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* 1/2 (/ (* a (+ b (* -1 b))) (* b (sqrt -1)))) (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* a (+ (* 1/2 (/ (* a (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2)))))) (* b (sqrt -1)))) (* 1/2 (/ (+ b (* -1 b)) (* b (sqrt -1)))))) (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* a (+ (* 1/2 (/ (+ b (* -1 b)) (* b (sqrt -1)))) (* a (+ (* -1/4 (/ (* a (* (+ b (* -1 b)) (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2))))))) (* (pow b 3) (pow (sqrt -1) 3)))) (* 1/2 (/ (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2))))) (* b (sqrt -1)))))))) (* b (sqrt -1)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (sqrt b)))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* 1/2 (* a (sqrt (/ 1 b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* a (+ (* -1/8 (* a (sqrt (/ 1 (pow b 3))))) (* 1/2 (sqrt (/ 1 b))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* a (+ (* 1/2 (sqrt (/ 1 b))) (* a (+ (* -1/8 (sqrt (/ 1 (pow b 3)))) (* 1/16 (* a (sqrt (/ 1 (pow b 5))))))))))))
#s(approx (+ b a) #s(hole binary64 b))
#s(approx (+ b a) #s(hole binary64 (+ a b)))
#s(approx a #s(hole binary64 a))
#s(approx (sqrt (- a b)) #s(hole binary64 (* (sqrt b) (sqrt -1))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 b)))) (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* a (+ (* -1/8 (* (/ a (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 3))))) (* 1/2 (* (sqrt (/ 1 b)) (/ 1 (sqrt -1)))))) (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* a (+ (* 1/2 (* (sqrt (/ 1 b)) (/ 1 (sqrt -1)))) (* a (+ (* -1/8 (* (sqrt (/ 1 (pow b 3))) (/ 1 (pow (sqrt -1) 3)))) (* 1/16 (* (/ a (pow (sqrt -1) 5)) (sqrt (/ 1 (pow b 5))))))))) (* (sqrt b) (sqrt -1)))))
#s(approx (- a b) #s(hole binary64 (* -1 b)))
#s(approx (- a b) #s(hole binary64 (- a b)))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 a))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (* 1/2 (/ (+ b (* -1 b)) a))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (+ (* 1/2 (/ (+ b (* -1 b)) a)) (* 1/2 (/ (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2))) (pow a 2))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (+ (* -1/4 (/ (* (+ b (* -1 b)) (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2)))) (pow a 3))) (+ (* 1/2 (/ (+ b (* -1 b)) a)) (* 1/2 (/ (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2))) (pow a 2)))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (sqrt a)))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (pow b 3))) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b))))))))
#s(approx (+ b a) #s(hole binary64 a))
#s(approx (+ b a) #s(hole binary64 (* a (+ 1 (/ b a)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (sqrt a)))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (* -1/2 (* (sqrt (/ 1 (pow a 3))) b))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/2 (* (sqrt (/ 1 (pow a 3))) b)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (* -1/16 (* (sqrt (/ 1 (pow a 7))) (pow b 3)))))))))
#s(approx (- a b) #s(hole binary64 a))
#s(approx (- a b) #s(hole binary64 (* a (+ 1 (* -1 (/ b a))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (pow (sqrt -1) 2)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1/2 (/ (+ b (* -1 b)) a)) (pow (sqrt -1) 2))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1 (/ (+ (* -1/2 (/ (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2)))) a)) (* 1/2 (+ b (* -1 b)))) a)) (pow (sqrt -1) 2))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1 (/ (+ (* -1 (/ (+ (* 1/4 (/ (* (+ b (* -1 b)) (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2))))) (* a (pow (sqrt -1) 2)))) (* 1/2 (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2)))))) a)) (* 1/2 (+ b (* -1 b)))) a)) (pow (sqrt -1) 2))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* (sqrt a) (pow (sqrt -1) 2)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (/ (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (/ (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (/ (pow b 3) (pow (sqrt -1) 4)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))))
#s(approx (+ b a) #s(hole binary64 (* -1 (* a (- (* -1 (/ b a)) 1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* (sqrt a) (pow (sqrt -1) 2)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (* (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (* (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (* (pow b 3) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))))
#s(approx (- a b) #s(hole binary64 (* -1 (* a (- (/ b a) 1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* 1/2 (/ (* b (+ a (* -1 a))) a)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* b (+ (* -1/2 (/ (* b (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2))))) a)) (* 1/2 (/ (+ a (* -1 a)) a)))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* b (+ (* 1/2 (/ (+ a (* -1 a)) a)) (* b (+ (* -1/2 (/ (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2)))) a)) (* 1/4 (/ (* b (* (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2)))) (+ a (* -1 a)))) (pow a 3))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* 1/2 (* (sqrt (/ 1 a)) b)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/8 (* (sqrt (/ 1 (pow a 3))) b)) (* 1/2 (sqrt (/ 1 a))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* b (+ (* 1/2 (sqrt (/ 1 a))) (* b (+ (* -1/8 (sqrt (/ 1 (pow a 3)))) (* 1/16 (* (sqrt (/ 1 (pow a 5))) b)))))))))
#s(approx b #s(hole binary64 b))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* -1/2 (* (sqrt (/ 1 a)) b)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/2 (sqrt (/ 1 a))) (* -1/8 (* (sqrt (/ 1 (pow a 3))) b)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/2 (sqrt (/ 1 a))) (* b (+ (* -1/8 (sqrt (/ 1 (pow a 3)))) (* -1/16 (* (sqrt (/ 1 (pow a 5))) b)))))))))
#s(approx (- a b) #s(hole binary64 (+ a (* -1 b))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* (pow b 2) (sqrt -1)))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* -1/4 (/ (* (+ a (* -1 a)) (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2))))) (* (pow b 3) (pow (sqrt -1) 3)))) (+ (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* (pow b 2) (sqrt -1))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (* 1/2 (* a (sqrt (/ 1 (pow b 3)))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (+ (* -1/8 (* (pow a 2) (sqrt (/ 1 (pow b 5))))) (* 1/2 (* a (sqrt (/ 1 (pow b 3))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (+ (* -1/8 (* (pow a 2) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (pow a 3) (sqrt (/ 1 (pow b 7))))) (* 1/2 (* a (sqrt (/ 1 (pow b 3)))))))))))
#s(approx (+ b a) #s(hole binary64 (* b (+ 1 (/ a b)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (/ (pow a 3) (pow (sqrt -1) 5)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))))
#s(approx (- a b) #s(hole binary64 (* b (- (/ a b) 1))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1 (/ (+ (* -1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* b (sqrt -1)))) (* 1/2 (/ (+ a (* -1 a)) (sqrt -1)))) b)))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1 (/ (+ (* -1 (/ (+ (* 1/4 (/ (* (+ a (* -1 a)) (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2))))) (* b (pow (sqrt -1) 3)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (sqrt -1)))) b)) (* 1/2 (/ (+ a (* -1 a)) (sqrt -1)))) b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* (sqrt b) (pow (sqrt -1) 2)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 2)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 2)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (/ (pow a 3) (pow (sqrt -1) 4)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2)))))))))
#s(approx (+ b a) #s(hole binary64 (* -1 (* b (- (* -1 (/ a b)) 1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (* (pow a 2) (sqrt -1)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (* (pow a 2) (sqrt -1)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (* (pow a 3) (sqrt -1)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))))
#s(approx (- a b) #s(hole binary64 (* -1 (* b (+ 1 (* -1 (/ a b)))))))

Calls

6 calls:

Time	Variable		Point	Expression
2.0ms	b	@	-inf	((* (sqrt (+ b a)) (sqrt (- a b))) (sqrt (+ b a)) (+ b a) b a (sqrt (- a b)) (- a b))
2.0ms	a	@	-inf	((* (sqrt (+ b a)) (sqrt (- a b))) (sqrt (+ b a)) (+ b a) b a (sqrt (- a b)) (- a b))
2.0ms	a	@	0	((* (sqrt (+ b a)) (sqrt (- a b))) (sqrt (+ b a)) (+ b a) b a (sqrt (- a b)) (- a b))
2.0ms	b	@	inf	((* (sqrt (+ b a)) (sqrt (- a b))) (sqrt (+ b a)) (+ b a) b a (sqrt (- a b)) (- a b))
1.0ms	a	@	inf	((* (sqrt (+ b a)) (sqrt (- a b))) (sqrt (+ b a)) (+ b a) b a (sqrt (- a b)) (- a b))

rewrite145.0ms (8.4%)

Memory

5.9MiB live, 95.2MiB allocated; 11ms collecting garbage

Iterations

Useful iterations: 0 (0.0ms)

Iter	Nodes	Cost
0	526	2709
0	532	2268
0	3177	2164

Stop Event

1×	iter limit
1×	iter limit
1×	node limit
1×	iter limit

Counts

90 → 101

Calls

Call 1

Inputs
(*.f64 (sqrt.f64 (+.f64 b a)) (sqrt.f64 (-.f64 a b)))
(sqrt.f64 (+.f64 b a))
(+.f64 b a)
b
a
(sqrt.f64 (-.f64 a b))
(-.f64 a b)
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (sqrt -1))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* 1/2 (/ (* a (+ b (* -1 b))) (* b (sqrt -1)))) (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* a (+ (* 1/2 (/ (* a (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2)))))) (* b (sqrt -1)))) (* 1/2 (/ (+ b (* -1 b)) (* b (sqrt -1)))))) (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* a (+ (* 1/2 (/ (+ b (* -1 b)) (* b (sqrt -1)))) (* a (+ (* -1/4 (/ (* a (* (+ b (* -1 b)) (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2))))))) (* (pow b 3) (pow (sqrt -1) 3)))) (* 1/2 (/ (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2))))) (* b (sqrt -1)))))))) (* b (sqrt -1)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (sqrt b)))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* 1/2 (* a (sqrt (/ 1 b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* a (+ (* -1/8 (* a (sqrt (/ 1 (pow b 3))))) (* 1/2 (sqrt (/ 1 b))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* a (+ (* 1/2 (sqrt (/ 1 b))) (* a (+ (* -1/8 (sqrt (/ 1 (pow b 3)))) (* 1/16 (* a (sqrt (/ 1 (pow b 5))))))))))))
#s(approx (+ b a) #s(hole binary64 b))
#s(approx (+ b a) #s(hole binary64 (+ a b)))
#s(approx a #s(hole binary64 a))
#s(approx (sqrt (- a b)) #s(hole binary64 (* (sqrt b) (sqrt -1))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 b)))) (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* a (+ (* -1/8 (* (/ a (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 3))))) (* 1/2 (* (sqrt (/ 1 b)) (/ 1 (sqrt -1)))))) (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* a (+ (* 1/2 (* (sqrt (/ 1 b)) (/ 1 (sqrt -1)))) (* a (+ (* -1/8 (* (sqrt (/ 1 (pow b 3))) (/ 1 (pow (sqrt -1) 3)))) (* 1/16 (* (/ a (pow (sqrt -1) 5)) (sqrt (/ 1 (pow b 5))))))))) (* (sqrt b) (sqrt -1)))))
#s(approx (- a b) #s(hole binary64 (* -1 b)))
#s(approx (- a b) #s(hole binary64 (- a b)))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 a))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (* 1/2 (/ (+ b (* -1 b)) a))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (+ (* 1/2 (/ (+ b (* -1 b)) a)) (* 1/2 (/ (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2))) (pow a 2))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (+ (* -1/4 (/ (* (+ b (* -1 b)) (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2)))) (pow a 3))) (+ (* 1/2 (/ (+ b (* -1 b)) a)) (* 1/2 (/ (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2))) (pow a 2)))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (sqrt a)))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (pow b 3))) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b))))))))
#s(approx (+ b a) #s(hole binary64 a))
#s(approx (+ b a) #s(hole binary64 (* a (+ 1 (/ b a)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (sqrt a)))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (* -1/2 (* (sqrt (/ 1 (pow a 3))) b))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/2 (* (sqrt (/ 1 (pow a 3))) b)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (* -1/16 (* (sqrt (/ 1 (pow a 7))) (pow b 3)))))))))
#s(approx (- a b) #s(hole binary64 a))
#s(approx (- a b) #s(hole binary64 (* a (+ 1 (* -1 (/ b a))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (pow (sqrt -1) 2)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1/2 (/ (+ b (* -1 b)) a)) (pow (sqrt -1) 2))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1 (/ (+ (* -1/2 (/ (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2)))) a)) (* 1/2 (+ b (* -1 b)))) a)) (pow (sqrt -1) 2))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1 (/ (+ (* -1 (/ (+ (* 1/4 (/ (* (+ b (* -1 b)) (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2))))) (* a (pow (sqrt -1) 2)))) (* 1/2 (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2)))))) a)) (* 1/2 (+ b (* -1 b)))) a)) (pow (sqrt -1) 2))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* (sqrt a) (pow (sqrt -1) 2)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (/ (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (/ (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (/ (pow b 3) (pow (sqrt -1) 4)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))))
#s(approx (+ b a) #s(hole binary64 (* -1 (* a (- (* -1 (/ b a)) 1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* (sqrt a) (pow (sqrt -1) 2)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (* (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (* (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (* (pow b 3) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))))
#s(approx (- a b) #s(hole binary64 (* -1 (* a (- (/ b a) 1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* 1/2 (/ (* b (+ a (* -1 a))) a)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* b (+ (* -1/2 (/ (* b (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2))))) a)) (* 1/2 (/ (+ a (* -1 a)) a)))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* b (+ (* 1/2 (/ (+ a (* -1 a)) a)) (* b (+ (* -1/2 (/ (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2)))) a)) (* 1/4 (/ (* b (* (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2)))) (+ a (* -1 a)))) (pow a 3))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* 1/2 (* (sqrt (/ 1 a)) b)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/8 (* (sqrt (/ 1 (pow a 3))) b)) (* 1/2 (sqrt (/ 1 a))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* b (+ (* 1/2 (sqrt (/ 1 a))) (* b (+ (* -1/8 (sqrt (/ 1 (pow a 3)))) (* 1/16 (* (sqrt (/ 1 (pow a 5))) b)))))))))
#s(approx b #s(hole binary64 b))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* -1/2 (* (sqrt (/ 1 a)) b)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/2 (sqrt (/ 1 a))) (* -1/8 (* (sqrt (/ 1 (pow a 3))) b)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/2 (sqrt (/ 1 a))) (* b (+ (* -1/8 (sqrt (/ 1 (pow a 3)))) (* -1/16 (* (sqrt (/ 1 (pow a 5))) b)))))))))
#s(approx (- a b) #s(hole binary64 (+ a (* -1 b))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* (pow b 2) (sqrt -1)))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* -1/4 (/ (* (+ a (* -1 a)) (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2))))) (* (pow b 3) (pow (sqrt -1) 3)))) (+ (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* (pow b 2) (sqrt -1))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (* 1/2 (* a (sqrt (/ 1 (pow b 3)))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (+ (* -1/8 (* (pow a 2) (sqrt (/ 1 (pow b 5))))) (* 1/2 (* a (sqrt (/ 1 (pow b 3))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (+ (* -1/8 (* (pow a 2) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (pow a 3) (sqrt (/ 1 (pow b 7))))) (* 1/2 (* a (sqrt (/ 1 (pow b 3)))))))))))
#s(approx (+ b a) #s(hole binary64 (* b (+ 1 (/ a b)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (/ (pow a 3) (pow (sqrt -1) 5)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))))
#s(approx (- a b) #s(hole binary64 (* b (- (/ a b) 1))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1 (/ (+ (* -1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* b (sqrt -1)))) (* 1/2 (/ (+ a (* -1 a)) (sqrt -1)))) b)))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1 (/ (+ (* -1 (/ (+ (* 1/4 (/ (* (+ a (* -1 a)) (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2))))) (* b (pow (sqrt -1) 3)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (sqrt -1)))) b)) (* 1/2 (/ (+ a (* -1 a)) (sqrt -1)))) b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* (sqrt b) (pow (sqrt -1) 2)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 2)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 2)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (/ (pow a 3) (pow (sqrt -1) 4)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2)))))))))
#s(approx (+ b a) #s(hole binary64 (* -1 (* b (- (* -1 (/ a b)) 1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (* (pow a 2) (sqrt -1)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (* (pow a 2) (sqrt -1)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (* (pow a 3) (sqrt -1)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))))
#s(approx (- a b) #s(hole binary64 (* -1 (* b (+ 1 (* -1 (/ a b)))))))

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

eval21.0ms (1.2%)

Memory: -28.9MiB live, 26.4MiB allocated; 17ms collecting garbage
Compiler: Compiled 2 759 to 560 computations (79.7% saved)

prune3.0ms (0.2%)

Memory

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

Pruning

5 alts after pruning (1 fresh and 4 done)

	Pruned	Kept	Total
New	104	1	105
Fresh	0	0	0
Picked	0	1	1
Done	0	3	3
Total	104	5	109

Accuracy

100.0%

Counts

109 → 5

Alt Table

Click to see full alt table

Status	Accuracy	Program
✓	99.3%	(*.f64 (sqrt.f64 (+.f64 b a)) (sqrt.f64 (-.f64 a b)))
▶	97.5%	(*.f64 (sqrt.f64 #s(approx (+ b a) a)) (sqrt.f64 (-.f64 a b)))
✓	51.5%	(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))
✓	51.9%	(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
✓	98.7%	#s(approx (sqrt (- (* a a) (* b b))) a)

Compiler

Compiled 46 to 34 computations (26.1% saved)

series5.0ms (0.3%)

Memory

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

Counts

7 → 83

Calls

Call 1

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

Outputs
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (sqrt -1))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* 1/2 (/ (* a (+ b (* -1 b))) (* b (sqrt -1)))) (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* a (+ (* 1/2 (/ (* a (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2)))))) (* b (sqrt -1)))) (* 1/2 (/ (+ b (* -1 b)) (* b (sqrt -1)))))) (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* a (+ (* 1/2 (/ (+ b (* -1 b)) (* b (sqrt -1)))) (* a (+ (* -1/4 (/ (* a (* (+ b (* -1 b)) (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2))))))) (* (pow b 3) (pow (sqrt -1) 3)))) (* 1/2 (/ (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2))))) (* b (sqrt -1)))))))) (* b (sqrt -1)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (sqrt b)))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* 1/2 (* a (sqrt (/ 1 b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* a (+ (* -1/8 (* a (sqrt (/ 1 (pow b 3))))) (* 1/2 (sqrt (/ 1 b))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* a (+ (* 1/2 (sqrt (/ 1 b))) (* a (+ (* -1/8 (sqrt (/ 1 (pow b 3)))) (* 1/16 (* a (sqrt (/ 1 (pow b 5))))))))))))
#s(approx (+ b a) #s(hole binary64 b))
#s(approx (+ b a) #s(hole binary64 (+ a b)))
#s(approx a #s(hole binary64 a))
#s(approx (sqrt (- a b)) #s(hole binary64 (* (sqrt b) (sqrt -1))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 b)))) (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* a (+ (* -1/8 (* (/ a (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 3))))) (* 1/2 (* (sqrt (/ 1 b)) (/ 1 (sqrt -1)))))) (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* a (+ (* 1/2 (* (sqrt (/ 1 b)) (/ 1 (sqrt -1)))) (* a (+ (* -1/8 (* (sqrt (/ 1 (pow b 3))) (/ 1 (pow (sqrt -1) 3)))) (* 1/16 (* (/ a (pow (sqrt -1) 5)) (sqrt (/ 1 (pow b 5))))))))) (* (sqrt b) (sqrt -1)))))
#s(approx (- a b) #s(hole binary64 (* -1 b)))
#s(approx (- a b) #s(hole binary64 (- a b)))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 a))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (* 1/2 (/ (+ b (* -1 b)) a))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (+ (* 1/2 (/ (+ b (* -1 b)) a)) (* 1/2 (/ (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2))) (pow a 2))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (+ (* -1/4 (/ (* (+ b (* -1 b)) (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2)))) (pow a 3))) (+ (* 1/2 (/ (+ b (* -1 b)) a)) (* 1/2 (/ (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2))) (pow a 2)))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (sqrt a)))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (pow b 3))) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b))))))))
#s(approx (+ b a) #s(hole binary64 a))
#s(approx (+ b a) #s(hole binary64 (* a (+ 1 (/ b a)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (sqrt a)))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (* -1/2 (* (sqrt (/ 1 (pow a 3))) b))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/2 (* (sqrt (/ 1 (pow a 3))) b)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (* -1/16 (* (sqrt (/ 1 (pow a 7))) (pow b 3)))))))))
#s(approx (- a b) #s(hole binary64 a))
#s(approx (- a b) #s(hole binary64 (* a (+ 1 (* -1 (/ b a))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (pow (sqrt -1) 2)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1/2 (/ (+ b (* -1 b)) a)) (pow (sqrt -1) 2))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1 (/ (+ (* -1/2 (/ (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2)))) a)) (* 1/2 (+ b (* -1 b)))) a)) (pow (sqrt -1) 2))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1 (/ (+ (* -1 (/ (+ (* 1/4 (/ (* (+ b (* -1 b)) (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2))))) (* a (pow (sqrt -1) 2)))) (* 1/2 (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2)))))) a)) (* 1/2 (+ b (* -1 b)))) a)) (pow (sqrt -1) 2))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* (sqrt a) (pow (sqrt -1) 2)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (/ (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (/ (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (/ (pow b 3) (pow (sqrt -1) 4)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))))
#s(approx (+ b a) #s(hole binary64 (* -1 (* a (- (* -1 (/ b a)) 1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* (sqrt a) (pow (sqrt -1) 2)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (* (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (* (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (* (pow b 3) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))))
#s(approx (- a b) #s(hole binary64 (* -1 (* a (- (/ b a) 1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* 1/2 (/ (* b (+ a (* -1 a))) a)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* b (+ (* -1/2 (/ (* b (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2))))) a)) (* 1/2 (/ (+ a (* -1 a)) a)))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* b (+ (* 1/2 (/ (+ a (* -1 a)) a)) (* b (+ (* -1/2 (/ (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2)))) a)) (* 1/4 (/ (* b (* (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2)))) (+ a (* -1 a)))) (pow a 3))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* 1/2 (* (sqrt (/ 1 a)) b)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/8 (* (sqrt (/ 1 (pow a 3))) b)) (* 1/2 (sqrt (/ 1 a))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* b (+ (* 1/2 (sqrt (/ 1 a))) (* b (+ (* -1/8 (sqrt (/ 1 (pow a 3)))) (* 1/16 (* (sqrt (/ 1 (pow a 5))) b)))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* -1/2 (* (sqrt (/ 1 a)) b)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/2 (sqrt (/ 1 a))) (* -1/8 (* (sqrt (/ 1 (pow a 3))) b)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/2 (sqrt (/ 1 a))) (* b (+ (* -1/8 (sqrt (/ 1 (pow a 3)))) (* -1/16 (* (sqrt (/ 1 (pow a 5))) b)))))))))
#s(approx (- a b) #s(hole binary64 (+ a (* -1 b))))
#s(approx b #s(hole binary64 b))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* (pow b 2) (sqrt -1)))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* -1/4 (/ (* (+ a (* -1 a)) (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2))))) (* (pow b 3) (pow (sqrt -1) 3)))) (+ (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* (pow b 2) (sqrt -1))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (* 1/2 (* a (sqrt (/ 1 (pow b 3)))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (+ (* -1/8 (* (pow a 2) (sqrt (/ 1 (pow b 5))))) (* 1/2 (* a (sqrt (/ 1 (pow b 3))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (+ (* -1/8 (* (pow a 2) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (pow a 3) (sqrt (/ 1 (pow b 7))))) (* 1/2 (* a (sqrt (/ 1 (pow b 3)))))))))))
#s(approx (+ b a) #s(hole binary64 (* b (+ 1 (/ a b)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (/ (pow a 3) (pow (sqrt -1) 5)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))))
#s(approx (- a b) #s(hole binary64 (* b (- (/ a b) 1))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1 (/ (+ (* -1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* b (sqrt -1)))) (* 1/2 (/ (+ a (* -1 a)) (sqrt -1)))) b)))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1 (/ (+ (* -1 (/ (+ (* 1/4 (/ (* (+ a (* -1 a)) (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2))))) (* b (pow (sqrt -1) 3)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (sqrt -1)))) b)) (* 1/2 (/ (+ a (* -1 a)) (sqrt -1)))) b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* (sqrt b) (pow (sqrt -1) 2)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 2)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 2)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (/ (pow a 3) (pow (sqrt -1) 4)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2)))))))))
#s(approx (+ b a) #s(hole binary64 (* -1 (* b (- (* -1 (/ a b)) 1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (* (pow a 2) (sqrt -1)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (* (pow a 2) (sqrt -1)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (* (pow a 3) (sqrt -1)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))))
#s(approx (- a b) #s(hole binary64 (* -1 (* b (+ 1 (* -1 (/ a b)))))))

Calls

6 calls:

Time	Variable		Point	Expression
1.0ms	a	@	-inf	((* (sqrt (+ b a)) (sqrt (- a b))) (sqrt (+ b a)) (+ b a) a (sqrt (- a b)) (- a b) b)
1.0ms	b	@	-inf	((* (sqrt (+ b a)) (sqrt (- a b))) (sqrt (+ b a)) (+ b a) a (sqrt (- a b)) (- a b) b)
1.0ms	b	@	inf	((* (sqrt (+ b a)) (sqrt (- a b))) (sqrt (+ b a)) (+ b a) a (sqrt (- a b)) (- a b) b)
0.0ms	a	@	inf	((* (sqrt (+ b a)) (sqrt (- a b))) (sqrt (+ b a)) (+ b a) a (sqrt (- a b)) (- a b) b)
0.0ms	a	@	0	((* (sqrt (+ b a)) (sqrt (- a b))) (sqrt (+ b a)) (+ b a) a (sqrt (- a b)) (- a b) b)

rewrite117.0ms (6.8%)

Memory

8.6MiB live, 142.0MiB allocated; 14ms collecting garbage

Iterations

Useful iterations: 0 (0.0ms)

Iter	Nodes	Cost
0	525	2715
0	530	2274
0	3179	2170

Stop Event

1×	iter limit
1×	iter limit
1×	node limit
1×	iter limit

Counts

90 → 96

Calls

Call 1

Inputs
(*.f64 (sqrt.f64 #s(approx (+ b a) a)) (sqrt.f64 (-.f64 a b)))
(sqrt.f64 #s(approx (+ b a) a))
#s(approx (+ b a) a)
a
(sqrt.f64 (-.f64 a b))
(-.f64 a b)
b
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (sqrt -1))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* 1/2 (/ (* a (+ b (* -1 b))) (* b (sqrt -1)))) (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* a (+ (* 1/2 (/ (* a (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2)))))) (* b (sqrt -1)))) (* 1/2 (/ (+ b (* -1 b)) (* b (sqrt -1)))))) (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ (* a (+ (* 1/2 (/ (+ b (* -1 b)) (* b (sqrt -1)))) (* a (+ (* -1/4 (/ (* a (* (+ b (* -1 b)) (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2))))))) (* (pow b 3) (pow (sqrt -1) 3)))) (* 1/2 (/ (- 1 (* 1/4 (/ (pow (+ b (* -1 b)) 2) (* (pow b 2) (pow (sqrt -1) 2))))) (* b (sqrt -1)))))))) (* b (sqrt -1)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (sqrt b)))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* 1/2 (* a (sqrt (/ 1 b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* a (+ (* -1/8 (* a (sqrt (/ 1 (pow b 3))))) (* 1/2 (sqrt (/ 1 b))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt b) (* a (+ (* 1/2 (sqrt (/ 1 b))) (* a (+ (* -1/8 (sqrt (/ 1 (pow b 3)))) (* 1/16 (* a (sqrt (/ 1 (pow b 5))))))))))))
#s(approx (+ b a) #s(hole binary64 b))
#s(approx (+ b a) #s(hole binary64 (+ a b)))
#s(approx a #s(hole binary64 a))
#s(approx (sqrt (- a b)) #s(hole binary64 (* (sqrt b) (sqrt -1))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 b)))) (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* a (+ (* -1/8 (* (/ a (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 3))))) (* 1/2 (* (sqrt (/ 1 b)) (/ 1 (sqrt -1)))))) (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (* a (+ (* 1/2 (* (sqrt (/ 1 b)) (/ 1 (sqrt -1)))) (* a (+ (* -1/8 (* (sqrt (/ 1 (pow b 3))) (/ 1 (pow (sqrt -1) 3)))) (* 1/16 (* (/ a (pow (sqrt -1) 5)) (sqrt (/ 1 (pow b 5))))))))) (* (sqrt b) (sqrt -1)))))
#s(approx (- a b) #s(hole binary64 (* -1 b)))
#s(approx (- a b) #s(hole binary64 (- a b)))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 a))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (* 1/2 (/ (+ b (* -1 b)) a))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (+ (* 1/2 (/ (+ b (* -1 b)) a)) (* 1/2 (/ (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2))) (pow a 2))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* a (+ 1 (+ (* -1/4 (/ (* (+ b (* -1 b)) (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2)))) (pow a 3))) (+ (* 1/2 (/ (+ b (* -1 b)) a)) (* 1/2 (/ (- (* -1 (pow b 2)) (* 1/4 (pow (+ b (* -1 b)) 2))) (pow a 2)))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (sqrt a)))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (pow b 3))) (* 1/2 (* (sqrt (/ 1 (pow a 3))) b))))))))
#s(approx (+ b a) #s(hole binary64 a))
#s(approx (+ b a) #s(hole binary64 (* a (+ 1 (/ b a)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (sqrt a)))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (* -1/2 (* (sqrt (/ 1 (pow a 3))) b))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* a (+ (sqrt (/ 1 a)) (+ (* -1/2 (* (sqrt (/ 1 (pow a 3))) b)) (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (pow b 2))) (* -1/16 (* (sqrt (/ 1 (pow a 7))) (pow b 3)))))))))
#s(approx (- a b) #s(hole binary64 a))
#s(approx (- a b) #s(hole binary64 (* a (+ 1 (* -1 (/ b a))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (pow (sqrt -1) 2)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1/2 (/ (+ b (* -1 b)) a)) (pow (sqrt -1) 2))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1 (/ (+ (* -1/2 (/ (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2)))) a)) (* 1/2 (+ b (* -1 b)))) a)) (pow (sqrt -1) 2))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* a (+ (* -1 (/ (+ (* -1 (/ (+ (* 1/4 (/ (* (+ b (* -1 b)) (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2))))) (* a (pow (sqrt -1) 2)))) (* 1/2 (- (pow b 2) (* 1/4 (/ (pow (+ b (* -1 b)) 2) (pow (sqrt -1) 2)))))) a)) (* 1/2 (+ b (* -1 b)))) a)) (pow (sqrt -1) 2))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* (sqrt a) (pow (sqrt -1) 2)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (/ (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (/ (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (/ (pow b 3) (pow (sqrt -1) 4)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) b)) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))))
#s(approx (+ b a) #s(hole binary64 (* -1 (* a (- (* -1 (/ b a)) 1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* (sqrt a) (pow (sqrt -1) 2)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (* (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* a (+ (* -1/8 (* (sqrt (/ 1 (pow a 5))) (* (pow b 2) (pow (sqrt -1) 2)))) (+ (* 1/16 (* (sqrt (/ 1 (pow a 7))) (* (pow b 3) (pow (sqrt -1) 2)))) (+ (* 1/2 (* (sqrt (/ 1 (pow a 3))) (* b (pow (sqrt -1) 2)))) (* (sqrt (/ 1 a)) (pow (sqrt -1) 2)))))))))
#s(approx (- a b) #s(hole binary64 (* -1 (* a (- (/ b a) 1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* 1/2 (/ (* b (+ a (* -1 a))) a)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* b (+ (* -1/2 (/ (* b (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2))))) a)) (* 1/2 (/ (+ a (* -1 a)) a)))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (+ a (* b (+ (* 1/2 (/ (+ a (* -1 a)) a)) (* b (+ (* -1/2 (/ (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2)))) a)) (* 1/4 (/ (* b (* (+ 1 (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow a 2)))) (+ a (* -1 a)))) (pow a 3))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* 1/2 (* (sqrt (/ 1 a)) b)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/8 (* (sqrt (/ 1 (pow a 3))) b)) (* 1/2 (sqrt (/ 1 a))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (+ (sqrt a) (* b (+ (* 1/2 (sqrt (/ 1 a))) (* b (+ (* -1/8 (sqrt (/ 1 (pow a 3)))) (* 1/16 (* (sqrt (/ 1 (pow a 5))) b)))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* -1/2 (* (sqrt (/ 1 a)) b)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/2 (sqrt (/ 1 a))) (* -1/8 (* (sqrt (/ 1 (pow a 3))) b)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (+ (sqrt a) (* b (+ (* -1/2 (sqrt (/ 1 a))) (* b (+ (* -1/8 (sqrt (/ 1 (pow a 3)))) (* -1/16 (* (sqrt (/ 1 (pow a 5))) b)))))))))
#s(approx (- a b) #s(hole binary64 (+ a (* -1 b))))
#s(approx b #s(hole binary64 b))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* (pow b 2) (sqrt -1)))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* b (+ (sqrt -1) (+ (* -1/4 (/ (* (+ a (* -1 a)) (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2))))) (* (pow b 3) (pow (sqrt -1) 3)))) (+ (* 1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* (pow b 2) (sqrt -1))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (* 1/2 (* a (sqrt (/ 1 (pow b 3)))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (+ (* -1/8 (* (pow a 2) (sqrt (/ 1 (pow b 5))))) (* 1/2 (* a (sqrt (/ 1 (pow b 3))))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* b (+ (sqrt (/ 1 b)) (+ (* -1/8 (* (pow a 2) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (pow a 3) (sqrt (/ 1 (pow b 7))))) (* 1/2 (* a (sqrt (/ 1 (pow b 3)))))))))))
#s(approx (+ b a) #s(hole binary64 (* b (+ 1 (/ a b)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 3)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (/ (pow a 3) (pow (sqrt -1) 5)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* (/ a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))))
#s(approx (- a b) #s(hole binary64 (* b (- (/ a b) 1))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (sqrt -1)))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1/2 (/ (+ a (* -1 a)) (* b (sqrt -1)))))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1 (/ (+ (* -1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (* b (sqrt -1)))) (* 1/2 (/ (+ a (* -1 a)) (sqrt -1)))) b)))))))
#s(approx (* (sqrt (+ b a)) (sqrt (- a b))) #s(hole binary64 (* -1 (* b (+ (sqrt -1) (* -1 (/ (+ (* -1 (/ (+ (* 1/4 (/ (* (+ a (* -1 a)) (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2))))) (* b (pow (sqrt -1) 3)))) (* 1/2 (/ (- (pow a 2) (* 1/4 (/ (pow (+ a (* -1 a)) 2) (pow (sqrt -1) 2)))) (sqrt -1)))) b)) (* 1/2 (/ (+ a (* -1 a)) (sqrt -1)))) b)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* (sqrt b) (pow (sqrt -1) 2)))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2)))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 2)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2))))))))
#s(approx (sqrt (+ b a)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (/ (pow a 2) (pow (sqrt -1) 2)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (/ (pow a 3) (pow (sqrt -1) 4)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* a (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (pow (sqrt -1) 2)))))))))
#s(approx (+ b a) #s(hole binary64 (* -1 (* b (- (* -1 (/ a b)) 1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* (sqrt b) (sqrt -1)))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (* (pow a 2) (sqrt -1)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1))))))))
#s(approx (sqrt (- a b)) #s(hole binary64 (* -1 (* b (+ (* -1/8 (* (* (pow a 2) (sqrt -1)) (sqrt (/ 1 (pow b 5))))) (+ (* 1/16 (* (* (pow a 3) (sqrt -1)) (sqrt (/ 1 (pow b 7))))) (+ (* 1/2 (* (* a (sqrt -1)) (sqrt (/ 1 (pow b 3))))) (* (sqrt (/ 1 b)) (sqrt -1)))))))))
#s(approx (- a b) #s(hole binary64 (* -1 (* b (+ 1 (* -1 (/ a b)))))))

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

eval24.0ms (1.4%)

Memory: -29.4MiB live, 31.1MiB allocated; 7ms collecting garbage
Compiler: Compiled 2 731 to 530 computations (80.6% saved)

prune3.0ms (0.2%)

Memory

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

Pruning

5 alts after pruning (0 fresh and 5 done)

	Pruned	Kept	Total
New	97	0	97
Fresh	0	0	0
Picked	0	1	1
Done	0	4	4
Total	97	5	102

Accuracy

100.0%

Counts

102 → 5

Alt Table

Click to see full alt table

Status	Accuracy	Program
✓	99.3%	(*.f64 (sqrt.f64 (+.f64 b a)) (sqrt.f64 (-.f64 a b)))
✓	97.5%	(*.f64 (sqrt.f64 #s(approx (+ b a) a)) (sqrt.f64 (-.f64 a b)))
✓	51.5%	(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))
✓	51.9%	(sqrt.f64 (*.f64 (+.f64 a b) (-.f64 a b)))
✓	98.7%	#s(approx (sqrt (- (* a a) (* b b))) a)

Compiler

Compiled 92 to 57 computations (38% saved)

regimes11.0ms (0.7%)

Memory

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

Counts

5 → 1

Calls

Call 1

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

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

Calls

6 calls:

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

Results

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

Compiler

Compiled 23 to 21 computations (8.7% saved)

regimes26.0ms (1.5%)

Memory

-18.8MiB live, 27.8MiB allocated; 2ms collecting garbage

Accuracy

Total -31.2b remaining (-3777.3%)

Threshold costs -31.2b (-3777.3%)

Counts

4 → 1

Calls

Call 1

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

Outputs
#s(approx (sqrt (- (* a a) (* b b))) a)

Calls

6 calls:

18.0ms	b
2.0ms	a
2.0ms	(*.f64 a a)
2.0ms	(sqrt.f64 (-.f64 (.f64 a a) (.f64 b b)))
2.0ms	(-.f64 (.f64 a a) (.f64 b b))

Results

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

Compiler

Compiled 23 to 21 computations (8.7% saved)

simplify57.0ms (3.3%)

Memory

-13.5MiB live, 30.7MiB allocated; 7ms collecting garbage

Iterations

Useful iterations: 0 (0.0ms)

Iter	Nodes	Cost
0	12	19
1	27	19
2	92	19
3	235	19
4	413	19
5	850	19
6	1403	19
7	2087	19
8	2916	19
9	3226	19
10	3288	19
11	3288	19
12	3366	19

Stop Event

1×

saturated

Calls

Call 1

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

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

derivations9.0ms (0.5%)

Memory

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

Stop Event

1×

done

Compiler

Compiled 19 to 12 computations (36.8% saved)

preprocess20.0ms (1.2%)

Memory: 26.8MiB live, 26.8MiB allocated; 0ms collecting garbage
Remove: (abs b)
Compiler: Compiled 134 to 96 computations (28.4% saved)

end0.0ms (0%)

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

start0.0ms (0%)

analyze34.0ms (2%)

sample798.0ms (46%)

explain106.0ms (6.1%)

preprocess119.0ms (6.9%)

eval0.0ms (0%)

prune0.0ms (0%)

series7.0ms (0.4%)

rewrite50.0ms (2.9%)

eval6.0ms (0.4%)

prune2.0ms (0.1%)

series9.0ms (0.5%)

rewrite131.0ms (7.6%)

eval14.0ms (0.8%)

prune2.0ms (0.1%)

series11.0ms (0.7%)

rewrite145.0ms (8.4%)

eval21.0ms (1.2%)

prune3.0ms (0.2%)

series5.0ms (0.3%)

rewrite117.0ms (6.8%)

eval24.0ms (1.4%)

prune3.0ms (0.2%)

regimes11.0ms (0.7%)

regimes26.0ms (1.5%)

simplify57.0ms (3.3%)

derivations9.0ms (0.5%)

preprocess20.0ms (1.2%)

end0.0ms (0%)

Profiling