From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <gsl-discuss-return-5162-listarch-gsl-discuss=sources.redhat.com@sourceware.org>
Received: (qmail 18922 invoked by alias); 25 Sep 2007 02:58:09 -0000
Received: (qmail 18914 invoked by uid 22791); 25 Sep 2007 02:58:08 -0000
X-Spam-Check-By: sourceware.org
Received: from mx02.em-net.ne.jp (HELO mx02.em-net.ne.jp) (210.131.159.52)
    by sourceware.org (qpsmtpd/0.31) with ESMTP; Tue, 25 Sep 2007 02:58:06 +0000
Received: from smtp02.em-net.ne.jp (smtp02.em-net.ne.jp [210.131.159.14])
	by mx02.em-net.ne.jp (Switch-3.2.5/Switch-3.2.5) with ESMTP id l8P2w0ql032055
	for <gsl-discuss@sourceware.org>; Tue, 25 Sep 2007 11:58:00 +0900
Received: from mm.em-net.ne.jp (static-200-105-171-95.acelerate.net [200.105.171.95])
	by smtp02.em-net.ne.jp (Postfix) with SMTP id C709C50
	for <gsl-discuss@sourceware.org>; Tue, 25 Sep 2007 11:57:59 +0900 (JST)
From: Yoshiki TSUNESADA <y-tsunesada@mm.em-net.ne.jp>
To: gsl-discuss@sourceware.org
Subject: test_bessel.c: bessel_Jn_e(0, 20000)
Date: Tue, 25 Sep 2007 02:58:00 -0000
MIME-Version: 1.0
Content-Type: text/plain; charset=us-ascii
Content-Transfer-Encoding: 7bit
X-Mailer: HidemaruMail 4.81 (WinNT,501)
Message-Id: <82C7FF1FE827DBy-tsunesada@mm.em-net.ne.jp>
Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm
Precedence: bulk
List-Id: <gsl-discuss.sourceware.org>
List-Subscribe: <mailto:gsl-discuss-subscribe@sourceware.org>
List-Archive: <http://sourceware.org/ml/gsl-discuss/>
List-Post: <mailto:gsl-discuss@sourceware.org>
List-Help: <mailto:gsl-discuss-help@sourceware.org>, <http://sourceware.org/ml/#faqs>
Sender: gsl-discuss-owner@sourceware.org
X-SW-Source: 2007-q3/txt/msg00039.txt.bz2

 I encountered a glitch: A test code fails, in "specfunc/test_bessel.c", 
line 61,

  TEST_SF(s, gsl_sf_bessel_Jn_e, (0, 20000, &r), 
          -7.676508175684157103e-06, TEST_TOL4, GSL_SUCCESS);
          
The "expected" value -7.67651e-06 is exactly same as the one for
the test of line 60, the preceding line,

  TEST_SF(s, gsl_sf_bessel_Jn_e, (1, 1.0e+10, &r), 
          -7.676508175684157103e-06, TEST_TOL4, GSL_SUCCESS);

and this test passes. I guess the expected value in the line 61 
is mis-pasted. The following is an example session using Ruby/GSL:

   yoshiki@poirot[71] irb
   irb(main):001:0> require("rbgsl")
   => true
   irb(main):002:0> GSL::Sf::bessel_Jn(0, 20000)
   => 0.00556597490495494
   irb(main):003:0> GSL::Sf::bessel_Jn(1, 1.0e10)
   => -7.67650817568416e-06
 
 I don't have any other tools to compare and verify these results.
 
 Best regards,
 Yoshiki