25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/*
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * CDDL HEADER START
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis *
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 *
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 *
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 *
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * CDDL HEADER END
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/*
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Use is subject to license terms.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/*
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Copyright 2011, Richard Lowe.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/* Functions in this file are duplicated in locallibm.il. Keep them in sync */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#ifndef _LIBM_INLINES_H
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#define _LIBM_INLINES_H
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#ifdef __GNUC__
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#include <sys/types.h>
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#include <sys/ieeefp.h>
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#ifdef __cplusplus
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtisextern "C" {
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#endif
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE enum fp_class_type
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtisfp_classf(float f)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis enum fp_class_type ret;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis int fint; /* scratch for f as int */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis uint64_t tmp;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__(
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "fabss %3,%3\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "st %3,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ld %1,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "orcc %%g0,%0,%%g0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "be,pn %%icc,2f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "1:\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "sethi %%hi(0x7f800000),%2\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andcc %0,%2,%%g0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "bne,pt %%icc,1f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %%g0,1,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ba 2f\n\t" /* subnormal */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "1:\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "subcc %0,%2,%%g0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "bge,pn %%icc,1f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %%g0,2,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ba 2f\n\t" /* normal */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "1:\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "bg,pn %%icc,1f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %%g0,3,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ba 2f\n\t" /* infinity */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "1:\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "sethi %%hi(0x00400000),%2\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andcc %0,%2,%%g0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %%g0,4,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "bne,pt %%icc,2f\n\t" /* quiet NaN */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %%g0,5,%0\n\t" /* signalling NaN */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "2:\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "=r" (ret), "=m" (fint), "=r" (tmp), "+f" (f)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis :
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (ret);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE enum fp_class_type
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtisfp_class(double d)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis enum fp_class_type ret;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis uint64_t dint; /* Scratch for d-as-long */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis uint64_t tmp;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__(
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "fabsd %3,%3\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "std %3,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ldx %1,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "orcc %%g0,%0,%%g0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "be,pn %%xcc,2f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "sethi %%hi(0x7ff00000),%2\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "sllx %2,32,%2\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andcc %0,%2,%%g0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "bne,pt %%xcc,1f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %%g0,1,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ba 2f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "1:\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "subcc %0,%2,%%g0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "bge,pn %%xcc,1f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %%g0,2,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ba 2f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "1:\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andncc %0,%2,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "bne,pn %%xcc,1f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %%g0,3,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ba 2f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "1:\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "sethi %%hi(0x00080000),%2\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "sllx %2,32,%2\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andcc %0,%2,%%g0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %%g0,4,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "bne,pt %%xcc,2f\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "nop\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %%g0,5,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "2:\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "=r" (ret), "=m" (dint), "=r" (tmp), "+e" (d)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis :
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (ret);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE float
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis__inline_sqrtf(float f)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis float ret;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fsqrts %1,%0\n\t" : "=f" (ret) : "f" (f));
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (ret);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE double
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis__inline_sqrt(double d)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis double ret;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fsqrtd %1,%0\n\t" : "=f" (ret) : "f" (d));
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (ret);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE int
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis__swapEX(int i)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis int ret;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis uint32_t fsr;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis uint64_t tmp1, tmp2;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__(
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "and %4,0x1f,%2\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "sll %2,5,%2\n\t" /* shift input to aexc bit location */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ".volatile\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "st %%fsr,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ld %1,%0\n\t" /* %0 = fsr */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andn %0,0x3e0,%3\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %2,%3,%2\n\t" /* %2 = new fsr */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "st %2,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ld %1,%%fsr\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "srl %0,5,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "and %0,0x1f,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ".nonvolatile\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "=r" (ret), "=m" (fsr), "=r" (tmp1), "=r" (tmp2)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "r" (i)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (ret);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/*
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * On the SPARC, __swapRP is a no-op; always return 0 for backward
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * compatibility
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/* ARGSUSED */
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE enum fp_precision_type
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis__swapRP(enum fp_precision_type i)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (0);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE enum fp_direction_type
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis__swapRD(enum fp_direction_type d)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis enum fp_direction_type ret;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis uint32_t fsr;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis uint64_t tmp1, tmp2, tmp3;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__(
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "and %5,0x3,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "sll %0,30,%2\n\t" /* shift input to RD bit location */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ".volatile\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "st %%fsr,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ld %1,%0\n\t" /* %0 = fsr */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* mask of rounding direction bits */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "sethi %%hi(0xc0000000),%4\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andn %0,%4,%3\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %2,%3,%2\n\t" /* %2 = new fsr */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "st %2,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ld %1,%%fsr\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "srl %0,30,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "and %0,0x3,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ".nonvolatile\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "=r" (ret), "=m" (fsr), "=r" (tmp1), "=r" (tmp2), "=r" (tmp3)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "r" (d)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (ret);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE int
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis__swapTE(int i)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis int ret;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis uint32_t fsr;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis uint64_t tmp1, tmp2, tmp3;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__(
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "and %5,0x1f,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "sll %0,23,%2\n\t" /* shift input to TEM bit location */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ".volatile\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "st %%fsr,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ld %1,%0\n\t" /* %0 = fsr */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* mask of TEM (Trap Enable Mode bits) */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "sethi %%hi(0x0f800000),%4\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "andn %0,%4,%3\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "or %2,%3,%2\n\t" /* %2 = new fsr */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "st %2,%1\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "ld %1,%%fsr\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "srl %0,23,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis "and %0,0x1f,%0\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis ".nonvolatile\n\t"
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "=r" (ret), "=m" (fsr), "=r" (tmp1), "=r" (tmp2), "=r" (tmp3)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "r" (i)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis : "cc");
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (ret);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE double
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtissqrt(double d)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (__inline_sqrt(d));
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE float
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtissqrtf(float f)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (__inline_sqrtf(f));
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE double
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtisfabs(double d)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis double ret;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fabsd %1,%0\n\t" : "=e" (ret) : "e" (d));
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (ret);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
7f11fd00fc23e2af7ae21cc8837a2b86380dcfa7Richard Loweextern __GNU_INLINE float
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtisfabsf(float f)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis{
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis float ret;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __asm__ __volatile__("fabss %1,%0\n\t" : "=f" (ret) : "f" (f));
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (ret);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#ifdef __cplusplus
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis}
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#endif
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#endif /* __GNUC__ */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis#endif /* _LIBM_INLINES_H */