1 /* $OpenBSD: s_cprojf.c,v 1.3 2011/05/31 21:13:44 martynas Exp $ */
2 /*
3 * Copyright (c) 2008 Martynas Venckus <martynas@openbsd.org>
4 *
5 * Permission to use, copy, modify, and distribute this software for any
6 * purpose with or without fee is hereby granted, provided that the above
7 * copyright notice and this permission notice appear in all copies.
8 *
9 * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
10 * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
11 * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
12 * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
13 * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
14 * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
15 * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
16 */
17
18 #include <complex.h>
19 #include <math.h>
20
21 float complex
cprojf(float complex z)22 cprojf(float complex z)
23 {
24 float complex res;
25
26 if (isinf(__real__ z) || isinf(__imag__ z)) {
27 __real__ res = INFINITY;
28 __imag__ res = copysignf(0.0, __imag__ z);
29 } else {
30 res = z;
31 }
32
33 return res;
34 }
35