OLD | NEW |
| (Empty) |
1 /* mpf_abs -- Compute the absolute value of a float. | |
2 | |
3 Copyright 1993, 1994, 1995, 2001 Free Software Foundation, Inc. | |
4 | |
5 This file is part of the GNU MP Library. | |
6 | |
7 The GNU MP Library is free software; you can redistribute it and/or modify | |
8 it under the terms of the GNU Lesser General Public License as published by | |
9 the Free Software Foundation; either version 3 of the License, or (at your | |
10 option) any later version. | |
11 | |
12 The GNU MP Library is distributed in the hope that it will be useful, but | |
13 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY | |
14 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public | |
15 License for more details. | |
16 | |
17 You should have received a copy of the GNU Lesser General Public License | |
18 along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */ | |
19 | |
20 #include "gmp.h" | |
21 #include "gmp-impl.h" | |
22 | |
23 void | |
24 mpf_abs (mpf_ptr r, mpf_srcptr u) | |
25 { | |
26 mp_size_t size; | |
27 | |
28 size = ABS (u->_mp_size); | |
29 if (r != u) | |
30 { | |
31 mp_size_t prec; | |
32 mp_ptr rp, up; | |
33 | |
34 prec = r->_mp_prec + 1; /* lie not to lose precision in assignment */ | |
35 rp = r->_mp_d; | |
36 up = u->_mp_d; | |
37 | |
38 if (size > prec) | |
39 { | |
40 up += size - prec; | |
41 size = prec; | |
42 } | |
43 | |
44 MPN_COPY (rp, up, size); | |
45 r->_mp_exp = u->_mp_exp; | |
46 } | |
47 r->_mp_size = size; | |
48 } | |
OLD | NEW |