1/* Compute projection of complex float type value to Riemann sphere.
2 Copyright (C) 1997-2022 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
4
5 The GNU C Library is free software; you can redistribute it and/or
6 modify it under the terms of the GNU Lesser General Public
7 License as published by the Free Software Foundation; either
8 version 2.1 of the License, or (at your option) any later version.
9
10 The GNU C Library is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 Lesser General Public License for more details.
14
15 You should have received a copy of the GNU Lesser General Public
16 License along with the GNU C Library; if not, see
17 <https://www.gnu.org/licenses/>. */
18
19#include <complex.h>
20#include <math.h>
21
22
23CFLOAT
24M_DECL_FUNC (__cproj) (CFLOAT x)
25{
26 if (isinf (__real__ x) || isinf (__imag__ x))
27 {
28 CFLOAT res;
29
30 __real__ res = INFINITY;
31 __imag__ res = M_COPYSIGN (0, __imag__ x);
32
33 return res;
34 }
35
36 return x;
37}
38
39declare_mgen_alias (__cproj, cproj)
40

source code of glibc/math/s_cproj_template.c