From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 30639 invoked by alias); 14 Jun 2005 21:09:07 -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 30588 invoked by uid 22791); 14 Jun 2005 21:09:04 -0000 Received: from network-theory.com (HELO mail.network-theory.co.uk) (66.199.228.187) by sourceware.org (qpsmtpd/0.30-dev) with ESMTP; Tue, 14 Jun 2005 21:09:04 +0000 Received: from dsl-80-42-160-45.access.as9105.com ([80.42.160.45]:35131 helo=hp.network-theory.co.uk) by mail.network-theory.co.uk with esmtpa (Exim 4.44) id 1DiIcx-00072z-Eh; Tue, 14 Jun 2005 22:07:19 +0100 Received: by hp.network-theory.co.uk with local (Exim 3.35 #1 (Debian)) id 1DiIeT-0004wl-00; Tue, 14 Jun 2005 22:08:53 +0100 MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Message-ID: <17071.18277.391105.281999@network-theory.co.uk> Date: Tue, 14 Jun 2005 21:09:00 -0000 From: Brian Gough To: stevenj@math.mit.edu Cc: gsl-discuss@sources.redhat.com Subject: Re: contribution: code for multi-dimensional adaptive integration In-Reply-To: References: X-SW-Source: 2005-q2/txt/msg00021.txt.bz2 Steven G. Johnson writes: > For some time now, I've felt that GSL has cried out for a multidimensional > cubature routine, similar to its 1d adaptive quadrature routines. Now, > I've written one (~700 lines in C) based on the well-known Genz-Malik > embedded cubature rule (with help from a GPL'ed library called HIntLib, > see below), and I hope that you're interested to have it included in GSL. > I've attached the unmodified code; I'm willing to do the work to modify it > for your coding style, etcetera. Thanks.. that does look useful. You're right that cubature is an obvious gap in GSL, and I would like to see that filled. The first step would be make a package that builds outside of gsl, there are some instructions at http://sources.redhat.com/gsl/devel.html (notes on the coding conventions can also be found in the design document linked there) -- Brian Gough p.s. Sorry for the delay in replying.