xref: /openbsd-src/lib/libm/src/s_ccoshf.c (revision 60d908e137a0315099596747e9403ad9e8c76849)
1*60d908e1Sguenther /*	$OpenBSD: s_ccoshf.c,v 1.2 2010/07/18 18:42:26 guenther Exp $	*/
27b36286aSmartynas /*
37b36286aSmartynas  * Copyright (c) 2008 Stephen L. Moshier <steve@moshier.net>
47b36286aSmartynas  *
57b36286aSmartynas  * Permission to use, copy, modify, and distribute this software for any
67b36286aSmartynas  * purpose with or without fee is hereby granted, provided that the above
77b36286aSmartynas  * copyright notice and this permission notice appear in all copies.
87b36286aSmartynas  *
97b36286aSmartynas  * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
107b36286aSmartynas  * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
117b36286aSmartynas  * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
127b36286aSmartynas  * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
137b36286aSmartynas  * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
147b36286aSmartynas  * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
157b36286aSmartynas  * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
167b36286aSmartynas  */
177b36286aSmartynas 
187b36286aSmartynas /*							ccoshf
197b36286aSmartynas  *
207b36286aSmartynas  *	Complex hyperbolic cosine
217b36286aSmartynas  *
227b36286aSmartynas  *
237b36286aSmartynas  *
247b36286aSmartynas  * SYNOPSIS:
257b36286aSmartynas  *
267b36286aSmartynas  * float complex ccoshf();
277b36286aSmartynas  * float complex z, w;
287b36286aSmartynas  *
297b36286aSmartynas  * w = ccoshf (z);
307b36286aSmartynas  *
317b36286aSmartynas  *
327b36286aSmartynas  *
337b36286aSmartynas  * DESCRIPTION:
347b36286aSmartynas  *
357b36286aSmartynas  * ccosh(z) = cosh x  cos y + i sinh x sin y .
367b36286aSmartynas  *
377b36286aSmartynas  * ACCURACY:
387b36286aSmartynas  *
397b36286aSmartynas  *                      Relative error:
407b36286aSmartynas  * arithmetic   domain     # trials      peak         rms
417b36286aSmartynas  *    IEEE      -10,+10     30000       2.9e-16     8.1e-17
427b36286aSmartynas  *
437b36286aSmartynas  */
447b36286aSmartynas 
457b36286aSmartynas #include <complex.h>
467b36286aSmartynas #include <math.h>
477b36286aSmartynas 
487b36286aSmartynas float complex
ccoshf(float complex z)497b36286aSmartynas ccoshf(float complex z)
507b36286aSmartynas {
517b36286aSmartynas 	float complex w;
527b36286aSmartynas 	float x, y;
537b36286aSmartynas 
54*60d908e1Sguenther 	x = crealf(z);
55*60d908e1Sguenther 	y = cimagf(z);
567b36286aSmartynas 	w = coshf (x) * cosf (y)  +  (sinhf (x) * sinf (y)) * I;
577b36286aSmartynas 	return (w);
587b36286aSmartynas }
59