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