[Bug-gsl] [bug #29773] gsl_complex_pow_real(gsl_complex_rect(0, 0), 0) returns (0, 0)

[Rhys Ulerich]

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/