Am Dienstag, den 12.08.2008, 08:52 +0000 schrieb chadia kanaan:
>
> --- En date de : Mar 12.8.08, chadia kanaan <[EMAIL PROTECTED]> a écrit :
> De: chadia kanaan <[EMAIL PROTECTED]>
> Objet: Re: [Help-gsl] segmentation fault using qag.h with GSL
> À: "Thomas Weber" <[EMAIL PROTECTED]>
> 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
_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl
