sysdeps/i386/fpu/e_scalbf.S

   1 /*
   2  * Written by J.T. Conklin <jtc@netbsd.org>.
   3  * Public domain.
   4  * Adapted for float type by Ulrich Drepper <drepper@cygnus.com>.
   5  *
   6  * Correct handling of y==-inf <drepper@gnu>
   7  */
   8
   9 #include <machine/asm.h>
  10
  11 RCSID("$NetBSD: $")
  12
  13 #ifdef __ELF__
  14         .section .rodata
  15 #else
  16         .text
  17 #endif
  18
  19         .align ALIGNARG(4)
  20         ASM_TYPE_DIRECTIVE(zero_nan,@object)
  21 zero_nan:
  22         .double 0.0
  23 nan:    .byte 0, 0, 0, 0, 0, 0, 0xff, 0x7f
  24 minus_zero:
  25         .byte 0, 0, 0, 0, 0, 0, 0, 0x80
  26         .byte 0, 0, 0, 0, 0, 0, 0xff, 0x7f
  27         ASM_SIZE_DIRECTIVE(zero_nan)
  28
  29
  30 #ifdef PIC
  31 #define MO(op) op##@GOTOFF(%ecx)
  32 #define MOX(op,x,f) op##@GOTOFF(%ecx,x,f)
  33 #else
  34 #define MO(op) op
  35 #define MOX(op,x,f) op(,x,f)
  36 #endif
  37
  38
  39         .text
  40 ENTRY(__ieee754_scalbf)
  41         flds    8(%esp)
  42         fxam
  43         fnstsw
  44         flds    4(%esp)
  45         andl    $0x4700, %eax
  46         cmpl    $0x0700, %eax
  47         je      1f
  48         andl    $0x4500, %eax
  49         cmpl    $0x0100, %eax
  50         je      2f
  51         fxam
  52         fnstsw
  53         andl    $0x4500, %eax
  54         cmpl    $0x0100, %eax
  55         je      3f
  56         fld     %st(1)
  57         frndint
  58         fcomp   %st(2)
  59         fnstsw
  60         sahf
  61         jne     4f
  62         fscale
  63         fstp    %st(1)
  64         ret
  65
  66         /* y is -inf */
  67 1:      fxam
  68 #ifdef  PIC
  69         LOAD_PIC_REG (cx)
  70 #endif
  71         fnstsw
  72         movl    4(%esp), %edx
  73         shrl    $5, %eax
  74         fstp    %st
  75         fstp    %st
  76         andl    $0x80000000, %edx
  77         andl    $8, %eax
  78         jnz     4f
  79         shrl    $27, %edx
  80         addl    %edx, %eax
  81         fldl    MOX(zero_nan, %eax, 1)
  82         ret
  83
  84         /* The result is NaN, but we must not raise an exception.
  85            So use a variable.  */
  86 2:      fstp    %st
  87         fstp    %st
  88 #ifdef  PIC
  89         LOAD_PIC_REG (cx)
  90 #endif
  91         fldl    MO(nan)
  92         ret
  93
  94         /* The first parameter is a NaN.  Return it.  */
  95 3:      fstp    %st(1)
  96         ret
  97
  98         /* Return NaN and raise the invalid exception.  */
  99 4:      fstp    %st
 100         fstp    %st
 101         fldz
 102         fdiv    %st
 103         ret
 104 END(__ieee754_scalbf)