1/* Copy a sign bit between floating-point values.
2   IBM extended format long double version.
3   Copyright (C) 2004-2021 Free Software Foundation, Inc.
4   This file is part of the GNU C Library.
5
6   The GNU C Library is free software; you can redistribute it and/or
7   modify it under the terms of the GNU Lesser General Public
8   License as published by the Free Software Foundation; either
9   version 2.1 of the License, or (at your option) any later version.
10
11   The GNU C Library is distributed in the hope that it will be useful,
12   but WITHOUT ANY WARRANTY; without even the implied warranty of
13   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14   Lesser General Public License for more details.
15
16   You should have received a copy of the GNU Lesser General Public
17   License along with the GNU C Library; if not, see
18   <https://www.gnu.org/licenses/>.  */
19
20#include <sysdep.h>
21#include <math_ldbl_opt.h>
22
23ENTRY_TOCLESS (__copysignl)
24/* long double [f1,f2] copysign (long double [f1,f2] x, long double [f3,f4] y);
25   copysign(x,y) returns a value with the magnitude of x and
26   with the sign bit of y.  */
27	stfd	fp3,-16(r1)
28	ld	r3,-16(r1)
29	cmpdi   r3,0
30	blt     L(0)
31	fmr	fp0,fp1
32	fabs	fp1,fp1
33	fneg	fp3,fp2
34	fsel	fp2,fp0,fp2,fp3
35	blr
36L(0):
37	fmr	fp0,fp1
38	fnabs   fp1,fp1
39	fneg	fp3,fp2
40	fsel	fp2,fp0,fp3,fp2
41	blr
42END (__copysignl)
43
44#if IS_IN (libm)
45long_double_symbol (libm, __copysignl, copysignl)
46#else
47long_double_symbol (libc, __copysignl, copysignl)
48#endif
49