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.005 * * * * [points]: Setting MPFR precision to 320 0.006 * * * * [points]: Computing exacts on every 8 of 256 points to ramp up precision 0.008 * * * * [points]: Setting MPFR precision to 64 0.010 * * * * [points]: Setting MPFR precision to 320 0.012 * * * * [points]: Computing exacts on every 4 of 256 points to ramp up precision 0.013 * * * * [points]: Setting MPFR precision to 64 0.016 * * * * [points]: Setting MPFR precision to 320 0.019 * * * * [points]: Computing exacts on every 2 of 256 points to ramp up precision 0.021 * * * * [points]: Setting MPFR precision to 64 0.026 * * * * [points]: Setting MPFR precision to 320 0.032 * * * * [points]: Computing exacts for 256 points 0.034 * * * * [points]: Setting MPFR precision to 64 0.050 * * * * [points]: Setting MPFR precision to 320 0.065 * * * * [points]: Filtering points with unrepresentable outputs 0.066 * * * * [points]: Sampled 256 points with exact outputs 0.066 * * * [progress]: [2/2] Setting up program. 0.080 * [progress]: [Phase 2 of 3] Improving. 0.080 * * * * [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.081 * [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.081 * * [simplify]: iteration 1: (14 enodes) 0.084 * * [simplify]: iteration 2: (29 enodes) 0.090 * * [simplify]: iteration 3: (49 enodes) 0.099 * * [simplify]: iteration 4: (121 enodes) 0.145 * * [simplify]: iteration 5: (476 enodes) 0.440 * * [simplify]: Extracting #0: cost 1 inf + 0 0.440 * * [simplify]: Extracting #1: cost 3 inf + 0 0.440 * * [simplify]: Extracting #2: cost 5 inf + 0 0.444 * * [simplify]: Extracting #3: cost 89 inf + 1 0.445 * * [simplify]: Extracting #4: cost 456 inf + 2 0.448 * * [simplify]: Extracting #5: cost 590 inf + 29084 0.460 * * [simplify]: Extracting #6: cost 512 inf + 285351 0.495 * * [simplify]: Extracting #7: cost 159 inf + 842146 0.548 * * [simplify]: Extracting #8: cost 11 inf + 1120113 0.605 * * [simplify]: Extracting #9: cost 0 inf + 1147578 0.668 * [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)) 0.684 * * [progress]: iteration 1 / 4 0.684 * * * [progress]: picking best candidate 0.698 * * * * [pick]: Picked #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> 0.698 * * * [progress]: localizing error 0.899 * * * [progress]: generating rewritten candidates 0.899 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2) 0.901 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2 1) 0.901 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2) 0.901 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2 1 2 1 1) 0.905 * * * [progress]: generating series expansions 0.905 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2) 0.905 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2 1) 0.906 * * * * [progress]: [ 3 / 4 ] generating series at (2 2) 0.906 * * * * [progress]: [ 4 / 4 ] generating series at (2 2 1 2 1 1) 0.906 * * * [progress]: simplifying candidates 0.906 * * * * [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))))))> 0.906 * * * * [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))))))> 0.906 * * * * [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)))))> 0.906 * * * * [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)))))> 0.906 * * * * [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)))))> 0.906 * * * * [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)))))> 0.906 * * * * [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)))))> 0.906 * [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)))) 0.906 * * [simplify]: iteration 1: (18 enodes) 0.911 * * [simplify]: iteration 2: (38 enodes) 0.918 * * [simplify]: iteration 3: (98 enodes) 0.943 * * [simplify]: iteration 4: (439 enodes) 1.270 * * [simplify]: Extracting #0: cost 4 inf + 0 1.271 * * [simplify]: Extracting #1: cost 74 inf + 0 1.279 * * [simplify]: Extracting #2: cost 416 inf + 404 1.286 * * [simplify]: Extracting #3: cost 611 inf + 60189 1.323 * * [simplify]: Extracting #4: cost 345 inf + 574431 1.384 * * [simplify]: Extracting #5: cost 70 inf + 1032645 1.460 * * [simplify]: Extracting #6: cost 0 inf + 1142765 1.517 * [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)))) 1.517 * * * [progress]: adding candidates to table 1.682 * * [progress]: iteration 2 / 4 1.682 * * * [progress]: picking best candidate 1.696 * * * * [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))))))> 1.696 * * * [progress]: localizing error 1.902 * * * [progress]: generating rewritten candidates 1.902 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2 2) 1.905 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2) 1.909 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2 1 2 2 1) 1.910 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2) 1.910 * * * [progress]: generating series expansions 1.910 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2 2) 1.910 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2) 1.910 * * * * [progress]: [ 3 / 4 ] generating series at (2 2 1 2 2 1) 1.910 * * * * [progress]: [ 4 / 4 ] generating series at (2 2) 1.910 * * * [progress]: simplifying candidates 1.910 * * * * [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)))))))))> 1.910 * * * * [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))))))> 1.910 * * * * [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))))))> 1.910 * * * * [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))))))> 1.910 * * * * [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))))))> 1.910 * * * * [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))))))> 1.910 * [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))))) 1.911 * * [simplify]: iteration 1: (16 enodes) 1.914 * * [simplify]: iteration 2: (29 enodes) 1.919 * * [simplify]: iteration 3: (40 enodes) 1.929 * * [simplify]: iteration 4: (52 enodes) 1.943 * * [simplify]: iteration 5: (68 enodes) 1.959 * * [simplify]: iteration 6: (73 enodes) 1.979 * * [simplify]: Extracting #0: cost 2 inf + 0 1.979 * * [simplify]: Extracting #1: cost 6 inf + 0 1.979 * * [simplify]: Extracting #2: cost 9 inf + 1 1.979 * * [simplify]: Extracting #3: cost 18 inf + 4 1.979 * * [simplify]: Extracting #4: cost 19 inf + 325 1.980 * * [simplify]: Extracting #5: cost 14 inf + 1612 1.980 * * [simplify]: Extracting #6: cost 12 inf + 3016 1.981 * * [simplify]: Extracting #7: cost 3 inf + 14359 1.983 * * [simplify]: Extracting #8: cost 0 inf + 19008 1.985 * [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)) 1.986 * * * [progress]: adding candidates to table 2.129 * * [progress]: iteration 3 / 4 2.129 * * * [progress]: picking best candidate 2.149 * * * * [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.149 * * * [progress]: localizing error 2.509 * * * [progress]: generating rewritten candidates 2.509 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2 2 1) 2.515 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2 2 1 1) 2.515 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2 1 2 2) 2.522 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2) 2.522 * * * [progress]: generating series expansions 2.522 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2 2 1) 2.522 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2 2 1 1) 2.522 * * * * [progress]: [ 3 / 4 ] generating series at (2 2 1 2 2) 2.523 * * * * [progress]: [ 4 / 4 ] generating series at (2 2) 2.523 * * * [progress]: simplifying candidates 2.523 * * * * [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))))))> 2.523 * * * * [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))))))> 2.523 * * * * [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))))))> 2.523 * * * * [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))))))> 2.523 * * * * [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))))))> 2.523 * [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) 2.523 * * [simplify]: iteration 1: (8 enodes) 2.527 * * [simplify]: iteration 2: (14 enodes) 2.531 * * [simplify]: Extracting #0: cost 1 inf + 0 2.531 * * [simplify]: Extracting #1: cost 3 inf + 0 2.531 * * [simplify]: Extracting #2: cost 4 inf + 1 2.531 * * [simplify]: Extracting #3: cost 4 inf + 2 2.531 * * [simplify]: Extracting #4: cost 6 inf + 2 2.531 * * [simplify]: Extracting #5: cost 0 inf + 4012 2.531 * [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) 2.532 * * * [progress]: adding candidates to table 2.710 * [progress]: [Phase 3 of 3] Extracting. 2.710 * * [regime]: Finding splitpoints for: (#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))))))> #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 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))>) 2.712 * * * [regime-changes]: Trying 2 branch expressions: (im re) 2.712 * * * * [regimes]: Trying to branch on im from (#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))))))> #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 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))>) 2.894 * * * * [regimes]: Trying to branch on re from (#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))))))> #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 im (/.p16 (+.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re) im))))))>) 3.086 * * * [regime]: Found split indices: #