Details

analyze4.8s (18.2%)

Algorithm

1×

Search

True	Other	False	Iter
0%	99.8%	0.2%	0
0%	99.8%	0.2%	1
0%	99.8%	0.2%	2
0%	99.8%	0.2%	3
0%	99.8%	0.2%	4
0%	99.8%	0.2%	5
0%	99.8%	0.2%	6
0%	99.8%	0.2%	7
0%	99.8%	0.2%	8
0%	99.8%	0.2%	9
0%	99.8%	0.2%	10
0.8%	99%	0.2%	11
1.2%	98.6%	0.2%	12
9.9%	89.8%	0.2%	13
9.9%	89.8%	0.2%	14

Compiler

Compiled 39 to 24 computations (38.5% saved)

sample11.2s (42.2%)

Results

8.3s	4048×	body	1024	valid
1.5s	1362×	body	512	valid
813.0ms	2288×	body	128	valid
383.0ms	544×	body	256	valid
67.0ms	14×	body	2048	valid

Compiler

Compiled 116 to 71 computations (38.8% saved)

simplify52.0ms (0.2%)

Algorithm

1×	egg-herbie

Rules

571×	`fma-def_binary64`
244×	`fma-neg_binary64`
103×	`unsub-neg_binary64`
90×	`associate-l_binary64`
88×	`distribute-rgt-in_binary64`

Counts

1 → 2

Iterations

Useful iterations: 0 (0.0ms)

Iter	Nodes	Cost
0	23	53
1	40	53
2	66	53
3	85	53
4	99	53
5	133	53
6	173	53
7	217	53
8	281	53
9	305	53
10	362	53
11	407	53
12	468	53
13	604	53
14	728	53
15	836	53
16	943	53
17	891	53
18	893	53
19	899	53
20	904	53
21	841	53

prune4.0ms (0%)

Pruning

1 alts after pruning (1 fresh and 0 done)

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

Error

0.2b

Counts

3 → 1

Alt Table

Status	Error	Program
▶	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))

Status

Error

Program

▶

0.2b

(+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))

Compiler

Compiled 144 to 76 computations (47.2% saved)

localize23.0ms (0.1%)

Local error

Found 4 expressions with local error:

New	Error	Program
✓	0.2b	(*.f64 (sin.f64 theta) (sin.f64 delta))
✓	0.3b	(.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))
✓	0.3b	(asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))
✓	0.5b	(-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))

series118.0ms (0.4%)

Counts

4 → 132

Calls

4 calls:

56.0ms	(-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))
46.0ms	(.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))
9.0ms	(*.f64 (sin.f64 theta) (sin.f64 delta))
7.0ms	(asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))

rewrite13.0ms (0%)

Algorithm

1×	rewrite-expression-head

Rules

14×	`add-sqr-sqrt_binary64`
10×	`*-un-lft-identity_binary64`
10×	`add-cube-cbrt_binary64`
8×	`add-cbrt-cube_binary64`
8×	`add-exp-log_binary64`

Counts

4 → 72

Calls

4 calls:

3.0ms	(.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))
3.0ms	(*.f64 (sin.f64 theta) (sin.f64 delta))
3.0ms	(-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))
0.0ms	(asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))

simplify202.0ms (0.8%)

Algorithm

1×	egg-herbie

Rules

447×	`fma-neg_binary64`
376×	`distribute-rgt-neg-in_binary64`
325×	`distribute-lft-neg-in_binary64`
257×	`associate-r_binary64`
241×	`neg-mul-1_binary64`

Counts

204 → 185

Iterations

Useful iterations: 2 (0.0ms)

Iter	Nodes	Cost
0	262	6113
1	696	5704
2	2227	5493
3	4796	5493
4	4989	5493
5	4924	5493

prune295.0ms (1.1%)

Pruning

9 alts after pruning (9 fresh and 0 done)

	Pruned	Kept	Total
New	176	9	185
Fresh	0	0	0
Picked	1	0	1
Done	0	0	0
Total	177	9	186

Error

0.0b

Counts

186 → 9

Alt Table

