From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 3781 invoked by alias); 17 Sep 2009 19:14:19 -0000 Received: (qmail 3758 invoked by uid 22791); 17 Sep 2009 19:14:17 -0000 X-SWARE-Spam-Status: No, hits=-0.7 required=5.0 tests=AWL,BAYES_40,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; Thu, 17 Sep 2009 19:14:14 +0000 Date: Thu, 17 Sep 2009 19:14:00 -0000 Message-ID: From: Brian Gough To: Jonathan Underwood Cc: gsl-discuss@sourceware.org Subject: Re: GSL 2.0 roadmap (one man's view) In-Reply-To: <645d17210909090818u474f32f0q19a6334578b9f02c@mail.gmail.com> References: <48E25CA9.6080306@iki.fi> <498727E5.6080407@iki.fi> <49AA9DB5.6030908@iki.fi> <49FB01D1.30000@iki.fi> <4A7ADFDC.9080408@iki.fi> <1251414774.23092.80.camel@manticore.lanl.gov> <1251414939.23092.82.camel@manticore.lanl.gov> <645d17210909090818u474f32f0q19a6334578b9f02c@mail.gmail.com> User-Agent: Wanderlust/2.14.0 (Africa) Emacs/22.2 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: c32ac735f17e8ec234ffa7c1f72bbe9b 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-q3/txt/msg00052.txt.bz2 At Wed, 9 Sep 2009 16:18:54 +0100, Jonathan Underwood wrote: > It'd be nice if it was possible to implement an approach such as > CESTAC/CADNA in GSL ..... however, I don't see how to do it in C, as > it seems to require overloading of the basic operators (+/-* etc). Yes that would be the ideal way to do it, with the error estimate determined automatically, but as you suggest C is not capable of handling that--it would have to be done by code generation or a parser that preprocesses the source (which is how this sort of thing used to be done in FORTRAN 77 I believe). I found the SIAM book "Accuracy and reliability in scientific computing" a useful reference on this subject, although it is a few years old now so it may not cover the latest developments. -- Brian Gough