From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 15749 invoked by alias); 16 Dec 2009 17:21:27 -0000 Received: (qmail 15739 invoked by uid 22791); 16 Dec 2009 17:21:26 -0000 X-SWARE-Spam-Status: No, hits=-1.8 required=5.0 tests=AWL,BAYES_00,SPF_PASS X-Spam-Check-By: sourceware.org Received: from mail.network-theory.co.uk (HELO mail.network-theory.co.uk) (66.199.228.187) by sourceware.org (qpsmtpd/0.43rc1) with ESMTP; Wed, 16 Dec 2009 17:21:23 +0000 Date: Wed, 16 Dec 2009 17:21:00 -0000 Message-ID: <87zl5ig7i0.wl%bjg@network-theory.co.uk> From: Brian Gough To: Rhys Ulerich Cc: gsl-discuss@sourceware.org Subject: Re: Where a generalized Richardson extrapolation routine would fit in GSL? In-Reply-To: <4a00655d0912150931q4b34fd24p94594ed08857254f@mail.gmail.com> References: <4a00655d0908201247g7d7bd9a1t466f4a66f08df4@mail.gmail.com> <4a00655d0911291536t5a11752fp27ab9c274148f822@mail.gmail.com> <4a00655d0911291538y9f29830v984d1a796fdd5d1@mail.gmail.com> <4a00655d0912131549w19638273nb51d723e9ddd9273@mail.gmail.com> <87d42ggnsv.wl%bjg@network-theory.co.uk> <4a00655d0912150931q4b34fd24p94594ed08857254f@mail.gmail.com> User-Agent: Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI) MIME-Version: 1.0 (generated by SEMI 1.14.6 - "Maruoka") Content-Type: text/plain; charset=US-ASCII X-Message-Mac: b813c3d3e40e759ba0246f61a2525a48 Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sourceware.org X-SW-Source: 2009-q4/txt/msg00065.txt.bz2 At Tue, 15 Dec 2009 11:31:10 -0600, Rhys Ulerich wrote: > You're right that each component could be extrapolated independently. > I wrote it using vectors because doing so allows using BLAS calls for > the linear algebra and speeds up many component use cases. I'd prefer > to keep the code vector-capable under the covers under the theory that > people extrapolating only a scalar at a time aren't all that worried > about speed. Ideally I'd like to start with a scalar version using normal C arrays, similar to the gsl_sum functions, for implementation simplicity. Also, are the "exact/normtable" arguments essential or just for convenience? If they can be computed easily by the user it would be good to have only a minimal number of arguments in the base function. -- Brian Gough