rename api functions...
[blender.git] / source / blender / blenlib / intern / math_base_inline.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_inline.c
27  *  \ingroup bli
28  */
29
30
31 #include <float.h>
32 #include <stdio.h>
33 #include <stdlib.h>
34 #include <string.h>
35
36 #include "BLI_math.h"
37
38 #ifndef __MATH_BASE_INLINE_C__
39 #define __MATH_BASE_INLINE_C__
40
41 /* A few small defines. Keep'em local! */
42 #define SMALL_NUMBER  1.e-8f
43
44 MINLINE float sqrt3f(float f)
45 {
46         if (f == 0.0f) return 0.0f;
47         if (f < 0) return (float)(-exp(log(-f) / 3));
48         else return (float)(exp(log(f) / 3));
49 }
50
51 MINLINE double sqrt3d(double d)
52 {
53         if (d == 0.0) return 0;
54         if (d < 0) return -exp(log(-d) / 3);
55         else return exp(log(d) / 3);
56 }
57
58 MINLINE float saacos(float fac)
59 {
60         if (fac <= -1.0f) return (float)M_PI;
61         else if (fac >= 1.0f) return 0.0;
62         else return acosf(fac);
63 }
64
65 MINLINE float saasin(float fac)
66 {
67         if (fac <= -1.0f) return (float)-M_PI / 2.0f;
68         else if (fac >= 1.0f) return (float)M_PI / 2.0f;
69         else return asinf(fac);
70 }
71
72 MINLINE float sasqrt(float fac)
73 {
74         if (fac <= 0.0f) return 0.0f;
75         return sqrtf(fac);
76 }
77
78 MINLINE float saacosf(float fac)
79 {
80         if (fac <= -1.0f) return (float)M_PI;
81         else if (fac >= 1.0f) return 0.0f;
82         else return acosf(fac);
83 }
84
85 MINLINE float saasinf(float fac)
86 {
87         if (fac <= -1.0f) return (float)-M_PI / 2.0f;
88         else if (fac >= 1.0f) return (float)M_PI / 2.0f;
89         else return asinf(fac);
90 }
91
92 MINLINE float sasqrtf(float fac)
93 {
94         if (fac <= 0.0f) return 0.0f;
95         return sqrtf(fac);
96 }
97
98 MINLINE float interpf(float target, float origin, float fac)
99 {
100         return (fac * target) + (1.0f - fac) * origin;
101 }
102
103 /* useful to calculate an even width shell, by taking the angle between 2 planes.
104  * The return value is a scale on the offset.
105  * no angle between planes is 1.0, as the angle between the 2 planes approaches 180d
106  * the distance gets very high, 180d would be inf, but this case isn't valid */
107 MINLINE float shell_angle_to_dist(const float angle)
108 {
109         return (angle < SMALL_NUMBER) ? 1.0f : fabsf(1.0f / cosf(angle));
110 }
111
112 /* used for zoom values*/
113 MINLINE float power_of_2(float val)
114 {
115         return (float)pow(2.0, ceil(log((double)val) / M_LN2));
116 }
117
118 MINLINE int is_power_of_2_i(int n)
119 {
120         return (n & (n - 1)) == 0;
121 }
122
123 MINLINE int power_of_2_max_i(int n)
124 {
125         if (is_power_of_2_i(n))
126                 return n;
127
128         while (!is_power_of_2_i(n))
129                 n = n & (n - 1);
130
131         return n * 2;
132 }
133
134 MINLINE int power_of_2_min_i(int n)
135 {
136         while (!is_power_of_2_i(n))
137                 n = n & (n - 1);
138
139         return n;
140 }
141
142 MINLINE float min_ff(float a, float b)
143 {
144         return (a < b) ? a : b;
145 }
146 MINLINE float max_ff(float a, float b)
147 {
148         return (a > b) ? a : b;
149 }
150
151 MINLINE int min_ii(int a, int b)
152 {
153         return (a < b) ? a : b;
154 }
155 MINLINE int max_ii(int a, int b)
156 {
157         return (b < a) ? a : b;
158 }
159
160 MINLINE float signf(float f)
161 {
162         return (f < 0.f) ? -1.f : 1.f;
163 }
164
165
166 #endif /* __MATH_BASE_INLINE_C__ */