Details

analyze1.5s (9.4%)

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
49.9%	49.9%	0.2%	3
49.9%	49.9%	0.2%	4
49.9%	49.9%	0.2%	5
49.9%	49.9%	0.2%	6
74.9%	25%	0.2%	7
74.9%	25%	0.2%	8
74.9%	25%	0.2%	9
74.9%	25%	0.2%	10
87.3%	12.5%	0.2%	11
87.3%	12.5%	0.2%	12
87.3%	12.5%	0.2%	13
87.3%	12.5%	0.2%	14

Compiler

Compiled 34 to 30 computations (11.8% saved)

sample6.8s (42%)

Symmetry

(sort kx ky)

Results

3.4s	6572×	body	128	valid
2.0s	802×	body	1024	valid
886.0ms	556×	body	512	valid
340.0ms	326×	body	256	valid

Compiler

Compiled 101 to 89 computations (11.9% saved)

simplify85.0ms (0.5%)

Algorithm

1×	egg-herbie

Rules

665×	`fma-def_binary64`
575×	`associate-l_binary64`
512×	`unswap-sqr_binary64`
437×	`associate-r_binary64`
369×	`associate-/r/_binary64`

Counts

1 → 2

Iterations

Useful iterations: 1 (0.0ms)

Iter	Nodes	Cost
0	24	42
1	56	37
2	113	37
3	465	37
4	2149	37
5	4982	37
6	4847	37
7	4964	37

prune6.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

1.0b

Counts

3 → 1

Alt Table

Status	Error	Program
▶	1.0b	(sqrt.f64 (+.f64 1/2 (/.f64 1/2 (sqrt.f64 (fma.f64 (pow.f64 (/.f64 (*.f64 2 l) Om) 2) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)) 1)))))

Compiler

Compiled 113 to 95 computations (15.9% saved)

localize18.0ms (0.1%)

Local error

Found 4 expressions with local error:

New	Error	Program
✓	0.2b	(pow.f64 (sin.f64 kx) 2)
✓	0.3b	(pow.f64 (sin.f64 ky) 2)
✓	2.4b	(fma.f64 (pow.f64 (/.f64 (*.f64 2 l) Om) 2) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)) 1)
✓	8.4b	(sqrt.f64 (fma.f64 (pow.f64 (/.f64 (*.f64 2 l) Om) 2) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)) 1))

series647.0ms (4%)

Counts

4 → 104

Calls

4 calls:

404.0ms	(sqrt.f64 (fma.f64 (pow.f64 (/.f64 (*.f64 2 l) Om) 2) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)) 1))
229.0ms	(fma.f64 (pow.f64 (/.f64 (*.f64 2 l) Om) 2) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)) 1)
6.0ms	(pow.f64 (sin.f64 kx) 2)
6.0ms	(pow.f64 (sin.f64 ky) 2)

rewrite9.0ms (0.1%)

Algorithm

1×	rewrite-expression-head

Rules

10×	`add-sqr-sqrt_binary64`
9×	`*-un-lft-identity_binary64`
9×	`add-cube-cbrt_binary64`
7×	`pow1_binary64`
6×	`pow-unpow_binary64`

Counts

4 → 65

Calls

4 calls:

4.0ms	(pow.f64 (sin.f64 kx) 2)
2.0ms	(pow.f64 (sin.f64 ky) 2)
1.0ms	(sqrt.f64 (fma.f64 (pow.f64 (/.f64 (*.f64 2 l) Om) 2) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)) 1))
0.0ms	(fma.f64 (pow.f64 (/.f64 (*.f64 2 l) Om) 2) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)) 1)

simplify288.0ms (1.8%)

Algorithm

1×	egg-herbie

Rules

527×	`cancel-sign-sub-inv_binary64`
299×	`unswap-sqr_binary64`
244×	`fma-neg_binary64`
216×	`fma-def_binary64`
206×	`associate-r_binary64`

Counts

169 → 161

Iterations

Useful iterations: 2 (0.0ms)

Iter	Nodes	Cost
0	348	7820
1	970	7080
2	3540	7073
3	5110	7073

prune399.0ms (2.5%)

Pruning

3 alts after pruning (2 fresh and 1 done)

	Pruned	Kept	Total
New	159	2	161
Fresh	0	0	0
Picked	0	1	1
Done	0	0	0
Total	159	3	162

Error

Counts

162 → 3

Alt Table

Status	Error	Program
	29.7b	(sqrt.f64 (+.f64 1/2 (/.f64 1/2 (fma.f64 1/2 (.f64 Om (sqrt.f64 (/.f64 1/4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))) (/.f64 (sqrt.f64 (.f64 4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2))))) Om)))))
✓	1.0b	(sqrt.f64 (+.f64 1/2 (/.f64 1/2 (sqrt.f64 (fma.f64 (pow.f64 (/.f64 (*.f64 2 l) Om) 2) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)) 1)))))
▶	22.0b	(sqrt.f64 (+.f64 1/2 (/.f64 1/2 1)))

Compiler

Compiled 7831 to 5595 computations (28.6% saved)

localize8.0ms (0%)

prune8.0ms (0.1%)

Pruning

3 alts after pruning (1 fresh and 2 done)

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

Error

Counts

3 → 3

Alt Table

