1 | /* Implement modf for m68k. |
2 | Copyright (C) 1996-2024 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 Lesser General Public |
7 | License as published by the Free Software Foundation; either |
8 | version 2.1 of the 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 | Lesser General Public License for more details. |
14 | |
15 | You should have received a copy of the GNU Lesser General Public |
16 | License along with the GNU C Library. If not, see |
17 | <https://www.gnu.org/licenses/>. */ |
18 | |
19 | #include <math.h> |
20 | #include "mathimpl.h" |
21 | |
22 | FLOAT |
23 | M_DECL_FUNC (__modf) (FLOAT x, FLOAT *iptr) |
24 | { |
25 | FLOAT x_int, result; |
26 | unsigned long x_cond; |
27 | |
28 | __asm ("fintrz%.x %1, %0" : "=f" (x_int) : "f" (x)); |
29 | *iptr = x_int; |
30 | x_cond = __m81_test (val: x); |
31 | if (x_cond & __M81_COND_INF) |
32 | { |
33 | result = 0; |
34 | if (x_cond & __M81_COND_NEG) |
35 | result = -result; |
36 | } |
37 | else if (x_cond & __M81_COND_ZERO) |
38 | result = x; |
39 | else |
40 | result = x - x_int; |
41 | return result; |
42 | } |
43 | declare_mgen_alias (__modf, modf) |
44 | |