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