Merged changes in the trunk up to revision 45619.
[blender.git] / source / blender / blenlib / intern / math_base.c
1 /*
2  * ***** BEGIN GPL LICENSE BLOCK *****
3  *
4  * This program is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU General Public License
6  * as published by the Free Software Foundation; either version 2
7  * of the License, or (at your option) any later version.
8  *
9  * This program is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12  * GNU General Public License for more details.
13  *
14  * You should have received a copy of the GNU General Public License
15  * along with this program; if not, write to the Free Software Foundation,
16  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  *
18  * The Original Code is Copyright (C) 2001-2002 by NaN Holding BV.
19  * All rights reserved.
20  *
21  * The Original Code is: some of this file.
22  *
23  * ***** END GPL LICENSE BLOCK *****
24  * */
25
26 /** \file blender/blenlib/intern/math_base.c
27  *  \ingroup bli
28  */
29
30
31
32 #include "BLI_math.h"
33
34 /* WARNING: MSVC compiling hack for double_round() */
35 #if (defined(WIN32) || defined(WIN64)) && !(defined(FREE_WINDOWS))
36
37 /* from python 3.1 pymath.c */
38 double copysign(double x, double y)
39 {
40         /* use atan2 to distinguish -0.0 from 0.0 */
41         if (y > 0.0 || (y == 0.0 && atan2(y, -1.0) > 0.0)) {
42                 return fabs(x);
43         }
44         else {
45                 return -fabs(x);
46         }
47 }
48
49 /* from python 3.1 pymath.c */
50 double round(double x)
51 {
52         double absx, y;
53         absx = fabs(x);
54         y = floor(absx);
55         if (absx - y >= 0.5)
56                 y += 1.0;
57         return copysign(y, x);
58 }
59 #else /* OpenSuse 11.1 seems to need this. */
60 double round(double x);
61 #endif
62
63 /* from python 3.1 floatobject.c
64  * ndigits must be between 0 and 21 */
65 double double_round(double x, int ndigits)
66 {
67         double pow1, pow2, y, z;
68         if (ndigits >= 0) {
69                 pow1 = pow(10.0, (double)ndigits);
70                 pow2 = 1.0;
71                 y = (x * pow1) * pow2;
72                 /* if y overflows, then rounded value is exactly x */
73                 if (!finite(y))
74                         return x;
75         }
76         else {
77                 pow1 = pow(10.0, (double)-ndigits);
78                 pow2 = 1.0; /* unused; silences a gcc compiler warning */
79                 y = x / pow1;
80         }
81
82         z = round(y);
83         if (fabs(y - z) == 0.5)
84                 /* halfway between two integers; use round-half-even */
85                 z = 2.0 * round(y / 2.0);
86
87         if (ndigits >= 0)
88                 z = (z / pow2) / pow1;
89         else
90                 z *= pow1;
91
92         /* if computation resulted in overflow, raise OverflowError */
93         return z;
94 }