0.002 * [progress]: [Phase 1 of 3] Setting up. 0.003 * * * [progress]: [1/2] Preparing points 0.003 * * * * [points]: Sampling 256 additional inputs, on iter 0 have 0 / 256 0.004 * * * * [points]: Computing exacts on every 16 of 256 points to ramp up precision 0.008 * * * * [points]: Setting MPFR precision to 64 0.010 * * * * [points]: Setting MPFR precision to 320 0.011 * * * * [points]: Computing exacts on every 8 of 256 points to ramp up precision 0.015 * * * * [points]: Setting MPFR precision to 64 0.019 * * * * [points]: Setting MPFR precision to 320 0.022 * * * * [points]: Computing exacts on every 4 of 256 points to ramp up precision 0.026 * * * * [points]: Setting MPFR precision to 64 0.031 * * * * [points]: Setting MPFR precision to 320 0.037 * * * * [points]: Computing exacts on every 2 of 256 points to ramp up precision 0.042 * * * * [points]: Setting MPFR precision to 64 0.051 * * * * [points]: Setting MPFR precision to 320 0.060 * * * * [points]: Computing exacts for 256 points 0.065 * * * * [points]: Setting MPFR precision to 64 0.091 * * * * [points]: Setting MPFR precision to 320 0.118 * * * * [points]: Filtering points with unrepresentable outputs 0.119 * * * * [points]: Sampled 256 points with exact outputs 0.119 * * * [progress]: [2/2] Setting up program. 0.143 * [progress]: [Phase 2 of 3] Improving. 0.143 * * * * [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.143 * [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.144 * * [simplify]: iteration 1: (14 enodes) 0.150 * * [simplify]: iteration 2: (29 enodes) 0.182 * * [simplify]: iteration 3: (49 enodes) 0.203 * * [simplify]: iteration 4: (121 enodes) 0.234 * * [simplify]: iteration 5: (476 enodes) 0.633 * * [simplify]: Extracting #0: cost 1 inf + 0 0.633 * * [simplify]: Extracting #1: cost 3 inf + 0 0.633 * * [simplify]: Extracting #2: cost 5 inf + 0 0.634 * * [simplify]: Extracting #3: cost 89 inf + 1 0.636 * * [simplify]: Extracting #4: cost 456 inf + 2 0.641 * * [simplify]: Extracting #5: cost 590 inf + 29084 0.661 * * [simplify]: Extracting #6: cost 512 inf + 285351 0.696 * * [simplify]: Extracting #7: cost 159 inf + 842146 0.777 * * [simplify]: Extracting #8: cost 11 inf + 1120113 0.841 * * [simplify]: Extracting #9: cost 0 inf + 1147578 0.911 * [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.933 * * [progress]: iteration 1 / 4 0.933 * * * [progress]: picking best candidate 0.958 * * * * [pick]: Picked #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> 0.958 * * * [progress]: localizing error 1.141 * * * [progress]: generating rewritten candidates 1.141 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2) 1.142 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2 1) 1.143 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2) 1.143 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2 1 2 1 1) 1.147 * * * [progress]: generating series expansions 1.147 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2) 1.147 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2 1) 1.147 * * * * [progress]: [ 3 / 4 ] generating series at (2 2) 1.147 * * * * [progress]: [ 4 / 4 ] generating series at (2 2 1 2 1 1) 1.147 * * * [progress]: simplifying candidates 1.147 * * * * [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.147 * * * * [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.147 * * * * [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.147 * * * * [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.147 * * * * [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.147 * * * * [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.147 * * * * [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.148 * [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.148 * * [simplify]: iteration 1: (18 enodes) 1.152 * * [simplify]: iteration 2: (38 enodes) 1.158 * * [simplify]: iteration 3: (98 enodes) 1.184 * * [simplify]: iteration 4: (439 enodes) 1.433 * * [simplify]: Extracting #0: cost 4 inf + 0 1.433 * * [simplify]: Extracting #1: cost 74 inf + 0 1.434 * * [simplify]: Extracting #2: cost 416 inf + 404 1.437 * * [simplify]: Extracting #3: cost 611 inf + 60189 1.460 * * [simplify]: Extracting #4: cost 345 inf + 574431 1.508 * * [simplify]: Extracting #5: cost 70 inf + 1032645 1.566 * * [simplify]: Extracting #6: cost 0 inf + 1142765 1.623 * [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.623 * * * [progress]: adding candidates to table 1.784 * * [progress]: iteration 2 / 4 1.784 * * * [progress]: picking best candidate 1.810 * * * * [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.810 * * * [progress]: localizing error 2.056 * * * [progress]: generating rewritten candidates 2.056 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2 2) 2.072 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2) 2.081 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2 1 2 2 1) 2.082 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2) 2.082 * * * [progress]: generating series expansions 2.082 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2 2) 2.082 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2) 2.082 * * * * [progress]: [ 3 / 4 ] generating series at (2 2 1 2 2 1) 2.082 * * * * [progress]: [ 4 / 4 ] generating series at (2 2) 2.082 * * * [progress]: simplifying candidates 2.082 * * * * [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.083 * * * * [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.083 * * * * [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.083 * * * * [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.083 * * * * [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.083 * * * * [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.083 * [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.083 * * [simplify]: iteration 1: (16 enodes) 2.090 * * [simplify]: iteration 2: (29 enodes) 2.099 * * [simplify]: iteration 3: (40 enodes) 2.110 * * [simplify]: iteration 4: (52 enodes) 2.124 * * [simplify]: iteration 5: (68 enodes) 2.143 * * [simplify]: iteration 6: (73 enodes) 2.161 * * [simplify]: Extracting #0: cost 2 inf + 0 2.161 * * [simplify]: Extracting #1: cost 6 inf + 0 2.161 * * [simplify]: Extracting #2: cost 9 inf + 1 2.161 * * [simplify]: Extracting #3: cost 18 inf + 4 2.162 * * [simplify]: Extracting #4: cost 19 inf + 325 2.162 * * [simplify]: Extracting #5: cost 14 inf + 1612 2.162 * * [simplify]: Extracting #6: cost 12 inf + 3016 2.164 * * [simplify]: Extracting #7: cost 3 inf + 14359 2.166 * * [simplify]: Extracting #8: cost 0 inf + 19008 2.168 * [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.168 * * * [progress]: adding candidates to table 2.288 * * [progress]: iteration 3 / 4 2.288 * * * [progress]: picking best candidate 2.308 * * * * [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.308 * * * [progress]: localizing error 2.640 * * * [progress]: generating rewritten candidates 2.641 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2 2 1) 2.645 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2 2 1 1) 2.646 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2 1 2 2) 2.652 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2) 2.653 * * * [progress]: generating series expansions 2.653 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2 2 1) 2.653 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2 2 1 1) 2.653 * * * * [progress]: [ 3 / 4 ] generating series at (2 2 1 2 2) 2.653 * * * * [progress]: [ 4 / 4 ] generating series at (2 2) 2.653 * * * [progress]: simplifying candidates 2.653 * * * * [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.654 * * * * [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.654 * * * * [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.654 * * * * [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.654 * * * * [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.654 * [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.654 * * [simplify]: iteration 1: (8 enodes) 2.657 * * [simplify]: iteration 2: (14 enodes) 2.661 * * [simplify]: Extracting #0: cost 1 inf + 0 2.661 * * [simplify]: Extracting #1: cost 3 inf + 0 2.661 * * [simplify]: Extracting #2: cost 4 inf + 1 2.661 * * [simplify]: Extracting #3: cost 4 inf + 2 2.661 * * [simplify]: Extracting #4: cost 6 inf + 2 2.662 * * [simplify]: Extracting #5: cost 0 inf + 4012 2.662 * [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.662 * * * [progress]: adding candidates to table 2.828 * [progress]: [Phase 3 of 3] Extracting. 2.828 * * [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))))))>) 2.830 * * * [regime-changes]: Trying 2 branch expressions: (im re) 2.830 * * * * [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))))))>) 2.973 * * * * [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.112 * * * [regime]: Found split indices: #