5 #ifndef SHINOZUKAPP_LAYEREDCOMP_HPP 6 #define SHINOZUKAPP_LAYEREDCOMP_HPP 9 #include <boost/random/mersenne_twister.hpp> 10 #include <boost/random/uniform_real_distribution.hpp> 11 #include <boost/math/special_functions/erf.hpp> 12 #include <boost/random/normal_distribution.hpp> 13 #include <boost/math/special_functions/gamma.hpp> 14 #include <boost/math/special_functions/beta.hpp> 24 template<
typename typearg>
26 double s_tau(
double & tau);
31 void icdf_gamma(Epetra_Vector & V, Epetra_Vector & G,
double & alpha,
double & beta);
32 void icdf_beta(Epetra_Vector & V, Epetra_Vector & B,
double & tau1,
double & tau2);
38 boost::random::mt19937
rng;
39 boost::random::uniform_real_distribution<>
phi_;
40 boost::random::uniform_real_distribution<>
psi_;
boost::random::mt19937 rng
double s_tau(double &tau)
double tau_beta(typearg &beta)
boost::random::uniform_real_distribution psi_
void icdf_beta(Epetra_Vector &V, Epetra_Vector &B, double &tau1, double &tau2)
boost::random::uniform_real_distribution phi_
void icdf_gamma(Epetra_Vector &V, Epetra_Vector &G, double &alpha, double &beta)
void generator_gauss_points(Epetra_SerialDenseVector &v, mesh &Mesh, std::vector< int > &phase)
void generator_one_gauss_point(Epetra_SerialDenseVector &v, mesh &Mesh, std::vector< int > &phase, double &xi, double &eta, double &zeta)