[Help-gsl] gsl: ..\permutation\init.c:51: ERROR

[Mohamed Ali]

Hello,

I am using GSL on  CodeBlocks with Mingw. I have a program working in
debug mode but not release !!


gsl: ..\permutation\init.c:51: ERROR: failed to allocate space for
permutation data
Default GSL error handler invoked.

This application has requested the Runtime to terminate it in an unusual way.
Please contact the application's support team for more information.




void Snake::InitInternalEnergyMatrix()
{
    double a = 2 * Alpha + 6 * Beta;
    double b = - Alpha -4 * Beta;
    double c = Beta;
    try{
    gsl_matrix *InternalEnergyMatrix = gsl_matrix_calloc( Size , Size );
    for (int i=0;i<Size;i++)
    {
        gsl_matrix_set(InternalEnergyMatrix, i, i, a + 1); // May be 1 * Gamma
        if (i<Size-1)gsl_matrix_set(InternalEnergyMatrix, i, i + 1, b);
        if (i<Size-2)gsl_matrix_set(InternalEnergyMatrix, i, i + 2, c);
        if (i>0)gsl_matrix_set(InternalEnergyMatrix, i, i - 1, b);
        if (i>1)gsl_matrix_set(InternalEnergyMatrix, i, i - 2, c);
    }
    gsl_matrix_set(InternalEnergyMatrix, 0, Size-1, b);
    gsl_matrix_set(InternalEnergyMatrix, 0, Size-2, c);
    gsl_matrix_set(InternalEnergyMatrix, 1, Size-1, c);
    gsl_matrix_set(InternalEnergyMatrix, Size-2, 0, c);
    gsl_matrix_set(InternalEnergyMatrix, Size-1, 0, b);
    gsl_matrix_set(InternalEnergyMatrix, Size-1, 1, c);
    InternalEnergyMatrixInv = gsl_matrix_calloc( Size , Size );
    gsl_permutation *p = gsl_permutation_calloc( Size );
    gsl_linalg_LU_invert ( InternalEnergyMatrix , p , InternalEnergyMatrixInv );
    gsl_permutation_free( p );
    gsl_matrix_free( InternalEnergyMatrix );
    }
    catch(...){
        Log::Write("Snake::InitInternalEnergyMatrix",Config::FLog, 0);
    }
}