0.001 * [progress]: [Phase 1 of 3] Setting up. 0.001 * * * [progress]: [1/2] Preparing points 0.001 * * * * [points]: Sampling 256 additional inputs, on iter 0 have 0 / 256 0.002 * * * * [points]: Computing exacts on every 16 of 256 points to ramp up precision 0.004 * * * * [points]: Setting MPFR precision to 64 0.006 * * * * [points]: Setting MPFR precision to 320 0.007 * * * * [points]: Computing exacts on every 8 of 256 points to ramp up precision 0.009 * * * * [points]: Setting MPFR precision to 64 0.011 * * * * [points]: Setting MPFR precision to 320 0.012 * * * * [points]: Computing exacts on every 4 of 256 points to ramp up precision 0.015 * * * * [points]: Setting MPFR precision to 64 0.018 * * * * [points]: Setting MPFR precision to 320 0.022 * * * * [points]: Computing exacts on every 2 of 256 points to ramp up precision 0.024 * * * * [points]: Setting MPFR precision to 64 0.029 * * * * [points]: Setting MPFR precision to 320 0.034 * * * * [points]: Computing exacts for 256 points 0.036 * * * * [points]: Setting MPFR precision to 64 0.056 * * * * [points]: Setting MPFR precision to 320 0.086 * * * * [points]: Filtering points with unrepresentable outputs 0.087 * * * * [points]: Sampled 256 points with exact outputs 0.087 * * * [progress]: [2/2] Setting up program. 0.108 * [progress]: [Phase 2 of 3] Improving. 0.108 * * * * [progress]: [ 1 / 1 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> 0.108 * [simplify]: Simplifying (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))) 0.108 * * [simplify]: iteration 1: (14 enodes) 0.112 * * [simplify]: iteration 2: (29 enodes) 0.117 * * [simplify]: iteration 3: (49 enodes) 0.144 * * [simplify]: iteration 4: (121 enodes) 0.185 * * [simplify]: iteration 5: (476 enodes) 0.583 * * [simplify]: Extracting #0: cost 1 inf + 0 0.583 * * [simplify]: Extracting #1: cost 3 inf + 0 0.583 * * [simplify]: Extracting #2: cost 5 inf + 0 0.583 * * [simplify]: Extracting #3: cost 89 inf + 1 0.585 * * [simplify]: Extracting #4: cost 456 inf + 2 0.587 * * [simplify]: Extracting #5: cost 590 inf + 29084 0.601 * * [simplify]: Extracting #6: cost 512 inf + 285351 0.662 * * [simplify]: Extracting #7: cost 159 inf + 842146 0.750 * * [simplify]: Extracting #8: cost 11 inf + 1120113 0.863 * * [simplify]: Extracting #9: cost 0 inf + 1147578 0.978 * [simplify]: Simplified to (*.p16 (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 im im) (*.p16 re re))) re))) (real->posit16 0.5)) 1.006 * * [progress]: iteration 1 / 4 1.006 * * * [progress]: picking best candidate 1.031 * * * * [pick]: Picked #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> 1.031 * * * [progress]: localizing error 1.307 * * * [progress]: generating rewritten candidates 1.307 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2) 1.311 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2 1) 1.311 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2) 1.312 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2 1 2 1 1) 1.320 * * * [progress]: generating series expansions 1.320 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2) 1.320 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2 1) 1.321 * * * * [progress]: [ 3 / 4 ] generating series at (2 2) 1.321 * * * * [progress]: [ 4 / 4 ] generating series at (2 2 1 2 1 1) 1.321 * * * [progress]: simplifying candidates 1.321 * * * * [progress]: [ 1 / 7 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) (neg.p16 re))))))> 1.321 * * * * [progress]: [ 2 / 7 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (-.p16 (*.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im)))) (*.p16 re re)) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))))> 1.321 * * * * [progress]: [ 3 / 7 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 im im) (*.p16 re re))) re)))))> 1.321 * * * * [progress]: [ 4 / 7 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> 1.321 * * * * [progress]: [ 5 / 7 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> 1.321 * * * * [progress]: [ 6 / 7 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> 1.321 * * * * [progress]: [ 7 / 7 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> 1.322 * [simplify]: Simplifying (neg.p16 re), (-.p16 (*.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im)))) (*.p16 re re)), (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))) 1.322 * * [simplify]: iteration 1: (18 enodes) 1.330 * * [simplify]: iteration 2: (38 enodes) 1.343 * * [simplify]: iteration 3: (98 enodes) 1.396 * * [simplify]: iteration 4: (439 enodes) 1.800 * * [simplify]: Extracting #0: cost 4 inf + 0 1.801 * * [simplify]: Extracting #1: cost 74 inf + 0 1.803 * * [simplify]: Extracting #2: cost 416 inf + 404 1.810 * * [simplify]: Extracting #3: cost 611 inf + 60189 1.855 * * [simplify]: Extracting #4: cost 345 inf + 574431 1.946 * * [simplify]: Extracting #5: cost 70 inf + 1032645 2.037 * * [simplify]: Extracting #6: cost 0 inf + 1142765 2.131 * [simplify]: Simplified to (neg.p16 re), (*.p16 im im), (+.p16 re (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im)))), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) (real->posit16 2.0)))), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) (real->posit16 2.0)))), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) (real->posit16 2.0)))), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) (real->posit16 2.0)))) 2.132 * * * [progress]: adding candidates to table 2.385 * * [progress]: iteration 2 / 4 2.385 * * * [progress]: picking best candidate 2.411 * * * * [pick]: Picked #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))))> 2.411 * * * [progress]: localizing error 2.709 * * * [progress]: generating rewritten candidates 2.709 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2 2) 2.712 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2) 2.718 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2 1 2 2 1) 2.719 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2) 2.719 * * * [progress]: generating series expansions 2.719 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2 2) 2.719 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2) 2.719 * * * * [progress]: [ 3 / 4 ] generating series at (2 2 1 2 2 1) 2.719 * * * * [progress]: [ 4 / 4 ] generating series at (2 2) 2.719 * * * [progress]: simplifying candidates 2.719 * * * * [progress]: [ 1 / 6 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 re (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im)))))))))> 2.719 * * * * [progress]: [ 2 / 6 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))> 2.719 * * * * [progress]: [ 3 / 6 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))))> 2.719 * * * * [progress]: [ 4 / 6 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))))> 2.719 * * * * [progress]: [ 5 / 6 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))))> 2.719 * * * * [progress]: [ 6 / 6 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))))> 2.720 * [simplify]: Simplifying (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))), (*.p16 (real->posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))) 2.720 * * [simplify]: iteration 1: (16 enodes) 2.724 * * [simplify]: iteration 2: (29 enodes) 2.729 * * [simplify]: iteration 3: (40 enodes) 2.735 * * [simplify]: iteration 4: (52 enodes) 2.743 * * [simplify]: iteration 5: (68 enodes) 2.765 * * [simplify]: iteration 6: (73 enodes) 2.785 * * [simplify]: Extracting #0: cost 2 inf + 0 2.785 * * [simplify]: Extracting #1: cost 6 inf + 0 2.785 * * [simplify]: Extracting #2: cost 9 inf + 1 2.785 * * [simplify]: Extracting #3: cost 18 inf + 4 2.786 * * [simplify]: Extracting #4: cost 19 inf + 325 2.786 * * [simplify]: Extracting #5: cost 14 inf + 1612 2.786 * * [simplify]: Extracting #6: cost 12 inf + 3016 2.788 * * [simplify]: Extracting #7: cost 3 inf + 14359 2.790 * * [simplify]: Extracting #8: cost 0 inf + 19008 2.791 * [simplify]: Simplified to (/.p16 (+.p16 re (sqrt.p16 (+.p16 (*.p16 im im) (*.p16 re re)))) im), (*.p16 (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 re (sqrt.p16 (+.p16 (*.p16 im im) (*.p16 re re))))))) (real->posit16 0.5)), (*.p16 (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 re (sqrt.p16 (+.p16 (*.p16 im im) (*.p16 re re))))))) (real->posit16 0.5)), (*.p16 (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 re (sqrt.p16 (+.p16 (*.p16 im im) (*.p16 re re))))))) (real->posit16 0.5)), (*.p16 (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 re (sqrt.p16 (+.p16 (*.p16 im im) (*.p16 re re))))))) (real->posit16 0.5)) 2.792 * * * [progress]: adding candidates to table 2.950 * * [progress]: iteration 3 / 4 2.950 * * * [progress]: picking best candidate 2.973 * * * * [pick]: Picked #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))> 2.973 * * * [progress]: localizing error 3.295 * * * [progress]: generating rewritten candidates 3.295 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2 2 1) 3.301 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2 2 1 1) 3.301 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2 1 2 2) 3.308 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2) 3.309 * * * [progress]: generating series expansions 3.309 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2 2 1) 3.309 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2 2 1 1) 3.309 * * * * [progress]: [ 3 / 4 ] generating series at (2 2 1 2 2) 3.309 * * * * [progress]: [ 4 / 4 ] generating series at (2 2) 3.309 * * * [progress]: simplifying candidates 3.309 * * * * [progress]: [ 1 / 5 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 im (/.p16 (+.p16 re (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im)))) im))))))> 3.309 * * * * [progress]: [ 2 / 5 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))> 3.309 * * * * [progress]: [ 3 / 5 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))> 3.309 * * * * [progress]: [ 4 / 5 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))> 3.309 * * * * [progress]: [ 5 / 5 ] simplifiying candidate #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))> 3.310 * [simplify]: Simplifying (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im), (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im), (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im), (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im) 3.310 * * [simplify]: iteration 1: (8 enodes) 3.313 * * [simplify]: iteration 2: (14 enodes) 3.318 * * [simplify]: Extracting #0: cost 1 inf + 0 3.318 * * [simplify]: Extracting #1: cost 3 inf + 0 3.318 * * [simplify]: Extracting #2: cost 4 inf + 1 3.318 * * [simplify]: Extracting #3: cost 4 inf + 2 3.318 * * [simplify]: Extracting #4: cost 6 inf + 2 3.319 * * [simplify]: Extracting #5: cost 0 inf + 4012 3.319 * [simplify]: Simplified to (/.p16 (+.p16 re (sqrt.p16 (+.p16 (*.p16 im im) (*.p16 re re)))) im), (/.p16 (+.p16 re (sqrt.p16 (+.p16 (*.p16 im im) (*.p16 re re)))) im), (/.p16 (+.p16 re (sqrt.p16 (+.p16 (*.p16 im im) (*.p16 re re)))) im), (/.p16 (+.p16 re (sqrt.p16 (+.p16 (*.p16 im im) (*.p16 re re)))) im) 3.320 * * * [progress]: adding candidates to table 3.485 * [progress]: [Phase 3 of 3] Extracting. 3.485 * * [regime]: Finding splitpoints for: (#posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))> #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))))>) 3.486 * * * [regime-changes]: Trying 2 branch expressions: (im re) 3.487 * * * * [regimes]: Trying to branch on im from (#posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))> #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))))>) 3.637 * * * * [regimes]: Trying to branch on re from (#posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))> #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (/.p16 (*.p16 im im) (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re))))))>) 3.786 * * * [regime]: Found split indices: #