25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * CDDL HEADER START
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * The contents of this file are subject to the terms of the
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Common Development and Distribution License (the "License").
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * You may not use this file except in compliance with the License.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * or http://www.opensolaris.org/os/licensing.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * See the License for the specific language governing permissions
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * and limitations under the License.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * When distributing Covered Code, include this CDDL HEADER in each
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * If applicable, add the following below this CDDL HEADER, with the
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * fields enclosed by brackets "[]" replaced with your own identifying
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * information: Portions Copyright [yyyy] [name of copyright owner]
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * CDDL HEADER END
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Use is subject to license terms.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Copyright 2011, Richard Lowe
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/* Functions in this file are duplicated in locallibm.il. Keep them in sync */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fsqrt\n\t" : "=t" (ret) : "0" (a) : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fsqrt\n\t" : "=t" (ret) : "0" (a) : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andl $0x7fffffff,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "cmpl $0x43300000,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "1: fwait\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * 00 - 24 bits
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * 01 - reserved
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * 10 - 53 bits
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * 11 - 64 bits
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fstcw %0\n\t" : "=m" (cw));
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fldcw %0\n\t" : : "m" (cw));
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * 00 - Round to nearest, with even preferred
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * 01 - Round down
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * 10 - Round up
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fstcw %0\n\t" : "=m" (cw));
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fldcw %0\n\t" : : "m" (cw));
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Let's set a Rounding Control (RC) bits from x87 FPU Control Word
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * to fp_positive and save old bits in rd.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * The FRNDINT instruction returns a floating-point value that is the
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * integral value closest to the source value in the direction of the
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * rounding mode specified in the RC field of the x87 FPU control word.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Rounds the source value in the ST(0) register to the nearest
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * integral value, depending on the current rounding mode
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * (setting of the RC field of the FPU control word),
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * and stores the result in ST(0).
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("frndint" : "+t" (d) : : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* restore old RC bits */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andl $0x7fffffff,%0\n\t" /* %0 <-- hi_32(abs(d)) */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andl $0x80000000,%1\n\t" /* %1[31] <-- sign_bit(d2) */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "orl %1,%0\n\t" /* %0 <-- hi_32(copysign(x,y)) */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "+&r" (_HI_WORD(d1)), "+r" (_HI_WORD(d2))
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fabs\n\t" : "+t" (d) : : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fabs\n\t" : "+t" (d) : : "cc");
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE long double
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fabs\n\t" : "+t" (d) : : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andl $0x7ff00000,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "shrl $31,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("frndint" : "+t" (d), "+r" (rd) : : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * branchless __isnan
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * ((0x7ff00000-[((lx|-lx)>>31)&1]|ahx)>>31)&1 = 1 iff x is NaN
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Lowe "movl %1,%%ecx\n\t"
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Lowe "orl %%ecx,%1\n\t"
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Lowe "andl $0x7fffffff,%2\n\t" /* ecx <-- hi_32(abs(x)) */
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Lowe "orl %2,%1\n\t"
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Lowe "subl $0x7ff00000,%1\n\t"
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Lowe "negl %1\n\t"
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Lowe "shrl $31,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andl $0x7fffffff,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "addl $0x7f800000,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "shrl $31,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "fildl %2\n\t" /* Convert N to extended */
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE long double
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fsqrt" : "+t" (ld) : : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andl $0x00007fff,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "xorl $0x00007fff,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "testl $0x80000000,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "jz 3f\n\t" /* jump if leading bit is 0 */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "xorl %0,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "2:\n\t" /* note that %0 = 0 from before */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "cmpl $0x80000000,%1\n\t" /* what is first half of significand? */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "jnz 3f\n\t" /* jump if not equal to 0x80000000 */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "testl $0xffffffff,%2\n\t" /* is second half of significand 0? */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "movl $1,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#endif /* __GNUC__ */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#endif /* _LIBM_INLINES_H */