\[1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|}\]
Test:
Jmat.Real.erf
Bits:
128 bits
Bits error versus x
Time: 24.3 s
Input Error: 13.6
Output Error: 10.2
Log:
Profile: 🕒
\(\begin{cases} \left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) - \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} + \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(\sqrt[3]{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}}\right)}^3} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right) & \text{when } x \le -3.878602931916025 \cdot 10^{-24} \\ \frac{\left({1}^2 - {\left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)}^2\right) \cdot \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} - \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right) - \left(1 - \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) \cdot \left({\left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|}\right)}^2 - {\left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)}^2\right)}{\left(1 - \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) \cdot \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} - \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right)} & \text{when } x \le 1.1454841341683684 \cdot 10^{-10} \\ 1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + {\left(\sqrt{\frac{1}{1 + 0.3275911 \cdot \left|x\right|}}\right)}^2 \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|} & \text{otherwise} \end{cases}\)

    if x < -3.878602931916025e-24

    1. Started with
      \[1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|}\]
      1.6
    2. Applied taylor to get
      \[1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|} \leadsto \left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) - \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} + \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right)\]
      1.6
    3. Taylor expanded around 0 to get
      \[\color{red}{\left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) - \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} + \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right)} \leadsto \color{blue}{\left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) - \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} + \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right)}\]
      1.6
    4. Using strategy rm
      1.6
    5. Applied add-cube-cbrt to get
      \[\left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) - \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} + \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{\color{red}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right) \leadsto \left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) - \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} + \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{\color{blue}{{\left(\sqrt[3]{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}}\right)}^3}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right)\]
      1.6

    if -3.878602931916025e-24 < x < 1.1454841341683684e-10

    1. Started with
      \[1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|}\]
      28.0
    2. Applied taylor to get
      \[1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|} \leadsto \left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) - \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} + \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right)\]
      28.0
    3. Taylor expanded around 0 to get
      \[\color{red}{\left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) - \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} + \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right)} \leadsto \color{blue}{\left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) - \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} + \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right)}\]
      28.0
    4. Using strategy rm
      28.0
    5. Applied flip-+ to get
      \[\left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) - \color{red}{\left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} + \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right)} \leadsto \left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) - \color{blue}{\frac{{\left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|}\right)}^2 - {\left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)}^2}{0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} - \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)}}\]
      28.0
    6. Applied flip-+ to get
      \[\color{red}{\left(1 + \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right)} - \frac{{\left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|}\right)}^2 - {\left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)}^2}{0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} - \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)} \leadsto \color{blue}{\frac{{1}^2 - {\left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)}^2}{1 - \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)}} - \frac{{\left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|}\right)}^2 - {\left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)}^2}{0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} - \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)}\]
      30.8
    7. Applied frac-sub to get
      \[\color{red}{\frac{{1}^2 - {\left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)}^2}{1 - \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)} - \frac{{\left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|}\right)}^2 - {\left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)}^2}{0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} - \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)}} \leadsto \color{blue}{\frac{\left({1}^2 - {\left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)}^2\right) \cdot \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} - \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right) - \left(1 - \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) \cdot \left({\left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|}\right)}^2 - {\left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)}^2\right)}{\left(1 - \left(1.453152027 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{4}} + 0.284496736 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^2}\right)\right) \cdot \left(0.254829592 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{1 + 0.3275911 \cdot \left|x\right|} - \left(1.061405429 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{5}} + 1.421413741 \cdot \frac{e^{-{\left(\left|x\right|\right)}^2}}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{3}}\right)\right)}}\]
      20.7

    if 1.1454841341683684e-10 < x

    1. Started with
      \[1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|}\]
      0.5
    2. Using strategy rm
      0.5
    3. Applied add-sqr-sqrt to get
      \[1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \color{red}{\frac{1}{1 + 0.3275911 \cdot \left|x\right|}} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|} \leadsto 1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \color{blue}{{\left(\sqrt{\frac{1}{1 + 0.3275911 \cdot \left|x\right|}}\right)}^2} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|}\]
      0.5

  1. Removed slow pow expressions

Original test:


(lambda ((x default))
  #:name "Jmat.Real.erf"
  (- 1 (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))