diff options
Diffstat (limited to 'sysdeps/ia64/fpu/e_remainderl.S')
| -rw-r--r-- | sysdeps/ia64/fpu/e_remainderl.S | 115 |
1 files changed, 65 insertions, 50 deletions
diff --git a/sysdeps/ia64/fpu/e_remainderl.S b/sysdeps/ia64/fpu/e_remainderl.S index 8c1630e3ab..5856861442 100644 --- a/sysdeps/ia64/fpu/e_remainderl.S +++ b/sysdeps/ia64/fpu/e_remainderl.S @@ -1,10 +1,10 @@ -.file "remainderl.s" - - -// Copyright (c) 2000 - 2003, Intel Corporation +.file "remainderl.asm" +// Copyright (C) 2000, 2001, Intel Corporation // All rights reserved. // -// Contributed 2000 by the Intel Numerics Group, Intel Corporation +// Contributed 2/2/2000 by John Harrison, Cristina Iordache, Ted Kubaska, +// Bob Norin, Shane Story, and Ping Tak Peter Tang of the Computational +// Software Lab, Intel Corporation. // // Redistribution and use in source and binary forms, with or without // modification, are permitted provided that the following conditions are @@ -20,7 +20,7 @@ // * The name of Intel Corporation may not be used to endorse or promote // products derived from this software without specific prior written // permission. - +// // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS // "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT // LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR @@ -35,19 +35,17 @@ // // Intel Corporation is the author of this code, and requests that all // problem reports or change requests be submitted to it directly at -// http://www.intel.com/software/products/opensource/libraries/num.htm. +// http://developer.intel.com/opensource. // // History //==================================================================== -// 02/02/00 Initial version -// 03/02/00 New algorithm -// 04/04/00 Unwind support added -// 07/21/00 Fixed quotient=2^{24*m+23}*1.q1...q23 1 bug -// 08/15/00 Bundle added after call to __libm_error_support to properly +// 2/02/00 Initial version +// 3/02/00 New algorithm +// 4/04/00 Unwind support added +// 7/21/00 Fixed quotient=2^{24*m+23}*1.q1...q23 1 bug +// 8/15/00 Bundle added after call to __libm_error_support to properly // set [the previously overwritten] GR_Parameter_RESULT. -// 11/29/00 Set FR_Y to f9 -// 05/20/02 Cleaned up namespace and sf0 syntax -// 02/10/03 Reordered header: .section, .global, .proc, .align +//11/29/00 Set FR_Y to f9 // // API //==================================================================== @@ -79,6 +77,9 @@ //==================================================================== // a=+/- Inf, or b=+/-0: return NaN, call libm_error_support // a=NaN or b=NaN: return NaN + +#include "libm_support.h" + // // Registers used //==================================================================== @@ -86,6 +87,8 @@ // General registers: r2,r3,r28,r29,r32 (ar.pfs), r33-r39 // Floating point registers: f6-f15,f32 // +.section .text + GR_SAVE_B0 = r33 GR_SAVE_PFS = r34 @@ -102,9 +105,19 @@ FR_Y = f9 FR_RESULT = f8 -.section .text -GLOBAL_IEEE754_ENTRY(remainderl) + + .proc remainderl# + .align 32 + .global remainderl# + .align 32 + +remainderl: +#ifdef _LIBC +.global __remainderl +.type __remainderl,@function +__remainderl: +#endif // inputs in f8, f9 // result in f8 @@ -146,7 +159,7 @@ cmp.eq p11,p10=r29,r0;; // X +-NAN, +-inf, ? p9 { .mfi nop.m 999 - fclass.m.unc p9,p8 = f8, 0xe3 +(p0) fclass.m.unc p9,p8 = f8, 0xe3 nop.i 999;; } @@ -183,8 +196,8 @@ cmp.eq p11,p10=r29,r0;; } {.bbb - (p9) br.cond.spnt FREM_X_NAN_INF - (p11) br.cond.spnt FREM_Y_NAN_INF_ZERO + (p9) br.cond.spnt L(FREM_X_NAN_INF) + (p11) br.cond.spnt L(FREM_Y_NAN_INF_ZERO) nop.b 0 } {.mfi nop.m 0 @@ -193,7 +206,7 @@ cmp.eq p11,p10=r29,r0;; nop.i 0;; } -remloop24: +L(remloop24): { .mfi nop.m 0 // Step (2) @@ -215,7 +228,7 @@ remloop24: {.mfi nop.m 0 // q1=q0*(1+e0) - (p6) fma.s1 f15=f12,f7,f12 + fma.s1 f15=f12,f7,f12 nop.i 0 } { .mfi @@ -345,7 +358,7 @@ remloop24: // (p9) set r=r2 (new a, if not last iteration) // (p10) new a =r (p10) mov f13=f6 - (p12) br.cond.sptk remloop24;; + (p12) br.cond.sptk L(remloop24);; } // last iteration @@ -403,7 +416,7 @@ remloop24: -FREM_X_NAN_INF: +L(FREM_X_NAN_INF): // Y zero ? {.mfi @@ -420,19 +433,19 @@ FREM_X_NAN_INF: nop.m 0 nop.i 0 // if Y zero - (p11) br.cond.spnt FREM_Y_ZERO;; + (p11) br.cond.spnt L(FREM_Y_ZERO);; } // X infinity? Return QNAN indefinite { .mfi nop.m 999 - fclass.m.unc p8,p0 = f8, 0x23 +(p0) fclass.m.unc p8,p0 = f8, 0x23 nop.i 999 } // X infinity? Return QNAN indefinite { .mfi nop.m 999 - fclass.m.unc p11,p0 = f8, 0x23 +(p0) fclass.m.unc p11,p0 = f8, 0x23 nop.i 999;; } // Y NaN ? @@ -460,14 +473,14 @@ FREM_X_NAN_INF: } { .mfi nop.m 999 -(p8) fma.s0 f8=f8,f1,f0 +(p8) fma f8=f8,f1,f0 nop.i 0 ;; } { .mfb nop.m 999 frcpa.s0 f8,p7=f8,f9 - (p11) br.cond.spnt EXP_ERROR_RETURN;; + (p11) br.cond.spnt L(EXP_ERROR_RETURN);; } { .mib nop.m 0 @@ -476,24 +489,24 @@ FREM_X_NAN_INF: } -FREM_Y_NAN_INF_ZERO: +L(FREM_Y_NAN_INF_ZERO): // Y INF { .mfi nop.m 999 - fclass.m.unc p7,p0 = f9, 0x23 +(p0) fclass.m.unc p7,p0 = f9, 0x23 nop.i 999 ;; } { .mfb nop.m 999 -(p7) fma.s0 f8=f8,f1,f0 +(p7) fma f8=f8,f1,f0 (p7) br.ret.spnt b0 ;; } // Y NAN? { .mfi nop.m 999 - fclass.m.unc p9,p10 = f9, 0xc3 +(p0) fclass.m.unc p9,p10 = f9, 0xc3 nop.i 999 ;; } { .mfi @@ -504,11 +517,11 @@ FREM_Y_NAN_INF_ZERO: { .mfb nop.m 999 -(p9) fma.s0 f8=f9,f1,f0 +(p9) fma f8=f9,f1,f0 (p9) br.ret.spnt b0 ;; } -FREM_Y_ZERO: +L(FREM_Y_ZERO): // Y zero? Must be zero at this point // because it is the only choice left. // Return QNAN indefinite @@ -516,7 +529,7 @@ FREM_Y_ZERO: // X NAN? { .mfi nop.m 999 - fclass.m.unc p9,p10 = f8, 0xc3 +(p0) fclass.m.unc p9,p10 = f8, 0xc3 nop.i 999 ;; } { .mfi @@ -527,38 +540,43 @@ FREM_Y_ZERO: {.mfi nop.m 999 - (p9) frcpa.s0 f11,p7=f8,f0 + (p9) frcpa f11,p7=f8,f0 nop.i 0;; } { .mfi nop.m 999 -(p10) frcpa.s0 f11,p7 = f0,f0 +(p10) frcpa f11,p7 = f0,f0 nop.i 999;; } { .mfi nop.m 999 - fmerge.s f10 = f8, f8 +(p0) fmerge.s f10 = f8, f8 nop.i 999 } { .mfi nop.m 999 - fma.s0 f8=f11,f1,f0 +(p0) fma f8=f11,f1,f0 nop.i 999;; } -EXP_ERROR_RETURN: +L(EXP_ERROR_RETURN): { .mib - mov GR_Parameter_TAG = 123 +(p0) mov GR_Parameter_TAG = 123 nop.i 999 - br.sptk __libm_error_region;; +(p0) br.sptk __libm_error_region;; } -GLOBAL_IEEE754_END(remainderl) +.endp remainderl +ASM_SIZE_DIRECTIVE(remainderl) +#ifdef _LIBC +ASM_SIZE_DIRECTIVE(__remainderl) +#endif -LOCAL_LIBM_ENTRY(__libm_error_region) +.proc __libm_error_region +__libm_error_region: .prologue { .mfi add GR_Parameter_Y=-32,sp // Parameter 2 value @@ -606,12 +624,9 @@ LOCAL_LIBM_ENTRY(__libm_error_region) br.ret.sptk b0 // Return };; -LOCAL_LIBM_END(__libm_error_region) +.endp __libm_error_region +ASM_SIZE_DIRECTIVE(__libm_error_region) .type __libm_error_support#,@function .global __libm_error_support# - - - - |