Re: Laguerre Polynomial overflow...

[Gert Van den Eynde <gvdeynde@sckcen.be>]

[-- Attachment #1: Type: text/plain, Size: 197 bytes --]


Dear Robert,

I cannot reproduce this using gsl-1.2 (my own compilation) on a SuSE 8.0 linux machine. Could you try to compile and run the attached program and see if that fails too?

Bye,

Gert


[-- Attachment #2: lag.c --]
[-- Type: text/plain, Size: 729 bytes --]

#include <stdlib.h>
#include <stdio.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_sf_laguerre.h>

int main () {

    double x = 0.0, step = 0.2;
    int    k = 0, l = 0;

    gsl_sf_result res;

    for (k = 0; k < 4; k ++) {
        for (x = 4.0; x < 6.1; x += step)  {
            gsl_sf_laguerre_n_e(k,0.5,x,&res);
            printf("L_%d(%.1f) = %.8f        error = %.8f\n",k,x,res.val,res.err);
        }
        printf("\n");
    }


    gsl_sf_laguerre_n_e(2,0.5,4.97290e+00,&res);
    printf("L_%d(%.5f) = %.8f        error = %.8f\n",2,4.97290e+00,res.val,res.err);

    gsl_sf_laguerre_n_e(3,0.5,4.97290e+00,&res);
    printf("L_%d(%.5f) = %.8f        error = %.8f\n",3,4.97290e+00,res.val,res.err);

    return 0;
}