0

私はここでこのコードを確認しています：

#include <stdio.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_permutation.h>



int main ( ) {
    const long N = 10;
    const gsl_rng_type *T = NULL;
    gsl_rng *r = NULL;

    gsl_permutation *p = NULL;
    p = gsl_permutation_alloc ( N );
    gsl_permutation *q = NULL;
    q = gsl_permutation_alloc ( N );

    gsl_rng_env_setup ();
    T = gsl_rng_default;
    r = gsl_rng_alloc ( T );

    printf ( "initial permutation: " );
    gsl_permutation_init ( p );
    gsl_permutation_fprintf ( stdout, p, " %u" );
    printf ( "\n" );

    printf (" random permutation:");  
    gsl_ran_shuffle (r, p->data, N, sizeof(size_t));
    gsl_permutation_fprintf (stdout, p, " %u");
    printf ("\n");

    printf ( "inverse permutation: " );
    gsl_permutation_inverse ( q, p );
    gsl_permutation_fprintf ( stdout, q, "%u" );
    printf ( "\n" );

    gsl_permutation_free ( p );
    gsl_permutation_free ( q );

    return 0;
}

ただし、コンパイルしようとすると、未処理の例外が発生します。アクセス違反の書き込み場所エラー。私はすべてのポインタにメモリを割り当てたと思いますが、それでもエラーは解決しません。私はCにあまり詳しくないので、どんな入力でも歓迎されます。PSデバッガーを使用すると、次の行を指します。gsl_permutation_fprintf ( stdout, p, " %u" ); ありがとう

編集：私はコードをいじってみて、問題がそこのstdoutステートメントにある可能性があることに気づきました。

4

1 に答える 1

0

gsl_rng_env_setup ();最初に呼ばれるべきではありませんか？あなたはgsl_permutation_alloc前に2、3をします。

于 2011-10-13T05:59:35.543 に答える