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.007 * * * * [points]: Setting MPFR precision to 64 0.009 * * * * [points]: Setting MPFR precision to 320 0.011 * * * * [points]: Computing exacts on every 8 of 256 points to ramp up precision 0.016 * * * * [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.032 * * * * [points]: Setting MPFR precision to 320 0.037 * * * * [points]: Computing exacts on every 2 of 256 points to ramp up precision 0.041 * * * * [points]: Setting MPFR precision to 64 0.050 * * * * [points]: Setting MPFR precision to 320 0.059 * * * * [points]: Computing exacts for 256 points 0.063 * * * * [points]: Setting MPFR precision to 64 0.090 * * * * [points]: Setting MPFR precision to 320 0.119 * * * * [points]: Filtering points with unrepresentable outputs 0.120 * * * * [points]: Sampled 256 points with exact outputs 0.120 * * * [progress]: [2/2] Setting up program. 0.147 * [progress]: [Phase 2 of 3] Improving. 0.147 * * * * [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.147 * [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.148 * * [simplify]: iteration 1: (14 enodes) 0.176 * * [simplify]: iteration 2: (29 enodes) 0.188 * * [simplify]: iteration 3: (49 enodes) 0.206 * * [simplify]: iteration 4: (121 enodes) 0.265 * * [simplify]: iteration 5: (476 enodes) 0.700 * * [simplify]: Extracting #0: cost 1 inf + 0 0.700 * * [simplify]: Extracting #1: cost 3 inf + 0 0.700 * * [simplify]: Extracting #2: cost 5 inf + 0 0.700 * * [simplify]: Extracting #3: cost 89 inf + 1 0.701 * * [simplify]: Extracting #4: cost 456 inf + 2 0.704 * * [simplify]: Extracting #5: cost 590 inf + 29084 0.716 * * [simplify]: Extracting #6: cost 512 inf + 285351 0.756 * * [simplify]: Extracting #7: cost 159 inf + 842146 0.808 * * [simplify]: Extracting #8: cost 11 inf + 1120113 0.884 * * [simplify]: Extracting #9: cost 0 inf + 1147578 0.945 * [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.973 * * [progress]: iteration 1 / 4 0.973 * * * [progress]: picking best candidate 1.000 * * * * [pick]: Picked #posit16 0.5) (sqrt.p16 (*.p16 (real->posit16 2.0) (-.p16 (sqrt.p16 (+.p16 (*.p16 re re) (*.p16 im im))) re)))))> 1.000 * * * [progress]: localizing error 1.283 * * * [progress]: generating rewritten candidates 1.283 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2) 1.285 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2 1) 1.285 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2) 1.286 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2 1 2 1 1) 1.290 * * * [progress]: generating series expansions 1.290 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2) 1.290 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2 1) 1.290 * * * * [progress]: [ 3 / 4 ] generating series at (2 2) 1.290 * * * * [progress]: [ 4 / 4 ] generating series at (2 2 1 2 1 1) 1.290 * * * [progress]: simplifying candidates 1.290 * * * * [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.290 * * * * [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.290 * * * * [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.290 * * * * [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.290 * * * * [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.290 * * * * [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.290 * * * * [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.291 * [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.291 * * [simplify]: iteration 1: (18 enodes) 1.295 * * [simplify]: iteration 2: (38 enodes) 1.302 * * [simplify]: iteration 3: (98 enodes) 1.343 * * [simplify]: iteration 4: (439 enodes) 1.675 * * [simplify]: Extracting #0: cost 4 inf + 0 1.675 * * [simplify]: Extracting #1: cost 74 inf + 0 1.676 * * [simplify]: Extracting #2: cost 416 inf + 404 1.680 * * [simplify]: Extracting #3: cost 611 inf + 60189 1.701 * * [simplify]: Extracting #4: cost 345 inf + 574431 1.754 * * [simplify]: Extracting #5: cost 70 inf + 1032645 1.822 * * [simplify]: Extracting #6: cost 0 inf + 1142765 1.879 * [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.880 * * * [progress]: adding candidates to table 2.032 * * [progress]: iteration 2 / 4 2.032 * * * [progress]: picking best candidate 2.046 * * * * [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.046 * * * [progress]: localizing error 2.248 * * * [progress]: generating rewritten candidates 2.248 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2 2) 2.251 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2) 2.255 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2 1 2 2 1) 2.255 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2) 2.256 * * * [progress]: generating series expansions 2.256 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2 2) 2.256 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2) 2.256 * * * * [progress]: [ 3 / 4 ] generating series at (2 2 1 2 2 1) 2.256 * * * * [progress]: [ 4 / 4 ] generating series at (2 2) 2.256 * * * [progress]: simplifying candidates 2.256 * * * * [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.256 * * * * [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.256 * * * * [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.256 * * * * [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.256 * * * * [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.256 * * * * [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.256 * [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.257 * * [simplify]: iteration 1: (16 enodes) 2.260 * * [simplify]: iteration 2: (29 enodes) 2.265 * * [simplify]: iteration 3: (40 enodes) 2.271 * * [simplify]: iteration 4: (52 enodes) 2.278 * * [simplify]: iteration 5: (68 enodes) 2.287 * * [simplify]: iteration 6: (73 enodes) 2.301 * * [simplify]: Extracting #0: cost 2 inf + 0 2.301 * * [simplify]: Extracting #1: cost 6 inf + 0 2.301 * * [simplify]: Extracting #2: cost 9 inf + 1 2.301 * * [simplify]: Extracting #3: cost 18 inf + 4 2.301 * * [simplify]: Extracting #4: cost 19 inf + 325 2.301 * * [simplify]: Extracting #5: cost 14 inf + 1612 2.301 * * [simplify]: Extracting #6: cost 12 inf + 3016 2.302 * * [simplify]: Extracting #7: cost 3 inf + 14359 2.303 * * [simplify]: Extracting #8: cost 0 inf + 19008 2.304 * [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.304 * * * [progress]: adding candidates to table 2.518 * * [progress]: iteration 3 / 4 2.518 * * * [progress]: picking best candidate 2.546 * * * * [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.546 * * * [progress]: localizing error 2.962 * * * [progress]: generating rewritten candidates 2.962 * * * * [progress]: [ 1 / 4 ] rewriting at (2 2 1 2 2 1) 2.969 * * * * [progress]: [ 2 / 4 ] rewriting at (2 2 1 2 2 1 1) 2.969 * * * * [progress]: [ 3 / 4 ] rewriting at (2 2 1 2 2) 2.986 * * * * [progress]: [ 4 / 4 ] rewriting at (2 2) 2.987 * * * [progress]: generating series expansions 2.987 * * * * [progress]: [ 1 / 4 ] generating series at (2 2 1 2 2 1) 2.987 * * * * [progress]: [ 2 / 4 ] generating series at (2 2 1 2 2 1 1) 2.987 * * * * [progress]: [ 3 / 4 ] generating series at (2 2 1 2 2) 2.987 * * * * [progress]: [ 4 / 4 ] generating series at (2 2) 2.988 * * * [progress]: simplifying candidates 2.988 * * * * [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.988 * * * * [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.988 * * * * [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.988 * * * * [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.988 * * * * [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.988 * [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.988 * * [simplify]: iteration 1: (8 enodes) 2.991 * * [simplify]: iteration 2: (14 enodes) 2.996 * * [simplify]: Extracting #0: cost 1 inf + 0 2.996 * * [simplify]: Extracting #1: cost 3 inf + 0 2.996 * * [simplify]: Extracting #2: cost 4 inf + 1 2.996 * * [simplify]: Extracting #3: cost 4 inf + 2 2.996 * * [simplify]: Extracting #4: cost 6 inf + 2 2.996 * * [simplify]: Extracting #5: cost 0 inf + 4012 2.996 * [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.998 * * * [progress]: adding candidates to table 3.167 * [progress]: [Phase 3 of 3] Extracting. 3.167 * * [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 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)))))>) 3.168 * * * [regime-changes]: Trying 2 branch expressions: (im re) 3.168 * * * * [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 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)))))>) 3.332 * * * * [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 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)))))>) 3.467 * * * [regime]: Found split indices: #