gsl_complex_pow_real(gsl_complex_rect(0,0),0) returns (0,0)

Brian Gough bjg@gnu.org
Tue May 4 11:52:00 GMT 2010


At Mon, 3 May 2010 13:40:46 -0500,
Rhys Ulerich wrote:
> I've not seen any references to this on the mailing list or in
> savannah.  Let me know if you'd like me to file a bug.

Thanks for finding that.  Please do file a bug in the tracker at
http://savannah.gnu.org/bugs/?group=gsl.  (A patch and test case would
also be handy!)

Looking at the source, there's an explicit check for 0^0 in
gsl_complex_pow but it is missing from gsl_complex_pow_real, so it
looks like a simple oversight.

-- 
Brian Gough

GNU Scientific Library -
http://www.gnu.org/software/gsl/



More information about the Gsl-discuss mailing list