Toniolo and Linder, Equation (13)

Average Error: 34.7 → 29.7

Time: 46.1s

Precision: 64

Math

\[\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\]

↓

\[\begin{array}{l} \mathbf{if}\;t \le -3.75422355315464714 \cdot 10^{71}:\\ \;\;\;\;\sqrt{\left(\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)\right)\right) \cdot \left(2 \cdot n\right)\right) \cdot U}\\ \mathbf{elif}\;t \le 9.3445529337517344 \cdot 10^{-110}:\\ \;\;\;\;\sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\\ \mathbf{else}:\\ \;\;\;\;\sqrt{t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)\right)} \cdot \sqrt{\left(2 \cdot n\right) \cdot U}\\ \end{array}\]

C

double f(double n, double U, double t, double l, double Om, double U_) {
        double r83931 = 2.0;
        double r83932 = n;
        double r83933 = r83931 * r83932;
        double r83934 = U;
        double r83935 = r83933 * r83934;
        double r83936 = t;
        double r83937 = l;
        double r83938 = r83937 * r83937;
        double r83939 = Om;
        double r83940 = r83938 / r83939;
        double r83941 = r83931 * r83940;
        double r83942 = r83936 - r83941;
        double r83943 = r83937 / r83939;
        double r83944 = pow(r83943, r83931);
        double r83945 = r83932 * r83944;
        double r83946 = U_;
        double r83947 = r83934 - r83946;
        double r83948 = r83945 * r83947;
        double r83949 = r83942 - r83948;
        double r83950 = r83935 * r83949;
        double r83951 = sqrt(r83950);
        return r83951;
}

↓

double f(double n, double U, double t, double l, double Om, double U_) {
        double r83952 = t;
        double r83953 = -3.754223553154647e+71;
        bool r83954 = r83952 <= r83953;
        double r83955 = 2.0;
        double r83956 = l;
        double r83957 = Om;
        double r83958 = r83956 / r83957;
        double r83959 = r83956 * r83958;
        double r83960 = n;
        double r83961 = 2.0;
        double r83962 = r83955 / r83961;
        double r83963 = pow(r83958, r83962);
        double r83964 = r83960 * r83963;
        double r83965 = U;
        double r83966 = U_;
        double r83967 = r83965 - r83966;
        double r83968 = r83963 * r83967;
        double r83969 = r83964 * r83968;
        double r83970 = fma(r83955, r83959, r83969);
        double r83971 = r83952 - r83970;
        double r83972 = r83955 * r83960;
        double r83973 = r83971 * r83972;
        double r83974 = r83973 * r83965;
        double r83975 = sqrt(r83974);
        double r83976 = 9.344552933751734e-110;
        bool r83977 = r83952 <= r83976;
        double r83978 = r83972 * r83965;
        double r83979 = r83971 * r83978;
        double r83980 = sqrt(r83979);
        double r83981 = sqrt(r83971);
        double r83982 = sqrt(r83978);
        double r83983 = r83981 * r83982;
        double r83984 = r83977 ? r83980 : r83983;
        double r83985 = r83954 ? r83975 : r83984;
        return r83985;
}

Error

Bits error versus n

Bits error versus U

Bits error versus t

Bits error versus l

Bits error versus Om

Bits error versus U*

Derivation

1. Split input into 3 regimes

2. if t < -3.754223553154647e+71

   1. Initial program 35.8

      \[\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\]

   2. Simplified 35.8

      \[\leadsto \color{blue}{\sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}}\]
   3. Using strategy rm

   4. Applied *-un-lft-identity 35.8

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{\color{blue}{1 \cdot Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   5. Applied times-frac 33.1

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\frac{\ell}{1} \cdot \frac{\ell}{Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   6. Simplified 33.1

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\ell} \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   7. Using strategy rm

   8. Applied sqr-pow 33.1

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot \color{blue}{\left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right)}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   9. Applied associate-*r* 32.8

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \color{blue}{\left(\left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right)} \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   10. Using strategy rm

   11. Applied associate-*l* 33.1

       \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \color{blue}{\left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)}\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   12. Using strategy rm

   13. Applied associate-*r* 32.0

       \[\leadsto \sqrt{\color{blue}{\left(\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)\right)\right) \cdot \left(2 \cdot n\right)\right) \cdot U}}\]

   if -3.754223553154647e+71 < t < 9.344552933751734e-110

   1. Initial program 34.4

      \[\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\]

   2. Simplified 34.4

      \[\leadsto \color{blue}{\sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}}\]
   3. Using strategy rm

   4. Applied *-un-lft-identity 34.4

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{\color{blue}{1 \cdot Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   5. Applied times-frac 31.6

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\frac{\ell}{1} \cdot \frac{\ell}{Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   6. Simplified 31.6

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\ell} \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   7. Using strategy rm

   8. Applied sqr-pow 31.6

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot \color{blue}{\left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right)}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   9. Applied associate-*r* 30.4

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \color{blue}{\left(\left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right)} \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   10. Using strategy rm

   11. Applied associate-*l* 30.1

       \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \color{blue}{\left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)}\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   if 9.344552933751734e-110 < t

   1. Initial program 34.5

      \[\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\]

   2. Simplified 34.5

      \[\leadsto \color{blue}{\sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}}\]
   3. Using strategy rm

   4. Applied *-un-lft-identity 34.5

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{\color{blue}{1 \cdot Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   5. Applied times-frac 31.2

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\frac{\ell}{1} \cdot \frac{\ell}{Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   6. Simplified 31.2

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\ell} \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   7. Using strategy rm

   8. Applied sqr-pow 31.2

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot \color{blue}{\left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right)}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   9. Applied associate-*r* 30.7

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \color{blue}{\left(\left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right)} \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   10. Using strategy rm

   11. Applied associate-*l* 31.0

       \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \color{blue}{\left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)}\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   12. Using strategy rm

   13. Applied sqrt-prod 27.8

       \[\leadsto \color{blue}{\sqrt{t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)\right)} \cdot \sqrt{\left(2 \cdot n\right) \cdot U}}\]
3. Recombined 3 regimes into one program.

4. Final simplification 29.7

   \[\leadsto \begin{array}{l} \mathbf{if}\;t \le -3.75422355315464714 \cdot 10^{71}:\\ \;\;\;\;\sqrt{\left(\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)\right)\right) \cdot \left(2 \cdot n\right)\right) \cdot U}\\ \mathbf{elif}\;t \le 9.3445529337517344 \cdot 10^{-110}:\\ \;\;\;\;\sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\\ \mathbf{else}:\\ \;\;\;\;\sqrt{t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)\right)} \cdot \sqrt{\left(2 \cdot n\right) \cdot U}\\ \end{array}\]

Reproduce

herbie shell --seed 2019198 +o rules:numerics
(FPCore (n U t l Om U*)
  :name "Toniolo and Linder, Equation (13)"
  (sqrt (* (* (* 2.0 n) U) (- (- t (* 2.0 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2.0)) (- U U*))))))