From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <gsl-discuss-return-2711-listarch-gsl-discuss=sourceware.cygnus.com@sources.redhat.com>
Received: (qmail 15596 invoked by alias); 2 Oct 2002 07:03:53 -0000
Mailing-List: contact gsl-discuss-help@sources.redhat.com; run by ezmlm
Precedence: bulk
List-Subscribe: <mailto:gsl-discuss-subscribe@sources.redhat.com>
List-Archive: <http://sources.redhat.com/ml/gsl-discuss/>
List-Post: <mailto:gsl-discuss@sources.redhat.com>
List-Help: <mailto:gsl-discuss-help@sources.redhat.com>, <http://sources.redhat.com/ml/#faqs>
Sender: gsl-discuss-owner@sources.redhat.com
Received: (qmail 15587 invoked from network); 2 Oct 2002 07:03:52 -0000
Received: from unknown (HELO relay.sckcen.be) (193.190.187.211)
  by sources.redhat.com with SMTP; 2 Oct 2002 07:03:52 -0000
Received: from 193.190.143.188 by relay.sckcen.be (InterScan E-Mail VirusWall NT); Wed, 02 Oct 2002 09:03:47 +0200
Received: from pc1727.sck.be ([193.190.184.23]) by mailsrv3.sck.be over TLS secured channel with Microsoft SMTPSVC(5.0.2195.5329);
	 Wed, 2 Oct 2002 09:03:47 +0200
Received: from pc1727.sck.be (localhost [127.0.0.1])
	by pc1727.sck.be (8.12.3/8.12.3/SuSE Linux 0.6) with SMTP id g9273fut029390;
	Wed, 2 Oct 2002 09:03:41 +0200
Date: Wed, 02 Oct 2002 00:03:00 -0000
From: Gert Van den Eynde <gvdeynde@sckcen.be>
To: "Robert G. Brown" <rgb@phy.duke.edu>
Cc: gsl-discuss@sources.redhat.com
Subject: Re: Laguerre Polynomial overflow...
Message-Id: <20021002090340.413f6a83.gvdeynde@sckcen.be>
In-Reply-To: <Pine.LNX.4.44.0210011421180.25903-100000@ganesh.phy.duke.edu>
References: <Pine.LNX.4.44.0210011421180.25903-100000@ganesh.phy.duke.edu>
Organization: SCK-CEN
Mime-Version: 1.0
Content-Type: multipart/mixed;
 boundary="Multipart_Wed__2_Oct_2002_09:03:40_+0200_082b9010"
X-OriginalArrivalTime: 02 Oct 2002 07:03:47.0545 (UTC) FILETIME=[DA937490:01C269E1]
X-SW-Source: 2002-q4/txt/msg00001.txt.bz2

This is a multi-part message in MIME format.

--Multipart_Wed__2_Oct_2002_09:03:40_+0200_082b9010
Content-Type: text/plain; charset=US-ASCII
Content-Transfer-Encoding: 7bit
Content-length: 197


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


--Multipart_Wed__2_Oct_2002_09:03:40_+0200_082b9010
Content-Type: text/plain;
 name="lag.c"
Content-Disposition: attachment;
 filename="lag.c"
Content-Transfer-Encoding: 7bit
Content-length: 729

#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;
}


--Multipart_Wed__2_Oct_2002_09:03:40_+0200_082b9010--