-\frac{1}{\frac{\pi}{4}} \cdot \log \left(\frac{e^{\frac{\pi}{4} \cdot f} + e^{-\frac{\pi}{4} \cdot f}}{e^{\frac{\pi}{4} \cdot f} - e^{-\frac{\pi}{4} \cdot f}}\right)\left(\left(8 \cdot \frac{\left(f \cdot f\right) \cdot {\left(\pi \cdot 0.25\right)}^{5}}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}} + \left(16 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left(\left(f \cdot f\right) \cdot \left(0.00390625 \cdot {\pi}^{4}\right)\right)}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}} + \left(16 \cdot \frac{{\left(\pi \cdot 0.25\right)}^{5} \cdot \left(\left(f \cdot f\right) \cdot \left(\pi \cdot -0.25\right)\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{4}} + \left(0.046875 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left(\left(f \cdot f\right) \cdot {\pi}^{2}\right)}{\pi \cdot 0.25 - \pi \cdot -0.25} + \left(12 \cdot \frac{{\left(\pi \cdot 0.25\right)}^{2} \cdot \left(\left(f \cdot f\right) \cdot \left(-0.015625 \cdot {\pi}^{3}\right)\right)}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}} + \left(2 \cdot \frac{{\left(\pi \cdot 0.25\right)}^{4} \cdot \left(\left(f \cdot f\right) \cdot \pi\right)}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}} + \left(16 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left(\left(f \cdot f\right) \cdot \left(-0.0009765625 \cdot {\pi}^{5}\right)\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{4}} + \left(4 \cdot \frac{{\left(\pi \cdot 0.25\right)}^{2} \cdot \left(f \cdot \left(\pi \cdot -0.25\right)\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(0.25 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left(\left(f \cdot f\right) \cdot \left(-0.25 \cdot {\pi}^{3}\right)\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(0.75 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left(\left(f \cdot f\right) \cdot \left({\pi}^{2} \cdot 0.0625\right)\right)}{\pi \cdot 0.25 - \pi \cdot -0.25} + \left(0.125 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left({\pi}^{2} \cdot \left(\left(f \cdot f\right) \cdot -0.25\right)\right)}{\pi \cdot 0.25 - \pi \cdot -0.25} + \left(8 \cdot \frac{{\left(\pi \cdot 0.25\right)}^{4} \cdot \left(\left(f \cdot f\right) \cdot \left({\pi}^{2} \cdot 0.0625\right)\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{4}} + \left(8 \cdot \frac{\left(\left(f \cdot f\right) \cdot \left(0.00390625 \cdot {\pi}^{4}\right)\right) \cdot {\left(\pi \cdot 0.25\right)}^{2}}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{4}} + \left(4 \cdot \frac{{\left(\pi \cdot 0.25\right)}^{2} \cdot \left(\left(f \cdot f\right) \cdot \left({\pi}^{2} \cdot 0.0625\right)\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(4 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left(\left(f \cdot f\right) \cdot \left({\pi}^{4} \cdot -0.015625\right)\right)}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}} + \left(\frac{\left(\pi \cdot 0.25\right) \cdot \left(f \cdot \left(\pi \cdot -0.25\right)\right)}{\pi \cdot 0.25 - \pi \cdot -0.25} + \left(4 \cdot \frac{{\left(\pi \cdot 0.25\right)}^{2} \cdot \left(\left(f \cdot f\right) \cdot \left(-0.015625 \cdot {\pi}^{3}\right)\right)}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}} + \left(10.666666666666666 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left(\left(f \cdot f\right) \cdot \left(-0.015625 \cdot {\pi}^{3}\right)\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(8 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left(\left(f \cdot f\right) \cdot \left(0.00390625 \cdot {\pi}^{4}\right)\right)}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}} + \left(0.25 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left(f \cdot \pi\right)}{\pi \cdot 0.25 - \pi \cdot -0.25} + \left(\frac{\left(\pi \cdot 0.25\right) \cdot \left({\pi}^{3} \cdot \left(\left(f \cdot f\right) \cdot 0.0625\right)\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(2.6666666666666665 \cdot \frac{\left(\left(f \cdot f\right) \cdot \left(\pi \cdot -0.25\right)\right) \cdot {\left(\pi \cdot 0.25\right)}^{3}}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(2 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left({\pi}^{3} \cdot \left(\left(f \cdot f\right) \cdot 0.0625\right)\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(\log \left(\frac{4}{\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2}\right) + 4 \cdot \frac{\left(\pi \cdot 0.25\right) \cdot \left(f \cdot \left({\pi}^{2} \cdot 0.0625\right)\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}}\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right) - \left(32 \cdot \frac{\left(\left(f \cdot f\right) \cdot \left(-0.015625 \cdot {\pi}^{3}\right)\right) \cdot {\left(\pi \cdot 0.25\right)}^{3}}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{4}} + \left(16 \cdot \frac{\left(\left(f \cdot f\right) \cdot \left({\pi}^{2} \cdot 0.0625\right)\right) \cdot {\left(\pi \cdot 0.25\right)}^{3}}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}} + \left(4 \cdot \frac{f \cdot {\left(\pi \cdot 0.25\right)}^{3}}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(\frac{f \cdot \left({\pi}^{2} \cdot 0.0625\right)}{\pi \cdot 0.25 - \pi \cdot -0.25} + \left(8 \cdot \frac{\left(f \cdot f\right) \cdot {\left(\pi \cdot 0.25\right)}^{6}}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{4}} + \left(4 \cdot \frac{\left(\left(f \cdot f\right) \cdot \left(\pi \cdot -0.25\right)\right) \cdot {\left(\pi \cdot 0.25\right)}^{3}}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(0.046875 \cdot \frac{\left(f \cdot f\right) \cdot \left(-0.25 \cdot {\pi}^{3}\right)}{\pi \cdot 0.25 - \pi \cdot -0.25} + \left(0.125 \cdot \frac{{\pi}^{3} \cdot \left(\left(f \cdot f\right) \cdot 0.0625\right)}{\pi \cdot 0.25 - \pi \cdot -0.25} + \left(0.75 \cdot \frac{\left(f \cdot f\right) \cdot \left(-0.015625 \cdot {\pi}^{3}\right)}{\pi \cdot 0.25 - \pi \cdot -0.25} + \left(8 \cdot \frac{{\pi}^{6} \cdot \left(\left(f \cdot f\right) \cdot 0.000244140625\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{4}} + \left(2.6666666666666665 \cdot \frac{\left(f \cdot f\right) \cdot {\left(\pi \cdot 0.25\right)}^{4}}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(\log f + \left(0.25 \cdot \frac{f \cdot \left(-0.25 \cdot {\pi}^{2}\right)}{\pi \cdot 0.25 - \pi \cdot -0.25} + \left(8.666666666666666 \cdot \frac{\left(f \cdot f\right) \cdot \left(0.00390625 \cdot {\pi}^{4}\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(2 \cdot \frac{\left(f \cdot f\right) \cdot \left({\pi}^{4} \cdot -0.015625\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(4 \cdot \frac{\left({\pi}^{2} \cdot \left(\left(f \cdot f\right) \cdot -0.25\right)\right) \cdot {\left(\pi \cdot 0.25\right)}^{3}}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}} + \left(16 \cdot \frac{\left(\left(f \cdot f\right) \cdot \left({\pi}^{2} \cdot 0.0625\right)\right) \cdot {\left(\pi \cdot 0.25\right)}^{3}}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}} + \left(\frac{\left(\left(f \cdot f\right) \cdot \pi\right) \cdot {\left(\pi \cdot 0.25\right)}^{3}}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(2 \cdot \frac{{\pi}^{5} \cdot \left(\left(f \cdot f\right) \cdot 0.00390625\right)}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}} + \left(0.125 \cdot \frac{{\pi}^{4} \cdot \left(\left(f \cdot f\right) \cdot 0.0625\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(0.125 \cdot \frac{\left(\left(f \cdot f\right) \cdot {\pi}^{2}\right) \cdot {\left(\pi \cdot 0.25\right)}^{2}}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(4 \cdot \frac{f \cdot \left(-0.015625 \cdot {\pi}^{3}\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + \left(2 \cdot \frac{{\left(\pi \cdot 0.25\right)}^{2} \cdot \left(\left(f \cdot f\right) \cdot \left({\pi}^{2} \cdot 0.0625\right)\right)}{{\left(\left(\pi \cdot 0.25 - \pi \cdot -0.25\right) \cdot 2\right)}^{2}} + 16 \cdot \frac{\left(f \cdot f\right) \cdot \left(-0.0009765625 \cdot {\pi}^{5}\right)}{8 \cdot {\left(\pi \cdot 0.25 - \pi \cdot -0.25\right)}^{3}}\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right) \cdot \frac{-4}{\pi}(FPCore (f)
:precision binary64
(-
(*
(/ 1.0 (/ PI 4.0))
(log
(/
(+ (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f))))
(- (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f)))))))))(FPCore (f)
:precision binary64
(*
(-
(+
(*
8.0
(/
(* (* f f) (pow (* PI 0.25) 5.0))
(* 8.0 (pow (- (* PI 0.25) (* PI -0.25)) 3.0))))
(+
(*
16.0
(/
(* (* PI 0.25) (* (* f f) (* 0.00390625 (pow PI 4.0))))
(* 8.0 (pow (- (* PI 0.25) (* PI -0.25)) 3.0))))
(+
(*
16.0
(/
(* (pow (* PI 0.25) 5.0) (* (* f f) (* PI -0.25)))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 4.0)))
(+
(*
0.046875
(/
(* (* PI 0.25) (* (* f f) (pow PI 2.0)))
(- (* PI 0.25) (* PI -0.25))))
(+
(*
12.0
(/
(* (pow (* PI 0.25) 2.0) (* (* f f) (* -0.015625 (pow PI 3.0))))
(* 8.0 (pow (- (* PI 0.25) (* PI -0.25)) 3.0))))
(+
(*
2.0
(/
(* (pow (* PI 0.25) 4.0) (* (* f f) PI))
(* 8.0 (pow (- (* PI 0.25) (* PI -0.25)) 3.0))))
(+
(*
16.0
(/
(* (* PI 0.25) (* (* f f) (* -0.0009765625 (pow PI 5.0))))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 4.0)))
(+
(*
4.0
(/
(* (pow (* PI 0.25) 2.0) (* f (* PI -0.25)))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(*
0.25
(/
(* (* PI 0.25) (* (* f f) (* -0.25 (pow PI 3.0))))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(*
0.75
(/
(* (* PI 0.25) (* (* f f) (* (pow PI 2.0) 0.0625)))
(- (* PI 0.25) (* PI -0.25))))
(+
(*
0.125
(/
(* (* PI 0.25) (* (pow PI 2.0) (* (* f f) -0.25)))
(- (* PI 0.25) (* PI -0.25))))
(+
(*
8.0
(/
(* (pow (* PI 0.25) 4.0) (* (* f f) (* (pow PI 2.0) 0.0625)))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 4.0)))
(+
(*
8.0
(/
(*
(* (* f f) (* 0.00390625 (pow PI 4.0)))
(pow (* PI 0.25) 2.0))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 4.0)))
(+
(*
4.0
(/
(*
(pow (* PI 0.25) 2.0)
(* (* f f) (* (pow PI 2.0) 0.0625)))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(*
4.0
(/
(* (* PI 0.25) (* (* f f) (* (pow PI 4.0) -0.015625)))
(* 8.0 (pow (- (* PI 0.25) (* PI -0.25)) 3.0))))
(+
(/
(* (* PI 0.25) (* f (* PI -0.25)))
(- (* PI 0.25) (* PI -0.25)))
(+
(*
4.0
(/
(*
(pow (* PI 0.25) 2.0)
(* (* f f) (* -0.015625 (pow PI 3.0))))
(* 8.0 (pow (- (* PI 0.25) (* PI -0.25)) 3.0))))
(+
(*
10.666666666666666
(/
(* (* PI 0.25) (* (* f f) (* -0.015625 (pow PI 3.0))))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(*
8.0
(/
(* (* PI 0.25) (* (* f f) (* 0.00390625 (pow PI 4.0))))
(* 8.0 (pow (- (* PI 0.25) (* PI -0.25)) 3.0))))
(+
(*
0.25
(/
(* (* PI 0.25) (* f PI))
(- (* PI 0.25) (* PI -0.25))))
(+
(/
(* (* PI 0.25) (* (pow PI 3.0) (* (* f f) 0.0625)))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0))
(+
(*
2.6666666666666665
(/
(* (* (* f f) (* PI -0.25)) (pow (* PI 0.25) 3.0))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(*
2.0
(/
(* (* PI 0.25) (* (pow PI 3.0) (* (* f f) 0.0625)))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(log (/ 4.0 (* (- (* PI 0.25) (* PI -0.25)) 2.0)))
(*
4.0
(/
(* (* PI 0.25) (* f (* (pow PI 2.0) 0.0625)))
(pow
(* (- (* PI 0.25) (* PI -0.25)) 2.0)
2.0)))))))))))))))))))))))))))
(+
(*
32.0
(/
(* (* (* f f) (* -0.015625 (pow PI 3.0))) (pow (* PI 0.25) 3.0))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 4.0)))
(+
(*
16.0
(/
(* (* (* f f) (* (pow PI 2.0) 0.0625)) (pow (* PI 0.25) 3.0))
(* 8.0 (pow (- (* PI 0.25) (* PI -0.25)) 3.0))))
(+
(*
4.0
(/
(* f (pow (* PI 0.25) 3.0))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(/ (* f (* (pow PI 2.0) 0.0625)) (- (* PI 0.25) (* PI -0.25)))
(+
(*
8.0
(/
(* (* f f) (pow (* PI 0.25) 6.0))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 4.0)))
(+
(*
4.0
(/
(* (* (* f f) (* PI -0.25)) (pow (* PI 0.25) 3.0))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(*
0.046875
(/ (* (* f f) (* -0.25 (pow PI 3.0))) (- (* PI 0.25) (* PI -0.25))))
(+
(*
0.125
(/
(* (pow PI 3.0) (* (* f f) 0.0625))
(- (* PI 0.25) (* PI -0.25))))
(+
(*
0.75
(/
(* (* f f) (* -0.015625 (pow PI 3.0)))
(- (* PI 0.25) (* PI -0.25))))
(+
(*
8.0
(/
(* (pow PI 6.0) (* (* f f) 0.000244140625))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 4.0)))
(+
(*
2.6666666666666665
(/
(* (* f f) (pow (* PI 0.25) 4.0))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(log f)
(+
(*
0.25
(/ (* f (* -0.25 (pow PI 2.0))) (- (* PI 0.25) (* PI -0.25))))
(+
(*
8.666666666666666
(/
(* (* f f) (* 0.00390625 (pow PI 4.0)))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(*
2.0
(/
(* (* f f) (* (pow PI 4.0) -0.015625))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(*
4.0
(/
(*
(* (pow PI 2.0) (* (* f f) -0.25))
(pow (* PI 0.25) 3.0))
(* 8.0 (pow (- (* PI 0.25) (* PI -0.25)) 3.0))))
(+
(*
16.0
(/
(*
(* (* f f) (* (pow PI 2.0) 0.0625))
(pow (* PI 0.25) 3.0))
(* 8.0 (pow (- (* PI 0.25) (* PI -0.25)) 3.0))))
(+
(/
(* (* (* f f) PI) (pow (* PI 0.25) 3.0))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0))
(+
(*
2.0
(/
(* (pow PI 5.0) (* (* f f) 0.00390625))
(* 8.0 (pow (- (* PI 0.25) (* PI -0.25)) 3.0))))
(+
(*
0.125
(/
(* (pow PI 4.0) (* (* f f) 0.0625))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(*
0.125
(/
(* (* (* f f) (pow PI 2.0)) (pow (* PI 0.25) 2.0))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(*
4.0
(/
(* f (* -0.015625 (pow PI 3.0)))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(+
(*
2.0
(/
(*
(pow (* PI 0.25) 2.0)
(* (* f f) (* (pow PI 2.0) 0.0625)))
(pow (* (- (* PI 0.25) (* PI -0.25)) 2.0) 2.0)))
(*
16.0
(/
(* (* f f) (* -0.0009765625 (pow PI 5.0)))
(*
8.0
(pow
(- (* PI 0.25) (* PI -0.25))
3.0))))))))))))))))))))))))))))
(/ -4.0 PI)))double code(double f) {
return -((1.0 / (((double) M_PI) / 4.0)) * log((exp((((double) M_PI) / 4.0) * f) + exp(-((((double) M_PI) / 4.0) * f))) / (exp((((double) M_PI) / 4.0) * f) - exp(-((((double) M_PI) / 4.0) * f)))));
}
double code(double f) {
return (((8.0 * (((f * f) * pow((((double) M_PI) * 0.25), 5.0)) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))) + ((16.0 * (((((double) M_PI) * 0.25) * ((f * f) * (0.00390625 * pow(((double) M_PI), 4.0)))) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))) + ((16.0 * ((pow((((double) M_PI) * 0.25), 5.0) * ((f * f) * (((double) M_PI) * -0.25))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 4.0))) + ((0.046875 * (((((double) M_PI) * 0.25) * ((f * f) * pow(((double) M_PI), 2.0))) / ((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)))) + ((12.0 * ((pow((((double) M_PI) * 0.25), 2.0) * ((f * f) * (-0.015625 * pow(((double) M_PI), 3.0)))) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))) + ((2.0 * ((pow((((double) M_PI) * 0.25), 4.0) * ((f * f) * ((double) M_PI))) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))) + ((16.0 * (((((double) M_PI) * 0.25) * ((f * f) * (-0.0009765625 * pow(((double) M_PI), 5.0)))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 4.0))) + ((4.0 * ((pow((((double) M_PI) * 0.25), 2.0) * (f * (((double) M_PI) * -0.25))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + ((0.25 * (((((double) M_PI) * 0.25) * ((f * f) * (-0.25 * pow(((double) M_PI), 3.0)))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + ((0.75 * (((((double) M_PI) * 0.25) * ((f * f) * (pow(((double) M_PI), 2.0) * 0.0625))) / ((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)))) + ((0.125 * (((((double) M_PI) * 0.25) * (pow(((double) M_PI), 2.0) * ((f * f) * -0.25))) / ((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)))) + ((8.0 * ((pow((((double) M_PI) * 0.25), 4.0) * ((f * f) * (pow(((double) M_PI), 2.0) * 0.0625))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 4.0))) + ((8.0 * ((((f * f) * (0.00390625 * pow(((double) M_PI), 4.0))) * pow((((double) M_PI) * 0.25), 2.0)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 4.0))) + ((4.0 * ((pow((((double) M_PI) * 0.25), 2.0) * ((f * f) * (pow(((double) M_PI), 2.0) * 0.0625))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + ((4.0 * (((((double) M_PI) * 0.25) * ((f * f) * (pow(((double) M_PI), 4.0) * -0.015625))) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))) + ((((((double) M_PI) * 0.25) * (f * (((double) M_PI) * -0.25))) / ((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25))) + ((4.0 * ((pow((((double) M_PI) * 0.25), 2.0) * ((f * f) * (-0.015625 * pow(((double) M_PI), 3.0)))) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))) + ((10.666666666666666 * (((((double) M_PI) * 0.25) * ((f * f) * (-0.015625 * pow(((double) M_PI), 3.0)))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + ((8.0 * (((((double) M_PI) * 0.25) * ((f * f) * (0.00390625 * pow(((double) M_PI), 4.0)))) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))) + ((0.25 * (((((double) M_PI) * 0.25) * (f * ((double) M_PI))) / ((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)))) + ((((((double) M_PI) * 0.25) * (pow(((double) M_PI), 3.0) * ((f * f) * 0.0625))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0)) + ((2.6666666666666665 * ((((f * f) * (((double) M_PI) * -0.25)) * pow((((double) M_PI) * 0.25), 3.0)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + ((2.0 * (((((double) M_PI) * 0.25) * (pow(((double) M_PI), 3.0) * ((f * f) * 0.0625))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + (log(4.0 / (((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0)) + (4.0 * (((((double) M_PI) * 0.25) * (f * (pow(((double) M_PI), 2.0) * 0.0625))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))))))))))))))))))))))))))) - ((32.0 * ((((f * f) * (-0.015625 * pow(((double) M_PI), 3.0))) * pow((((double) M_PI) * 0.25), 3.0)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 4.0))) + ((16.0 * ((((f * f) * (pow(((double) M_PI), 2.0) * 0.0625)) * pow((((double) M_PI) * 0.25), 3.0)) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))) + ((4.0 * ((f * pow((((double) M_PI) * 0.25), 3.0)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + (((f * (pow(((double) M_PI), 2.0) * 0.0625)) / ((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25))) + ((8.0 * (((f * f) * pow((((double) M_PI) * 0.25), 6.0)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 4.0))) + ((4.0 * ((((f * f) * (((double) M_PI) * -0.25)) * pow((((double) M_PI) * 0.25), 3.0)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + ((0.046875 * (((f * f) * (-0.25 * pow(((double) M_PI), 3.0))) / ((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)))) + ((0.125 * ((pow(((double) M_PI), 3.0) * ((f * f) * 0.0625)) / ((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)))) + ((0.75 * (((f * f) * (-0.015625 * pow(((double) M_PI), 3.0))) / ((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)))) + ((8.0 * ((pow(((double) M_PI), 6.0) * ((f * f) * 0.000244140625)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 4.0))) + ((2.6666666666666665 * (((f * f) * pow((((double) M_PI) * 0.25), 4.0)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + (log(f) + ((0.25 * ((f * (-0.25 * pow(((double) M_PI), 2.0))) / ((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)))) + ((8.666666666666666 * (((f * f) * (0.00390625 * pow(((double) M_PI), 4.0))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + ((2.0 * (((f * f) * (pow(((double) M_PI), 4.0) * -0.015625)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + ((4.0 * (((pow(((double) M_PI), 2.0) * ((f * f) * -0.25)) * pow((((double) M_PI) * 0.25), 3.0)) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))) + ((16.0 * ((((f * f) * (pow(((double) M_PI), 2.0) * 0.0625)) * pow((((double) M_PI) * 0.25), 3.0)) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))) + (((((f * f) * ((double) M_PI)) * pow((((double) M_PI) * 0.25), 3.0)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0)) + ((2.0 * ((pow(((double) M_PI), 5.0) * ((f * f) * 0.00390625)) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))) + ((0.125 * ((pow(((double) M_PI), 4.0) * ((f * f) * 0.0625)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + ((0.125 * ((((f * f) * pow(((double) M_PI), 2.0)) * pow((((double) M_PI) * 0.25), 2.0)) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + ((4.0 * ((f * (-0.015625 * pow(((double) M_PI), 3.0))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + ((2.0 * ((pow((((double) M_PI) * 0.25), 2.0) * ((f * f) * (pow(((double) M_PI), 2.0) * 0.0625))) / pow((((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)) * 2.0), 2.0))) + (16.0 * (((f * f) * (-0.0009765625 * pow(((double) M_PI), 5.0))) / (8.0 * pow(((((double) M_PI) * 0.25) - (((double) M_PI) * -0.25)), 3.0)))))))))))))))))))))))))))) * (-4.0 / ((double) M_PI));
}



Bits error versus f
Results
Initial program 61.5
Simplified61.5
rmApplied flip--_binary6461.6
Simplified61.6
Taylor expanded around 0 2.3
Simplified2.3
Final simplification2.3
herbie shell --seed 2020342
(FPCore (f)
:name "VandenBroeck and Keller, Equation (20)"
:precision binary64
(- (* (/ 1.0 (/ PI 4.0)) (log (/ (+ (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f)))) (- (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f)))))))))