Erreur de segmentation

Question

C'est probablement une question stupide, mais je ne vois pas ce que je fais de mal ici. J'ai la classe:

#include <sys/time.h> 
#include <gsl/gsl_cdf.h> 
#include <gsl/gsl_randist.h> 
#include <cmath> 
#include "randomnumbergenerator.h" 

class RandomNumberGenerator 
{ 
    gsl_rng * rn; 
public: 
    RandomNumberGenerator(); 
    ~RandomNumberGenerator(); 
    double univariate(); 
    void bivariateGaussian(double rho, double &x, double &y); 
}; 

long currentMicroseconds() 
{ 
    struct timeval now; 
    gettimeofday(&now, NULL); 
    return now.tv_usec; 
} 

RandomNumberGenerator::RandomNumberGenerator() 
{ 
    const gsl_rng_type * T; 


    gsl_rng_env_setup(); 

    //T = gsl_rng_default; 
    T = gsl_rng_mt19937; 
    rn = gsl_rng_alloc (T); 
    gsl_rng_set(rn,currentMicroseconds()); 
} 

double RandomNumberGenerator::univariate() 
{ 
    return gsl_rng_uniform(rn); 
} 

void RandomNumberGenerator::bivariateGaussian(double rho, double &x, double &y) 
{ 
    gsl_ran_bivariate_gaussian (rn, 1.0, 1.0, rho, &x, &y); 
} 

RandomNumberGenerator::~RandomNumberGenerator() 
{ 
    gsl_rng_free (rn); 
} 

que j'appelle d'ici:

double x; 
double y; 
rng.bivariateGaussian(rho, x, y); 

mais je reçois une erreur de segmentation sur gsl_ran_bivariate_gaussian (rn, 1.0, 1.0, rho, &x, &y);

Toute idée?

Answer 1

Vérifiez si rn a bien été affecté. C'est probablement la seule chose qui peut causer une erreur de segmentation.

J'ai testé votre code sur mon ordinateur, il fonctionne correctement pour autant qu'il sache. Peut être vérifier l'installation de GSL, ils ont une suite de test, vous pouvez utiliser

Answer 2

Quel compilateur? Je suppose que rn est une variable membre de RandomNumberGenerator. L'initialisez-vous à 0 dans le constructeur? Vous ne semblez pas vérifier un retour d'erreur de gsl_rng_alloc, vous devriez probablement être parce que la seule chose que je peux voir tout de suite qui peut causer un problème est si rn ne pointe pas vers quelque chose de valide à l'appel qui segfaulting.

En regardant le manual for gsl_rng_alloc vous pouvez vérifier pour voir si elle renvoie NULL ou 0 et ensuite jeter une exception si ce n'est pas le cas. Par exemple:

#include <stdexcept> 

RandomNumberGenerator::RandomNumberGenerator() 
{ 
    const gsl_rng_type * T; 


    gsl_rng_env_setup(); 

    //T = gsl_rng_default; 
    T = gsl_rng_mt19937; 
    rn = gsl_rng_alloc (T); 
    if (rn == 0) { 
     throw ::std::runtime_error("Failed to allocation a random number generator."); 
    } 
    gsl_rng_set(rn,currentMicroseconds()); 
} 

, vous avez également essayé de compiler avec -O0 tourner de toute optimisation?

Answer 3

Dans:

double x; 
double y; 
rng.bivariateGaussian(rho, x, y); 

sont x et y peut-être censé être des tableaux plutôt que des variables simples? Je m'attendrais à ce qu'une distribution produise des valeurs N plutôt qu'une (ou deux).

Answer 4

Je ne programme pas C++, mais C. J'espère que cela s'appliquera aussi à vous. Mais sur C, j'utilise parfois un débogueur comme GDB ou le débogueur dans Eclipse. J'utilise aussi valgrind (j'aime vraiment beaucoup cet outil) pour réparer les fuites de mémoire/défauts de segmentation. Je vous conseille d'aimer à cette tutorial pour avoir une meilleure compréhension de ce que valgrind peut faire pour vous. Valgrind peut en faire beaucoup plus, je vous conseille donc de lire sur valgrind/helgrind.

Answer 5

Merci à tous pour vos réponses. Le bug était dans le morceau de code que je ne collais pas :(Je passais une instance de RandomNumberGenerator comme un paramètre normal.Quand je l'ai changé en passant comme référence il a commencé à fonctionner magiquement