Status	Error	Program
▶	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (cbrt.f64 (.f64 (.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))))
	0.3b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (.f64 (.f64 (cbrt.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))) (cbrt.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))))) (cbrt.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))))
	14.7b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (.f64 (cos.f64 phi1) (-.f64 (cos.f64 phi1) (.f64 (*.f64 (sin.f64 phi1) delta) (cos.f64 theta))))))
	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (log.f64 (exp.f64 (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))))
	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (cbrt.f64 (pow.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (sin.f64 phi1) (cos.f64 delta) (.f64 (*.f64 (cos.f64 phi1) (cos.f64 theta)) (sin.f64 delta))))))) 3))))
	0.3b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (.f64 (cbrt.f64 (sin.f64 theta)) (cbrt.f64 (sin.f64 theta))) (.f64 (cbrt.f64 (sin.f64 theta)) (sin.f64 delta))) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))
	5.2b	(+.f64 lambda1 (atan2.f64 (.f64 (cbrt.f64 (.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 theta)) (sin.f64 theta)) (.f64 (.f64 (sin.f64 delta) (sin.f64 delta)) (sin.f64 delta)))) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))
	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (log.f64 (/.f64 (exp.f64 (cos.f64 delta)) (exp.f64 (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))))
	24.3b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (+.f64 (.f64 (.f64 (sin.f64 phi1) (sin.f64 delta)) (+.f64 (.f64 (.f64 (cos.f64 phi1) (.f64 theta theta)) 1/2) (.f64 (.f64 (cos.f64 phi1) (pow.f64 theta 6)) 1/720))) (-.f64 (.f64 (.f64 (cos.f64 phi1) (cos.f64 phi1)) (cos.f64 delta)) (.f64 (.f64 (sin.f64 phi1) (sin.f64 delta)) (+.f64 (cos.f64 phi1) (.f64 (*.f64 (cos.f64 phi1) (pow.f64 theta 4)) 1/24)))))))

Compiler

Compiled 8698 to 1950 computations (77.6% saved)

localize27.0ms (0.1%)

Local error

Found 4 expressions with local error:

New	Error	Program
	0.3b	(asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))
✓	0.4b	(.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))
✓	0.5b	(-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (cbrt.f64 (.f64 (.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))
✓	7.0b	(cbrt.f64 (.f64 (.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))

series87.0ms (0.3%)

Counts

3 → 108

Calls

3 calls:

58.0ms

(-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (cbrt.f64 (*.f64 (*.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))

22.0ms

(*.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))

7.0ms

(cbrt.f64 (*.f64 (*.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))

rewrite17.0ms (0.1%)

Algorithm

1×	rewrite-expression-head

Rules

11×	`add-sqr-sqrt_binary64`
10×	`pow1_binary64`
9×	`*-un-lft-identity_binary64`
9×	`add-cube-cbrt_binary64`
5×	`add-cbrt-cube_binary64`

Counts

3 → 59

Calls

3 calls:

6.0ms

(cbrt.f64 (*.f64 (*.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))

3.0ms

(-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (cbrt.f64 (*.f64 (*.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))

3.0ms

(*.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))

simplify210.0ms (0.8%)

Algorithm

1×	egg-herbie

Rules

661×	`fma-def_binary64`
411×	`cancel-sign-sub-inv_binary64`
387×	`sub-neg_binary64`
346×	`*-commutative_binary64`
269×	`distribute-rgt-neg-in_binary64`

Counts

167 → 116

Iterations

Useful iterations: 4 (0.0ms)

Iter	Nodes	Cost
0	171	7110
1	434	5461
2	1313	4997
3	3463	4957
4	4776	4771
5	4904	4771
6	4904	4771

prune438.0ms (1.6%)

Pruning

8 alts after pruning (8 fresh and 0 done)

	Pruned	Kept	Total
New	172	4	176
Fresh	4	4	8
Picked	1	0	1
Done	0	0	0
Total	177	8	185

Error

0.0b

Counts

185 → 8

Alt Table

Status	Error	Program
	0.3b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (.f64 (cbrt.f64 (sin.f64 theta)) (cbrt.f64 (sin.f64 theta))) (.f64 (cbrt.f64 (sin.f64 theta)) (sin.f64 delta))) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))
▶	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (.f64 (cbrt.f64 (.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))) (cbrt.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))))
	5.2b	(+.f64 lambda1 (atan2.f64 (.f64 (cbrt.f64 (.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 theta)) (sin.f64 theta)) (.f64 (.f64 (sin.f64 delta) (sin.f64 delta)) (sin.f64 delta)))) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))
	14.7b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (.f64 (cos.f64 phi1) (-.f64 (cos.f64 phi1) (.f64 (*.f64 (sin.f64 phi1) delta) (cos.f64 theta))))))
	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (cbrt.f64 (pow.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (sin.f64 phi1) (cos.f64 delta) (.f64 (*.f64 (cos.f64 phi1) (cos.f64 theta)) (sin.f64 delta))))))) 3))))
	10.5b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (.f64 (.f64 (sin.f64 phi1) (sin.f64 delta)) (-.f64 (.f64 1/2 (.f64 (cos.f64 phi1) (.f64 theta theta))) (cos.f64 phi1))) (.f64 (cos.f64 delta) (+.f64 (pow.f64 (sin.f64 phi1) 2) -1)))))
	15.9b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (fma.f64 (sqrt.f64 (cos.f64 delta)) (sqrt.f64 (cos.f64 delta)) (neg.f64 (.f64 (sin.f64 phi1) (sin.f64 (cbrt.f64 (.f64 (.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))))))))
	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (expm1.f64 (log1p.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta))))))))))))

