Fix for expression template use in chebyshev.hpp.

jzmaddock · jzmaddock · commit 0a014fd5ca01 · 2023-06-21T08:56:02.000+01:00
diff --git a/include/boost/math/special_functions/chebyshev.hpp b/include/boost/math/special_functions/chebyshev.hpp
@@ -77,7 +77,7 @@ inline Real chebyshev_imp(unsigned n, Real const & x, const Policy&)
         if (x > 1 || x < -1)
         {
             Real t = sqrt(x*x -1);
-            return static_cast<Real>((expt(x+t, n+1) - expt(x-t, n+1))/(2*t));
+            return static_cast<Real>((expt(static_cast<Real>(x+t), n+1) - expt(static_cast<Real>(x-t), n+1))/(2*t));
         }
         T1 = 2*x;
     }

Original file line number	Diff line number	Diff line change
`@@ -77,7 +77,7 @@ inline Real chebyshev_imp(unsigned n, Real const & x, const Policy&)`
`77`	`77`	`if (x > 1 \|\| x < -1)`
`78`	`78`	`{`
`79`	`79`	`Real t = sqrt(x*x -1);`
`80`		`- return static_cast<Real>((expt(x+t, n+1) - expt(x-t, n+1))/(2*t));`
	`80`	`+ return static_cast<Real>((expt(static_cast<Real>(x+t), n+1) - expt(static_cast<Real>(x-t), n+1))/(2*t));`
`81`	`81`	`}`
`82`	`82`	`T1 = 2*x;`
`83`	`83`	`}`