From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 15405 invoked by alias); 29 Jun 2005 18:31:01 -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 15360 invoked by uid 22791); 29 Jun 2005 18:30:54 -0000 Received: from mail.phy.duke.edu (HELO mail.phy.duke.edu) (152.3.182.2) by sourceware.org (qpsmtpd/0.30-dev) with ESMTP; Wed, 29 Jun 2005 18:30:54 +0000 Received: from lilith (client212-5.dsl.intrex.net [209.42.212.5]) by mail.phy.duke.edu (Postfix) with ESMTP id 131A993289 for ; Wed, 29 Jun 2005 14:30:49 -0400 (EDT) Message-ID: From: "Robert G. Brown" To: gsl-discuss@sources.redhat.com Subject: multidimensional integration Date: Wed, 29 Jun 2005 18:31:00 -0000 Mime-Version: 1.0 Content-Type: multipart/signed; boundary="=_mimegpg-lilith-2075-1120069503-0013"; micalg=pgp-sha1; protocol="application/pgp-signature" X-SW-Source: 2005-q2/txt/msg00029.txt.bz2 This is a MIME GnuPG-signed message. If you see this text, it means that your E-mail or Usenet software does not support MIME signed messages. --=_mimegpg-lilith-2075-1120069503-0013 Content-Type: text/plain; format=flowed; charset="US-ASCII" Content-Disposition: inline Content-Transfer-Encoding: 7bit Content-length: 2333 Hi fellow GSL'ers. We have a postdoc in our department who is preparing to integrate something. In his previous position at another place, he used NAG to do this, and has the requisite code already in place. He requested that we buy and install a single copy of NAG just for him and a student to be able to use this one routine to do this one integral on just one computer, at a cost of many hundreds of dollars. I suggested that he look into using the GSL instead, since it is a very high-quality library to my own direct experience and of course is both free and universally installed in our department. GSL and NAG both use QUADPACK as the basis for their 1D integrals (and have nearly identical call structure) so I figured that the transition would actually be painless. However, the integrand he has to integrate is actually defined and integrated over somewhere between 5 to 7 dimensions (with rectangular limits). The routine he used from NAG was actually d01fcc, which is NOT from QUADPACK but rather implements the multidimensional adaptive routine HALF with a custom interval rule. When I looked at GSL's online manual (version 1.6 as of this last December) I didn't see a multidimensional integration routine equivalent to d01fcc. SO, questions: a) Is a multidimensional integration routine equivalent to d01fcc implemented or under development, and if so, where is it and/or how do I get a version that has it? I looked at the CVS tree and didn't immediately see one. In principle I could probably use e.g. a multidimentional ODE solver but I'd think that having a d01fcc equivalent would be much more efficient. b) If not, does anybody have any suggestions on the "best" way to attack this sort of integral using existing tools? At five dimensions I suspect that just calling 1 dim integrations five levels deep would result in an awful lot of wasted energy and time. Framing it as an ODE set also seems like it would work but likely not be terribly efficient or terribly easy to control error-wise. c) On a related note, has anybody done a head-to-head performance comparison of GSL with NAG -- either time/efficiency performance or numerical accuracy type performance? This isn't a significant issue on this particular project but is an issue that I expect to see come up in the future. rgb --=_mimegpg-lilith-2075-1120069503-0013 Content-Type: application/pgp-signature Content-Transfer-Encoding: 7bit Content-length: 189 -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.2.6 (GNU/Linux) iD8DBQBCwud/UfPgIH3s/SIRArQmAKCv9m7nbCtRezl6XxIUEDd8UPW8ygCfeQNZ +Rnzr3VGf28ngHK8AT/BHTw= =KQTp -----END PGP SIGNATURE----- --=_mimegpg-lilith-2075-1120069503-0013--