sysdeps/hppa/fpu/fesetround.c

   1 /* Set current rounding direction.
   2    Copyright (C) 2000-2015 Free Software Foundation, Inc.
   3    This file is part of the GNU C Library.
   4    Contributed by David Huggins-Daines <dhd@debian.org>, 2000
   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 <fenv.h>
  21
  22 int
  23 fesetround (int round)
  24 {
  25   union { unsigned long long l; unsigned int sw[2]; } s;
  26
  27   if (round & ~FE_DOWNWARD)
  28     /* round is not a valid rounding mode. */
  29     return 1;
  30
  31   /* Get the current status word. */
  32   __asm__ ("fstd %%fr0,0(%1)" : "=m" (s.l) : "r" (&s.l) : "%r0");
  33   s.sw[0] &= ~FE_DOWNWARD;
  34   s.sw[0] |= round & FE_DOWNWARD;
  35   __asm__ ("fldd 0(%0),%%fr0" : : "r" (&s.l), "m" (s.l) : "%r0");
  36
  37   return 0;
  38 }
  39
  40 libm_hidden_def (fesetround)