1 /* Set current rounding direction.
2 Copyright (C) 2000 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
4 Contributed by David Huggins-Daines <dhd@debian.org>, 2000
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.
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.
16 You should have received a copy of the GNU Lesser General Public
17 License along with the GNU C Library; if not, write to the Free
18 Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
24 fesetround (int round
)
26 union { unsigned long long l
; unsigned int sw
[2]; } s
;
28 if (round
& ~FE_DOWNWARD
)
29 /* round is not a valid rounding mode. */
32 /* Get the current status word. */
33 __asm__ ("fstd %%fr0,0(%1)" : "=m" (s
.l
) : "r" (&s
.l
) : "%r0");
34 s
.sw
[0] &= ~FE_DOWNWARD
;
35 s
.sw
[0] |= round
& FE_DOWNWARD
;
36 __asm__ ("fldd 0(%0),%%fr0" : : "r" (&s
.l
), "m" (s
.l
) : "%r0");
41 libm_hidden_def (fesetround
)