using boost in OpenFOAM

服务器上已经安装了 boost，在 /usr/include/boost

在 OpenFOAM 中，不需要修改 options ，直接

#include <boost/lambda/lambda.hpp>

就可以用了：

Info<<"beta in boost ="<<boost::math::beta(1,1,1)<<endl;

boost 中的函数

beta(a,b)	$\int_0^1t^{a-1}(1-t)^{b-1}dt=\dfrac{\Gamma(a)\Gamma(b)}{\Gamma(a+b)}$
beta(a,b,x)	$\int_0^xt^{a-1}(1-t)^{b-1}dt$
ibeta(a,b,x)	$\dfrac{\int_0^xt^{a-1}(1-t)^{b-1}dt}{\int_0^1t^{a-1}(1-t)^{b-1}dt}=\dfrac{beta(a,b,x)}{beta(a,b)}$

GSL 中的函数

gsl_sf_beta(a,b)	$\int_0^1t^{a-1}(1-t)^{b-1}dt=\dfrac{\Gamma(a)\Gamma(b)}{\Gamma(a+b)}$
gsl_sf_beta_inc(a,b,x)	$\dfrac{\int_0^xt^{a-1}(1-t)^{b-1}dt}{\int_0^1t^{a-1}(1-t)^{b-1}dt}$