01c00ddb48e626325370cc678692b5a565387c52
[blender.git] / source / blender / blenlib / BLI_math_vector.h
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 #ifndef __BLI_MATH_VECTOR_H__
27 #define __BLI_MATH_VECTOR_H__
28
29 /** \file BLI_math_vector.h
30  *  \ingroup bli
31  */
32
33 #ifdef __cplusplus
34 extern "C" {
35 #endif
36
37 #include "BLI_compiler_attrs.h"
38 #include "BLI_math_inline.h"
39
40 /************************************* Init ***********************************/
41
42 #ifdef BLI_MATH_GCC_WARN_PRAGMA
43 #  pragma GCC diagnostic push
44 #  pragma GCC diagnostic ignored "-Wredundant-decls"
45 #endif
46
47
48 MINLINE void zero_v2(float r[2]);
49 MINLINE void zero_v3(float r[3]);
50 MINLINE void zero_v4(float r[4]);
51
52 MINLINE void copy_v2_v2(float r[2], const float a[2]);
53 MINLINE void copy_v3_v3(float r[3], const float a[3]);
54 MINLINE void copy_v4_v4(float r[4], const float a[4]);
55
56 MINLINE void copy_v2_fl(float r[2], float f);
57 MINLINE void copy_v3_fl(float r[3], float f);
58 MINLINE void copy_v4_fl(float r[4], float f);
59
60 MINLINE void swap_v2_v2(float a[2], float b[2]);
61 MINLINE void swap_v3_v3(float a[3], float b[3]);
62 MINLINE void swap_v4_v4(float a[4], float b[4]);
63
64 /* unsigned char */
65 MINLINE void copy_v2_v2_uchar(unsigned char r[2], const unsigned char a[2]);
66 MINLINE void copy_v3_v3_uchar(unsigned char r[3], const unsigned char a[3]);
67 MINLINE void copy_v4_v4_uchar(unsigned char r[4], const unsigned char a[4]);
68 /* char */
69 MINLINE void copy_v2_v2_char(char r[2], const char a[2]);
70 MINLINE void copy_v3_v3_char(char r[3], const char a[3]);
71 MINLINE void copy_v4_v4_char(char r[4], const char a[4]);
72 /* short */
73 MINLINE void copy_v2_v2_short(short r[2], const short a[2]);
74 MINLINE void copy_v3_v3_short(short r[3], const short a[3]);
75 MINLINE void copy_v4_v4_short(short r[4], const short a[4]);
76 /* int */
77 MINLINE void zero_v3_int(int r[3]);
78 MINLINE void copy_v2_v2_int(int r[2], const int a[2]);
79 MINLINE void copy_v3_v3_int(int r[3], const int a[3]);
80 MINLINE void copy_v4_v4_int(int r[4], const int a[4]);
81 /* double -> float */
82 MINLINE void copy_v2fl_v2db(float r[2], const double a[2]);
83 MINLINE void copy_v3fl_v3db(float r[3], const double a[3]);
84 MINLINE void copy_v4fl_v4db(float r[4], const double a[4]);
85 /* float -> double */
86 MINLINE void copy_v2db_v2fl(double r[2], const float a[2]);
87 MINLINE void copy_v3db_v3fl(double r[3], const float a[3]);
88 MINLINE void copy_v4db_v4fl(double r[4], const float a[4]);
89 /* float args -> vec */
90 MINLINE void copy_v2_fl2(float v[2], float x, float y);
91 MINLINE void copy_v3_fl3(float v[3], float x, float y, float z);
92 MINLINE void copy_v4_fl4(float v[4], float x, float y, float z, float w);
93
94 /********************************* Arithmetic ********************************/
95
96 MINLINE void add_v2_fl(float r[2], float f);
97 MINLINE void add_v3_fl(float r[3], float f);
98 MINLINE void add_v4_fl(float r[4], float f);
99 MINLINE void add_v2_v2(float r[2], const float a[2]);
100 MINLINE void add_v2_v2v2(float r[2], const float a[2], const float b[2]);
101 MINLINE void add_v2_v2v2_int(int r[2], const int a[2], const int b[2]);
102 MINLINE void add_v3_v3(float r[3], const float a[3]);
103 MINLINE void add_v3_v3v3(float r[3], const float a[3], const float b[3]);
104 MINLINE void add_v4_v4(float r[4], const float a[4]);
105 MINLINE void add_v4_v4v4(float r[4], const float a[4], const float b[4]);
106
107 MINLINE void sub_v2_v2(float r[2], const float a[2]);
108 MINLINE void sub_v2_v2v2(float r[2], const float a[2], const float b[2]);
109 MINLINE void sub_v2_v2v2_int(int r[2], const int a[2], const int b[2]);
110 MINLINE void sub_v3_v3(float r[3], const float a[3]);
111 MINLINE void sub_v3_v3v3(float r[3], const float a[3], const float b[3]);
112 MINLINE void sub_v4_v4(float r[4], const float a[4]);
113 MINLINE void sub_v4_v4v4(float r[4], const float a[4], const float b[4]);
114
115 MINLINE void mul_v2_fl(float r[2], float f);
116 MINLINE void mul_v2_v2fl(float r[2], const float a[2], float f);
117 MINLINE void mul_v3_fl(float r[3], float f);
118 MINLINE void mul_v3_v3fl(float r[3], const float a[3], float f);
119 MINLINE void mul_v2_v2(float r[2], const float a[2]);
120 MINLINE void mul_v3_v3(float r[3], const float a[3]);
121 MINLINE void mul_v3_v3v3(float r[3], const float a[3], const float b[3]);
122 MINLINE void mul_v4_fl(float r[4], float f);
123 MINLINE void mul_v4_v4fl(float r[3], const float a[3], float f);
124 MINLINE void mul_v2_v2_cw(float r[2], const float mat[2], const float vec[2]);
125 MINLINE void mul_v2_v2_ccw(float r[2], const float mat[2], const float vec[2]);
126 MINLINE float mul_project_m4_v3_zfac(float mat[4][4], const float co[3]) ATTR_WARN_UNUSED_RESULT;
127 MINLINE float dot_m3_v3_row_x(float M[3][3], const float a[3]) ATTR_WARN_UNUSED_RESULT;
128 MINLINE float dot_m3_v3_row_y(float M[3][3], const float a[3]) ATTR_WARN_UNUSED_RESULT;
129 MINLINE float dot_m3_v3_row_z(float M[3][3], const float a[3]) ATTR_WARN_UNUSED_RESULT;
130 MINLINE float dot_m4_v3_row_x(float M[4][4], const float a[3]) ATTR_WARN_UNUSED_RESULT;
131 MINLINE float dot_m4_v3_row_y(float M[4][4], const float a[3]) ATTR_WARN_UNUSED_RESULT;
132 MINLINE float dot_m4_v3_row_z(float M[4][4], const float a[3]) ATTR_WARN_UNUSED_RESULT;
133
134 MINLINE void madd_v2_v2fl(float r[2], const float a[2], float f);
135 MINLINE void madd_v3_v3fl(float r[3], const float a[3], float f);
136 MINLINE void madd_v3_v3v3(float r[3], const float a[3], const float b[3]);
137 MINLINE void madd_v2_v2v2fl(float r[2], const float a[2], const float b[2], float f);
138 MINLINE void madd_v3_v3v3fl(float r[3], const float a[3], const float b[3], float f);
139 MINLINE void madd_v3_v3v3v3(float r[3], const float a[3], const float b[3], const float c[3]);
140 MINLINE void madd_v4_v4fl(float r[4], const float a[4], float f);
141 MINLINE void madd_v4_v4v4(float r[4], const float a[4], const float b[4]);
142
143 MINLINE void negate_v2(float r[2]);
144 MINLINE void negate_v2_v2(float r[2], const float a[2]);
145 MINLINE void negate_v3(float r[3]);
146 MINLINE void negate_v3_v3(float r[3], const float a[3]);
147 MINLINE void negate_v4(float r[4]);
148 MINLINE void negate_v4_v4(float r[4], const float a[3]);
149
150 MINLINE void negate_v3_short(short r[3]);
151
152 MINLINE void abs_v2(float r[2]);
153 MINLINE void abs_v2_v2(float r[2], const float a[2]);
154 MINLINE void abs_v3(float r[3]);
155 MINLINE void abs_v3_v3(float r[3], const float a[3]);
156 MINLINE void abs_v4(float r[4]);
157 MINLINE void abs_v4_v4(float r[4], const float a[4]);
158
159 MINLINE float dot_v2v2(const float a[2], const float b[2]) ATTR_WARN_UNUSED_RESULT;
160 MINLINE float dot_v3v3(const float a[3], const float b[3]) ATTR_WARN_UNUSED_RESULT;
161 MINLINE float dot_v3v3v3(const float p[3], const float a[3], const float b[3]) ATTR_WARN_UNUSED_RESULT;
162 MINLINE float dot_v4v4(const float a[4], const float b[4]) ATTR_WARN_UNUSED_RESULT;
163
164 MINLINE double dot_v3db_v3fl(const double a[3], const float b[3]) ATTR_WARN_UNUSED_RESULT;
165
166 MINLINE float cross_v2v2(const float a[2], const float b[2]) ATTR_WARN_UNUSED_RESULT;
167 MINLINE void cross_v3_v3v3(float r[3], const float a[3], const float b[3]);
168
169 MINLINE void add_newell_cross_v3_v3v3(float n[3], const float v_prev[3], const float v_curr[3]);
170
171 MINLINE void star_m3_v3(float rmat[3][3], float a[3]);
172
173 /*********************************** Length **********************************/
174
175 MINLINE float len_squared_v2(const float v[2]) ATTR_WARN_UNUSED_RESULT;
176 MINLINE float len_squared_v3(const float v[3]) ATTR_WARN_UNUSED_RESULT;
177 MINLINE float len_manhattan_v2(const float v[2]) ATTR_WARN_UNUSED_RESULT;
178 MINLINE int   len_manhattan_v2_int(const int v[2]) ATTR_WARN_UNUSED_RESULT;
179 MINLINE float len_manhattan_v3(const float v[3]) ATTR_WARN_UNUSED_RESULT;
180 MINLINE float len_v2(const float a[2]) ATTR_WARN_UNUSED_RESULT;
181 MINLINE float len_v2v2(const float a[2], const float b[2]) ATTR_WARN_UNUSED_RESULT;
182 MINLINE float len_v2v2_int(const int v1[2], const int v2[2]);
183 MINLINE float len_squared_v2v2(const float a[2], const float b[2]) ATTR_WARN_UNUSED_RESULT;
184 MINLINE float len_squared_v3v3(const float a[3], const float b[3]) ATTR_WARN_UNUSED_RESULT;
185 MINLINE float len_manhattan_v2v2(const float a[2], const float b[2]) ATTR_WARN_UNUSED_RESULT;
186 MINLINE int   len_manhattan_v2v2_int(const int a[2], const int b[2]) ATTR_WARN_UNUSED_RESULT;
187 MINLINE float len_manhattan_v3v3(const float a[3], const float b[3]) ATTR_WARN_UNUSED_RESULT;
188 MINLINE float len_v3(const float a[3]) ATTR_WARN_UNUSED_RESULT;
189 MINLINE float len_v3v3(const float a[3], const float b[3]) ATTR_WARN_UNUSED_RESULT;
190
191 MINLINE float normalize_v2(float r[2]);
192 MINLINE float normalize_v2_v2(float r[2], const float a[2]);
193 MINLINE float normalize_v3(float r[3]);
194 MINLINE float normalize_v3_v3(float r[3], const float a[3]);
195 MINLINE double normalize_v3_d(double n[3]);
196
197 /******************************* Interpolation *******************************/
198
199 void interp_v2_v2v2(float r[2], const float a[2], const float b[2], const float t);
200 void interp_v2_v2v2v2(float r[2], const float a[2], const float b[2], const float c[2], const float t[3]);
201 void interp_v3_v3v3(float r[3], const float a[3], const float b[3], const float t);
202 void interp_v3_v3v3v3(float p[3], const float v1[3], const float v2[3], const float v3[3], const float w[3]);
203 void interp_v3_v3v3v3v3(float p[3], const float v1[3], const float v2[3], const float v3[3], const float v4[3], const float w[4]);
204 void interp_v4_v4v4(float r[4], const float a[4], const float b[4], const float t);
205 void interp_v4_v4v4v4(float p[4], const float v1[4], const float v2[4], const float v3[4], const float w[3]);
206 void interp_v4_v4v4v4v4(float p[4], const float v1[4], const float v2[4], const float v3[4], const float v4[4], const float w[4]);
207 void interp_v3_v3v3v3_uv(float p[3], const float v1[3], const float v2[3], const float v3[3], const float uv[2]);
208
209 bool interp_v3_v3v3_slerp(float target[3], const float a[3], const float b[3], const float t)  ATTR_WARN_UNUSED_RESULT;
210 bool interp_v2_v2v2_slerp(float target[2], const float a[2], const float b[2], const float t)  ATTR_WARN_UNUSED_RESULT;
211
212 void interp_v3_v3v3_slerp_safe(float target[3], const float a[3], const float b[3], const float t);
213 void interp_v2_v2v2_slerp_safe(float target[2], const float a[2], const float b[2], const float t);
214
215 void interp_v3_v3v3_char(char target[3], const char a[3], const char b[3], const float t);
216 void interp_v3_v3v3_uchar(unsigned char target[3], const unsigned char a[3], const unsigned char b[3], const float t);
217 void interp_v4_v4v4_char(char target[4], const char a[4], const char b[4], const float t);
218 void interp_v4_v4v4_uchar(unsigned char target[4], const unsigned char a[4], const unsigned char b[4], const float t);
219
220 void mid_v3_v3v3(float r[3], const float a[3], const float b[3]);
221 void mid_v2_v2v2(float r[2], const float a[2], const float b[2]);
222 void mid_v3_v3v3v3(float v[3], const float v1[3], const float v2[3], const float v3[3]);
223 void mid_v3_v3v3v3v3(float v[3], const float v1[3], const float v2[3], const float v3[3], const float v4[3]);
224
225 void mid_v3_v3v3_angle_weighted(float r[3], const float a[3], const float b[3]);
226 void mid_v3_angle_weighted(float r[3]);
227
228 void flip_v4_v4v4(float v[4], const float v1[4], const float v2[4]);
229 void flip_v3_v3v3(float v[3], const float v1[3], const float v2[3]);
230 void flip_v2_v2v2(float v[2], const float v1[2], const float v2[2]);
231
232 /********************************* Comparison ********************************/
233
234 MINLINE bool is_zero_v2(const float a[3])  ATTR_WARN_UNUSED_RESULT;
235 MINLINE bool is_zero_v3(const float a[3])  ATTR_WARN_UNUSED_RESULT;
236 MINLINE bool is_zero_v4(const float a[4])  ATTR_WARN_UNUSED_RESULT;
237
238 MINLINE bool is_finite_v2(const float a[3])  ATTR_WARN_UNUSED_RESULT;
239 MINLINE bool is_finite_v3(const float a[3])  ATTR_WARN_UNUSED_RESULT;
240 MINLINE bool is_finite_v4(const float a[4])  ATTR_WARN_UNUSED_RESULT;
241
242 MINLINE bool is_one_v3(const float a[3])  ATTR_WARN_UNUSED_RESULT;
243
244 MINLINE bool equals_v2v2(const float v1[2], const float v2[2])  ATTR_WARN_UNUSED_RESULT;
245 MINLINE bool equals_v3v3(const float a[3], const float b[3])  ATTR_WARN_UNUSED_RESULT;
246 MINLINE bool equals_v4v4(const float a[4], const float b[4])  ATTR_WARN_UNUSED_RESULT;
247
248 MINLINE bool compare_v2v2(const float a[2], const float b[2], const float limit)  ATTR_WARN_UNUSED_RESULT;
249 MINLINE bool compare_v3v3(const float a[3], const float b[3], const float limit)  ATTR_WARN_UNUSED_RESULT;
250 MINLINE bool compare_v4v4(const float a[4], const float b[4], const float limit)  ATTR_WARN_UNUSED_RESULT;
251
252 MINLINE bool compare_v2v2_relative(const float a[2], const float b[2], const float limit, const int max_ulps)  ATTR_WARN_UNUSED_RESULT;
253 MINLINE bool compare_v3v3_relative(const float a[3], const float b[3], const float limit, const int max_ulps)  ATTR_WARN_UNUSED_RESULT;
254 MINLINE bool compare_v4v4_relative(const float a[4], const float b[4], const float limit, const int max_ulps)  ATTR_WARN_UNUSED_RESULT;
255
256 MINLINE bool compare_len_v3v3(const float a[3], const float b[3], const float limit)  ATTR_WARN_UNUSED_RESULT;
257 MINLINE bool compare_len_squared_v3v3(const float a[3], const float b[3], const float limit)  ATTR_WARN_UNUSED_RESULT;
258
259 MINLINE float line_point_side_v2(const float l1[2], const float l2[2], const float pt[2]) ATTR_WARN_UNUSED_RESULT;
260
261 /********************************** Angles ***********************************/
262 /* - angle with 2 arguments is angle between vector                          */
263 /* - angle with 3 arguments is angle between 3 points at the middle point    */
264 /* - angle_normalized_* is faster equivalent if vectors are normalized       */
265
266 float angle_v2v2(const float a[2], const float b[2]) ATTR_WARN_UNUSED_RESULT;
267 float angle_signed_v2v2(const float v1[2], const float v2[2]) ATTR_WARN_UNUSED_RESULT;
268 float angle_v2v2v2(const float a[2], const float b[2], const float c[2]) ATTR_WARN_UNUSED_RESULT;
269 float angle_normalized_v2v2(const float a[2], const float b[2]) ATTR_WARN_UNUSED_RESULT;
270 float angle_v3v3(const float a[3], const float b[3]) ATTR_WARN_UNUSED_RESULT;
271 float angle_v3v3v3(const float a[3], const float b[3], const float c[3]) ATTR_WARN_UNUSED_RESULT;
272 float cos_v3v3v3(const float p1[3], const float p2[3], const float p3[3]) ATTR_WARN_UNUSED_RESULT;
273 float angle_normalized_v3v3(const float v1[3], const float v2[3]) ATTR_WARN_UNUSED_RESULT;
274 float angle_on_axis_v3v3v3_v3(const float v1[3], const float v2[3], const float v3[3], const float axis[3]) ATTR_WARN_UNUSED_RESULT;
275 float angle_signed_on_axis_v3v3v3_v3(const float v1[3], const float v2[3], const float v3[3], const float axis[3]) ATTR_WARN_UNUSED_RESULT;
276 void angle_tri_v3(float angles[3], const float v1[3], const float v2[3], const float v3[3]);
277 void angle_quad_v3(float angles[4], const float v1[3], const float v2[3], const float v3[3], const float v4[3]);
278 void angle_poly_v3(float *angles, const float *verts[3], int len);
279
280 /********************************* Geometry **********************************/
281
282 void project_v2_v2v2(float c[2], const float v1[2], const float v2[2]);
283 void project_v3_v3v3(float r[3], const float p[3], const float n[3]);
284 void project_plane_v3_v3v3(float c[3], const float v[3], const float v_plane[3]);
285 void project_plane_v2_v2v2(float c[2], const float v[2], const float v_plane[2]);
286 void project_v3_plane(float v[3], const float n[3], const float p[3]);
287 void reflect_v3_v3v3(float r[3], const float v[3], const float n[3]);
288 void ortho_basis_v3v3_v3(float r_n1[3], float r_n2[3], const float n[3]);
289 void ortho_v3_v3(float p[3], const float v[3]);
290 void ortho_v2_v2(float p[2], const float v[2]);
291 void bisect_v3_v3v3v3(float r[3], const float a[3], const float b[3], const float c[3]);
292 void rotate_v3_v3v3fl(float v[3], const float p[3], const float axis[3], const float angle);
293 void rotate_normalized_v3_v3v3fl(float v[3], const float p[3], const float axis[3], const float angle);
294
295 /*********************************** Other ***********************************/
296
297 void print_v2(const char *str, const float a[2]);
298 void print_v3(const char *str, const float a[3]);
299 void print_v4(const char *str, const float a[4]);
300 void print_vn(const char *str, const float v[], const int n);
301
302 #define print_v2_id(v) print_v2(STRINGIFY(v), v)
303 #define print_v3_id(v) print_v3(STRINGIFY(v), v)
304 #define print_v4_id(v) print_v4(STRINGIFY(v), v)
305 #define print_vn_id(v, n) print_vn(STRINGIFY(v), v, n)
306
307 MINLINE void normal_short_to_float_v3(float r[3], const short n[3]);
308 MINLINE void normal_float_to_short_v3(short r[3], const float n[3]);
309
310 void minmax_v3v3_v3(float min[3], float max[3], const float vec[3]);
311 void minmax_v2v2_v2(float min[2], float max[2], const float vec[2]);
312
313 void minmax_v3v3_v3_array(float r_min[3], float r_max[3], const float (*vec_arr)[3], int nbr);
314
315 void dist_ensure_v3_v3fl(float v1[3], const float v2[3], const float dist);
316 void dist_ensure_v2_v2fl(float v1[2], const float v2[2], const float dist);
317
318 void axis_sort_v3(const float axis_values[3], int r_axis_order[3]);
319
320 /***************************** Array Functions *******************************/
321 /* follow fixed length vector function conventions. */
322 double dot_vn_vn(const float *array_src_a, const float *array_src_b, const int size) ATTR_WARN_UNUSED_RESULT;
323 double len_squared_vn(const float *array, const int size) ATTR_WARN_UNUSED_RESULT;
324 float normalize_vn_vn(float *array_tar, const float *array_src, const int size);
325 float normalize_vn(float *array_tar, const int size);
326 void range_vn_i(int *array_tar, const int size, const int start);
327 void range_vn_u(unsigned int *array_tar, const int size, const unsigned int start);
328 void range_vn_fl(float *array_tar, const int size, const float start, const float step);
329 void negate_vn(float *array_tar, const int size);
330 void negate_vn_vn(float *array_tar, const float *array_src, const int size);
331 void mul_vn_fl(float *array_tar, const int size, const float f);
332 void mul_vn_vn_fl(float *array_tar, const float *array_src, const int size, const float f);
333 void add_vn_vn(float *array_tar, const float *array_src, const int size);
334 void add_vn_vnvn(float *array_tar, const float *array_src_a, const float *array_src_b, const int size);
335 void madd_vn_vn(float *array_tar, const float *array_src, const float f, const int size);
336 void madd_vn_vnvn(float *array_tar, const float *array_src_a, const float *array_src_b, const float f, const int size);
337 void sub_vn_vn(float *array_tar, const float *array_src, const int size);
338 void sub_vn_vnvn(float *array_tar, const float *array_src_a, const float *array_src_b, const int size);
339 void msub_vn_vn(float *array_tar, const float *array_src, const float f, const int size);
340 void msub_vn_vnvn(float *array_tar, const float *array_src_a, const float *array_src_b, const float f, const int size);
341 void interp_vn_vn(float *array_tar, const float *array_src, const float t, const int size);
342 void copy_vn_i(int *array_tar, const int size, const int val);
343 void copy_vn_short(short *array_tar, const int size, const short val);
344 void copy_vn_ushort(unsigned short *array_tar, const int size, const unsigned short val);
345 void copy_vn_uchar(unsigned char *array_tar, const int size, const unsigned char val);
346 void copy_vn_fl(float *array_tar, const int size, const float val);
347
348 void add_vn_vn_d(double *array_tar, const double *array_src, const int size);
349 void add_vn_vnvn_d(double *array_tar, const double *array_src_a, const double *array_src_b, const int size);
350 void mul_vn_db(double *array_tar, const int size, const double f);
351
352 /**************************** Inline Definitions ******************************/
353
354 #if BLI_MATH_DO_INLINE
355 #include "intern/math_vector_inline.c"
356 #endif
357
358 #ifdef BLI_MATH_GCC_WARN_PRAGMA
359 #  pragma GCC diagnostic pop
360 #endif
361
362 #ifdef __cplusplus
363 }
364 #endif
365
366 #endif /* __BLI_MATH_VECTOR_H__ */