x
Yes
No
Do you want to visit DriveHQ English website?
首页
产品服务
价格
免费试用
下载客户端
关于我们
云文件服务
|
云备份服务
|
FTP服务
|
企业邮箱服务
|
网站托管
|
客户端软件
云文件服务
云备份服务
FTP服务
企业级邮箱服务
网站托管
客户端软件
laguerre.hpp - Hosted on DriveHQ Cloud IT Platform
返回上层目录
上传
下载
共享
发布
新建文件夹
新建文件
复制
剪切
删除
粘贴
评论
升级服务
路径: \\game3dprogramming\materials\GameFactory\GameFactoryDemo\references\boost_1_35_0\boost\math\special_functions\laguerre.hpp
旋转
特效
属性
历史版本
// (C) Copyright John Maddock 2006. // Use, modification and distribution are subject to the // Boost Software License, Version 1.0. (See accompanying file // LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt) #ifndef BOOST_MATH_SPECIAL_LAGUERRE_HPP #define BOOST_MATH_SPECIAL_LAGUERRE_HPP #include
#include
#include
namespace boost{ namespace math{ // Recurrance relation for Laguerre polynomials: template
inline typename tools::promote_args
::type laguerre_next(unsigned n, T1 x, T2 Ln, T3 Lnm1) { typedef typename tools::promote_args
::type result_type; return ((2 * n + 1 - result_type(x)) * result_type(Ln) - n * result_type(Lnm1)) / (n + 1); } namespace detail{ // Implement Laguerre polynomials via recurrance: template
T laguerre_imp(unsigned n, T x) { T p0 = 1; T p1 = 1 - x; if(n == 0) return p0; unsigned c = 1; while(c < n) { std::swap(p0, p1); p1 = laguerre_next(c, x, p0, p1); ++c; } return p1; } template
inline typename tools::promote_args
::type laguerre(unsigned n, T x, const Policy&, const mpl::true_&) { typedef typename tools::promote_args
::type result_type; typedef typename policies::evaluation
::type value_type; return policies::checked_narrowing_cast
(detail::laguerre_imp(n, static_cast
(x)), "boost::math::laguerre<%1%>(unsigned, %1%)"); } template
inline typename tools::promote_args
::type laguerre(unsigned n, unsigned m, T x, const mpl::false_&) { return boost::math::laguerre(n, m, x, policies::policy<>()); } } // namespace detail template
inline typename tools::promote_args
::type laguerre(unsigned n, T x) { return laguerre(n, x, policies::policy<>()); } // Recurrence for associated polynomials: template
inline typename tools::promote_args
::type laguerre_next(unsigned n, unsigned l, T1 x, T2 Pl, T3 Plm1) { typedef typename tools::promote_args
::type result_type; return ((2 * n + l + 1 - result_type(x)) * result_type(Pl) - (n + l) * result_type(Plm1)) / (n+1); } namespace detail{ // Laguerre Associated Polynomial: template
T laguerre_imp(unsigned n, unsigned m, T x, const Policy& pol) { // Special cases: if(m == 0) return boost::math::laguerre(n, x, pol); T p0 = 1; if(n == 0) return p0; T p1 = m + 1 - x; unsigned c = 1; while(c < n) { std::swap(p0, p1); p1 = laguerre_next(c, m, x, p0, p1); ++c; } return p1; } } template
inline typename tools::promote_args
::type laguerre(unsigned n, unsigned m, T x, const Policy& pol) { typedef typename tools::promote_args
::type result_type; typedef typename policies::evaluation
::type value_type; return policies::checked_narrowing_cast
(detail::laguerre_imp(n, m, static_cast
(x), pol), "boost::math::laguerre<%1%>(unsigned, unsigned, %1%)"); } template
inline typename laguerre_result
::type laguerre(unsigned n, T1 m, T2 x) { typedef typename policies::is_policy
::type tag_type; return detail::laguerre(n, m, x, tag_type()); } } // namespace math } // namespace boost #endif // BOOST_MATH_SPECIAL_LAGUERRE_HPP
laguerre.hpp
网页地址
文件地址
上一页
24/35
下一页
下载
( 3 KB )
Comments
Total ratings:
0
Average rating:
无评论
of 10
Would you like to comment?
Join now
, or
Logon
if you are already a member.