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.014 * * * * [points]: Setting MPFR precision to 64 0.017 * * * * [points]: Setting MPFR precision to 320 0.020 * * * * [points]: Computing exacts on every 2 of 256 points to ramp up precision 0.022 * * * * [points]: Setting MPFR precision to 64 0.031 * * * * [points]: Setting MPFR precision to 320 0.040 * * * * [points]: Computing exacts for 256 points 0.045 * * * * [points]: Setting MPFR precision to 64 0.071 * * * * [points]: Setting MPFR precision to 320 0.099 * * * * [points]: Filtering points with unrepresentable outputs 0.100 * * * * [points]: Sampled 256 points with exact outputs 0.100 * * * [progress]: [2/2] Setting up program. 0.126 * [progress]: [Phase 2 of 3] Improving. 0.126 * * * * [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.127 * [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.127 * * [simplify]: iteration 1: (14 enodes) 0.153 * * [simplify]: iteration 2: (29 enodes) 0.158 * * [simplify]: iteration 3: (49 enodes) 0.168 * * [simplify]: iteration 4: (121 enodes) 0.199 * * [simplify]: iteration 5: (476 enodes) 0.517 * * [simplify]: Extracting #0: cost 1 inf + 0 0.517 * * [simplify]: Extracting #1: cost 3 inf + 0 0.517 * * [simplify]: Extracting #2: cost 5 inf + 0 0.517 * * [simplify]: Extracting #3: cost 89 inf + 1 0.518 * * [simplify]: Extracting #4: cost 456 inf + 2 0.521 * * [simplify]: Extracting #5: cost 590 inf + 29084 0.533 * * [simplify]: Extracting #6: cost 512 inf + 285351 0.569 * * [simplify]: Extracting #7: cost 159 inf + 842146 0.664 * * [simplify]: Extracting #8: cost 11 inf + 1120113 0.775 * * [simplify]: Extracting #9: cost 0 inf + 1147578 0.885 * [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.910 * * [progress]: iteration 1 / 4 0.910 * * * [progress]: picking best candidate 0.933 * * * * [pick]: Picked #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> 0.933 * * * [progress]: localizing error 1.213 * * * [progress]: generating rewritten candidates 1.213 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2) 1.216 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2 1) 1.217 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2) 1.217 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2 1 2 1 1) 1.226 * * * [progress]: generating series expansions 1.226 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2) 1.226 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2 1) 1.226 * * * * [progress]: [ 3 / 4 ] generating series at (2 2) 1.226 * * * * [progress]: [ 4 / 4 ] generating series at (2 2 1 2 1 1) 1.226 * * * [progress]: simplifying candidates 1.226 * * * * [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.226 * * * * [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.226 * * * * [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.226 * * * * [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.227 * * * * [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.227 * * * * [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.227 * * * * [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.227 * [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.227 * * [simplify]: iteration 1: (18 enodes) 1.235 * * [simplify]: iteration 2: (38 enodes) 1.241 * * [simplify]: iteration 3: (98 enodes) 1.271 * * [simplify]: iteration 4: (439 enodes) 1.643 * * [simplify]: Extracting #0: cost 4 inf + 0 1.643 * * [simplify]: Extracting #1: cost 74 inf + 0 1.645 * * [simplify]: Extracting #2: cost 416 inf + 404 1.652 * * [simplify]: Extracting #3: cost 611 inf + 60189 1.689 * * [simplify]: Extracting #4: cost 345 inf + 574431 1.745 * * [simplify]: Extracting #5: cost 70 inf + 1032645 1.845 * * [simplify]: Extracting #6: cost 0 inf + 1142765 1.951 * [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.951 * * * [progress]: adding candidates to table 2.178 * * [progress]: iteration 2 / 4 2.178 * * * [progress]: picking best candidate 2.203 * * * * [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.204 * * * [progress]: localizing error 2.461 * * * [progress]: generating rewritten candidates 2.461 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2 2) 2.467 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2) 2.476 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2 1 2 2 1) 2.476 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2) 2.477 * * * [progress]: generating series expansions 2.477 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2 2) 2.477 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2) 2.477 * * * * [progress]: [ 3 / 4 ] generating series at (2 2 1 2 2 1) 2.477 * * * * [progress]: [ 4 / 4 ] generating series at (2 2) 2.477 * * * [progress]: simplifying candidates 2.477 * * * * [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.477 * * * * [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.477 * * * * [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.477 * * * * [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.477 * * * * [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.477 * * * * [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.478 * [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.478 * * [simplify]: iteration 1: (16 enodes) 2.484 * * [simplify]: iteration 2: (29 enodes) 2.488 * * [simplify]: iteration 3: (40 enodes) 2.494 * * [simplify]: iteration 4: (52 enodes) 2.501 * * [simplify]: iteration 5: (68 enodes) 2.510 * * [simplify]: iteration 6: (73 enodes) 2.519 * * [simplify]: Extracting #0: cost 2 inf + 0 2.519 * * [simplify]: Extracting #1: cost 6 inf + 0 2.519 * * [simplify]: Extracting #2: cost 9 inf + 1 2.519 * * [simplify]: Extracting #3: cost 18 inf + 4 2.519 * * [simplify]: Extracting #4: cost 19 inf + 325 2.519 * * [simplify]: Extracting #5: cost 14 inf + 1612 2.520 * * [simplify]: Extracting #6: cost 12 inf + 3016 2.520 * * [simplify]: Extracting #7: cost 3 inf + 14359 2.521 * * [simplify]: Extracting #8: cost 0 inf + 19008 2.522 * [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.522 * * * [progress]: adding candidates to table 2.655 * * [progress]: iteration 3 / 4 2.655 * * * [progress]: picking best candidate 2.668 * * * * [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.668 * * * [progress]: localizing error 3.008 * * * [progress]: generating rewritten candidates 3.008 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2 2 1) 3.013 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2 2 1 1) 3.014 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2 1 2 2) 3.021 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2) 3.022 * * * [progress]: generating series expansions 3.022 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2 2 1) 3.022 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2 2 1 1) 3.022 * * * * [progress]: [ 3 / 4 ] generating series at (2 2 1 2 2) 3.022 * * * * [progress]: [ 4 / 4 ] generating series at (2 2) 3.022 * * * [progress]: simplifying candidates 3.022 * * * * [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.022 * * * * [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.022 * * * * [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.022 * * * * [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.022 * * * * [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.023 * [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.023 * * [simplify]: iteration 1: (8 enodes) 3.026 * * [simplify]: iteration 2: (14 enodes) 3.030 * * [simplify]: Extracting #0: cost 1 inf + 0 3.030 * * [simplify]: Extracting #1: cost 3 inf + 0 3.030 * * [simplify]: Extracting #2: cost 4 inf + 1 3.030 * * [simplify]: Extracting #3: cost 4 inf + 2 3.030 * * [simplify]: Extracting #4: cost 6 inf + 2 3.030 * * [simplify]: Extracting #5: cost 0 inf + 4012 3.031 * [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.031 * * * [progress]: adding candidates to table 3.201 * [progress]: [Phase 3 of 3] Extracting. 3.201 * * [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.203 * * * [regime-changes]: Trying 2 branch expressions: (im re) 3.203 * * * * [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.311 * * * * [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.428 * * * [regime]: Found split indices: #