22#error "This file must not be included by compiled programs."
25#include <plugins/include/clad/Differentiator/Differentiator.h>
36namespace custom_derivatives {
104#ifdef R__HAS_MATHMORE
126template <
typename T,
typename U>
163template <
typename T,
typename U>
183template <
typename T,
typename U>
202template <
typename T,
typename U>
206 *
d_x += t.pushforward *
p;
208 *
d_y += t.pushforward *
p;
253 return {3.1415926535897931, 0.};
258 return {2.3025850929940459, 0.};
272 static double p1[5] = {0.4259894875,-0.1249762550, 0.03984243700, -0.006298287635, 0.001511162253};
273 static double q1[5] = {1.0 ,-0.3388260629, 0.09594393323, -0.01608042283, 0.003778942063};
275 static double p2[5] = {0.1788541609, 0.1173957403, 0.01488850518, -0.001394989411, 0.0001283617211};
276 static double q2[5] = {1.0 , 0.7428795082, 0.3153932961, 0.06694219548, 0.008790609714};
278 static double p3[5] = {0.1788544503, 0.09359161662,0.006325387654, 0.00006611667319,-0.000002031049101};
279 static double q3[5] = {1.0 , 0.6097809921, 0.2560616665, 0.04746722384, 0.006957301675};
281 static double p4[5] = {0.9874054407, 118.6723273, 849.2794360, -743.7792444, 427.0262186};
282 static double q4[5] = {1.0 , 106.8615961, 337.6496214, 2016.712389, 1597.063511};
284 static double p5[5] = {1.003675074, 167.5702434, 4789.711289, 21217.86767, -22324.94910};
285 static double q5[5] = {1.0 , 156.9424537, 3745.310488, 9834.698876, 66924.28357};
287 static double p6[5] = {1.000827619, 664.9143136, 62972.92665, 475554.6998, -5743609.109};
288 static double q6[5] = {1.0 , 651.4101098, 56974.73333, 165917.4725, -2815759.939};
290 static double a1[3] = {0.04166666667,-0.01996527778, 0.02709538966};
292 static double a2[2] = {-1.845568670,-4.284640743};
294 const double _const0 = 0.3989422803;
295 double v = (
x - x0) / xi;
299 double u = ::std::exp(
v + 1.);
302 const double ue = ::std::exp(-1 /
u);
303 const double us = ::std::sqrt(
u);
322 _r4 +=
_r_d2 * clad::custom_derivatives::sqrt_pushforward(
u, 1.).pushforward;
338 double u = ::std::exp(-
v - 1);
341 double _t8 = ::std::exp(-
u);
342 double _t7 = ::std::sqrt(
u);
355 clad::custom_derivatives::sqrt_pushforward(
u, 1.).pushforward;
460 }
else if (
v < 300) {
490 double _t25 = ::std::log(
v);
491 double _t24 = (
v + 1);
520 double _r0 =
_d_v * -((
x - x0) / (xi * xi));
527 static double p1[5] = {0.2514091491e+0,-0.6250580444e-1, 0.1458381230e-1,-0.2108817737e-2, 0.7411247290e-3};
528 static double q1[5] = {1.0 ,-0.5571175625e-2, 0.6225310236e-1,-0.3137378427e-2, 0.1931496439e-2};
530 static double p2[4] = {0.2868328584e+0, 0.3564363231e+0, 0.1523518695e+0, 0.2251304883e-1};
531 static double q2[4] = {1.0 , 0.6191136137e+0, 0.1720721448e+0, 0.2278594771e-1};
533 static double p3[4] = {0.2868329066e+0, 0.3003828436e+0, 0.9950951941e-1, 0.8733827185e-2};
534 static double q3[4] = {1.0 , 0.4237190502e+0, 0.1095631512e+0, 0.8693851567e-2};
536 static double p4[4] = {0.1000351630e+1, 0.4503592498e+1, 0.1085883880e+2, 0.7536052269e+1};
537 static double q4[4] = {1.0 , 0.5539969678e+1, 0.1933581111e+2, 0.2721321508e+2};
539 static double p5[4] = {0.1000006517e+1, 0.4909414111e+2, 0.8505544753e+2, 0.1532153455e+3};
540 static double q5[4] = {1.0 , 0.5009928881e+2, 0.1399819104e+3, 0.4200002909e+3};
542 static double p6[4] = {0.1000000983e+1, 0.1329868456e+3, 0.9162149244e+3,-0.9605054274e+3};
543 static double q6[4] = {1.0 , 0.1339887843e+3, 0.1055990413e+4, 0.5532224619e+3};
545 static double a1[4] = {0 ,-0.4583333333e+0, 0.6675347222e+0,-0.1641741416e+1};
546 static double a2[4] = {0 , 1.0 ,-0.4227843351e+0,-0.2043403138e+1};
549 const double v = (
x - x0) / xi;
553 const double _const0 = 0.3989422803;
554 double u = ::std::exp(
v + 1);
555 double _t3 = ::std::exp(-1. /
u);
556 double _t2 = ::std::sqrt(
u);
559 double _r3 =
_r2 * -(-1. / (
u *
u));
563 clad::custom_derivatives::sqrt_pushforward(
u, 1.).pushforward;
571 double u = ::std::exp(-
v - 1);
572 double _t8 = ::std::exp(-
u);
573 double _t7 = ::std::sqrt(
u);
580 _r8 +=
_r7 * clad::custom_derivatives::sqrt_pushforward(
u, 1.).pushforward;
636 }
else if (
v < 300) {
651 double _t25 = ::std::log(
v);
652 double _t24 = (
v + 1);
654 double u = 1. /
_t23;
674#ifdef R__HAS_MATHMORE
681 constexpr double kMACHEP = 1.11022302462515654042363166809e-16;
682 constexpr double kMAXLOG = 709.782712893383973096206318587;
683 constexpr double kMINLOG = -708.396418532264078748994506896;
684 constexpr double kMAXSTIR = 108.116855767857671821730036754;
685 constexpr double kMAXLGM = 2.556348e305;
686 constexpr double kBig = 4.503599627370496e15;
687 constexpr double kBiginv = 2.22044604925031308085e-16;
695 clad::tape<double>
_t7 = {};
696 clad::tape<double>
_t8 = {};
697 clad::tape<double>
_t9 = {};
703 if ((
x > 1.) && (
x >
a)) {
712 ax =
a *
_t1 -
x - ::std::lgamma(
a);
727 unsigned long _t6 = 0;
800 constexpr double kMACHEP = 1.11022302462515654042363166809e-16;
801 constexpr double kMAXLOG = 709.782712893383973096206318587;
802 constexpr double kMINLOG = -708.396418532264078748994506896;
803 constexpr double kMAXSTIR = 108.116855767857671821730036754;
804 constexpr double kMAXLGM = 2.556348e305;
805 constexpr double kBig = 4.503599627370496e15;
806 constexpr double kBiginv = 2.22044604925031308085e-16;
821 clad::tape<double>
_t12 = {};
822 clad::tape<double>
_t13 = {};
823 clad::tape<double>
_t14 = {};
824 clad::tape<double>
_t15 = {};
825 clad::tape<double>
_t16 = {};
826 clad::tape<double>
_t17 = {};
827 clad::tape<double>
_t19 = {};
828 clad::tape<double>
_t20 = {};
829 clad::tape<double>
_t21 = {};
830 clad::tape<double>
_t22 = {};
831 clad::tape<double>
_t23 = {};
832 clad::tape<double>
_t24 = {};
833 clad::tape<double>
_t25 = {};
834 clad::tape<double>
_t26 = {};
835 clad::tape<double>
_t27 = {};
836 clad::tape<bool>
_t29 = {};
837 clad::tape<double>
_t30 = {};
838 clad::tape<double>
_t31 = {};
839 clad::tape<double>
_t32 = {};
840 clad::tape<double>
_t33 = {};
847 if ((
x < 1.) || (
x <
a)) {
856 ax =
a *
_t1 -
x - ::std::lgamma(
a);
902 t = ::std::abs((
ans -
r) /
r);
919 bool _t28 = ::std::abs(
pk) > kBig;
940 if (clad::pop(
_t29)) {
990 if (clad::pop(
_t19)) {
1002 _r7 +=
_r_d17 * clad::custom_derivatives::std::abs_pushforward((
ans -
r) /
r, 1.).pushforward;
1009 r = clad::pop(
_t20);
1017 t = clad::pop(
_t23);
1047 z = clad::pop(
_t14);
1051 y = clad::pop(
_t13);
1055 c = clad::pop(
_t12);
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
winID h TVirtualViewer3D TVirtualGLPainter p
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t r
Namespace for new Math classes and functions.
tbb::task_arena is an alias of tbb::interface7::task_arena, which doesn't allow to forward declare tb...
Double_t CosH(Double_t)
Returns the hyperbolic cosine of x.
Double_t ACos(Double_t)
Returns the principal value of the arc cosine of x, expressed in radians.
Short_t Max(Short_t a, Short_t b)
Returns the largest of a and b.
Double_t ASin(Double_t)
Returns the principal value of the arc sine of x, expressed in radians.
Double_t Log2(Double_t x)
Returns the binary (base-2) logarithm of x.
Double_t Exp(Double_t x)
Returns the base-e exponential function of x, which is e raised to the power x.
Double_t Erf(Double_t x)
Computation of the error function erf(x).
Double_t ATan(Double_t)
Returns the principal value of the arc tangent of x, expressed in radians.
Double_t ASinH(Double_t)
Returns the area hyperbolic sine of x.
Double_t TanH(Double_t)
Returns the hyperbolic tangent of x.
Double_t ACosH(Double_t)
Returns the nonnegative area hyperbolic cosine of x.
Double_t Log(Double_t x)
Returns the natural logarithm of x.
Double_t Erfc(Double_t x)
Computes the complementary error function erfc(x).
Double_t Sq(Double_t x)
Returns x*x.
Double_t Sqrt(Double_t x)
Returns the square root of x.
LongDouble_t Power(LongDouble_t x, LongDouble_t y)
Returns x raised to the power y.
Short_t Min(Short_t a, Short_t b)
Returns the smallest of a and b.
constexpr Double_t Ln10()
Natural log of 10 (to convert log to ln)
Double_t Hypot(Double_t x, Double_t y)
Returns sqrt(x*x + y*y)
Double_t Cos(Double_t)
Returns the cosine of an angle of x radians.
Double_t LnGamma(Double_t z)
Computation of ln[gamma(z)] for all z.
Double_t Sin(Double_t)
Returns the sine of an angle of x radians.
Double_t Tan(Double_t)
Returns the tangent of an angle of x radians.
Double_t ATanH(Double_t)
Returns the area hyperbolic tangent of x.
Double_t Log10(Double_t x)
Returns the common (base-10) logarithm of x.
Short_t Abs(Short_t d)
Returns the absolute value of parameter Short_t d.
Double_t SinH(Double_t)
Returns the hyperbolic sine of `x.
void landau_pdf_pullback(double x, double xi, double x0, double d_out, double *d_x, double *d_xi, double *d_x0)
void landau_cdf_pullback(double x, double xi, double x0, double d_out, double *d_x, double *d_xi, double *d_x0)
ValueAndPushforward< T, T > CosH_pushforward(T x, T d_x)
ValueAndPushforward< T, T > Abs_pushforward(T x, T d_x)
void Min_pullback(T a, T b, U p, clad::array_ref< T > d_a, clad::array_ref< T > d_b)
ValueAndPushforward< T, T > Sq_pushforward(T x, T d_x)
void Max_pullback(T a, T b, U p, clad::array_ref< T > d_a, clad::array_ref< T > d_b)
ValueAndPushforward< T, T > Erf_pushforward(T x, T d_x)
ValueAndPushforward< T, T > Erfc_pushforward(T x, T d_x)
ValueAndPushforward< T, T > Sin_pushforward(T x, T d_x)
ValueAndPushforward< T, T > Max_pushforward(T x, T y, T d_x, T d_y)
ValueAndPushforward< T, T > Hypot_pushforward(T x, T y, T d_x, T d_y)
ValueAndPushforward< T, T > ASinH_pushforward(T x, T d_x)
ValueAndPushforward< T, T > ACosH_pushforward(T x, T d_x)
ValueAndPushforward< T, T > ASin_pushforward(T x, T d_x)
ValueAndPushforward< T, T > Cos_pushforward(T x, T d_x)
ValueAndPushforward< T, T > Sqrt_pushforward(T x, T d_x)
ValueAndPushforward< T, T > Tan_pushforward(T x, T d_x)
void Hypot_pullback(T x, T y, U p, clad::array_ref< T > d_x, clad::array_ref< T > d_y)
ValueAndPushforward< T, T > Power_pushforward(T x, T y, T d_x, T d_y)
void Power_pullback(T x, T y, U p, clad::array_ref< T > d_x, clad::array_ref< T > d_y)
ValueAndPushforward< T, T > Min_pushforward(T x, T y, T d_x, T d_y)
ValueAndPushforward< T, T > Log_pushforward(T x, T d_x)
ValueAndPushforward< T, T > Log10_pushforward(T x, T d_x)
ValueAndPushforward< T, T > TanH_pushforward(T x, T d_x)
ValueAndPushforward< T, T > ACos_pushforward(T x, T d_x)
ValueAndPushforward< T, T > SinH_pushforward(T x, T d_x)
ValueAndPushforward< T, T > Exp_pushforward(T x, T d_x)
ValueAndPushforward< T, T > Log2_pushforward(T x, T d_x)
ValueAndPushforward< T, T > ATanH_pushforward(T x, T d_x)
ValueAndPushforward< T, T > ATan_pushforward(T x, T d_x)