Re: [Help-gsl] QAWS & conditional IF / ELSE definition in integrand function

[jeremy theler]

yes, although it would be a good idea to avoid generating NaNs by evaluating f 
after checking the sign of x (and also the sign of the product alpha*x)

On Wednesday 17 July 2013 14:02:45 Mátyás Benke wrote:
> Hi guys,
> 
> I have a quick question regarding the QAWS algorithm.
> Does QAWS allow you to define the integrand function using IF / ELSE
> conditions. For example, will it work in the following way?
> 
> double f (double x, void * params)
> {
>        double alpha = *(double *) params;
>        double f = log(alpha*x) / sqrt(x);
>        if (x > 0.0)
>            return f;
>        else
>            return 0.0;
> }Cheers, Matyas