diff options
Diffstat (limited to 'sysdeps/ia64/fpu/e_remainderl.S')
| -rw-r--r-- | sysdeps/ia64/fpu/e_remainderl.S | 116 |
1 files changed, 50 insertions, 66 deletions
diff --git a/sysdeps/ia64/fpu/e_remainderl.S b/sysdeps/ia64/fpu/e_remainderl.S index 5856861442..1c1a3c3072 100644 --- a/sysdeps/ia64/fpu/e_remainderl.S +++ b/sysdeps/ia64/fpu/e_remainderl.S @@ -1,10 +1,10 @@ -.file "remainderl.asm" -// Copyright (C) 2000, 2001, Intel Corporation +.file "remainderl.s" + + +// Copyright (c) 2000 - 2003, Intel Corporation // All rights reserved. // -// 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. +// Contributed 2000 by the Intel Numerics Group, 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,17 +35,19 @@ // // Intel Corporation is the author of this code, and requests that all // problem reports or change requests be submitted to it directly at -// http://developer.intel.com/opensource. +// http://www.intel.com/software/products/opensource/libraries/num.htm. // // History //==================================================================== -// 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 +// 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 // set [the previously overwritten] GR_Parameter_RESULT. -//11/29/00 Set FR_Y to f9 +// 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 // // API //==================================================================== @@ -77,9 +79,6 @@ //==================================================================== // a=+/- Inf, or b=+/-0: return NaN, call libm_error_support // a=NaN or b=NaN: return NaN - -#include "libm_support.h" - // // Registers used //==================================================================== @@ -87,8 +86,6 @@ // 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 @@ -105,19 +102,9 @@ 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 @@ -159,7 +146,7 @@ cmp.eq p11,p10=r29,r0;; // X +-NAN, +-inf, ? p9 { .mfi nop.m 999 -(p0) fclass.m.unc p9,p8 = f8, 0xe3 + fclass.m.unc p9,p8 = f8, 0xe3 nop.i 999;; } @@ -196,8 +183,8 @@ cmp.eq p11,p10=r29,r0;; } {.bbb - (p9) br.cond.spnt L(FREM_X_NAN_INF) - (p11) br.cond.spnt L(FREM_Y_NAN_INF_ZERO) + (p9) br.cond.spnt FREM_X_NAN_INF + (p11) br.cond.spnt FREM_Y_NAN_INF_ZERO nop.b 0 } {.mfi nop.m 0 @@ -206,7 +193,7 @@ cmp.eq p11,p10=r29,r0;; nop.i 0;; } -L(remloop24): +remloop24: { .mfi nop.m 0 // Step (2) @@ -228,7 +215,7 @@ L(remloop24): {.mfi nop.m 0 // q1=q0*(1+e0) - fma.s1 f15=f12,f7,f12 + (p6) fma.s1 f15=f12,f7,f12 nop.i 0 } { .mfi @@ -358,7 +345,7 @@ L(remloop24): // (p9) set r=r2 (new a, if not last iteration) // (p10) new a =r (p10) mov f13=f6 - (p12) br.cond.sptk L(remloop24);; + (p12) br.cond.sptk remloop24;; } // last iteration @@ -416,7 +403,7 @@ L(remloop24): -L(FREM_X_NAN_INF): +FREM_X_NAN_INF: // Y zero ? {.mfi @@ -433,19 +420,19 @@ L(FREM_X_NAN_INF): nop.m 0 nop.i 0 // if Y zero - (p11) br.cond.spnt L(FREM_Y_ZERO);; + (p11) br.cond.spnt FREM_Y_ZERO;; } // X infinity? Return QNAN indefinite { .mfi nop.m 999 -(p0) fclass.m.unc p8,p0 = f8, 0x23 + fclass.m.unc p8,p0 = f8, 0x23 nop.i 999 } // X infinity? Return QNAN indefinite { .mfi nop.m 999 -(p0) fclass.m.unc p11,p0 = f8, 0x23 + fclass.m.unc p11,p0 = f8, 0x23 nop.i 999;; } // Y NaN ? @@ -473,14 +460,14 @@ L(FREM_X_NAN_INF): } { .mfi nop.m 999 -(p8) fma f8=f8,f1,f0 +(p8) fma.s0 f8=f8,f1,f0 nop.i 0 ;; } { .mfb nop.m 999 frcpa.s0 f8,p7=f8,f9 - (p11) br.cond.spnt L(EXP_ERROR_RETURN);; + (p11) br.cond.spnt EXP_ERROR_RETURN;; } { .mib nop.m 0 @@ -489,24 +476,24 @@ L(FREM_X_NAN_INF): } -L(FREM_Y_NAN_INF_ZERO): +FREM_Y_NAN_INF_ZERO: // Y INF { .mfi nop.m 999 -(p0) fclass.m.unc p7,p0 = f9, 0x23 + fclass.m.unc p7,p0 = f9, 0x23 nop.i 999 ;; } { .mfb nop.m 999 -(p7) fma f8=f8,f1,f0 +(p7) fma.s0 f8=f8,f1,f0 (p7) br.ret.spnt b0 ;; } // Y NAN? { .mfi nop.m 999 -(p0) fclass.m.unc p9,p10 = f9, 0xc3 + fclass.m.unc p9,p10 = f9, 0xc3 nop.i 999 ;; } { .mfi @@ -517,11 +504,11 @@ L(FREM_Y_NAN_INF_ZERO): { .mfb nop.m 999 -(p9) fma f8=f9,f1,f0 +(p9) fma.s0 f8=f9,f1,f0 (p9) br.ret.spnt b0 ;; } -L(FREM_Y_ZERO): +FREM_Y_ZERO: // Y zero? Must be zero at this point // because it is the only choice left. // Return QNAN indefinite @@ -529,7 +516,7 @@ L(FREM_Y_ZERO): // X NAN? { .mfi nop.m 999 -(p0) fclass.m.unc p9,p10 = f8, 0xc3 + fclass.m.unc p9,p10 = f8, 0xc3 nop.i 999 ;; } { .mfi @@ -540,43 +527,37 @@ L(FREM_Y_ZERO): {.mfi nop.m 999 - (p9) frcpa f11,p7=f8,f0 + (p9) frcpa.s0 f11,p7=f8,f0 nop.i 0;; } { .mfi nop.m 999 -(p10) frcpa f11,p7 = f0,f0 +(p10) frcpa.s0 f11,p7 = f0,f0 nop.i 999;; } { .mfi nop.m 999 -(p0) fmerge.s f10 = f8, f8 + fmerge.s f10 = f8, f8 nop.i 999 } { .mfi nop.m 999 -(p0) fma f8=f11,f1,f0 + fma.s0 f8=f11,f1,f0 nop.i 999;; } -L(EXP_ERROR_RETURN): +EXP_ERROR_RETURN: { .mib -(p0) mov GR_Parameter_TAG = 123 + mov GR_Parameter_TAG = 123 nop.i 999 -(p0) br.sptk __libm_error_region;; + br.sptk __libm_error_region;; } -.endp remainderl -ASM_SIZE_DIRECTIVE(remainderl) -#ifdef _LIBC -ASM_SIZE_DIRECTIVE(__remainderl) -#endif - -.proc __libm_error_region -__libm_error_region: +GLOBAL_IEEE754_END(remainderl) +LOCAL_LIBM_ENTRY(__libm_error_region) .prologue { .mfi add GR_Parameter_Y=-32,sp // Parameter 2 value @@ -624,9 +605,12 @@ __libm_error_region: br.ret.sptk b0 // Return };; -.endp __libm_error_region -ASM_SIZE_DIRECTIVE(__libm_error_region) +LOCAL_LIBM_END(__libm_error_region) .type __libm_error_support#,@function .global __libm_error_support# + + + + |