Compiler

Compiled 12753 to 1921 computations (84.9% saved)

localize27.0ms (0.1%)

Local error

Found 4 expressions with local error:

New	Error	Program
	0.4b	(.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))
✓	0.5b	(-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (.f64 (cbrt.f64 (.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))) (cbrt.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))
✓	0.6b	(cbrt.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))
✓	4.5b	(cbrt.f64 (.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))

series2.9s (11.1%)

Counts

3 → 108

Calls

3 calls:

1.7s

(cbrt.f64 (*.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))

1.1s

(cbrt.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))

58.0ms

(-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (*.f64 (cbrt.f64 (*.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))) (cbrt.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))

rewrite8.0ms (0%)

Algorithm

1×	rewrite-expression-head

Rules

6×	`*-un-lft-identity_binary64`
6×	`add-sqr-sqrt_binary64`
6×	`add-cube-cbrt_binary64`
5×	`add-log-exp_binary64`
4×	`cbrt-prod_binary64`

Counts

3 → 44

Calls

3 calls:

3.0ms

(-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (*.f64 (cbrt.f64 (*.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))) (cbrt.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))

2.0ms

(cbrt.f64 (*.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))

1.0ms

(cbrt.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (*.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))

simplify200.0ms (0.8%)

Algorithm

1×	egg-herbie

Rules

661×	`fma-def_binary64`
411×	`cancel-sign-sub-inv_binary64`
387×	`sub-neg_binary64`
344×	`*-commutative_binary64`
269×	`distribute-rgt-neg-in_binary64`

Counts

152 → 101

Iterations

Useful iterations: 4 (0.0ms)

Iter	Nodes	Cost
0	163	6236
1	413	5852
2	1264	5538
3	3480	4736
4	4690	4732
5	4919	4732
6	4954	4732

prune351.0ms (1.3%)

Pruning

8 alts after pruning (7 fresh and 1 done)

	Pruned	Kept	Total
New	139	0	139
Fresh	0	7	7
Picked	0	1	1
Done	0	0	0
Total	139	8	147

Error

0.0b

Counts

147 → 8

Alt Table

Status	Error	Program
	0.3b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (.f64 (cbrt.f64 (sin.f64 theta)) (cbrt.f64 (sin.f64 theta))) (.f64 (cbrt.f64 (sin.f64 theta)) (sin.f64 delta))) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))
