bon j'ai mis ça il me sort un segmentation fault
double determ_tps_inf(double beta, int nbi,int N,int nbs)
{
const gsl_rng_type * T1;
gsl_rng * r1;
double mu,tps_infection,Nd,inter;
int i;
time_t * s;
Nd = (double) N;
// c'est la moyenne de la loi exponentielle
mu = (double) beta*nbi/Nd;
/* create a generator chosen by the
environment variable GSL_RNG_TYPE */
gsl_rng_env_setup();
s = malloc(sizeof(time_t));
T1 = gsl_rng_default;
r1 = gsl_rng_alloc (T1);
time(s);
gsl_rng_set (r1,(*s)); // segmentation fault sur cette ligne
// on tire un nombre dans la loi gamma, et on le place dans le vecteur tps_period
inter = gsl_ran_exponential(r1,mu);
// on libère la mémoire allouée pour r1
gsl_rng_free (r1);
return inter;
}
je pige pas
Message édité par Profil supprimé le 15-05-2007 à 11:15:11