From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 5099 invoked by alias); 20 Aug 2003 11:52:47 -0000 Mailing-List: contact gsl-discuss-help@sources.redhat.com; run by ezmlm Precedence: bulk List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sources.redhat.com Received: (qmail 5091 invoked from network); 20 Aug 2003 11:52:47 -0000 Received: from unknown (HELO relay.sckcen.be) (193.190.187.211) by sources.redhat.com with SMTP; 20 Aug 2003 11:52:47 -0000 Received: from 193.190.187.215 by relay.sckcen.be (InterScan E-Mail VirusWall NT); Wed, 20 Aug 2003 13:39:40 +0200 Received: from pc1727.sck.be ([193.190.184.23]) by mailsrv3.sck.be with Microsoft SMTPSVC(5.0.2195.5329); Wed, 20 Aug 2003 13:39:39 +0200 From: Gert Van den Eynde Reply-To: gvdeynde@sckcen.be To: Brian Gough , Axel Hutt Subject: Re: integration Date: Wed, 20 Aug 2003 11:52:00 -0000 User-Agent: KMail/1.5.3 Cc: gsl-discuss@sources.redhat.com References: <3F3CABD4.6060707@wias-berlin.de> <16195.22145.205768.677021@debian.local> In-Reply-To: <16195.22145.205768.677021@debian.local> MIME-Version: 1.0 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: 7bit Content-Disposition: inline Message-Id: <200308201339.19421.gvdeynde@sckcen.be> X-OriginalArrivalTime: 20 Aug 2003 11:39:39.0811 (UTC) FILETIME=[BD824F30:01C3670F] X-SW-Source: 2003-q3/txt/msg00162.txt.bz2 On Wednesday 20 August 2003 13:07, Brian Gough wrote: > The only thing I can contribute here is that all of the quadpack > routines assume smooth functions. Apart from that, I don't know much > about integro-differential equations. A small addendum: QUADPACK (the original FORTRAN code) allows functions with end-point singularities (routine DQAGS(E)). You would do that anyway when integrating a function with singularties: locate the singularities and split up the integration interval accordingly. With the routine DQAGP(E) you can even provide information on the kind of singularty.I suppose the GSL routine is the C translation of DQAGS. Of course, the function inside the integration boundaries should be smooth (numerical integration is for a big part based on the approximation of the function by orthogonal polynomials). For an overview of all original QUADPACK routines, have a look here http://netlib.org/quadpack/doc Hope this helps, gert