✓	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (.f64 (cbrt.f64 (.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))) (cbrt.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))))
	5.2b	(+.f64 lambda1 (atan2.f64 (.f64 (cbrt.f64 (.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 theta)) (sin.f64 theta)) (.f64 (.f64 (sin.f64 delta) (sin.f64 delta)) (sin.f64 delta)))) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))
	14.7b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (.f64 (cos.f64 phi1) (-.f64 (cos.f64 phi1) (.f64 (*.f64 (sin.f64 phi1) delta) (cos.f64 theta))))))
	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (cbrt.f64 (pow.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (sin.f64 phi1) (cos.f64 delta) (.f64 (*.f64 (cos.f64 phi1) (cos.f64 theta)) (sin.f64 delta))))))) 3))))
	10.5b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (.f64 (.f64 (sin.f64 phi1) (sin.f64 delta)) (-.f64 (.f64 1/2 (.f64 (cos.f64 phi1) (.f64 theta theta))) (cos.f64 phi1))) (.f64 (cos.f64 delta) (+.f64 (pow.f64 (sin.f64 phi1) 2) -1)))))
	15.9b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (fma.f64 (sqrt.f64 (cos.f64 delta)) (sqrt.f64 (cos.f64 delta)) (neg.f64 (.f64 (sin.f64 phi1) (sin.f64 (cbrt.f64 (.f64 (.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))))))))
▶	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (expm1.f64 (log1p.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta))))))))))))

Compiler

Compiled 10349 to 1604 computations (84.5% saved)

localize27.0ms (0.1%)

Local error

Found 4 expressions with local error:

New	Error	Program
✓	0.3b	(.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta)))))))
✓	0.3b	(log1p.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta)))))))))
✓	0.3b	(asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (.f64 (cos.f64 phi1) (cos.f64 theta)))))
✓	0.5b	(-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta))))))))

series2.0s (7.6%)

Counts

4 → 108

Calls

4 calls:

1.9s	(log1p.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta)))))))))
60.0ms	(-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta))))))))
48.0ms	(.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta)))))))
6.0ms	(asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (.f64 (cos.f64 phi1) (cos.f64 theta)))))

rewrite11.0ms (0%)

Algorithm

1×	rewrite-expression-head

Rules

10×	`add-sqr-sqrt_binary64`
8×	`*-un-lft-identity_binary64`
8×	`add-cube-cbrt_binary64`
6×	`add-cbrt-cube_binary64`
6×	`add-exp-log_binary64`

Counts

4 → 62

Calls

4 calls:

4.0ms	(.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta)))))))
3.0ms	(-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta))))))))
0.0ms	(log1p.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta)))))))))
0.0ms	(asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (.f64 (cos.f64 phi1) (cos.f64 theta)))))

simplify186.0ms (0.7%)

Algorithm

1×	egg-herbie

Rules

524×	`fma-neg_binary64`
364×	`associate-r_binary64`
307×	`sub-neg_binary64`
306×	`associate-l_binary64`
191×	`*-commutative_binary64`

Counts

170 → 178

Iterations

Useful iterations: 2 (0.0ms)

Iter	Nodes	Cost
0	319	7643
1	941	6975
2	3776	6949
3	5000	6949
4	5047	6949

prune804.0ms (3%)

Pruning

7 alts after pruning (6 fresh and 1 done)

	Pruned	Kept	Total
New	177	1	178
Fresh	1	5	6
Picked	1	0	1
Done	0	1	1
Total	179	7	186

Error

0.0b

Counts

186 → 7

Alt Table

Status	Error	Program
	0.3b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (.f64 (cbrt.f64 (sin.f64 theta)) (cbrt.f64 (sin.f64 theta))) (.f64 (cbrt.f64 (sin.f64 theta)) (sin.f64 delta))) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))
