FirstIndexLast

[bug-gsl] [bug #22651] arccosh for complex numbers - branch problem? (2008-11-19)

From: Brian Gough <INVALID.NOREPLY@gnu.org>
Subject: [bug #22651] arccosh for complex numbers - branch problem?
Date: Wed, 19 Nov 2008 09:17:19 +0000
To: bug-gsl@gnu.org


Update of bug #22651 (project gsl):

Status: Need Info => Invalid
Open/Closed: Open => Closed
Release: => 1.11

_______________________________________________________

Follow-up Comment #1:

The information given in the email was incorrect. GSL returns the correct
value +0.96 for both real and complex versions of the function.

#include <gsl/gsl_complex_math.h>
#include <math.h>
#include <stdio.h>

int main ()
{
gsl_complex z,f;
double x, y, zero = 0.0;
x = 1.5;
y = acosh(x);
GSL_SET_COMPLEX(&z, x, zero);
f = gsl_complex_arccosh(z);
printf("%.18en", x);
printf("%.18en", y);
printf("%.18e %.18en", GSL_REAL(z), GSL_IMAG(z));
printf("%.18e %.18en", GSL_REAL(f), GSL_IMAG(f));
}


$ ./a.out
1.500000000000000000e+00
9.624236501192069415e-01
1.500000000000000000e+00 0.000000000000000000e+00
9.624236501192069415e-01 -0.000000000000000000e+00


_______________________________________________________

Reply to this item at:

<http://savannah.gnu.org/bugs/?22651>

_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/