From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 122609 invoked by alias); 18 Nov 2015 01:50:10 -0000 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 Received: (qmail 122417 invoked by uid 89); 18 Nov 2015 01:50:10 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-2.6 required=5.0 tests=AWL,BAYES_05,FREEMAIL_FROM,RCVD_IN_DNSWL_LOW,RP_MATCHES_RCVD,SPF_PASS autolearn=ham version=3.3.2 X-HELO: BAY004-OMC1S13.hotmail.com Received: from bay004-omc1s13.hotmail.com (HELO BAY004-OMC1S13.hotmail.com) (65.54.190.24) by sourceware.org (qpsmtpd/0.93/v0.84-503-g423c35a) with (AES256-SHA256 encrypted) ESMTPS; Wed, 18 Nov 2015 01:50:08 +0000 Received: from BAY176-W4 ([65.54.190.59]) by BAY004-OMC1S13.hotmail.com over TLS secured channel with Microsoft SMTPSVC(7.5.7601.23008); Tue, 17 Nov 2015 17:50:07 -0800 X-TMN: [83/UOjz5Z076D51bye3HWaHBTGyq04m+] Message-ID: From: Bill Maier To: "gsl-discuss@sourceware.org" Subject: Bug #45926 Date: Wed, 18 Nov 2015 01:50:00 -0000 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable MIME-Version: 1.0 X-SW-Source: 2015-q4/txt/msg00009.txt.bz2 I've investigated bug #45929 "error in generalized hypergeometric function"= . The problem is in the file specfunc/hyperg_2F1.c in the first function the file, hyperg_2F1_serie= s(). There is a loop starting at line 59 which is intended to converge to the solution, and within this l= oop is a check to insure the loop is not executed more than 30000 times. This code is actually working a= nd computing the correct value, but the error terms del_pos and del_neg eventually get no smaller an= d so the loop executes the full 30000 iterations and then reports GSL_EMAXITER error, even though it d= oes have the correct value calculated. I would fix this by saving the previous error terms and insuring they are b= ecoming smaller with each iteration. If they are not getting smaller, machine precision has been reac= hed and the value can be returned to the user. I can supply code or could check in myself through gi= t if I can get access. I have 30+ years of experience with C code and want to begin helping maintain GSL. By the way, when I do a git clone to get the code and then run ./autogen, i= t fails with message "autoreconf: not found". Has this file been left out of the git repository? -Bill Maier=20=09=09=20=09=20=20=20=09=09=20=20