sysdeps/sparc/sparc32/sparcv9/fpu/multiarch/s_ceil-vis2.S

   1 /* ceil function, sparc32 v9 vis2 version.
   2    Copyright (C) 2013-2014 Free Software Foundation, Inc.
   3    This file is part of the GNU C Library.
   4    Contributed by David S. Miller <davem@davemloft.net>, 2013.
   5
   6    The GNU C Library is free software; you can redistribute it and/or
   7    modify it under the terms of the GNU Lesser General Public
   8    License as published by the Free Software Foundation; either
   9    version 2.1 of the License, or (at your option) any later version.
  10
  11    The GNU C Library is distributed in the hope that it will be useful,
  12    but WITHOUT ANY WARRANTY; without even the implied warranty of
  13    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14    Lesser General Public License for more details.
  15
  16    You should have received a copy of the GNU Lesser General Public
  17    License along with the GNU C Library; if not, see
  18    <http://www.gnu.org/licenses/>.  */
  19
  20 #include <sysdep.h>
  21
  22         /* 'siam' (Set Interval Arithmetic Mode) is used to quickly override
  23            the rounding mode during this routine.
  24
  25            We add then subtract (or subtract than add if the initial
  26            value was negative) 2**23 to the value, then subtract it
  27            back out.
  28
  29            This will clear out the fractional portion of the value and,
  30            with suitable 'siam' initiated rouding mode settings, round
  31            the final result in the proper direction.  */
  32
  33 #define TWO_FIFTYTWO    0x43300000              /* 2**52 */
  34
  35 #define ZERO            %f10                    /* 0.0 */
  36 #define SIGN_BIT        %f12                    /* -0.0 */
  37
  38 ENTRY (__ceil_vis2)
  39         sethi   %hi(TWO_FIFTYTWO), %o2
  40         sllx    %o0, 32, %o0
  41         or      %o0, %o1, %o0
  42         stx     %o0, [%sp + 72]
  43         sllx    %o2, 32, %o2
  44         fzero   ZERO
  45         ldd     [%sp + 72], %f0
  46         fnegd   ZERO, SIGN_BIT
  47         stx     %o2, [%sp + 72]
  48         fabsd   %f0, %f14
  49         ldd     [%sp + 72], %f16
  50         fcmpd   %fcc3, %f14, %f16
  51         fmovduge %fcc3, ZERO, %f16
  52         fand    %f0, SIGN_BIT, SIGN_BIT
  53         for     %f16, SIGN_BIT, %f16
  54         siam    (1 << 2) | 2
  55         faddd   %f0, %f16, %f18
  56         siam    (1 << 2) | 0
  57         fsubd   %f18, %f16, %f18
  58         siam    (0 << 2)
  59         retl
  60          for    %f18, SIGN_BIT, %f0
  61 END (__ceil_vis2)