From MAILER-DAEMON Fri Aug 06 09:09:16 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OhMfk-0004ID-4b for mharc-bug-gsl@gnu.org; Fri, 06 Aug 2010 09:09:16 -0400 Received: from [140.186.70.92] (port=38169 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OhMfh-0004Gw-62 for bug-gsl@gnu.org; Fri, 06 Aug 2010 09:09:14 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1OhMfg-0005ht-4y for bug-gsl@gnu.org; Fri, 06 Aug 2010 09:09:13 -0400 Received: from colonialone.fsf.org ([140.186.70.51]:51378 helo=internal.in.savannah.gnu.org) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OhMfg-0005ho-3k for bug-gsl@gnu.org; Fri, 06 Aug 2010 09:09:12 -0400 Received: from [10.1.0.103] (helo=frontend.in.savannah.gnu.org) by internal.in.savannah.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OhMff-00057H-RC; Fri, 06 Aug 2010 13:09:11 +0000 Received: from www-data by frontend.in.savannah.gnu.org with local (Exim 4.69) (envelope-from ) id 1OhMff-0004dd-Pt; Fri, 06 Aug 2010 13:09:11 +0000 To: Sergey B Kirpichev , bug-gsl@gnu.org From: Brian Gough X-Savane-Server: savannah.gnu.org:443 [10.1.0.103] MIME-Version: 1.0 Content-Type: text/plain;charset=UTF-8 X-Savane-Project: gsl X-Savane-Tracker: bugs X-Savane-Item-ID: 30540 User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-GB; rv:1.9.0.11) Gecko/2009061118 Fedora/3.0.11-1.fc9 Firefox/3.0.11 X-Apparently-From: 94.193.216.241 (Savane authenticated user bjg) Message-Id: <20100806-140911.sv1824.76026@savannah.gnu.org> References: <20100724-121507.sv79421.66062@savannah.gnu.org> In-Reply-To: <20100724-121507.sv79421.66062@savannah.gnu.org> Date: Fri, 06 Aug 2010 13:09:11 +0000 X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) Cc: Subject: [Bug-gsl] [bug #30540] please, add convergence checks in rk4imp/rk2imp X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Fri, 06 Aug 2010 13:09:14 -0000 Update of bug #30540 (project gsl): Status: None => Ready For Test _______________________________________________________ Follow-up Comment #1: Thanks for the bug report. I'm merging a completely new and upgraded version of the ode solvers "ode-initval2" written by Tuomo Keskitalo in the next release. It includes the new solver ode-initval2/imprk4.c for implicit RK4. You can get the branch with bzr checkout http://bzr.savannah.gnu.org/r/gsl/ode-initval2/ To compile it use autoreconf -i -f; ./configure --enable-maintainer-mode ; make Please try it and let us know if it solves this problem. _______________________________________________________ Reply to this item at: _______________________________________________ Message sent via/by Savannah http://savannah.gnu.org/ From MAILER-DAEMON Thu Aug 12 03:08:21 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OjRtl-0006wA-2X for mharc-bug-gsl@gnu.org; Thu, 12 Aug 2010 03:08:21 -0400 Received: from [140.186.70.92] (port=52662 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OjLRO-0004S9-1L for bug-gsl@gnu.org; Wed, 11 Aug 2010 20:14:39 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1OjLRM-0005PS-Ic for bug-gsl@gnu.org; Wed, 11 Aug 2010 20:14:37 -0400 Received: from mail-wy0-f169.google.com ([74.125.82.169]:44057) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OjLRM-0005P7-Ek for bug-gsl@gnu.org; Wed, 11 Aug 2010 20:14:36 -0400 Received: by wyg36 with SMTP id 36so830808wyg.0 for ; Wed, 11 Aug 2010 17:14:35 -0700 (PDT) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:sender:received:date :x-google-sender-auth:message-id:subject:from:to:cc:content-type; bh=KBI6svCHUM7NOsx8wg8xulYfl0CA814WzTTa7QrOOQY=; b=tEZsn4BFE5L2fOAb/sw9U+d2hTkfQbZFhDkgpHOY/+pQ4L6rBhEejkqo4y2/f9jtFF Vx6jsUUjCczI2EiKBDlK/MqyIy2n5PjT1eeGyoygjc/d7/Gd0s8cyLzYfp/lRJvTe9cR Ynj0BjBjvH7N+GJujdIVAW2p2EkzxN49M78hA= DomainKey-Signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:sender:date:x-google-sender-auth:message-id:subject :from:to:cc:content-type; b=NyfiQAYKf6RYTeUrfqoRB2O8iBGmAzjEgi8klcCstPhYoWqkgS2BHaB80uQfjcvFEx 5xvTVMU9iod46cG0f1/xIBA1msJg7gUYKD7oYrdJU+vYLimQXBXW4q9gKBDZygXIrGp9 F5YVbJZ2JrJpDSyV9iAS8um9th9EqpU33mkoU= MIME-Version: 1.0 Received: by 10.216.28.213 with SMTP id g63mr6008511wea.71.1281572075034; Wed, 11 Aug 2010 17:14:35 -0700 (PDT) Sender: thanhvv@gmail.com Received: by 10.216.70.135 with HTTP; Wed, 11 Aug 2010 17:14:34 -0700 (PDT) Date: Wed, 11 Aug 2010 17:14:34 -0700 X-Google-Sender-Auth: F-D2vJTpnlYaVP50naFPJMaY-wc Message-ID: From: Thanh Vo To: gsl-discuss@sourceware.org, bug-gsl@gnu.org Content-Type: text/plain; charset=UTF-8 X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) X-Mailman-Approved-At: Thu, 12 Aug 2010 03:08:19 -0400 Cc: "Earl T. Barr" Subject: [Bug-gsl] Potential Bugs? X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Thu, 12 Aug 2010 00:14:39 -0000 Hi, My name's Thanh Vo. My research group is developing a tool to detect floating point exceptions in numerical code. We've found inputs that cause some of the GSL 1.14 functions to throw floating-point exceptions. I've written this post to see if you consider any of them to be real bugs. For each exception, we describe the name of function, the type of exception thrown, the inputs that trigger the exception, and specific location of the code that causes the exception. 1. The function airy_aie throws a divide by zero exception when its input x = 0 at line 618 in file gsl/specfunc/airy.c : double z = 2.0/(x*sqx) - 1.0;. It throws an invalid exception when x = -1 at line 617: double sqx = sqrt(x);. 2. The function airy_bie throws a divide by zero exception when its input x = 0 at line 635 in file gsl/specfunc/airy.c: double z = ATR/(x*sqx) + BTR;. It throws an invalid exception when x=-1 at line 634: double sqx = sqrt(x);. 3. The function gsl_sf_asymp_thetanu_corr_e throws a divide by zero exception when its input x = 0 and nu != 0 at the line 181 in file gsl/specfunc/bessel_amp_phase.c: const double r = 2.0*nu/x;. It throws an invalid exception when x=0 and nu=0 at line 181: const double r = 2.0*nu/x;. 4. The function gsl_sf_bessel_asymp_Mnu_e throws a divide by zero exception when its input x=0 and nu !=0 at the line 167 in file gsl/specfunc/bessel_amp_phase.c: const double r = 2.0*nu/x;. It throws an invalid exception when x=0 and nu=0 at line 167: const double r = 2.0*nu/x;. 5. The function gsl_sf_bessel_J0 throws a overflow exception when its input x = 5.85177e+145 at the line 88 in the file gsl/specfunc/bessel_J0.c: const double z = 32.0/(y*y) - 1.0;. It throws an underflow exception when x = 6.8481e-155 at the line 81: result->err = y*y;. If anyone could confirm whether any of these exceptions is an actual bug, I would appreciate it. Looking forward to hearing from you, Best regards, Thanh Vo. From MAILER-DAEMON Tue Aug 17 16:38:39 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OlSvf-0001WZ-OV for mharc-bug-gsl@gnu.org; Tue, 17 Aug 2010 16:38:39 -0400 Received: from [199.232.76.173] (port=52036 helo=monty-python.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OlSvd-0001WQ-Rx for bug-gsl@gnu.org; Tue, 17 Aug 2010 16:38:37 -0400 Received: from Debian-exim by monty-python.gnu.org with spam-scanned (Exim 4.60) (envelope-from ) id 1OlSvb-0005Uv-PC for bug-gsl@gnu.org; Tue, 17 Aug 2010 16:38:37 -0400 Received: from fencepost.gnu.org ([140.186.70.10]:52193) by monty-python.gnu.org with esmtp (Exim 4.60) (envelope-from ) id 1OlSvb-0005Ur-Ho for bug-gsl@gnu.org; Tue, 17 Aug 2010 16:38:35 -0400 Received: from localhost ([127.0.0.1]:52810 helo=fencepost.gnu.org) by fencepost.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OlSva-00069S-Tb; Tue, 17 Aug 2010 16:38:34 -0400 Date: Tue, 17 Aug 2010 21:38:19 +0100 Message-ID: <43mxsl6klg.wl%bjg@gnu.org> From: Brian Gough To: Thanh Vo Subject: Re: [Bug-gsl] Potential Bugs? In-Reply-To: References: User-Agent: Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) MIME-Version: 1.0 (generated by SEMI 1.14.6 - "Maruoka") Content-Type: text/plain; charset=US-ASCII X-detected-operating-system: by monty-python.gnu.org: GNU/Linux 2.6 (newer, 3) Cc: gsl-discuss@sourceware.org, bug-gsl@gnu.org, "Earl T. Barr" X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 17 Aug 2010 20:38:38 -0000 At Wed, 11 Aug 2010 17:14:34 -0700, Thanh Vo wrote: > My name's Thanh Vo. My research group is developing a tool to detect > floating point exceptions in numerical code. We've found inputs that > cause some of > the GSL 1.14 functions to throw floating-point exceptions. I've written > this post to see if you consider any of them to be real bugs. Hi, Thanks for your email. I have made some comments below. This sounds like a useful tool, I hope you will be able to release it as free software. > > 1. The function airy_aie throws a divide by zero exception when its > input x = 0 at line 618 in file gsl/specfunc/airy.c : double z = > 2.0/(x*sqx) - 1.0;. It throws an invalid exception when x = -1 at > line 617: double sqx = sqrt(x);. > > 2. The function airy_bie throws a divide by zero exception when its > input x = 0 at line 635 in file gsl/specfunc/airy.c: double z = > ATR/(x*sqx) + BTR;. It throws an invalid exception when x=-1 at line > 634: double sqx = sqrt(x);. These would be bugs but the functions are static and only called with x > 1 or x > 2 so the exceptions cannot occur in practice. > 3. The function gsl_sf_asymp_thetanu_corr_e throws a divide by zero > exception when its input x = 0 and nu != 0 at the line 181 in file > gsl/specfunc/bessel_amp_phase.c: const double r = 2.0*nu/x;. It > throws an invalid exception when x=0 and nu=0 at line 181: const > double r = 2.0*nu/x;. > > 4. The function gsl_sf_bessel_asymp_Mnu_e throws a divide by zero > exception when its input x=0 and nu !=0 at the line 167 in file > gsl/specfunc/bessel_amp_phase.c: const double r = 2.0*nu/x;. It > throws an invalid exception when x=0 and nu=0 at line 167: const > double r = 2.0*nu/x;. It's a valid warning. The functions are publicly accessible, but not documented. Within the library they are only called with large values of x >= 1000 so the exception does not occur in practice. Ideally they should be static functions. > 5. The function gsl_sf_bessel_J0 throws a overflow exception when its > input x = 5.85177e+145 at the line 88 in the file > gsl/specfunc/bessel_J0.c: const double z = 32.0/(y*y) - 1.0;. It > throws an underflow exception when x = 6.8481e-155 at the line 81: > result->err = y*y;. The overflow exception is a valid warning. It's also a bug in the sense that it is avoidable (although the end result is correct, z=-1 since 1/inf=0). It would be better written as (32/y)/y. For the underflow it's a valid warning. There isn't really any alternative to returning zero in that case though. -- Brian Gough From MAILER-DAEMON Tue Aug 17 16:47:51 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OlT4Z-0003NJ-42 for mharc-bug-gsl@gnu.org; Tue, 17 Aug 2010 16:47:51 -0400 Received: from [199.232.76.173] (port=50854 helo=monty-python.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OlT4X-0003Mz-K1 for bug-gsl@gnu.org; Tue, 17 Aug 2010 16:47:49 -0400 Received: from Debian-exim by monty-python.gnu.org with spam-scanned (Exim 4.60) (envelope-from ) id 1OlT4W-0005ZB-EC for bug-gsl@gnu.org; Tue, 17 Aug 2010 16:47:49 -0400 Received: from fencepost.gnu.org ([140.186.70.10]:51007) by monty-python.gnu.org with esmtp (Exim 4.60) (envelope-from ) id 1OlT4W-0005Z6-4Q for bug-gsl@gnu.org; Tue, 17 Aug 2010 16:47:48 -0400 Received: from localhost ([127.0.0.1]:34164 helo=fencepost.gnu.org) by fencepost.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OlT4T-0006VJ-1U; Tue, 17 Aug 2010 16:47:46 -0400 Date: Tue, 17 Aug 2010 21:47:29 +0100 Message-ID: <43k4np6k66.wl%bjg@gnu.org> From: Brian Gough To: Thanh Vo Subject: Re: [Bug-gsl] Potential Bugs? In-Reply-To: <43mxsl6klg.wl%bjg@gnu.org> References: <43mxsl6klg.wl%bjg@gnu.org> User-Agent: Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) MIME-Version: 1.0 (generated by SEMI 1.14.6 - "Maruoka") Content-Type: text/plain; charset=US-ASCII X-detected-operating-system: by monty-python.gnu.org: GNU/Linux 2.6 (newer, 3) Cc: gsl-discuss@sourceware.org, bug-gsl@gnu.org, "Earl T. Barr" X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 17 Aug 2010 20:47:49 -0000 At Tue, 17 Aug 2010 21:38:19 +0100, Brian Gough wrote: > The overflow exception is a valid warning. It's also a bug in the > sense that it is avoidable (although the end result is correct, z=-1 > since 1/inf=0). It would be better written as (32/y)/y. Of course, apart from the extra factor of 2^5 that gives an underflow instead but it's probably preferable to an overflow. From MAILER-DAEMON Tue Aug 24 09:28:11 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OntXu-0007kz-CP for mharc-bug-gsl@gnu.org; Tue, 24 Aug 2010 09:28:10 -0400 Received: from [140.186.70.92] (port=53155 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OngEW-0007WO-Cp for bug-gsl@gnu.org; Mon, 23 Aug 2010 19:15:17 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1OngEV-00006V-DS for bug-gsl@gnu.org; Mon, 23 Aug 2010 19:15:16 -0400 Received: from mail-iw0-f169.google.com ([209.85.214.169]:54242) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OngEV-00006L-59 for bug-gsl@gnu.org; Mon, 23 Aug 2010 19:15:15 -0400 Received: by iwn33 with SMTP id 33so5353059iwn.0 for ; Mon, 23 Aug 2010 16:15:14 -0700 (PDT) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:received:date:message-id :subject:from:to:content-type; bh=tWjhHqi2Lek8w16JC3R+ogYYyQ7DkcZIXxvOJAu0bwQ=; b=q4S0p+uz4r2hhwupr5TwTFggZ2jATlzXAcuXM70D3izF1kv/YTWGJBqwC/ZGMXTZOd PKMippln4eUn7laaJ+DkyPN142jbrFAfAnilO2Er5AusqqMbl6A3w1bd2jKwCNtowNhM 9DWcqbjsJp7yUPddI7vLTE4R4bWf03KJCt14M= DomainKey-Signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:date:message-id:subject:from:to:content-type; b=jfPVEZtEfQeIfQWrnTete8yVlfAzOZqJvaRqi0/xv8deqkVvbaM5SwF+8I9coAVr4G Zumrgg+Bll7zE7L3aYCgpInaUONgXVMKPEqQGvRKWsxyDWK7dALGSpAWVH+w1vrrG+HH FMQIR+Ff5LIQvE2PtAzcMmBNwRquSXKp2JPNI= MIME-Version: 1.0 Received: by 10.231.177.25 with SMTP id bg25mr6412788ibb.154.1282605314074; Mon, 23 Aug 2010 16:15:14 -0700 (PDT) Received: by 10.231.182.18 with HTTP; Mon, 23 Aug 2010 16:15:13 -0700 (PDT) Date: Mon, 23 Aug 2010 17:15:13 -0600 Message-ID: From: shaohao chen To: bug-gsl@gnu.org X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) X-Mailman-Approved-At: Tue, 24 Aug 2010 09:28:09 -0400 Content-Type: text/plain; charset=ISO-8859-1 X-Content-Filtered-By: Mailman/MimeDel 2.1.5 Subject: [Bug-gsl] questions for Confluent Hypergeometric Function X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Mon, 23 Aug 2010 23:15:18 -0000 To gsl develop group: The gsl has offered a lot ot useful functions for scientific reseach. I usually calculate the special functions by gsl. But I have a problem now. I have to calculate the Confluent Hypergeometric Function with complex arguments in a C++ programme. I found the function "double gsl sf hyperg 1F1 (double a, double b, double x )" for double arguments, but is there any function for complex arguments? If not, could you please consider to develop such a function? The Confluent Hypergeometric Function with complex arguments is quite useful in many fields of physics. Thank you! Best, Shaohao From MAILER-DAEMON Tue Aug 24 23:59:41 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1Oo79J-0007Z2-JF for mharc-bug-gsl@gnu.org; Tue, 24 Aug 2010 23:59:41 -0400 Received: from [140.186.70.92] (port=39012 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1Oo79H-0007Xe-E1 for bug-gsl@gnu.org; Tue, 24 Aug 2010 23:59:40 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1Oo79G-0004Rj-96 for bug-gsl@gnu.org; Tue, 24 Aug 2010 23:59:39 -0400 Received: from outgoing-mail.its.caltech.edu ([131.215.239.19]:38736) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1Oo79G-0004Qv-44 for bug-gsl@gnu.org; Tue, 24 Aug 2010 23:59:38 -0400 Received: from fire-doxen.imss.caltech.edu (localhost [127.0.0.1]) by fire-doxen-postvirus (Postfix) with ESMTP id 3B58E32831E for ; Tue, 24 Aug 2010 20:59:34 -0700 (PDT) Received: from ajb-laptop-2009.localnet (adsl-69-230-183-128.dsl.irvnca.pacbell.net [69.230.183.128]) (Authenticated sender: abenson) by fire-doxen-ssl (Postfix) with ESMTP id 3F3703283D1 for ; Tue, 24 Aug 2010 20:59:32 -0700 (PDT) To: bug-gsl@gnu.org From: Andrew Benson Organization: California Institute of Technology Date: Tue, 24 Aug 2010 21:00:48 -0700 MIME-Version: 1.0 Content-Type: Multipart/Mixed; boundary="Boundary-00=_xVJdMyJ77+UlE93" Message-Id: <201008242100.49225.abenson@caltech.edu> X-detected-operating-system: by eggs.gnu.org: Genre and OS details not recognized. Subject: [Bug-gsl] Bug in gsl_sf_hyperg_2F1 X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 25 Aug 2010 03:59:40 -0000 --Boundary-00=_xVJdMyJ77+UlE93 Content-Type: text/plain; charset="us-ascii" Content-Transfer-Encoding: 7bit Under GSL 1.14 (compiled with either gcc version 4.4.4 or gcc version 4.1.2) the gsl_sf_hyperg_2F1 function crashes with the following arguments: gsl_sf_hyperg_2F1(-10.34,2.05,3.05,0.1725) The expected answer (computed using Maple for example) is 0.3104735522 I've attached an example case: $ g++ test.cpp -o test.x -lgsl -lgslcblas $ test.x gsl: hyperg_2F1.c:755: ERROR: error Default GSL error handler invoked. Abort The crash seems to occur only when the first argument is less than or equal to -10. For example, gsl_sf_hyperg_2F1(-10.0,2.05,3.05,0.1725) crashes, but gsl_sf_hyperg_2F1(-9.99999999999,2.05,3.05,0.1725) evaluates correctly to 0.321419346301974773. --Boundary-00=_xVJdMyJ77+UlE93 Content-Type: text/x-c++src; charset="utf-8"; name="test.cpp" Content-Transfer-Encoding: 7bit Content-Disposition: attachment; filename="test.cpp" #include #include #include #include using namespace std; #include int main (void) { double x = 0.1725; double a = -10.34; double b = 2.05; double c = 3.05; double y = gsl_sf_hyperg_2F1(a,b,c,x); cout << "= " << setprecision(18) << setw(20) << y << endl; return 0; } --Boundary-00=_xVJdMyJ77+UlE93-- From MAILER-DAEMON Wed Aug 25 00:01:17 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1Oo7Ar-0007nF-QM for mharc-bug-gsl@gnu.org; Wed, 25 Aug 2010 00:01:17 -0400 Received: from [140.186.70.92] (port=40215 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1Oo7Ao-0007mz-CS for bug-gsl@gnu.org; Wed, 25 Aug 2010 00:01:15 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1Oo7An-0004y6-68 for bug-gsl@gnu.org; Wed, 25 Aug 2010 00:01:14 -0400 Received: from outgoing-mail.its.caltech.edu ([131.215.239.19]:52268) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1Oo7Am-0004xz-Sb for bug-gsl@gnu.org; Wed, 25 Aug 2010 00:01:13 -0400 Received: from fire-doxen.imss.caltech.edu (localhost [127.0.0.1]) by fire-doxen-postvirus (Postfix) with ESMTP id 51B1D328306 for ; Tue, 24 Aug 2010 21:01:12 -0700 (PDT) Received: from ajb-laptop-2009.localnet (adsl-69-230-183-128.dsl.irvnca.pacbell.net [69.230.183.128]) (Authenticated sender: abenson) by fire-doxen-ssl (Postfix) with ESMTP id E2A3E32831E for ; Tue, 24 Aug 2010 21:01:09 -0700 (PDT) To: bug-gsl@gnu.org From: Andrew Benson Organization: California Institute of Technology Date: Tue, 24 Aug 2010 21:02:26 -0700 MIME-Version: 1.0 Content-Type: Multipart/Mixed; boundary="Boundary-00=_SXJdMAE3hXyD+js" Message-Id: <201008242102.27057.abenson@caltech.edu> X-detected-operating-system: by eggs.gnu.org: Genre and OS details not recognized. Subject: [Bug-gsl] Re: Bug in gsl_sf_hyperg_2F1 X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 25 Aug 2010 04:01:16 -0000 --Boundary-00=_SXJdMAE3hXyD+js Content-Type: text/plain; charset="iso-8859-15" Content-Transfer-Encoding: 7bit On Tuesday, August 24, 2010, Andrew Benson wrote: > Under GSL 1.14 (compiled with either gcc version 4.4.4 or gcc version > 4.1.2) the gsl_sf_hyperg_2F1 function crashes with the following > arguments: > > gsl_sf_hyperg_2F1(-10.34,2.05,3.05,0.1725) > > The expected answer (computed using Maple for example) is 0.3104735522 > > I've attached an example case: > > $ g++ test.cpp -o test.x -lgsl -lgslcblas > $ test.x > gsl: hyperg_2F1.c:755: ERROR: error > Default GSL error handler invoked. > Abort > > The crash seems to occur only when the first argument is less than or equal > to -10. For example, > > gsl_sf_hyperg_2F1(-10.0,2.05,3.05,0.1725) > > crashes, but > > gsl_sf_hyperg_2F1(-9.99999999999,2.05,3.05,0.1725) > > evaluates correctly to 0.321419346301974773. Some further investigation shows that this combination of parameters is falling through all of the various methods for computing the result in gsl_sf_hyperg_2F1_e and hitting the "/* We give up. */" block at the end. This seems to be related to http://savannah.gnu.org/bugs/?24812 - if I make the changes suggested there (patch attached) then my test code works successfully: $ test.x = 0.310473552213130111 I haven't checked that this doesn't break anything else (beyond running "make check" on the newly compiled code). --Boundary-00=_SXJdMAE3hXyD+js Content-Type: text/x-patch; charset="utf-8"; name="hyperg_2F1.c.patch" Content-Transfer-Encoding: 7bit Content-Disposition: attachment; filename="hyperg_2F1.c.patch" --- specfunc/hyperg_2F1.c 2010-03-10 02:57:14.000000000 -0800 +++ specfunc/hyperg_2F1_new.c 2010-08-24 13:56:20.000000000 -0700 @@ -734,14 +734,14 @@ gsl_sf_hyperg_2F1_e(double a, double b, return hyperg_2F1_luke(a, b, c, x, result); } - if(GSL_MAX_DBL(fabs(a),1.0)*fabs(bp)*fabs(x) < 2.0*fabs(c)) { + if(GSL_MAX_DBL(fabs(ap),1.0)*fabs(bp)*fabs(x) < 2.0*fabs(c)) { /* If c is large enough or x is small enough, * we can attempt the series anyway. */ return hyperg_2F1_series(a, b, c, x, result); } - if(fabs(bp*bp*x*x) < 0.001*fabs(bp) && fabs(a) < 10.0) { + if(fabs(bp*bp*x*x) < 0.001*fabs(bp) && fabs(ap) < 10.0) { /* The famous but nearly worthless "large b" asymptotic. */ int stat = gsl_sf_hyperg_1F1_e(a, c, bp*x, result); --Boundary-00=_SXJdMAE3hXyD+js-- From MAILER-DAEMON Wed Aug 25 08:09:25 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OoEnE-0000nH-S6 for mharc-bug-gsl@gnu.org; Wed, 25 Aug 2010 08:09:24 -0400 Received: from [140.186.70.92] (port=42365 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1Onziw-0008Bz-LP for bug-gsl@gnu.org; Tue, 24 Aug 2010 16:03:59 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1Onziv-0001Fh-89 for bug-gsl@gnu.org; Tue, 24 Aug 2010 16:03:58 -0400 Received: from outgoing-mail.its.caltech.edu ([131.215.239.19]:56081) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1Onziu-0001FO-Vl for bug-gsl@gnu.org; Tue, 24 Aug 2010 16:03:57 -0400 Received: from fire-doxen.imss.caltech.edu (localhost [127.0.0.1]) by fire-doxen-postvirus (Postfix) with ESMTP id 4FA593283AE for ; Tue, 24 Aug 2010 13:03:53 -0700 (PDT) Received: from ajb-laptop-2009.localnet (adsl-69-230-183-128.dsl.irvnca.pacbell.net [69.230.183.128]) (Authenticated sender: abenson) by fire-doxen-ssl (Postfix) with ESMTP id 41F1632834E for ; Tue, 24 Aug 2010 13:03:51 -0700 (PDT) From: Andrew Benson Organization: California Institute of Technology To: bug-gsl@gnu.org Date: Tue, 24 Aug 2010 13:05:07 -0700 User-Agent: KMail/1.13.5 (Linux/2.6.33.6-147.2.4.fc13.i686.PAE; KDE/4.4.5; i686; ; ) MIME-Version: 1.0 Content-Type: Multipart/Mixed; boundary="Boundary-00=_0XCdMoaEYKrU8gG" Message-Id: <201008241305.08275.abenson@its.caltech.edu> X-detected-operating-system: by eggs.gnu.org: Genre and OS details not recognized. X-Mailman-Approved-At: Wed, 25 Aug 2010 08:09:23 -0400 Subject: [Bug-gsl] Bug in gsl_sf_hyperg_2F1 X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 24 Aug 2010 20:03:59 -0000 --Boundary-00=_0XCdMoaEYKrU8gG Content-Type: Text/Plain; charset="us-ascii" Content-Transfer-Encoding: 7bit Under GSL 1.14 (compiled with either gcc version 4.4.4 or gcc version 4.1.2) the gsl_sf_hyperg_2F1 function crashes with the following arguments: gsl_sf_hyperg_2F1(-10.34,2.05,3.05,0.1725) The expected answer (computed using Maple for example) is 0.3104735522 I've attached an example case: $ g++ test.cpp -o test.x -lgsl -lgslcblas $ test.x gsl: hyperg_2F1.c:755: ERROR: error Default GSL error handler invoked. Abort The crash seems to occur only when the first argument is less than or equal to -10. For example, gsl_sf_hyperg_2F1(-10.0,2.05,3.05,0.1725) crashes, but gsl_sf_hyperg_2F1(-9.99999999999,2.05,3.05,0.1725) evaluates correctly to 0.321419346301974773. --Boundary-00=_0XCdMoaEYKrU8gG Content-Type: text/x-c++src; charset="UTF-8"; name="test.cpp" Content-Transfer-Encoding: 7bit Content-Disposition: attachment; filename="test.cpp" #include #include #include #include using namespace std; #include int main (void) { double x = 0.1725; double a = -10.34; double b = 2.05; double c = 3.05; double y = gsl_sf_hyperg_2F1(a,b,c,x); cout << "= " << setprecision(18) << setw(20) << y << endl; return 0; } --Boundary-00=_0XCdMoaEYKrU8gG-- From MAILER-DAEMON Wed Aug 25 08:09:26 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OoEnF-0000ns-QA for mharc-bug-gsl@gnu.org; Wed, 25 Aug 2010 08:09:25 -0400 Received: from [140.186.70.92] (port=40994 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1Oo1gN-0006yh-0V for bug-gsl@gnu.org; Tue, 24 Aug 2010 18:09:27 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1Oo1gL-00028j-H1 for bug-gsl@gnu.org; Tue, 24 Aug 2010 18:09:26 -0400 Received: from outgoing-mail.its.caltech.edu ([131.215.239.19]:56246) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1Oo1gL-00028e-8I for bug-gsl@gnu.org; Tue, 24 Aug 2010 18:09:25 -0400 Received: from fire-doxen.imss.caltech.edu (localhost [127.0.0.1]) by fire-doxen-postvirus (Postfix) with ESMTP id C8BFE328258 for ; Tue, 24 Aug 2010 15:09:24 -0700 (PDT) Received: from ajb-laptop-2009.localnet (adsl-69-230-183-128.dsl.irvnca.pacbell.net [69.230.183.128]) (Authenticated sender: abenson) by fire-doxen-ssl (Postfix) with ESMTP id 93D593282EB for ; Tue, 24 Aug 2010 15:09:22 -0700 (PDT) From: Andrew Benson Organization: California Institute of Technology To: bug-gsl@gnu.org Date: Tue, 24 Aug 2010 15:10:39 -0700 User-Agent: KMail/1.13.5 (Linux/2.6.33.6-147.2.4.fc13.i686.PAE; KDE/4.4.5; i686; ; ) References: <201008241305.08275.abenson@its.caltech.edu> In-Reply-To: <201008241305.08275.abenson@its.caltech.edu> MIME-Version: 1.0 Content-Type: Multipart/Mixed; boundary="Boundary-00=_fNEdM1kxuoQGUl0" Message-Id: <201008241510.39528.abenson@its.caltech.edu> X-detected-operating-system: by eggs.gnu.org: Genre and OS details not recognized. X-Mailman-Approved-At: Wed, 25 Aug 2010 08:09:23 -0400 Subject: [Bug-gsl] Re: Bug in gsl_sf_hyperg_2F1 X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 24 Aug 2010 22:09:28 -0000 --Boundary-00=_fNEdM1kxuoQGUl0 Content-Type: Text/Plain; charset="iso-8859-15" Content-Transfer-Encoding: 7bit On Tuesday, August 24, 2010, Andrew Benson wrote: > Under GSL 1.14 (compiled with either gcc version 4.4.4 or gcc version > 4.1.2) the gsl_sf_hyperg_2F1 function crashes with the following > arguments: > > gsl_sf_hyperg_2F1(-10.34,2.05,3.05,0.1725) > > The expected answer (computed using Maple for example) is 0.3104735522 > > I've attached an example case: > > $ g++ test.cpp -o test.x -lgsl -lgslcblas > $ test.x > gsl: hyperg_2F1.c:755: ERROR: error > Default GSL error handler invoked. > Abort > > The crash seems to occur only when the first argument is less than or equal > to -10. For example, > > gsl_sf_hyperg_2F1(-10.0,2.05,3.05,0.1725) > > crashes, but > > gsl_sf_hyperg_2F1(-9.99999999999,2.05,3.05,0.1725) > > evaluates correctly to 0.321419346301974773. Some further investigation shows that this combination of parameters is falling through all of the various methods for computing the result in gsl_sf_hyperg_2F1_e and hitting the "/* We give up. */" block at the end. This seems to be related to http://savannah.gnu.org/bugs/?24812 - if I make the changes suggested there (patch attached) then my test code works successfully: $ test.x = 0.310473552213130111 I haven't checked that this doesn't break anything else (beyond running "make check" on the newly compiled code). --Boundary-00=_fNEdM1kxuoQGUl0 Content-Type: text/x-patch; charset="UTF-8"; name="hyperg_2F1.c.patch" Content-Transfer-Encoding: 7bit Content-Disposition: attachment; filename="hyperg_2F1.c.patch" --- specfunc/hyperg_2F1.c 2010-03-10 02:57:14.000000000 -0800 +++ specfunc/hyperg_2F1_new.c 2010-08-24 13:56:20.000000000 -0700 @@ -734,14 +734,14 @@ gsl_sf_hyperg_2F1_e(double a, double b, return hyperg_2F1_luke(a, b, c, x, result); } - if(GSL_MAX_DBL(fabs(a),1.0)*fabs(bp)*fabs(x) < 2.0*fabs(c)) { + if(GSL_MAX_DBL(fabs(ap),1.0)*fabs(bp)*fabs(x) < 2.0*fabs(c)) { /* If c is large enough or x is small enough, * we can attempt the series anyway. */ return hyperg_2F1_series(a, b, c, x, result); } - if(fabs(bp*bp*x*x) < 0.001*fabs(bp) && fabs(a) < 10.0) { + if(fabs(bp*bp*x*x) < 0.001*fabs(bp) && fabs(ap) < 10.0) { /* The famous but nearly worthless "large b" asymptotic. */ int stat = gsl_sf_hyperg_1F1_e(a, c, bp*x, result); --Boundary-00=_fNEdM1kxuoQGUl0-- From MAILER-DAEMON Wed Aug 25 08:09:25 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OoEnF-0000nW-AE for mharc-bug-gsl@gnu.org; Wed, 25 Aug 2010 08:09:25 -0400 Received: from [140.186.70.92] (port=40914 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1Oo1gB-0006x9-NM for bug-gsl@gnu.org; Tue, 24 Aug 2010 18:09:17 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1Oo1gA-00026A-PF for bug-gsl@gnu.org; Tue, 24 Aug 2010 18:09:15 -0400 Received: from outgoing-mail.its.caltech.edu ([131.215.239.19]:43839) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1Oo1gA-000264-JL for bug-gsl@gnu.org; Tue, 24 Aug 2010 18:09:14 -0400 Received: from fire-doxen.imss.caltech.edu (localhost [127.0.0.1]) by fire-doxen-postvirus (Postfix) with ESMTP id 31D4A328333 for ; Tue, 24 Aug 2010 15:09:13 -0700 (PDT) Received: from ajb-laptop-2009.localnet (adsl-69-230-183-128.dsl.irvnca.pacbell.net [69.230.183.128]) (Authenticated sender: abenson) by fire-doxen-ssl (Postfix) with ESMTP id 1C2433282BF for ; Tue, 24 Aug 2010 15:09:11 -0700 (PDT) To: bug-gsl@gnu.org From: Andrew Benson Organization: California Institute of Technology Date: Tue, 24 Aug 2010 15:10:26 -0700 MIME-Version: 1.0 Content-Type: Multipart/Mixed; boundary="Boundary-00=_TNEdMPQGvsBrbRk" Message-Id: <201008241510.27901.abenson@its.caltech.edu> X-detected-operating-system: by eggs.gnu.org: Genre and OS details not recognized. X-Mailman-Approved-At: Wed, 25 Aug 2010 08:09:23 -0400 Subject: [Bug-gsl] Bug in gsl_sf_hyperg_2F1 X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 24 Aug 2010 22:09:17 -0000 --Boundary-00=_TNEdMPQGvsBrbRk Content-Type: text/plain; charset="us-ascii" Content-Transfer-Encoding: 7bit Under GSL 1.14 (compiled with either gcc version 4.4.4 or gcc version 4.1.2) the gsl_sf_hyperg_2F1 function crashes with the following arguments: gsl_sf_hyperg_2F1(-10.34,2.05,3.05,0.1725) The expected answer (computed using Maple for example) is 0.3104735522 I've attached an example case: $ g++ test.cpp -o test.x -lgsl -lgslcblas $ test.x gsl: hyperg_2F1.c:755: ERROR: error Default GSL error handler invoked. Abort The crash seems to occur only when the first argument is less than or equal to -10. For example, gsl_sf_hyperg_2F1(-10.0,2.05,3.05,0.1725) crashes, but gsl_sf_hyperg_2F1(-9.99999999999,2.05,3.05,0.1725) evaluates correctly to 0.321419346301974773. --Boundary-00=_TNEdMPQGvsBrbRk Content-Type: text/x-c++src; charset="utf-8"; name="test.cpp" Content-Transfer-Encoding: 7bit Content-Disposition: attachment; filename="test.cpp" #include #include #include #include using namespace std; #include int main (void) { double x = 0.1725; double a = -10.34; double b = 2.05; double c = 3.05; double y = gsl_sf_hyperg_2F1(a,b,c,x); cout << "= " << setprecision(18) << setw(20) << y << endl; return 0; } --Boundary-00=_TNEdMPQGvsBrbRk-- From MAILER-DAEMON Wed Aug 25 23:39:46 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OoTJa-0004Zg-RE for mharc-bug-gsl@gnu.org; Wed, 25 Aug 2010 23:39:46 -0400 Received: from [140.186.70.92] (port=42430 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OoTJX-0004Yi-QH for bug-gsl@gnu.org; Wed, 25 Aug 2010 23:39:44 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1OoTJW-0004Av-Gs for bug-gsl@gnu.org; Wed, 25 Aug 2010 23:39:43 -0400 Received: from mail-iw0-f169.google.com ([209.85.214.169]:52651) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OoTJW-00049j-C6 for bug-gsl@gnu.org; Wed, 25 Aug 2010 23:39:42 -0400 Received: by iwn33 with SMTP id 33so1738552iwn.0 for ; Wed, 25 Aug 2010 20:39:41 -0700 (PDT) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from :user-agent:mime-version:to:subject:content-type; bh=cJwAP685lgW4WR/1XwaRWsqbbbj/ySIEHEU2/9++YuE=; b=E0Cuf+ygYRr/m8xVddkbv9cZvbH8L3iHgaLoHYtSgWAn0X0juno8WNnkn+/EM6u16k UmZ4ubXFKmLT7PN2jfJwaFw9vxbNTu80j4WNzgWkWliHcPBPImVuoId2GIElNmuSAn+r HJj+QDKiCs2lYkyUtr2zYStNiQT8QZLaXEL3M= DomainKey-Signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:user-agent:mime-version:to:subject :content-type; b=MCpgB+4Z7bKKzmFizVNkdTPrQnGkcB70YmYiQ7FKUvqR0ZamtH+4lZIGQZoDu9s34L 3VZRLKE5IfSoGFW9kDAh95ZcO53/LgVf5RFwJaALp6Gw2Rr5/W2fQ3At2sJwR+1AWZxb lwrw0SnCo7RoLR3+35jNRyl3Md+gRUDYOSEPY= Received: by 10.231.31.199 with SMTP id z7mr10995618ibc.50.1282793981146; Wed, 25 Aug 2010 20:39:41 -0700 (PDT) Received: from [192.168.0.101] (CPE0011950b723b-CM001225409cb6.cpe.net.cable.rogers.com [173.34.169.193]) by mx.google.com with ESMTPS id h8sm2004118ibk.3.2010.08.25.20.39.39 (version=TLSv1/SSLv3 cipher=RC4-MD5); Wed, 25 Aug 2010 20:39:40 -0700 (PDT) Message-ID: <4C75E25A.7050307@gmail.com> Date: Wed, 25 Aug 2010 23:41:14 -0400 From: "Alexey A. Illarionov" User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.9.2.4) Gecko/20100624 Lanikai/3.1 MIME-Version: 1.0 To: bug-gsl@gnu.org Content-Type: multipart/mixed; boundary="------------000907090308050908060603" X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) Subject: [Bug-gsl] bug in gsl_sf_coulomb_wave_FG_e (gsl 1.14 compiled with gcc 4.5.0) X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Thu, 26 Aug 2010 03:39:45 -0000 This is a multi-part message in MIME format. --------------000907090308050908060603 Content-Type: text/plain; charset=UTF-8; format=flowed Content-Transfer-Encoding: 7bit Hi guys, It seems that I find a bug in gsl_sf_coulomb_wave_FG_e. For large values of lambda and small values of x, with eta = 0 it produces nan values without raising of the flag. Here the sample cases (see attachment with example) gsl_sf_coulomb_wave_FG_e(0.,1.2693881947287221e-07, 37, 1, ...) gsl_sf_coulomb_wave_FG_e(0.,5.911135077240691e-07, 39, 1, ...) gsl_sf_coulomb_wave_FG_e(0.,6.118185507743884e-07, 40, 1, ...) and lots of others. It seems that the origin of the problem is at coulomb_F_recur function (call at coulomb.c:970), that return garbage-like big numbers but the error flag is not raised. Here since we are interested only in unnormalized values we can makes a workaround like checking if values greater than some number, renormalize them to 1 for some value. But if this function is used elsewhere with the need of normalized values, this workaround will produce problems there. Unfortunately I don't have time right now to dig in more, I'll try look at it at the weekend if the issue is not solved by that date. P.S. I'm not sure, but for me at this values of arguments (x << x_{turning_point}) using series approximation directly is more adequate (if we do not consider G and Gp which I'm not sure how to deal). -- Best regards, Alexey Alexandrovich Illarionov --------------000907090308050908060603 Content-Type: text/x-csrc; name="2.c" Content-Transfer-Encoding: 7bit Content-Disposition: attachment; filename="2.c" #include #include int main() { gsl_sf_result f, fp, g, gp; double exp_f, exp_g; int i; i = gsl_sf_coulomb_wave_FG_e(0.,1.2693881947287221e-07, 37, 1, &f,&fp, & g,&gp, &exp_f, &exp_g); printf(" f=(%g, %g)\n" " fp=(%g, %g)\n" " g=(%g, %g)\n" " gp=(%g, %g)\n" " exp_f=%g exp_g=%g\n", f.val, f.err, fp.val, fp.err, g.val, g.err, gp.val, gp.err, exp_f, exp_g ); if (i) printf("%i\n", i); return i; } --------------000907090308050908060603-- From MAILER-DAEMON Thu Aug 26 08:05:35 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OobD5-00010b-BC for mharc-bug-gsl@gnu.org; Thu, 26 Aug 2010 08:05:35 -0400 Received: from [140.186.70.92] (port=50423 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OobD1-0000zH-Td for bug-gsl@gnu.org; Thu, 26 Aug 2010 08:05:32 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1OobD0-000636-MF for bug-gsl@gnu.org; Thu, 26 Aug 2010 08:05:31 -0400 Received: from colonialone.fsf.org ([140.186.70.51]:50986 helo=internal.in.savannah.gnu.org) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OobD0-000630-Kl for bug-gsl@gnu.org; Thu, 26 Aug 2010 08:05:30 -0400 Received: from [10.1.0.103] (helo=frontend.in.savannah.gnu.org) by internal.in.savannah.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OobCy-0007Wg-GN; Thu, 26 Aug 2010 12:05:28 +0000 Received: from www-data by frontend.in.savannah.gnu.org with local (Exim 4.69) (envelope-from ) id 1OobCy-0000AM-E1; Thu, 26 Aug 2010 12:05:28 +0000 To: Brian Gough , Sergey B Kirpichev , bug-gsl@gnu.org From: Sergey B Kirpichev X-Savane-Server: savannah.gnu.org:443 [10.1.0.103] MIME-Version: 1.0 Content-Type: text/plain;charset=UTF-8 X-Savane-Project: gsl X-Savane-Tracker: bugs X-Savane-Item-ID: 30540 User-Agent: Mozilla/5.0 (X11; U; Linux x86_64; en; rv:1.9.0.19) Gecko/20080528 Epiphany/2.22 X-Apparently-From: 91.192.71.112 (Savane authenticated user skirpichev) Message-Id: <20100826-120528.sv79421.31356@savannah.gnu.org> References: <20100724-121507.sv79421.66062@savannah.gnu.org> <20100806-140911.sv1824.76026@savannah.gnu.org> In-Reply-To: <20100806-140911.sv1824.76026@savannah.gnu.org> Date: Thu, 26 Aug 2010 12:05:28 +0000 X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) Cc: Subject: [Bug-gsl] [bug #30540] please, add convergence checks in rk4imp/rk2imp X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Thu, 26 Aug 2010 12:05:33 -0000 Follow-up Comment #2, bug #30540 (project gsl): I can't figure out how to use this code for fixed step size. imprk2/imprk4 seems to be not working at all for this case: $ ./bug Segmentation fault Other steppers works (e.g. rk2) just fine. See bug-ode2.c. Second, Newton-type iterations could be inefficient in general. See: Ernst Hairer, Christian Lubich, Gerhard Wanner, Geometric Numerical Integration. Structure-Preserving Algorithms for Ordinary Differential Equations. Springer Series in Comput. Mathematics, Vol. 31, Springer-Verlag 2002. pp. 330-332. (file #21314) _______________________________________________________ Additional Item Attachment: File name: bug-ode2.c Size:2 KB _______________________________________________________ Reply to this item at: _______________________________________________ Message sent via/by Savannah http://savannah.gnu.org/ From MAILER-DAEMON Thu Aug 26 14:27:17 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OohAT-0006BO-68 for mharc-bug-gsl@gnu.org; Thu, 26 Aug 2010 14:27:17 -0400 Received: from [199.232.76.173] (port=37959 helo=monty-python.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OohAR-0006BA-9I for bug-gsl@gnu.org; Thu, 26 Aug 2010 14:27:15 -0400 Received: from Debian-exim by monty-python.gnu.org with spam-scanned (Exim 4.60) (envelope-from ) id 1OohAQ-0004U1-8i for bug-gsl@gnu.org; Thu, 26 Aug 2010 14:27:15 -0400 Received: from fencepost.gnu.org ([140.186.70.10]:51245) by monty-python.gnu.org with esmtp (Exim 4.60) (envelope-from ) id 1OohAO-0004Tw-JS for bug-gsl@gnu.org; Thu, 26 Aug 2010 14:27:14 -0400 Received: from localhost ([127.0.0.1]:34977 helo=fencepost.gnu.org) by fencepost.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1Ooh9s-0005kk-Ok; Thu, 26 Aug 2010 14:26:40 -0400 Date: Thu, 26 Aug 2010 19:26:32 +0100 Message-ID: <43zkw9nsbr.wl%bjg@gnu.org> From: Brian Gough To: shaohao chen Subject: Re: [Bug-gsl] questions for Confluent Hypergeometric Function In-Reply-To: References: User-Agent: Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) MIME-Version: 1.0 (generated by SEMI 1.14.6 - "Maruoka") Content-Type: text/plain; charset=US-ASCII X-detected-operating-system: by monty-python.gnu.org: GNU/Linux 2.6 (newer, 3) Cc: bug-gsl@gnu.org X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Thu, 26 Aug 2010 18:27:15 -0000 At Mon, 23 Aug 2010 17:15:13 -0600, shaohao chen wrote: > I have to calculate the Confluent Hypergeometric Function with > complex arguments in a C++ programme. I found the function "double > gsl sf hyperg 1F1 (double a, double b, double x )" for double > arguments, but is there any function for complex arguments? Sorry, we don't have that :-( > If not, could you please consider to develop such a function? The > Confluent Hypergeometric Function with complex arguments is quite > useful in many fields of physics. Thank you! I will note it in the TODO file. I am not able to work on it myself, but if you can work on it I encourage you to do that. -- Brian Gough From MAILER-DAEMON Fri Aug 27 04:42:23 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OouVy-0008OY-Sb for mharc-bug-gsl@gnu.org; Fri, 27 Aug 2010 04:42:22 -0400 Received: from [140.186.70.92] (port=54868 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OouVv-0008Ny-Gy for bug-gsl@gnu.org; Fri, 27 Aug 2010 04:42:20 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1OouVu-0007bW-8Q for bug-gsl@gnu.org; Fri, 27 Aug 2010 04:42:19 -0400 Received: from colonialone.fsf.org ([140.186.70.51]:60383 helo=internal.in.savannah.gnu.org) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OouVu-0007bQ-2O for bug-gsl@gnu.org; Fri, 27 Aug 2010 04:42:18 -0400 Received: from [10.1.0.103] (helo=frontend.in.savannah.gnu.org) by internal.in.savannah.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OouVt-0005og-DX for bug-gsl@gnu.org; Fri, 27 Aug 2010 08:42:17 +0000 Received: from www-data by frontend.in.savannah.gnu.org with local (Exim 4.69) (envelope-from ) id 1OouVt-0004ZK-9h for bug-gsl@gnu.org; Fri, 27 Aug 2010 08:42:17 +0000 To: bug-gsl@gnu.org From: Brian Gough X-Savane-Server: savannah.gnu.org:443 [10.1.0.103] MIME-Version: 1.0 Content-Type: text/plain;charset=UTF-8 X-Savane-Project: gsl X-Savane-Tracker: bugs X-Savane-Item-ID: 24812 User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-GB; rv:1.9.0.11) Gecko/2009061118 Fedora/3.0.11-1.fc9 Firefox/3.0.11 X-Apparently-From: 94.193.216.241 (Savane authenticated user bjg) Message-Id: <20100827-094216.sv1824.61592@savannah.gnu.org> References: <20081111-144308.sv1824.49604@savannah.gnu.org> <20081111-144501.sv1824.26308@savannah.gnu.org> In-Reply-To: <20081111-144501.sv1824.26308@savannah.gnu.org> Date: Fri, 27 Aug 2010 08:42:17 +0000 X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) Subject: [Bug-gsl] [bug #24812] gsl_sf_hyperg_2F1(11, -1 ; 11/2; 0.125) fails X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Fri, 27 Aug 2010 08:42:21 -0000 Follow-up Comment #2, bug #24812 (project gsl): From: Andrew Benson To: bug-gsl@gnu.org Subject: [Bug-gsl] Bug in gsl_sf_hyperg_2F1 Date: Tue, 24 Aug 2010 21:00:48 -0700 [1 ] Under GSL 1.14 (compiled with either gcc version 4.4.4 or gcc version 4.1.2) the gsl_sf_hyperg_2F1 function crashes with the following arguments: gsl_sf_hyperg_2F1(-10.34,2.05,3.05,0.1725) The expected answer (computed using Maple for example) is 0.3104735522 I've attached an example case: $ g++ test.cpp -o test.x -lgsl -lgslcblas $ test.x gsl: hyperg_2F1.c:755: ERROR: error Default GSL error handler invoked. Abort The crash seems to occur only when the first argument is less than or equal to -10. For example, gsl_sf_hyperg_2F1(-10.0,2.05,3.05,0.1725) crashes, but gsl_sf_hyperg_2F1(-9.99999999999,2.05,3.05,0.1725) evaluates correctly to 0.321419346301974773. _______________________________________________________ Reply to this item at: _______________________________________________ Message sent via/by Savannah http://savannah.gnu.org/ From MAILER-DAEMON Fri Aug 27 05:09:35 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OouwJ-0005Rg-5S for mharc-bug-gsl@gnu.org; Fri, 27 Aug 2010 05:09:35 -0400 Received: from [140.186.70.92] (port=47852 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OouwG-0005Qa-FF for bug-gsl@gnu.org; Fri, 27 Aug 2010 05:09:33 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1OouwF-0002Ny-Ak for bug-gsl@gnu.org; Fri, 27 Aug 2010 05:09:32 -0400 Received: from colonialone.fsf.org ([140.186.70.51]:44241 helo=internal.in.savannah.gnu.org) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OouwF-0002Nu-87 for bug-gsl@gnu.org; Fri, 27 Aug 2010 05:09:31 -0400 Received: from [10.1.0.103] (helo=frontend.in.savannah.gnu.org) by internal.in.savannah.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OouwE-00068B-Cj; Fri, 27 Aug 2010 09:09:30 +0000 Received: from www-data by frontend.in.savannah.gnu.org with local (Exim 4.69) (envelope-from ) id 1OouwE-0004ql-9s; Fri, 27 Aug 2010 09:09:30 +0000 To: pinchon.didier@orange.fr, bug-gsl@gnu.org From: Brian Gough X-Savane-Server: savannah.gnu.org:443 [10.1.0.103] MIME-Version: 1.0 Content-Type: text/plain;charset=UTF-8 X-Savane-Project: gsl X-Savane-Tracker: bugs X-Savane-Item-ID: 24812 User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-GB; rv:1.9.0.11) Gecko/2009061118 Fedora/3.0.11-1.fc9 Firefox/3.0.11 X-Apparently-From: 94.193.216.241 (Savane authenticated user bjg) Message-Id: <20100827-100929.sv1824.94467@savannah.gnu.org> References: <20081111-144308.sv1824.49604@savannah.gnu.org> <20081111-144501.sv1824.26308@savannah.gnu.org> <20100827-094216.sv1824.61592@savannah.gnu.org> In-Reply-To: <20100827-094216.sv1824.61592@savannah.gnu.org> Date: Fri, 27 Aug 2010 09:09:30 +0000 X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) Cc: Subject: [Bug-gsl] [bug #24812] gsl_sf_hyperg_2F1(11, -1 ; 11/2; 0.125) fails X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Fri, 27 Aug 2010 09:09:33 -0000 Update of bug #24812 (project gsl): Status: Confirmed => Fixed Open/Closed: Open => Closed _______________________________________________________ Follow-up Comment #3: fixed in revision 4676 _______________________________________________________ Reply to this item at: _______________________________________________ Message sent via/by Savannah http://savannah.gnu.org/ From MAILER-DAEMON Fri Aug 27 05:14:55 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1Oov1S-0006R3-TS for mharc-bug-gsl@gnu.org; Fri, 27 Aug 2010 05:14:54 -0400 Received: from [199.232.76.173] (port=40877 helo=monty-python.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1Oov1Q-0006Ph-Ac for bug-gsl@gnu.org; Fri, 27 Aug 2010 05:14:52 -0400 Received: from Debian-exim by monty-python.gnu.org with spam-scanned (Exim 4.60) (envelope-from ) id 1Oov1P-0000Do-3O for bug-gsl@gnu.org; Fri, 27 Aug 2010 05:14:52 -0400 Received: from fencepost.gnu.org ([140.186.70.10]:42161) by monty-python.gnu.org with esmtp (Exim 4.60) (envelope-from ) id 1Oov1O-0000Di-UT for bug-gsl@gnu.org; Fri, 27 Aug 2010 05:14:51 -0400 Received: from localhost ([127.0.0.1]:53418 helo=fencepost.gnu.org) by fencepost.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1Oov1O-0001vd-EL; Fri, 27 Aug 2010 05:14:50 -0400 Date: Fri, 27 Aug 2010 10:12:21 +0100 Message-ID: <43wrrco1vu.wl%bjg@gnu.org> From: Brian Gough To: Andrew Benson Subject: Re: [Bug-gsl] Re: Bug in gsl_sf_hyperg_2F1 In-Reply-To: <201008242102.27057.abenson@caltech.edu> References: <201008242102.27057.abenson@caltech.edu> User-Agent: Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) MIME-Version: 1.0 (generated by SEMI 1.14.6 - "Maruoka") Content-Type: text/plain; charset=US-ASCII X-detected-operating-system: by monty-python.gnu.org: GNU/Linux 2.6 (newer, 3) Cc: bug-gsl@gnu.org X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Fri, 27 Aug 2010 09:14:52 -0000 At Tue, 24 Aug 2010 21:02:26 -0700, Andrew Benson wrote: > Some further investigation shows that this combination of parameters is > falling through all of the various methods for computing the result in > gsl_sf_hyperg_2F1_e and hitting the "/* We give up. */" block at the end. > > This seems to be related to http://savannah.gnu.org/bugs/?24812 - if I make > the changes suggested there (patch attached) then my test code works > successfully: > > $ test.x > = 0.310473552213130111 > > I haven't checked that this doesn't break anything else (beyond running "make > check" on the newly compiled code). Thanks for the bug report. I've added a test case for it and applied the earlier patch (with a further correction because one of the a->ap changes in the function call arguments was also missing). It should work now from revision 4676 in the repository http://bzr.savannah.gnu.org/r/gsl/trunk. -- Brian Gough From MAILER-DAEMON Fri Aug 27 06:49:20 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OowUq-0005id-KQ for mharc-bug-gsl@gnu.org; Fri, 27 Aug 2010 06:49:20 -0400 Received: from [140.186.70.92] (port=52470 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OowUY-0005VV-JF for bug-gsl@gnu.org; Fri, 27 Aug 2010 06:49:19 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1OowU5-0006sh-3M for bug-gsl@gnu.org; Fri, 27 Aug 2010 06:48:33 -0400 Received: from colonialone.fsf.org ([140.186.70.51]:42321 helo=internal.in.savannah.gnu.org) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OowU5-0006sL-22 for bug-gsl@gnu.org; Fri, 27 Aug 2010 06:48:33 -0400 Received: from [10.1.0.103] (helo=frontend.in.savannah.gnu.org) by internal.in.savannah.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OowU3-0007Kv-UT for bug-gsl@gnu.org; Fri, 27 Aug 2010 10:48:32 +0000 Received: from www-data by frontend.in.savannah.gnu.org with local (Exim 4.69) (envelope-from ) id 1OowU3-0005wv-Ra for bug-gsl@gnu.org; Fri, 27 Aug 2010 10:48:31 +0000 To: bug-gsl@gnu.org From: Brian Gough X-Savane-Server: savannah.gnu.org:443 [10.1.0.103] MIME-Version: 1.0 Content-Type: text/plain;charset=UTF-8 X-Savane-Project: gsl X-Savane-Tracker: bugs X-Savane-Item-ID: 30885 User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-GB; rv:1.9.0.11) Gecko/2009061118 Fedora/3.0.11-1.fc9 Firefox/3.0.11 X-Apparently-From: 94.193.216.241 (Savane authenticated user bjg) Message-Id: <20100827-114831.sv1824.37784@savannah.gnu.org> References: In-Reply-To: Date: Fri, 27 Aug 2010 10:48:31 +0000 X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) Subject: [Bug-gsl] [bug #30885] nans from gsl_sf_coulomb_wave_FG_e(1.2693881947287221e-07, 0.0, lam_F=37, lam_G=36) X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Fri, 27 Aug 2010 10:49:19 -0000 URL: Summary: nans from gsl_sf_coulomb_wave_FG_e(1.2693881947287221e-07, 0.0, lam_F=37, lam_G=36) Project: GNU Scientific Library Submitted by: bjg Submitted on: Fri 27 Aug 2010 11:48:31 AM BST Category: Runtime error Severity: 3 - Normal Operating System: Status: Confirmed Assigned to: None Open/Closed: Open Release: 1.14 Discussion Lock: Any _______________________________________________________ Details: From: "Alexey A. Illarionov" To: bug-gsl@gnu.org Subject: [Bug-gsl] bug in gsl_sf_coulomb_wave_FG_e (gsl 1.14 compiled with \ gcc 4.5.0) Date: Wed, 25 Aug 2010 23:41:14 -0400 [1 ] Hi guys, It seems that I find a bug in gsl_sf_coulomb_wave_FG_e. For large values of lambda and small values of x, with eta = 0 it produces nan values without raising of the flag. Here the sample cases (see attachment with example) gsl_sf_coulomb_wave_FG_e(0.,1.2693881947287221e-07, 37, 1, ...) gsl_sf_coulomb_wave_FG_e(0.,5.911135077240691e-07, 39, 1, ...) gsl_sf_coulomb_wave_FG_e(0.,6.118185507743884e-07, 40, 1, ...) and lots of others. It seems that the origin of the problem is at coulomb_F_recur function (call at coulomb.c:970), that return garbage-like big numbers but the error flag is not raised. Here since we are interested only in unnormalized values we can makes a workaround like checking if values greater than some number, renormalize them to 1 for some value. But if this function is used elsewhere with the need of normalized values, this workaround will produce problems there. Unfortunately I don't have time right now to dig in more, I'll try look at it at the weekend if the issue is not solved by that date. P.S. I'm not sure, but for me at this values of arguments (x << x_{turning_point}) using series approximation directly is more adequate (if we do not consider G and Gp which I'm not sure how to deal). _______________________________________________________ Reply to this item at: _______________________________________________ Message sent via/by Savannah http://savannah.gnu.org/ From MAILER-DAEMON Fri Aug 27 09:32:02 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1Ooz2H-0007Bk-Tx for mharc-bug-gsl@gnu.org; Fri, 27 Aug 2010 09:32:01 -0400 Received: from [199.232.76.173] (port=49322 helo=monty-python.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1Ooz2E-0007Am-5w for bug-gsl@gnu.org; Fri, 27 Aug 2010 09:31:58 -0400 Received: from Debian-exim by monty-python.gnu.org with spam-scanned (Exim 4.60) (envelope-from ) id 1OoyzC-0001xQ-PB for bug-gsl@gnu.org; Fri, 27 Aug 2010 09:28:51 -0400 Received: from fencepost.gnu.org ([140.186.70.10]:38332) by monty-python.gnu.org with esmtp (Exim 4.60) (envelope-from ) id 1OoyzC-0001xL-Ie for bug-gsl@gnu.org; Fri, 27 Aug 2010 09:28:50 -0400 Received: from localhost ([127.0.0.1]:48302 helo=fencepost.gnu.org) by fencepost.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OoyzB-0005Tf-TW; Fri, 27 Aug 2010 09:28:49 -0400 Date: Fri, 27 Aug 2010 11:47:04 +0100 Message-ID: <43tymgnxhz.wl%bjg@gnu.org> From: Brian Gough To: "Alexey A. Illarionov" Subject: Re: [Bug-gsl] bug in gsl_sf_coulomb_wave_FG_e (gsl 1.14 compiled with gcc 4.5.0) In-Reply-To: <4C75E25A.7050307@gmail.com> References: <4C75E25A.7050307@gmail.com> User-Agent: Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) MIME-Version: 1.0 (generated by SEMI 1.14.6 - "Maruoka") Content-Type: text/plain; charset=US-ASCII X-detected-operating-system: by monty-python.gnu.org: GNU/Linux 2.6 (newer, 3) Cc: bug-gsl@gnu.org X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Fri, 27 Aug 2010 13:31:58 -0000 At Wed, 25 Aug 2010 23:41:14 -0400, Alexey A. Illarionov wrote: > It seems that I find a bug in gsl_sf_coulomb_wave_FG_e. For large values > of lambda and small values of x, with eta = 0 it produces nan values > without raising of the flag. > > Here the sample cases (see attachment with example) > gsl_sf_coulomb_wave_FG_e(0.,1.2693881947287221e-07, 37, 1, ...) > gsl_sf_coulomb_wave_FG_e(0.,5.911135077240691e-07, 39, 1, ...) > gsl_sf_coulomb_wave_FG_e(0.,6.118185507743884e-07, 40, 1, ...) > and lots of others. Thanks for the bug report. I'm adding a test case for this. Can you confirm if the following values are correct (I'm not working with these functions often): lam_F = 37.0; eta = 0.0; x = 1.2693881947287221e-07; k_G = 1; gsl_sf_coulomb_wave_FG_e(eta, x, lam_F, k_G, &F, &Fp, &G, &Gp, &Fe, &Ge); s = 0; message_buff[0] = 0; s += test_sf_check_result(message_buff, F, -2.020327525317343313380459069e299, TEST_TOL3); s += test_sf_check_result(message_buff, Fp, 5.729672862364037942289798904e307, TEST_TOL3); s += test_sf_check_result(message_buff, G, 4.397355593129492554472984282e299, TEST_TOL3); s += test_sf_check_result(message_buff, Gp, - 1.247095270068213211088454516e308, TEST_TOL3); I've computed them using Pari as c(L,n) = { 2^L * exp(-Pi*n/2) * abs(gamma(L+1+I*n)) / gamma(2*L+2)} cofac(L,n,r) = { I * exp(-I*r)*r^(-L)/(gamma(2*L+2)*c(L,n)) } FplusIG(L,n,r) = { cofac(L,n,r) * intnum(t=0,[[1],1],exp(-t)*t^(L-I*n)*(t+2*I*r)^(L+I*n)) } FplusIG(37,0,1.2693881947287221e-07+x) # for F and F' (real part) FplusIG(36,0,1.2693881947287221e-07+x) # for G and G' (imaginary part) -- Brian Gough From MAILER-DAEMON Fri Aug 27 13:34:54 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1Op2pK-0000ud-Dp for mharc-bug-gsl@gnu.org; Fri, 27 Aug 2010 13:34:54 -0400 Received: from [140.186.70.92] (port=47731 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1Op2pC-0000uP-V2 for bug-gsl@gnu.org; Fri, 27 Aug 2010 13:34:51 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1Op2p7-0001vt-Vl for bug-gsl@gnu.org; Fri, 27 Aug 2010 13:34:46 -0400 Received: from mail-gw0-f41.google.com ([74.125.83.41]:48774) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1Op2p1-0001uW-O3; Fri, 27 Aug 2010 13:34:35 -0400 Received: by gwj16 with SMTP id 16so1687023gwj.0 for ; Fri, 27 Aug 2010 10:34:34 -0700 (PDT) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from :user-agent:mime-version:to:cc:subject:references:in-reply-to :content-type:content-transfer-encoding; bh=B5Q4aS9ZysxQz8i3ogIJz+BeEXbB00HYZNVf8EbkgYE=; b=WlGn34Zv6I/axRZofuGPAQsV37+icFlAu3RmEYJ8E8gfHydJkjjCuTaRQfUJN60OVt hFJeyw1bI5ZEUAAIX98/jj7LCH7mAt3MXXcbzYwFUODvLvUNdOwqHBTLpSK8UgObeo4J k0zXK9tWp1I5bgdF9c5vGXMsicZUsAkkbX2QQ= DomainKey-Signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:user-agent:mime-version:to:cc:subject :references:in-reply-to:content-type:content-transfer-encoding; b=sPMfkaDxy3B9jfPgbiDK0AJAalBZoR8RV/XyA8qaKRcxce4jKYdnD92EzGZhwiPZEb JbMi94YIrT/Xvw+dfg6fSwV0c8Wj16BgVgdQnOQK/LPoT8gF6jAvQFZUrz3EOWvKL8sG S8YlMWMosCU7n6gfsfnNoOiLDadQ26814Hpxs= Received: by 10.151.92.4 with SMTP id u4mr1898941ybl.395.1282930474265; Fri, 27 Aug 2010 10:34:34 -0700 (PDT) Received: from [192.168.0.101] (CPE0011950b723b-CM001225409cb6.cpe.net.cable.rogers.com [173.34.169.193]) by mx.google.com with ESMTPS id t20sm9046500ybm.5.2010.08.27.10.34.32 (version=TLSv1/SSLv3 cipher=RC4-MD5); Fri, 27 Aug 2010 10:34:32 -0700 (PDT) Message-ID: <4C77F788.2090102@gmail.com> Date: Fri, 27 Aug 2010 13:36:08 -0400 From: "Alexey A. Illarionov" User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.9.2.4) Gecko/20100624 Lanikai/3.1 MIME-Version: 1.0 To: Brian Gough Subject: Re: [Bug-gsl] bug in gsl_sf_coulomb_wave_FG_e (gsl 1.14 compiled with gcc 4.5.0) References: <4C75E25A.7050307@gmail.com> <43tymgnxhz.wl%bjg@gnu.org> In-Reply-To: <43tymgnxhz.wl%bjg@gnu.org> Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) Cc: bug-gsl@gnu.org X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Fri, 27 Aug 2010 17:34:52 -0000 On 08/27/2010 06:47 AM, Brian Gough wrote: > At Wed, 25 Aug 2010 23:41:14 -0400, > Alexey A. Illarionov wrote: >> It seems that I find a bug in gsl_sf_coulomb_wave_FG_e. For large values >> of lambda and small values of x, with eta = 0 it produces nan values >> without raising of the flag. >> >> Here the sample cases (see attachment with example) >> gsl_sf_coulomb_wave_FG_e(0.,1.2693881947287221e-07, 37, 1, ...) >> gsl_sf_coulomb_wave_FG_e(0.,5.911135077240691e-07, 39, 1, ...) >> gsl_sf_coulomb_wave_FG_e(0.,6.118185507743884e-07, 40, 1, ...) >> and lots of others. > > Thanks for the bug report. I'm adding a test case for this. Can you > confirm if the following values are correct (I'm not working with > these functions often): > > lam_F = 37.0; > eta = 0.0; > x = 1.2693881947287221e-07; > k_G = 1; > gsl_sf_coulomb_wave_FG_e(eta, x, lam_F, k_G,&F,&Fp,&G,&Gp,&Fe,&Ge); > s = 0; > message_buff[0] = 0; > s += test_sf_check_result(message_buff, F, -2.020327525317343313380459069e299, TEST_TOL3); > s += test_sf_check_result(message_buff, Fp, 5.729672862364037942289798904e307, TEST_TOL3); > s += test_sf_check_result(message_buff, G, 4.397355593129492554472984282e299, TEST_TOL3); > s += test_sf_check_result(message_buff, Gp, - 1.247095270068213211088454516e308, TEST_TOL3); > > I've computed them using Pari as > > c(L,n) = { 2^L * exp(-Pi*n/2) * abs(gamma(L+1+I*n)) / gamma(2*L+2)} > cofac(L,n,r) = { I * exp(-I*r)*r^(-L)/(gamma(2*L+2)*c(L,n)) } > FplusIG(L,n,r) = { cofac(L,n,r) * intnum(t=0,[[1],1],exp(-t)*t^(L-I*n)*(t+2*I*r)^(L+I*n)) } > > FplusIG(37,0,1.2693881947287221e-07+x) # for F and F' (real part) > FplusIG(36,0,1.2693881947287221e-07+x) # for G and G' (imaginary part) > Hi Brian, It looks like that these values are incorrect since both F and Fp should be very small (correct me if I'm wrong since I'm not familiar with Pari and gsl test system). If you need verify the values than its better to use bessel function series. For eta == 0. and integer lambda coulomb function are just bessel-ricatti functions, i.e. spherical_bessel functions times argument F_l(x) = j_l (x) * x G_l(x) = n_l (x) * x Thus F_l(x) = 1/(2l+1)!! x^(l+1) [1 - x^2/(2l+3)/2 + ...] G_l(x) = (2l+1)!!/(2l+1) x^(-l) [1 + x^2/(2l-1)/2 + ...] I think Abramowitz, Stegun should have the complete expressions. -- Best regards, Alexey Alexandrovich Illarionov From MAILER-DAEMON Fri Aug 27 13:53:22 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1Op37C-0003yn-3r for mharc-bug-gsl@gnu.org; Fri, 27 Aug 2010 13:53:22 -0400 Received: from [140.186.70.92] (port=47782 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1Op379-0003yC-FA for bug-gsl@gnu.org; Fri, 27 Aug 2010 13:53:20 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1Op377-0004ku-GA for bug-gsl@gnu.org; Fri, 27 Aug 2010 13:53:19 -0400 Received: from mail-qw0-f41.google.com ([209.85.216.41]:43210) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1Op377-0004kn-Ch for bug-gsl@gnu.org; Fri, 27 Aug 2010 13:53:17 -0400 Received: by qwf7 with SMTP id 7so3909688qwf.0 for ; Fri, 27 Aug 2010 10:53:16 -0700 (PDT) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:sender:received :in-reply-to:references:date:x-google-sender-auth:message-id:subject :from:to:cc:content-type:content-transfer-encoding; bh=ltV8vnKPQ1ygR0EURQNQ9Ag7SWUNhxnO79HrXxeFL8M=; b=ZVVsSN05oHLWQaWRhj+dkqVjlx7+0NheMy2Yl6sZfgZg80aw4RAuexUUbeV9qf2ust lqhDF+kae6K9arKRKq67WPeq4CAWOtSLdb0h+DegzCMeUmbzpHLlpSuAuKpcU6YbVHJZ MdxaOynoCSlo4vo4VDC410m+Vn0WWdcsvjTCM= DomainKey-Signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:sender:in-reply-to:references:date :x-google-sender-auth:message-id:subject:from:to:cc:content-type :content-transfer-encoding; b=C622wB7EjmmX8geg6eI8cGquolpNKGwkxLuyhGf0yOALW1XVf5d+aarZqCV4PdmFJ3 yutkZAi1HsJUqUNaDi8+j3EUCus7thPrvW8E+EDC0wePYmC7MC9kr8J0+3pC6GGMbtzK Z2rdNjcdYS540INxSz2nYvFRVW1dRmr36NBMg= MIME-Version: 1.0 Received: by 10.229.214.13 with SMTP id gy13mr763747qcb.155.1282931596239; Fri, 27 Aug 2010 10:53:16 -0700 (PDT) Sender: tuomo.keskitalo@gmail.com Received: by 10.229.213.75 with HTTP; Fri, 27 Aug 2010 10:53:16 -0700 (PDT) In-Reply-To: <20100826-120528.sv79421.31356@savannah.gnu.org> References: <20100724-121507.sv79421.66062@savannah.gnu.org> <20100806-140911.sv1824.76026@savannah.gnu.org> <20100826-120528.sv79421.31356@savannah.gnu.org> Date: Fri, 27 Aug 2010 20:53:16 +0300 X-Google-Sender-Auth: hynsRCBDQNnlPCRCeGFq8Dc5ZJg Message-ID: Subject: Re: [Bug-gsl] [bug #30540] please, add convergence checks in rk4imp/rk2imp From: Tuomo Keskitalo To: Sergey B Kirpichev Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) Cc: Brian Gough , bug-gsl@gnu.org X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Fri, 27 Aug 2010 17:53:20 -0000 Hi, A quick comment: You need to specify the function for calculating the Jacobian for sys-henon, the new implicit codes require it. Please look at documentation in doc/ode-initval.texi. Use of fixed step size is similar to old ode-initval. BR, Tuomo On Thu, Aug 26, 2010 at 3:05 PM, Sergey B Kirpichev wrote: > > Follow-up Comment #2, bug #30540 (project gsl): > > I can't figure out how to use this code for fixed step size. =A0imprk2/im= prk4 > seems to be not working at all for this case: > $ ./bug > Segmentation fault > > Other steppers works (e.g. rk2) just fine. =A0See bug-ode2.c. > > Second, Newton-type iterations could be inefficient in general. =A0See: E= rnst > Hairer, Christian Lubich, Gerhard Wanner, Geometric Numerical Integration= . > Structure-Preserving Algorithms for Ordinary Differential Equations. Spri= nger > Series in Comput. Mathematics, Vol. 31, Springer-Verlag 2002. pp. 330-332= . > > > (file #21314) > =A0 =A0_______________________________________________________ > > Additional Item Attachment: > > File name: bug-ode2.c =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 Size:2 KB > > > =A0 =A0_______________________________________________________ > > Reply to this item at: > > =A0 > > _______________________________________________ > =A0Message sent via/by Savannah > =A0http://savannah.gnu.org/ > > > _______________________________________________ > Bug-gsl mailing list > Bug-gsl@gnu.org > http://lists.gnu.org/mailman/listinfo/bug-gsl > --=20 Tuomo.Keskitalo@iki.fi http://iki.fi/tuomo.keskitalo From MAILER-DAEMON Sun Aug 29 07:04:28 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1Opfga-0003p8-0m for mharc-bug-gsl@gnu.org; Sun, 29 Aug 2010 07:04:28 -0400 Received: from [140.186.70.92] (port=56719 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OpNXm-0002Fq-BQ for bug-gsl@gnu.org; Sat, 28 Aug 2010 11:42:13 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1OpNXf-0006Mv-AK for bug-gsl@gnu.org; Sat, 28 Aug 2010 11:42:04 -0400 Received: from mailout-de.gmx.net ([213.165.64.22]:53406 helo=mail.gmx.net) by eggs.gnu.org with smtp (Exim 4.69) (envelope-from ) id 1OpNXe-0006Mb-O3 for bug-gsl@gnu.org; Sat, 28 Aug 2010 11:42:03 -0400 Received: (qmail invoked by alias); 28 Aug 2010 15:41:57 -0000 Received: from g226140091.adsl.alicedsl.de (EHLO birne.lan) [92.226.140.91] by mail.gmx.net (mp058) with SMTP; 28 Aug 2010 17:41:57 +0200 X-Authenticated: #600002 X-Provags-ID: V01U2FsdGVkX19dPvQgqIpdzuOh32/Rqs1IH7cd7PlZ0wSIFwWUL6 8C6RNIGCZOIirl Message-ID: <4C792E41.3090007@gmx.de> Date: Sat, 28 Aug 2010 17:41:53 +0200 From: Thomas Tanner User-Agent: Mozilla/5.0 (Macintosh; U; Intel Mac OS X 10.6; de-DE; rv:1.9.2.8) Gecko/20100802 Thunderbird/3.1.2 MIME-Version: 1.0 To: bug-gsl@gnu.org X-Enigmail-Version: 1.1.1 Content-Type: text/plain; charset=ISO-8859-15 Content-Transfer-Encoding: 7bit X-Y-GMX-Trusted: 0 X-detected-operating-system: by eggs.gnu.org: Genre and OS details not recognized. X-Mailman-Approved-At: Sun, 29 Aug 2010 07:04:26 -0400 Subject: [Bug-gsl] gsl_sf_beta_inc underflow error X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Sat, 28 Aug 2010 15:42:13 -0000 Hello, calling gsl_sf_beta_inc(0.5,101.5,0.999457); with GSL version 1.14 results in gsl: exp.c:545: ERROR: underflow Default GSL error handler invoked. Abort trap Other libraries (e.g., cephes) correctly return 1.0 as result. cheers, -- Thomas Tanner ------ email: tanner@gmx.de GnuPG: 1024/5924D4DD From MAILER-DAEMON Mon Aug 30 10:33:28 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1Oq5QN-0004Je-EF for mharc-bug-gsl@gnu.org; Mon, 30 Aug 2010 10:33:27 -0400 Received: from [199.232.76.173] (port=53351 helo=monty-python.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1Oq5QK-0004JG-9s for bug-gsl@gnu.org; Mon, 30 Aug 2010 10:33:24 -0400 Received: from Debian-exim by monty-python.gnu.org with spam-scanned (Exim 4.60) (envelope-from ) id 1Oq5QJ-0004Bb-Jg for bug-gsl@gnu.org; Mon, 30 Aug 2010 10:33:24 -0400 Received: from fencepost.gnu.org ([140.186.70.10]:40753) by monty-python.gnu.org with esmtp (Exim 4.60) (envelope-from ) id 1Oq5Q9-0004BG-ET for bug-gsl@gnu.org; Mon, 30 Aug 2010 10:33:23 -0400 Received: from localhost ([127.0.0.1]:41481 helo=fencepost.gnu.org) by fencepost.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1Oq5PL-0007MD-T6; Mon, 30 Aug 2010 10:32:23 -0400 Date: Sun, 29 Aug 2010 14:12:28 +0100 Message-ID: <43pqx1o94z.wl%bjg@gnu.org> From: Brian Gough To: "Alexey A. Illarionov" Subject: Re: [Bug-gsl] bug in gsl_sf_coulomb_wave_FG_e (gsl 1.14 compiled with gcc 4.5.0) In-Reply-To: <4C77F788.2090102@gmail.com> References: <4C75E25A.7050307@gmail.com> <43tymgnxhz.wl%bjg@gnu.org> <4C77F788.2090102@gmail.com> User-Agent: Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) MIME-Version: 1.0 (generated by SEMI 1.14.6 - "Maruoka") Content-Type: text/plain; charset=US-ASCII X-detected-operating-system: by monty-python.gnu.org: GNU/Linux 2.6 (newer, 3) Cc: bug-gsl@gnu.org X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Mon, 30 Aug 2010 14:33:24 -0000 At Fri, 27 Aug 2010 13:36:08 -0400, Alexey A. Illarionov wrote: > It looks like that these values are incorrect since both F and Fp should > be very small (correct me if I'm wrong since I'm not familiar with Pari > and gsl test system). If you need verify the values than its better to > use bessel function series. > > For eta == 0. and integer lambda coulomb function are just > bessel-ricatti functions, i.e. spherical_bessel functions times argument > F_l(x) = j_l (x) * x > G_l(x) = n_l (x) * x > Thus > F_l(x) = 1/(2l+1)!! x^(l+1) [1 - x^2/(2l+3)/2 + ...] > G_l(x) = (2l+1)!!/(2l+1) x^(-l) [1 + x^2/(2l-1)/2 + ...] > I think Abramowitz, Stegun should have the complete expressions. Thanks that is much better, the calculation of the integral was not reliable. I've committed the new test values. From MAILER-DAEMON Tue Aug 31 14:04:12 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OqVBs-0002Cw-8l for mharc-bug-gsl@gnu.org; Tue, 31 Aug 2010 14:04:12 -0400 Received: from [140.186.70.92] (port=39829 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OqVBl-00027q-7V for bug-gsl@gnu.org; Tue, 31 Aug 2010 14:04:10 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1OqVBd-0001Av-Hj for bug-gsl@gnu.org; Tue, 31 Aug 2010 14:04:04 -0400 Received: from colonialone.fsf.org ([140.186.70.51]:36545 helo=internal.in.savannah.gnu.org) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OqVBd-0001Am-Fm for bug-gsl@gnu.org; Tue, 31 Aug 2010 14:03:57 -0400 Received: from [10.1.0.103] (helo=frontend.in.savannah.gnu.org) by internal.in.savannah.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OqVBd-0007bS-8C for bug-gsl@gnu.org; Tue, 31 Aug 2010 18:03:57 +0000 Received: from www-data by frontend.in.savannah.gnu.org with local (Exim 4.69) (envelope-from ) id 1OqVBb-0003M6-4C for bug-gsl@gnu.org; Tue, 31 Aug 2010 18:03:57 +0000 To: bug-gsl@gnu.org From: Brian Gough X-Savane-Server: savannah.gnu.org:443 [10.1.0.103] MIME-Version: 1.0 Content-Type: text/plain;charset=UTF-8 X-Savane-Project: gsl X-Savane-Tracker: bugs X-Savane-Item-ID: 30933 User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-GB; rv:1.9.0.11) Gecko/2009061118 Fedora/3.0.11-1.fc9 Firefox/3.0.11 X-Apparently-From: 94.193.216.241 (Savane authenticated user bjg) Message-Id: <20100831-190354.sv1824.33560@savannah.gnu.org> References: In-Reply-To: Date: Tue, 31 Aug 2010 18:03:55 +0000 X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) Subject: [Bug-gsl] [bug #30933] gsl_sf_beta_inc underflow error X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 31 Aug 2010 18:04:10 -0000 URL: Summary: gsl_sf_beta_inc underflow error Project: GNU Scientific Library Submitted by: bjg Submitted on: Tue 31 Aug 2010 07:03:54 PM BST Category: Runtime error Severity: 3 - Normal Operating System: Status: Confirmed Assigned to: None Open/Closed: Open Release: 1.14 Discussion Lock: Any _______________________________________________________ Details: From: Thomas Tanner To: bug-gsl@gnu.org Subject: [Bug-gsl] gsl_sf_beta_inc underflow error Date: Sat, 28 Aug 2010 17:41:53 +0200 Hello, calling gsl_sf_beta_inc(0.5,101.5,0.999457); with GSL version 1.14 results in gsl: exp.c:545: ERROR: underflow Default GSL error handler invoked. Abort trap Other libraries (e.g., cephes) correctly return 1.0 as result. cheers, -- Thomas Tanner ------ email: tanner@gmx.de GnuPG: 1024/5924D4DD _______________________________________________________ Reply to this item at: _______________________________________________ Message sent via/by Savannah http://savannah.gnu.org/ From MAILER-DAEMON Tue Aug 31 14:15:42 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OqVMz-0006A4-VQ for mharc-bug-gsl@gnu.org; Tue, 31 Aug 2010 14:15:41 -0400 Received: from [140.186.70.92] (port=36315 helo=eggs.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OqVMx-00069w-Oy for bug-gsl@gnu.org; Tue, 31 Aug 2010 14:15:40 -0400 Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69) (envelope-from ) id 1OqVMw-0003Qk-ME for bug-gsl@gnu.org; Tue, 31 Aug 2010 14:15:39 -0400 Received: from colonialone.fsf.org ([140.186.70.51]:48213 helo=internal.in.savannah.gnu.org) by eggs.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OqVMw-0003Qf-Ku for bug-gsl@gnu.org; Tue, 31 Aug 2010 14:15:38 -0400 Received: from [10.1.0.103] (helo=frontend.in.savannah.gnu.org) by internal.in.savannah.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OqVMw-0007m4-Af for bug-gsl@gnu.org; Tue, 31 Aug 2010 18:15:38 +0000 Received: from www-data by frontend.in.savannah.gnu.org with local (Exim 4.69) (envelope-from ) id 1OqVMu-0003W3-NF for bug-gsl@gnu.org; Tue, 31 Aug 2010 18:15:37 +0000 To: bug-gsl@gnu.org From: Brian Gough X-Savane-Server: savannah.gnu.org:443 [10.1.0.103] MIME-Version: 1.0 Content-Type: text/plain;charset=UTF-8 X-Savane-Project: gsl X-Savane-Tracker: bugs X-Savane-Item-ID: 30933 User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-GB; rv:1.9.0.11) Gecko/2009061118 Fedora/3.0.11-1.fc9 Firefox/3.0.11 X-Apparently-From: 94.193.216.241 (Savane authenticated user bjg) Message-Id: <20100831-191536.sv1824.64304@savannah.gnu.org> References: <20100831-190354.sv1824.33560@savannah.gnu.org> In-Reply-To: <20100831-190354.sv1824.33560@savannah.gnu.org> Date: Tue, 31 Aug 2010 18:15:36 +0000 X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2) Subject: [Bug-gsl] [bug #30933] gsl_sf_beta_inc underflow error X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 31 Aug 2010 18:15:40 -0000 Update of bug #30933 (project gsl): Status: Confirmed => Fixed Open/Closed: Open => Closed _______________________________________________________ Follow-up Comment #1: fixed by revision 4680. _______________________________________________________ Reply to this item at: _______________________________________________ Message sent via/by Savannah http://savannah.gnu.org/ From MAILER-DAEMON Tue Aug 31 14:16:59 2010 Received: from mailman by lists.gnu.org with archive (Exim 4.43) id 1OqVOF-0006Qs-5W for mharc-bug-gsl@gnu.org; Tue, 31 Aug 2010 14:16:59 -0400 Received: from [199.232.76.173] (port=59703 helo=monty-python.gnu.org) by lists.gnu.org with esmtp (Exim 4.43) id 1OqVOD-0006Qk-Bc for bug-gsl@gnu.org; Tue, 31 Aug 2010 14:16:57 -0400 Received: from Debian-exim by monty-python.gnu.org with spam-scanned (Exim 4.60) (envelope-from ) id 1OqVOC-0006J8-FH for bug-gsl@gnu.org; Tue, 31 Aug 2010 14:16:57 -0400 Received: from fencepost.gnu.org ([140.186.70.10]:56859) by monty-python.gnu.org with esmtp (Exim 4.60) (envelope-from ) id 1OqVOC-0006J3-8X for bug-gsl@gnu.org; Tue, 31 Aug 2010 14:16:56 -0400 Received: from localhost ([127.0.0.1]:59583 helo=fencepost.gnu.org) by fencepost.gnu.org with esmtp (Exim 4.69) (envelope-from ) id 1OqVOB-00073S-4y; Tue, 31 Aug 2010 14:16:55 -0400 Date: Tue, 31 Aug 2010 19:16:47 +0100 Message-ID: <43mxs2odf4.wl%bjg@gnu.org> From: Brian Gough To: Thomas Tanner Subject: Re: [Bug-gsl] gsl_sf_beta_inc underflow error In-Reply-To: <4C792E41.3090007@gmx.de> References: <4C792E41.3090007@gmx.de> User-Agent: Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) MIME-Version: 1.0 (generated by SEMI 1.14.6 - "Maruoka") Content-Type: text/plain; charset=US-ASCII X-detected-operating-system: by monty-python.gnu.org: GNU/Linux 2.6 (newer, 3) Cc: bug-gsl@gnu.org X-BeenThere: bug-gsl@gnu.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports for the GNU Scientific Library List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 31 Aug 2010 18:16:58 -0000 At Sat, 28 Aug 2010 17:41:53 +0200, Thomas Tanner wrote: > > Hello, > > calling gsl_sf_beta_inc(0.5,101.5,0.999457); > > with GSL version 1.14 results in > > gsl: exp.c:545: ERROR: underflow > Default GSL error handler invoked. > Abort trap > > Other libraries (e.g., cephes) correctly return 1.0 as result. > Thanks for the bug report. I've removed the incorrect error code, which came from an intermediate underflow that does not affect the final result. -- Brian Gough