diff options
Diffstat (limited to 'sysdeps/ia64/fpu/s_modfl.S')
| -rw-r--r-- | sysdeps/ia64/fpu/s_modfl.S | 49 |
1 files changed, 21 insertions, 28 deletions
diff --git a/sysdeps/ia64/fpu/s_modfl.S b/sysdeps/ia64/fpu/s_modfl.S index b5eb509..eaf410c 100644 --- a/sysdeps/ia64/fpu/s_modfl.S +++ b/sysdeps/ia64/fpu/s_modfl.S @@ -1,10 +1,10 @@ .file "modfl.s" -// Copyright (C) 2000, 2001, Intel Corporation + +// Copyright (c) 2000 - 2003, Intel Corporation // All rights reserved. -// -// Contributed 2/2/2000 by John Harrison, 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,15 +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://developer.intel.com/opensource. +// http://www.intel.com/software/products/opensource/libraries/num.htm. // // History //============================================================== -// 2/02/00: Initial version -// 4/04/00: Improved speed, corrected result for NaN input -// 5/30/00 Fixed bug for exponent 0x1003e +// 02/02/00 Initial version +// 04/04/00 Improved speed, corrected result for NaN input +// 05/30/00 Fixed bug for exponent 0x1003e // 12/22/00 Fixed so inexact flag is never set, and invalid is not set for // qnans nor for inputs larger than 2^63. +// 05/20/02 Cleaned up namespace and sf0 syntax +// 02/10/03 Reordered header: .section, .global, .proc, .align // // API //============================================================== @@ -92,8 +94,6 @@ // p13 --------------------------------------------------->| // -#include "libm_support.h" - // floating-point registers used: MODF_NORM_F8 = f9 MODF_FRACTION_PART = f10 @@ -110,23 +110,17 @@ modf_exp = r18 // r34 = iptr -.align 32 -.global modfl# - .section .text -.proc modfl# -.align 32 - +GLOBAL_LIBM_ENTRY(modfl) // Main path is p9, p11, p8 FALSE and p12 TRUE // Assume input is normalized and get signexp // Normalize input just in case // Form exponent bias -modfl: { .mfi getf.exp modf_signexp = f8 - fnorm MODF_NORM_F8 = f8 + fnorm.s0 MODF_NORM_F8 = f8 addl modf_GR_FFFF = 0xffff, r0 } // Get integer part of input @@ -171,10 +165,10 @@ modfl: { .mfb (p10) cmp.ge.unc p9,p12 = modf_exp, modf_GR_no_frac (p6) fclass.m.unc p6,p7 = f8, 0x23 -(p8) br.cond.spnt L(MODF_DENORM) ;; +(p8) br.cond.spnt MODF_DENORM ;; } -L(MODF_COMMON): +MODF_COMMON: // For HUGE set fraction to signed 0 { .mfi nop.m 999 @@ -184,7 +178,7 @@ L(MODF_COMMON): // For HUGE set integer part to normalized input { .mfi nop.m 999 -(p9) fnorm MODF_INTEGER_PART = MODF_NORM_F8 +(p9) fnorm.s0 MODF_INTEGER_PART = MODF_NORM_F8 nop.i 999 ;; } @@ -196,7 +190,7 @@ L(MODF_COMMON): } { .mfi nop.m 999 -(p11) fnorm f8 = MODF_NORM_F8 +(p11) fnorm.s0 f8 = MODF_NORM_F8 nop.i 999 ;; } @@ -237,7 +231,7 @@ L(MODF_COMMON): // For NORMAL test if fraction part is zero; if so append correct sign { .mfi nop.m 999 -(p12) fcmp.eq.unc p7,p0 = MODF_NORM_F8, MODF_INTEGER_PART +(p12) fcmp.eq.unc.s0 p7,p0 = MODF_NORM_F8, MODF_INTEGER_PART nop.i 999 ;; } @@ -254,7 +248,7 @@ L(MODF_COMMON): br.ret.sptk b0 ;; } -L(MODF_DENORM): +MODF_DENORM: // If x unorm get signexp from normalized input // If x unorm get integer part from normalized input { .mfi @@ -273,8 +267,7 @@ L(MODF_DENORM): { .mfb (p10) cmp.ge.unc p9,p12 = modf_exp, modf_GR_no_frac nop.f 999 - br.cond.spnt L(MODF_COMMON) ;; + br.cond.spnt MODF_COMMON ;; } -.endp modfl -ASM_SIZE_DIRECTIVE(modfl) +GLOBAL_LIBM_END(modfl) |