✓	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (.f64 (cbrt.f64 (.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))) (cbrt.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))))))))))
	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (expm1.f64 (log1p.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (log.f64 (exp.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (sin.f64 delta) (*.f64 (cos.f64 phi1) (cos.f64 theta))))))))))))))
	14.7b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (.f64 (cos.f64 phi1) (-.f64 (cos.f64 phi1) (.f64 (*.f64 (sin.f64 phi1) delta) (cos.f64 theta))))))
	0.2b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (cbrt.f64 (pow.f64 (-.f64 (cos.f64 delta) (.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (fma.f64 (sin.f64 phi1) (cos.f64 delta) (.f64 (*.f64 (cos.f64 phi1) (cos.f64 theta)) (sin.f64 delta))))))) 3))))
	10.5b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (.f64 (.f64 (sin.f64 phi1) (sin.f64 delta)) (-.f64 (.f64 1/2 (.f64 (cos.f64 phi1) (.f64 theta theta))) (cos.f64 phi1))) (.f64 (cos.f64 delta) (+.f64 (pow.f64 (sin.f64 phi1) 2) -1)))))
	15.9b	(+.f64 lambda1 (atan2.f64 (.f64 (.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (fma.f64 (sqrt.f64 (cos.f64 delta)) (sqrt.f64 (cos.f64 delta)) (neg.f64 (.f64 (sin.f64 phi1) (sin.f64 (cbrt.f64 (.f64 (.f64 (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta)))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))) (asin.f64 (fma.f64 (cos.f64 delta) (sin.f64 phi1) (.f64 (*.f64 (sin.f64 delta) (cos.f64 phi1)) (cos.f64 theta))))))))))))

Compiler

Compiled 9559 to 2554 computations (73.3% saved)

regimes2.2s (8.4%)

Accuracy

Total 0.1b remaining (57.7%)

Threshold costs 0.1b (57.7%)

Counts

177 → 1

Compiler

Compiled 54073 to 30247 computations (44.1% saved)

simplify10.0ms (0%)

Algorithm

1×	egg-herbie

Rules

9×	`*-commutative_binary64`
5×	`+-commutative_binary64`
3×	`sub-neg_binary64`
3×	`neg-sub0_binary64`
3×	`neg-mul-1_binary64`

Iterations

Useful iterations: 0 (0.0ms)

Iter	Nodes	Cost
0	30	276
1	48	276
2	56	276
3	60	276
4	61	276
5	60	276

end212.0ms (0.8%)

Compiler: Compiled 681 to 234 computations (65.6% saved)

Iter	Nodes	Cost
0	23	53
1	40	53
2	66	53
3	85	53
4	99	53
5	133	53
6	173	53
7	217	53
8	281	53
9	305	53
10	362	53
11	407	53
12	468	53
13	604	53
14	728	53
15	836	53
16	943	53
17	891	53
18	893	53
19	899	53
20	904	53
21	841	53

Iter	Nodes	Cost
0	23	53
1	40	53
2	66	53
3	85	53
4	99	53
5	133	53
6	173	53
7	217	53
8	281	53
9	305	53
10	362	53
11	407	53
12	468	53
13	604	53
14	728	53
15	836	53
16	943	53
17	891	53
18	893	53
19	899	53
20	904	53
21	841	53

Details

analyze4.8s (18.2%)

sample11.2s (42.2%)

simplify52.0ms (0.2%)

prune4.0ms (0%)

localize23.0ms (0.1%)

series118.0ms (0.4%)

rewrite13.0ms (0%)

simplify202.0ms (0.8%)

prune295.0ms (1.1%)

localize27.0ms (0.1%)

series87.0ms (0.3%)

rewrite17.0ms (0.1%)

simplify210.0ms (0.8%)

prune438.0ms (1.6%)

localize27.0ms (0.1%)

series2.9s (11.1%)

rewrite8.0ms (0%)

simplify200.0ms (0.8%)

prune351.0ms (1.3%)

localize27.0ms (0.1%)

series2.0s (7.6%)

rewrite11.0ms (0%)

simplify186.0ms (0.7%)

prune804.0ms (3%)

regimes2.2s (8.4%)

simplify10.0ms (0%)

end212.0ms (0.8%)

Profiling

Iter	Nodes	Cost
0	23	53
1	40	53
2	66	53
3	85	53
4	99	53
5	133	53
6	173	53
7	217	53
8	281	53
9	305	53
10	362	53
11	407	53
12	468	53
13	604	53
14	728	53
15	836	53
16	943	53
17	891	53
18	893	53
19	899	53
20	904	53
21	841	53