8eb6561003b060509976eb1c584a5e83130a8905
[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 #ifndef __MATH_BASE_INLINE_C__
31 #define __MATH_BASE_INLINE_C__
32
33 #include <float.h>
34 #include <stdio.h>
35 #include <stdlib.h>
36 #include <string.h>
37
38 #include "BLI_math.h"
39
40 /* A few small defines. Keep'em local! */
41 #define SMALL_NUMBER  1.e-8f
42
43 MINLINE float sqrt3f(float f)
44 {
45         if (f == 0.0f) return 0.0f;
46         if (f < 0) return (float)(-exp(log(-f) / 3));
47         else return (float)(exp(log(f) / 3));
48 }
49
50 MINLINE double sqrt3d(double d)
51 {
52         if (d == 0.0) return 0;
53         if (d < 0) return -exp(log(-d) / 3);
54         else return exp(log(d) / 3);
55 }
56
57 MINLINE float saacos(float fac)
58 {
59         if (fac <= -1.0f) return (float)M_PI;
60         else if (fac >= 1.0f) return 0.0;
61         else return acosf(fac);
62 }
63
64 MINLINE float saasin(float fac)
65 {
66         if (fac <= -1.0f) return (float)-M_PI / 2.0f;
67         else if (fac >= 1.0f) return (float)M_PI / 2.0f;
68         else return asinf(fac);
69 }
70
71 MINLINE float sasqrt(float fac)
72 {
73         if (fac <= 0.0f) return 0.0f;
74         return sqrtf(fac);
75 }
76
77 MINLINE float saacosf(float fac)
78 {
79         if (fac <= -1.0f) return (float)M_PI;
80         else if (fac >= 1.0f) return 0.0f;
81         else return acosf(fac);
82 }
83
84 MINLINE float saasinf(float fac)
85 {
86         if (fac <= -1.0f) return (float)-M_PI / 2.0f;
87         else if (fac >= 1.0f) return (float)M_PI / 2.0f;
88         else return asinf(fac);
89 }
90
91 MINLINE float sasqrtf(float fac)
92 {
93         if (fac <= 0.0f) return 0.0f;
94         return sqrtf(fac);
95 }
96
97 MINLINE float interpf(float target, float origin, float fac)
98 {
99         return (fac * target) + (1.0f - fac) * origin;
100 }
101
102 /* useful to calculate an even width shell, by taking the angle between 2 planes.
103  * The return value is a scale on the offset.
104  * no angle between planes is 1.0, as the angle between the 2 planes approaches 180d
105  * the distance gets very high, 180d would be inf, but this case isn't valid */
106 MINLINE float shell_angle_to_dist(const float angle)
107 {
108         return (angle < SMALL_NUMBER) ? 1.0f : fabsf(1.0f / cosf(angle));
109 }
110
111 /* used for zoom values*/
112 MINLINE float power_of_2(float val)
113 {
114         return (float)pow(2.0, ceil(log((double)val) / M_LN2));
115 }
116
117 MINLINE int is_power_of_2_i(int n)
118 {
119         return (n & (n - 1)) == 0;
120 }
121
122 MINLINE int power_of_2_max_i(int n)
123 {
124         if (is_power_of_2_i(n))
125                 return n;
126
127         while (!is_power_of_2_i(n))
128                 n = n & (n - 1);
129
130         return n * 2;
131 }
132
133 MINLINE int power_of_2_min_i(int n)
134 {
135         while (!is_power_of_2_i(n))
136                 n = n & (n - 1);
137
138         return n;
139 }
140
141 MINLINE unsigned int highest_order_bit_i(unsigned int n)
142 {
143         n |= (n >>  1);
144         n |= (n >>  2);
145         n |= (n >>  4);
146         n |= (n >>  8);
147         n |= (n >> 16);
148         return n - (n >> 1);
149 }
150
151 MINLINE unsigned short highest_order_bit_s(unsigned short n)
152 {
153         n |= (n >>  1);
154         n |= (n >>  2);
155         n |= (n >>  4);
156         n |= (n >>  8);
157         return n - (n >> 1);
158 }
159
160 MINLINE float min_ff(float a, float b)
161 {
162         return (a < b) ? a : b;
163 }
164 MINLINE float max_ff(float a, float b)
165 {
166         return (a > b) ? a : b;
167 }
168
169 MINLINE int min_ii(int a, int b)
170 {
171         return (a < b) ? a : b;
172 }
173 MINLINE int max_ii(int a, int b)
174 {
175         return (b < a) ? a : b;
176 }
177
178 MINLINE float min_fff(float a, float b, float c)
179 {
180         return min_ff(min_ff(a, b), c);
181 }
182 MINLINE float max_fff(float a, float b, float c)
183 {
184         return max_ff(max_ff(a, b), c);
185 }
186
187 MINLINE int min_iii(int a, int b, int c)
188 {
189         return min_ii(min_ii(a, b), c);
190 }
191 MINLINE int max_iii(int a, int b, int c)
192 {
193         return max_ii(max_ii(a, b), c);
194 }
195
196 MINLINE float min_ffff(float a, float b, float c, float d)
197 {
198         return min_ff(min_fff(a, b, c), d);
199 }
200 MINLINE float max_ffff(float a, float b, float c, float d)
201 {
202         return max_ff(max_fff(a, b, c), d);
203 }
204
205 MINLINE int min_iiii(int a, int b, int c, int d)
206 {
207         return min_ii(min_iii(a, b, c), d);
208 }
209 MINLINE int max_iiii(int a, int b, int c, int d)
210 {
211         return max_ii(max_iii(a, b, c), d);
212 }
213
214 MINLINE float signf(float f)
215 {
216         return (f < 0.f) ? -1.f : 1.f;
217 }
218
219
220 #endif /* __MATH_BASE_INLINE_C__ */