Re: Fw : Re: [Help-gsl] segmentation fault using qag.h with GSL

[Thomas Weber]

Am Dienstag, den 12.08.2008, 08:52 +0000 schrieb chadia kanaan:
> 
> --- En date de : Mar 12.8.08, chadia kanaan <address@hidden> a écrit :
> De: chadia kanaan <address@hidden>
> Objet: Re: [Help-gsl] segmentation fault using qag.h with GSL
> À: "Thomas Weber" <address@hidden>
> Date: Mardi 12 Août 2008, 11h49
> 
> Dear Thomas,
> 
> I re-enclose again the error file concerning the debugging in a plain
>  text format ( for the convenience of kind users that may help me ).
>  You'll find also attached the code file "code.c" . I really appreciate
>  your help, many thanks,  also for Joel :


Line 31 from gsl's integration/initialise.c:
          workspace->rlist[0] = 0.0;

Okay, let's look at your workspace w (code edited for the mail)
        double p_simfunc(double z, int SFR_index)
               {
                double p_sim ;
                double int_p_unor, err;
                /*normalize p_unor*/
        
                gsl_integration_workspace *w = 
                        gsl_integration_workspace_alloc(10000);
                gsl_integration_workspace_free(w);
                ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
                ...

If you want to use your workspace, you shouldn't free it right after
allocating :)

        Thomas