[Bug-gsl] [bug #22651] arccosh for complex numbers - branch problem?

[Brian Gough]

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/