sysdeps/powerpc/s_lrint.c

   1 /* Round floating-point to integer.  PowerPC version.
   2    Copyright (C) 1997 Free Software Foundation, Inc.
   3    This file is part of the GNU C Library.
   4
   5    The GNU C Library is free software; you can redistribute it and/or
   6    modify it under the terms of the GNU Library General Public License as
   7    published by the Free Software Foundation; either version 2 of the
   8    License, or (at your option) any later version.
   9
  10    The GNU C Library is distributed in the hope that it will be useful,
  11    but WITHOUT ANY WARRANTY; without even the implied warranty of
  12    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13    Library General Public License for more details.
  14
  15    You should have received a copy of the GNU Library General Public
  16    License along with the GNU C Library; see the file COPYING.LIB.  If not,
  17    write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
  18    Boston, MA 02111-1307, USA.  */
  19
  20 #include "math.h"
  21
  22 long int
  23 __lrint (double x)
  24 {
  25   union {
  26     double d;
  27     long int ll[2];
  28   } u;
  29   asm ("fctiw %0,%1" : "=f"(u.d) : "f"(x));
  30   return u.ll[1];
  31 }
  32 weak_alias (__lrint, lrint)
  33
  34 /* This code will also work for a 'float' argument.  */
  35 asm ("\
  36         .globl __lrintf
  37         .globl lrintf
  38         .weak lrintf
  39         .set __lrintf,__lrint
  40         .set lrintf,__lrint
  41 ");
  42
  43 #ifdef NO_LONG_DOUBLE
  44 strong_alias (__lrint, __lrintl)
  45 weak_alias (__lrint, lrintl)
  46 #endif