On Tue, May 19, 2009 at 8:52 PM, Alex Tsoi wrote:

Hi all,

I am trying to write a MCMC algorithm by using boost's math distribution library, what I want to do is to first declare an array of length 10 of gamma distribution, and then run a for loop to declare each gamma distribution's parameters. However, boost does not allow me to declare the gamma distribution array without having defining the parameter in advance:

gamma_distribution<double> mydist10[10];

error: no appropriate default constructor available.

because the library requires to know the parameter of the gamma distribution in advance.

Do you guys have any idea to overcome this ?

Use a gamma_distribution<double>* as contained element. For instance: gamma_distribution<double>*. mydist10_1[10]: mydist10_1[0] = new gamma_distribution<double>(1.2); mydist10_1[1] = new gamma_distribution<double>(1.3); //... or: std::vector mydist10_2(10); mydist10_2[0] = new gamma_distribution<double>(1.2); mydist10_2[1] = new gamma_distribution<double>(1.3); //... or (maybe better): std::vector< boost::shared_ptr< gamma_distribution<double> > > mydist10_3(10): mydist10_3[0] = boost::shared_ptr< gamma_distribution<double> >(new gamma_distribution<double>(1.2)); mydist10_3[1] = boost::shared_ptr< gamma_distribution<double> >(new gamma_distribution<double>(1.3)); //... Cheers, -- Marco

Thanks a lot for any comment or suggestion.

Alex-

-- Lam C. Tsoi (Alex) Medical University of South Carolina

_______________________________________________ Boost-users mailing list Boost-users@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users