1*7b36286aSmartynas /* $OpenBSD: s_catanhf.c,v 1.1 2008/09/07 20:36:09 martynas Exp $ */
2*7b36286aSmartynas /*
3*7b36286aSmartynas * Copyright (c) 2008 Stephen L. Moshier <steve@moshier.net>
4*7b36286aSmartynas *
5*7b36286aSmartynas * Permission to use, copy, modify, and distribute this software for any
6*7b36286aSmartynas * purpose with or without fee is hereby granted, provided that the above
7*7b36286aSmartynas * copyright notice and this permission notice appear in all copies.
8*7b36286aSmartynas *
9*7b36286aSmartynas * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
10*7b36286aSmartynas * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
11*7b36286aSmartynas * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
12*7b36286aSmartynas * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
13*7b36286aSmartynas * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
14*7b36286aSmartynas * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
15*7b36286aSmartynas * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
16*7b36286aSmartynas */
17*7b36286aSmartynas
18*7b36286aSmartynas /* catanhf
19*7b36286aSmartynas *
20*7b36286aSmartynas * Complex inverse hyperbolic tangent
21*7b36286aSmartynas *
22*7b36286aSmartynas *
23*7b36286aSmartynas *
24*7b36286aSmartynas * SYNOPSIS:
25*7b36286aSmartynas *
26*7b36286aSmartynas * float complex catanhf();
27*7b36286aSmartynas * float complex z, w;
28*7b36286aSmartynas *
29*7b36286aSmartynas * w = catanhf (z);
30*7b36286aSmartynas *
31*7b36286aSmartynas *
32*7b36286aSmartynas *
33*7b36286aSmartynas * DESCRIPTION:
34*7b36286aSmartynas *
35*7b36286aSmartynas * Inverse tanh, equal to -i catan (iz);
36*7b36286aSmartynas *
37*7b36286aSmartynas * ACCURACY:
38*7b36286aSmartynas *
39*7b36286aSmartynas * Relative error:
40*7b36286aSmartynas * arithmetic domain # trials peak rms
41*7b36286aSmartynas * IEEE -10,+10 30000 2.3e-16 6.2e-17
42*7b36286aSmartynas *
43*7b36286aSmartynas */
44*7b36286aSmartynas
45*7b36286aSmartynas #include <complex.h>
46*7b36286aSmartynas #include <math.h>
47*7b36286aSmartynas
48*7b36286aSmartynas float complex
catanhf(float complex z)49*7b36286aSmartynas catanhf(float complex z)
50*7b36286aSmartynas {
51*7b36286aSmartynas float complex w;
52*7b36286aSmartynas
53*7b36286aSmartynas w = -1.0f * I * catanf (z * I);
54*7b36286aSmartynas return (w);
55*7b36286aSmartynas }
56