Status	Error	Program
▶	29.7b	(sqrt.f64 (+.f64 1/2 (/.f64 1/2 (fma.f64 1/2 (.f64 Om (sqrt.f64 (/.f64 1/4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))) (/.f64 (sqrt.f64 (.f64 4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2))))) Om)))))
✓	1.0b	(sqrt.f64 (+.f64 1/2 (/.f64 1/2 (sqrt.f64 (fma.f64 (pow.f64 (/.f64 (*.f64 2 l) Om) 2) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)) 1)))))
✓	22.0b	(sqrt.f64 (+.f64 1/2 (/.f64 1/2 1)))

Compiler

Compiled 132 to 102 computations (22.7% saved)

localize52.0ms (0.3%)

Local error

Found 4 expressions with local error:

New	Error	Program
✓	0.7b	(/.f64 1/2 (fma.f64 1/2 (.f64 Om (sqrt.f64 (/.f64 1/4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))) (/.f64 (sqrt.f64 (.f64 4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2))))) Om)))
✓	3.5b	(.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))
✓	27.1b	(sqrt.f64 (.f64 4 (.f64 (*.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))
✓	27.2b	(sqrt.f64 (/.f64 1/4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))

series2.4s (14.8%)

Counts

4 → 156

Calls

4 calls:

1.8s	(/.f64 1/2 (fma.f64 1/2 (.f64 Om (sqrt.f64 (/.f64 1/4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))) (/.f64 (sqrt.f64 (.f64 4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2))))) Om)))
378.0ms	(sqrt.f64 (/.f64 1/4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))
153.0ms	(sqrt.f64 (.f64 4 (.f64 (*.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))
86.0ms	(.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))

rewrite134.0ms (0.8%)

Algorithm

1×	rewrite-expression-head

Rules

20×	`pow1_binary64`
19×	`*-un-lft-identity_binary64`
17×	`add-sqr-sqrt_binary64`
15×	`add-cube-cbrt_binary64`
12×	`times-frac_binary64`

Counts

4 → 95

Calls

4 calls:

35.0ms	(.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))
34.0ms	(sqrt.f64 (.f64 4 (.f64 (*.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))
31.0ms	(sqrt.f64 (/.f64 1/4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))
4.0ms	(/.f64 1/2 (fma.f64 1/2 (.f64 Om (sqrt.f64 (/.f64 1/4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))) (/.f64 (sqrt.f64 (.f64 4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2))))) Om)))

simplify711.0ms (4.4%)

Algorithm

1×	egg-herbie

Rules

895×	`times-frac_binary64`
515×	`fma-def_binary64`
223×	`*-commutative_binary64`
106×	`cancel-sign-sub-inv_binary64`
105×	`+-commutative_binary64`

Counts

251 → 248

Iterations

Useful iterations: 2 (0.0ms)

Iter	Nodes	Cost
0	684	14415
1	2199	13387
2	4902	13256
3	5011	13256

prune965.0ms (6%)

Pruning

3 alts after pruning (0 fresh and 3 done)

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

Error

Counts

306 → 3

Alt Table

Status	Error	Program
✓	29.7b	(sqrt.f64 (+.f64 1/2 (/.f64 1/2 (fma.f64 1/2 (.f64 Om (sqrt.f64 (/.f64 1/4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)))))) (/.f64 (sqrt.f64 (.f64 4 (.f64 (.f64 l l) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2))))) Om)))))
✓	1.0b	(sqrt.f64 (+.f64 1/2 (/.f64 1/2 (sqrt.f64 (fma.f64 (pow.f64 (/.f64 (*.f64 2 l) Om) 2) (+.f64 (pow.f64 (sin.f64 kx) 2) (pow.f64 (sin.f64 ky) 2)) 1)))))
✓	22.0b	(sqrt.f64 (+.f64 1/2 (/.f64 1/2 1)))

Compiler

Compiled 21880 to 13147 computations (39.9% saved)

regimes2.0s (12.3%)

Accuracy

Total 0.6b remaining (60.7%)

Threshold costs 0.6b (60.7%)

Counts

49 → 1

Compiler

Compiled 34722 to 26610 computations (23.4% saved)

simplify3.0ms (0%)

Algorithm

1×	egg-herbie

Rules

2×	`+-commutative_binary64`
1×	`*-commutative_binary64`
1×	`1-exp_binary64`

Iterations

Useful iterations: 0 (0.0ms)

Iter	Nodes	Cost
0	21	37
1	26	37
2	25	37

end111.0ms (0.7%)

Remove: (sort kx ky)
Compiler: Compiled 197 to 164 computations (16.8% saved)

Details

analyze1.5s (9.4%)

sample6.8s (42%)

simplify85.0ms (0.5%)

prune6.0ms (0%)

localize18.0ms (0.1%)

series647.0ms (4%)

rewrite9.0ms (0.1%)

simplify288.0ms (1.8%)

prune399.0ms (2.5%)

localize8.0ms (0%)

prune8.0ms (0.1%)

localize52.0ms (0.3%)

series2.4s (14.8%)

rewrite134.0ms (0.8%)

simplify711.0ms (4.4%)

prune965.0ms (6%)

regimes2.0s (12.3%)

simplify3.0ms (0%)

end111.0ms (0.7%)

Profiling