URL:
<http://savannah.gnu.org/bugs/?29773>
Summary: gsl_complex_pow_real(gsl_complex_rect(0,0),0)
returns (0,0)
Project: GNU Scientific Library
Submitted by: rhysu
Submitted on: Tue 04 May 2010 03:31:33 PM GMT
Category: Accuracy problem
Severity: 3 - Normal
Operating System: Ubuntu 8.04+
Status: None
Assigned to: None
Open/Closed: Open
Release: 1.14
Discussion Lock: Any
_______________________________________________________
Details:
gsl_complex_pow_real incorrectly returns that (0+0i)^0 = (0+0i). Based on
the consistency with gsl_sf_pow_int and gsl_complex_pow, it should return
(1+0i) for these inputs.
Simple standalone test case attached. Problem is reproducible from at least
1.12 onward and with both the GNU and Intel compilers.
Original discussion at
http://www.cygwin.com/ml/gsl-discuss/2010-q2/msg00004.html
_______________________________________________________
File Attachments:
-------------------------------------------------------
Date: Tue 04 May 2010 03:31:33 PM GMT Name: testcase.c Size: 256B By:
rhysu
Standalone test case
<http://savannah.gnu.org/bugs/download.php?file_id=20427>
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?29773>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl