Fix T60935: More numerically stable distance to ray computation
[blender.git] / source / blender / blenlib / BLI_math_geom.h
1 /*
2  * This program is free software; you can redistribute it and/or
3  * modify it under the terms of the GNU General Public License
4  * as published by the Free Software Foundation; either version 2
5  * of the License, or (at your option) any later version.
6  *
7  * This program is distributed in the hope that it will be useful,
8  * but WITHOUT ANY WARRANTY; without even the implied warranty of
9  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
10  * GNU General Public License for more details.
11  *
12  * You should have received a copy of the GNU General Public License
13  * along with this program; if not, write to the Free Software Foundation,
14  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
15  *
16  * The Original Code is Copyright (C) 2001-2002 by NaN Holding BV.
17  * All rights reserved.
18  *
19  * The Original Code is: some of this file.
20  *
21  * */
22
23 #ifndef __BLI_MATH_GEOM_H__
24 #define __BLI_MATH_GEOM_H__
25
26 /** \file \ingroup bli
27  */
28
29 #ifdef __cplusplus
30 extern "C" {
31 #endif
32
33 #include "BLI_compiler_attrs.h"
34 #include "BLI_math_inline.h"
35
36 #ifdef BLI_MATH_GCC_WARN_PRAGMA
37 #  pragma GCC diagnostic push
38 #  pragma GCC diagnostic ignored "-Wredundant-decls"
39 #endif
40
41 /********************************** Polygons *********************************/
42
43 float normal_tri_v3(float r[3], const float a[3], const float b[3], const float c[3]);
44 float normal_quad_v3(float r[3], const float a[3], const float b[3], const float c[3], const float d[3]);
45 float normal_poly_v3(float r[3], const float verts[][3], unsigned int nr);
46
47 MINLINE float area_tri_v2(const float a[2], const float b[2], const float c[2]);
48 MINLINE float area_squared_tri_v2(const float a[2], const float b[2], const float c[2]);
49 MINLINE float area_tri_signed_v2(const float v1[2], const float v2[2], const float v3[2]);
50 float area_tri_v3(const float a[3], const float b[3], const float c[3]);
51 float area_squared_tri_v3(const float a[3], const float b[3], const float c[3]);
52 float area_tri_signed_v3(const float v1[3], const float v2[3], const float v3[3], const float normal[3]);
53 float area_quad_v3(const float a[3], const float b[3], const float c[3], const float d[3]);
54 float area_squared_quad_v3(const float a[3], const float b[3], const float c[3], const float d[3]);
55 float area_poly_v3(const float verts[][3], unsigned int nr);
56 float area_poly_v2(const float verts[][2], unsigned int nr);
57 float area_squared_poly_v3(const float verts[][3], unsigned int nr);
58 float area_squared_poly_v2(const float verts[][2], unsigned int nr);
59 float area_poly_signed_v2(const float verts[][2], unsigned int nr);
60 float cotangent_tri_weight_v3(const float v1[3], const float v2[3], const float v3[3]);
61
62 void          cross_tri_v3(float n[3], const float v1[3], const float v2[3], const float v3[3]);
63 MINLINE float cross_tri_v2(const float v1[2], const float v2[2], const float v3[2]);
64 void cross_poly_v3(float n[3], const float verts[][3], unsigned int nr);
65 float cross_poly_v2(const float verts[][2], unsigned int nr);
66
67 /********************************* Planes **********************************/
68
69 void  plane_from_point_normal_v3(float r_plane[4], const float plane_co[3], const float plane_no[3]);
70 void  plane_to_point_vector_v3(const float plane[4], float r_plane_co[3], float r_plane_no[3]);
71 void  plane_to_point_vector_v3_normalized(const float plane[4], float r_plane_co[3], float r_plane_no[3]);
72
73 MINLINE float plane_point_side_v3(const float plane[4], const float co[3]);
74
75 /********************************* Volume **********************************/
76
77 float volume_tetrahedron_v3(const float v1[3], const float v2[3], const float v3[3], const float v4[3]);
78 float volume_tetrahedron_signed_v3(const float v1[3], const float v2[3], const float v3[3], const float v4[3]);
79
80 bool is_quad_convex_v3(const float v1[3], const float v2[3], const float v3[3], const float v4[3]);
81 bool is_quad_convex_v2(const float v1[2], const float v2[2], const float v3[2], const float v4[2]);
82 bool is_poly_convex_v2(const float verts[][2], unsigned int nr);
83 int  is_quad_flip_v3(const float v1[3], const float v2[3], const float v3[3], const float v4[3]);
84 bool is_quad_flip_v3_first_third_fast(const float v0[3], const float v1[3], const float v2[3], const float v3[3]);
85
86 /********************************* Distance **********************************/
87
88 float dist_squared_to_line_v2(const float p[2], const float l1[2], const float l2[2]);
89 float         dist_to_line_v2(const float p[2], const float l1[2], const float l2[2]);
90 float dist_squared_to_line_segment_v2(const float p[2], const float l1[2], const float l2[2]);
91 float         dist_to_line_segment_v2(const float p[2], const float l1[2], const float l2[2]);
92
93 float dist_signed_squared_to_plane_v3(const float p[3], const float plane[4]);
94 float        dist_squared_to_plane_v3(const float p[3], const float plane[4]);
95 float dist_signed_to_plane_v3(const float p[3], const float plane[4]);
96 float        dist_to_plane_v3(const float p[3], const float plane[4]);
97
98 /* plane3 versions */
99 float dist_signed_squared_to_plane3_v3(const float p[3], const float plane[4]);
100 float        dist_squared_to_plane3_v3(const float p[3], const float plane[4]);
101 float dist_signed_to_plane3_v3(const float p[3], const float plane[4]);
102 float        dist_to_plane3_v3(const float p[3], const float plane[4]);
103
104 float dist_squared_to_line_segment_v3(const float p[3], const float l1[3], const float l2[3]);
105 float         dist_to_line_segment_v3(const float p[3], const float l1[3], const float l2[3]);
106 float dist_squared_to_line_v3(const float p[3], const float l1[3], const float l2[3]);
107 float         dist_to_line_v3(const float p[3], const float l1[3], const float l2[3]);
108 float dist_signed_squared_to_corner_v3v3v3(
109         const float p[3],
110         const float v1[3], const float v2[3], const float v3[3],
111         const float axis_ref[3]);
112 float dist_squared_to_ray_v3_normalized(
113         const float ray_origin[3], const float ray_direction[3],
114         const float co[3]);
115 float dist_squared_ray_to_seg_v3(
116         const float ray_origin[3], const float ray_direction[3],
117         const float v0[3], const float v1[3],
118         float r_point[3], float *r_depth);
119
120 void aabb_get_near_far_from_plane(
121         const float plane_no[3], const float bbmin[3], const float bbmax[3],
122         float bb_near[3], float bb_afar[3]);
123
124 struct DistRayAABB_Precalc {
125         float ray_origin[3];
126         float ray_direction[3];
127         float ray_inv_dir[3];
128 };
129 void dist_squared_ray_to_aabb_v3_precalc(
130         struct DistRayAABB_Precalc *neasrest_precalc,
131         const float ray_origin[3], const float ray_direction[3]);
132 float dist_squared_ray_to_aabb_v3(
133         const struct DistRayAABB_Precalc *data,
134         const float bb_min[3], const float bb_max[3],
135         float r_point[3], float *r_depth);
136 /* when there is no advantage to precalc. */
137 float dist_squared_ray_to_aabb_v3_simple(
138         const float ray_origin[3], const float ray_direction[3],
139         const float bb_min[3], const float bb_max[3],
140         float r_point[3], float *r_depth);
141
142 struct DistProjectedAABBPrecalc {
143         float ray_origin[3];
144         float ray_direction[3];
145         float ray_inv_dir[3];
146         float pmat[4][4];
147         float mval[2];
148 };
149 void dist_squared_to_projected_aabb_precalc(
150         struct DistProjectedAABBPrecalc *precalc,
151         const float projmat[4][4], const float winsize[2], const float mval[2]);
152 float dist_squared_to_projected_aabb(
153         struct DistProjectedAABBPrecalc *data,
154         const float bbmin[3], const float bbmax[3],
155         bool r_axis_closest[3]);
156 float dist_squared_to_projected_aabb_simple(
157         const float projmat[4][4], const float winsize[2], const float mval[2],
158         const float bbmin[3], const float bbmax[3]);
159
160 float closest_to_line_v2(float r_close[2], const float p[2], const float l1[2], const float l2[2]);
161 float closest_to_line_v3(float r_close[3], const float p[3], const float l1[3], const float l2[3]);
162 void closest_to_line_segment_v2(float r_close[2], const float p[2], const float l1[2], const float l2[2]);
163 void closest_to_line_segment_v3(float r_close[3], const float p[3], const float l1[3], const float l2[3]);
164 void closest_to_plane_normalized_v3(float r_close[3], const float plane[4], const float pt[3]);
165 void closest_to_plane_v3(float r_close[3], const float plane[4], const float pt[3]);
166 void closest_to_plane3_normalized_v3(float r_close[3], const float plane[3], const float pt[3]);
167 void closest_to_plane3_v3(float r_close[3], const float plane[3], const float pt[3]);
168
169 /* Set 'r' to the point in triangle (t1, t2, t3) closest to point 'p' */
170 void closest_on_tri_to_point_v3(float r[3], const float p[3], const float t1[3], const float t2[3], const float t3[3]);
171
172 float ray_point_factor_v3_ex(
173         const float p[3], const float ray_origin[3], const float ray_direction[3],
174         const float epsilon, const float fallback);
175 float ray_point_factor_v3(
176         const float p[3], const float ray_origin[3], const float ray_direction[3]);
177
178 float line_point_factor_v3_ex(
179         const float p[3], const float l1[3], const float l2[3],
180         const float epsilon, const float fallback);
181 float line_point_factor_v3(
182         const float p[3], const float l1[3], const float l2[3]);
183
184 float line_point_factor_v2_ex(
185         const float p[2], const float l1[2], const float l2[2],
186         const float epsilon, const float fallback);
187 float line_point_factor_v2(
188         const float p[2], const float l1[2], const float l2[2]);
189
190 float line_plane_factor_v3(const float plane_co[3], const float plane_no[3],
191                            const float l1[3], const float l2[3]);
192
193 void limit_dist_v3(float v1[3], float v2[3], const float dist);
194
195 /******************************* Intersection ********************************/
196
197 /* TODO int return value consistency */
198
199 /* line-line */
200 #define ISECT_LINE_LINE_COLINEAR    -1
201 #define ISECT_LINE_LINE_NONE         0
202 #define ISECT_LINE_LINE_EXACT        1
203 #define ISECT_LINE_LINE_CROSS        2
204
205 int  isect_seg_seg_v2(const float a1[2], const float a2[2], const float b1[2], const float b2[2]);
206 void isect_seg_seg_v3(
207         const float a0[3], const float a1[3],
208         const float b0[3], const float b1[3],
209         float r_a[3], float r_b[3]);
210
211 int  isect_seg_seg_v2_int(const int a1[2], const int a2[2], const int b1[2], const int b2[2]);
212 int  isect_seg_seg_v2_point_ex(
213         const float v0[2], const float v1[2], const float v2[2], const float v3[2], const float endpoint_bias,
214         float vi[2]);
215 int  isect_seg_seg_v2_point(
216         const float v0[2], const float v1[2], const float v2[2], const float v3[2],
217         float vi[2]);
218 bool isect_seg_seg_v2_simple(
219         const float v1[2], const float v2[2], const float v3[2], const float v4[2]);
220
221 int isect_line_sphere_v3(const float l1[3], const float l2[3], const float sp[3], const float r, float r_p1[3], float r_p2[3]);
222 int isect_line_sphere_v2(const float l1[2], const float l2[2], const float sp[2], const float r, float r_p1[2], float r_p2[2]);
223
224 int isect_line_line_v2_point(
225         const float v0[2], const float v1[2],
226         const float v2[2], const float v3[2],
227         float r_vi[2]);
228 int isect_line_line_epsilon_v3(
229         const float v1[3], const float v2[3],
230         const float v3[3], const float v4[3],
231         float i1[3], float i2[3],
232         const float epsilon);
233 int isect_line_line_v3(
234         const float v1[3], const float v2[3],
235         const float v3[3], const float v4[3],
236         float r_i1[3], float r_i2[3]);
237 bool isect_line_line_strict_v3(
238         const float v1[3], const float v2[3],
239         const float v3[3], const float v4[3],
240         float vi[3], float *r_lambda);
241
242 bool isect_ray_plane_v3(
243         const float ray_origin[3], const float ray_direction[3],
244         const float plane[4],
245         float *r_lambda, const bool clip);
246
247 bool isect_point_planes_v3(float (*planes)[4], int totplane, const float p[3]);
248 bool isect_point_planes_v3_negated(
249         const float (*planes)[4], const int totplane, const float p[3]);
250
251 bool isect_line_plane_v3(
252         float r_isect_co[3], const float l1[3], const float l2[3],
253         const float plane_co[3], const float plane_no[3]) ATTR_WARN_UNUSED_RESULT;
254
255 bool isect_plane_plane_plane_v3(
256         const float plane_a[4], const float plane_b[4], const float plane_c[4],
257         float r_isect_co[3]) ATTR_WARN_UNUSED_RESULT;
258 bool isect_plane_plane_v3(
259         const float plane_a[4], const float plane_b[4],
260         float r_isect_co[3], float r_isect_no[3]) ATTR_WARN_UNUSED_RESULT;
261
262 /* line/ray triangle */
263 bool isect_line_segment_tri_v3(
264         const float p1[3], const float p2[3],
265         const float v0[3], const float v1[3], const float v2[3],
266         float *r_lambda, float r_uv[2]);
267 bool isect_line_segment_tri_epsilon_v3(
268         const float p1[3], const float p2[3],
269         const float v0[3], const float v1[3], const float v2[3],
270         float *r_lambda, float r_uv[2], const float epsilon);
271 bool isect_axial_line_segment_tri_v3(
272         const int axis, const float p1[3], const float p2[3],
273         const float v0[3], const float v1[3], const float v2[3], float *r_lambda);
274
275 bool isect_ray_tri_v3(
276         const float ray_origin[3], const float ray_direction[3],
277         const float v0[3], const float v1[3], const float v2[3],
278         float *r_lambda, float r_uv[2]);
279 bool isect_ray_tri_threshold_v3(
280         const float ray_origin[3], const float ray_direction[3],
281         const float v0[3], const float v1[3], const float v2[3],
282         float *r_lambda, float r_uv[2], const float threshold);
283 bool isect_ray_tri_epsilon_v3(
284         const float ray_origin[3], const float ray_direction[3],
285         const float v0[3], const float v1[3], const float v2[3],
286         float *r_lambda, float r_uv[2], const float epsilon);
287 bool isect_tri_tri_epsilon_v3(
288         const float t_a0[3], const float t_a1[3], const float t_a2[3],
289         const float t_b0[3], const float t_b1[3], const float t_b2[3],
290         float r_i1[3], float r_i2[3],
291         const float epsilon);
292
293 /* water-tight raycast (requires pre-calculation) */
294 struct IsectRayPrecalc {
295         /* Maximal dimension kz, and orthogonal dimensions. */
296         int kx, ky, kz;
297
298         /* Shear constants. */
299         float sx, sy, sz;
300 };
301
302 void isect_ray_tri_watertight_v3_precalc(
303         struct IsectRayPrecalc *isect_precalc, const float ray_direction[3]);
304 bool isect_ray_tri_watertight_v3(
305         const float ray_origin[3], const struct IsectRayPrecalc *isect_precalc,
306         const float v0[3], const float v1[3], const float v2[3],
307         float *r_dist, float r_uv[2]);
308 /* slower version which calculates IsectRayPrecalc each time */
309 bool isect_ray_tri_watertight_v3_simple(
310         const float ray_origin[3], const float ray_direction[3],
311         const float v0[3], const float v1[3], const float v2[3],
312         float *r_lambda, float r_uv[2]);
313
314 bool isect_ray_seg_v2(
315         const float ray_origin[2], const float ray_direction[2],
316         const float v0[2], const float v1[2],
317         float *r_lambda, float *r_u);
318
319 bool isect_ray_seg_v3(
320         const float ray_origin[3], const float ray_direction[3],
321         const float v0[3], const float v1[3],
322         float *r_lambda);
323
324 /* point in polygon */
325 bool isect_point_poly_v2(const float pt[2], const float verts[][2], const unsigned int nr, const bool use_holes);
326 bool isect_point_poly_v2_int(const int pt[2], const int verts[][2], const unsigned int nr, const bool use_holes);
327
328 int isect_point_quad_v2(const float p[2], const float a[2], const float b[2], const float c[2], const float d[2]);
329
330 int  isect_point_tri_v2(const float pt[2], const float v1[2], const float v2[2], const float v3[2]);
331 bool isect_point_tri_v2_cw(const float pt[2], const float v1[2], const float v2[2], const float v3[2]);
332 int  isect_point_tri_v2_int(const int x1, const int y1, const int x2, const int y2, const int a, const int b);
333 bool isect_point_tri_prism_v3(const float p[3], const float v1[3], const float v2[3], const float v3[3]);
334 bool isect_point_tri_v3(
335         const float p[3], const float v1[3], const float v2[3], const float v3[3],
336         float r_isect_co[3]);
337
338 /* axis-aligned bounding box */
339 bool isect_aabb_aabb_v3(const float min1[3], const float max1[3], const float min2[3], const float max2[3]);
340
341 struct IsectRayAABB_Precalc {
342         float ray_origin[3];
343         float ray_inv_dir[3];
344         int sign[3];
345 };
346
347 void isect_ray_aabb_v3_precalc(
348         struct IsectRayAABB_Precalc *data,
349         const float ray_origin[3], const float ray_direction[3]);
350 bool isect_ray_aabb_v3(
351         const struct IsectRayAABB_Precalc *data,
352         const float bb_min[3], const float bb_max[3], float *tmin);
353 bool isect_ray_aabb_v3_simple(
354         const float orig[3], const float dir[3],
355         const float bb_min[3], const float bb_max[3],
356         float *tmin, float *tmax);
357
358 /* other */
359 #define ISECT_AABB_PLANE_BEHIND_ANY   0
360 #define ISECT_AABB_PLANE_CROSS_ANY    1
361 #define ISECT_AABB_PLANE_IN_FRONT_ALL 2
362
363 int isect_aabb_planes_v3(
364         const float (*planes)[4], const int totplane,
365         const float bbmin[3], const float bbmax[3]);
366
367 bool isect_sweeping_sphere_tri_v3(const float p1[3], const float p2[3], const float radius,
368                                   const float v0[3], const float v1[3], const float v2[3], float *r_lambda, float ipoint[3]);
369
370 bool clip_segment_v3_plane(
371         const float p1[3], const float p2[3], const float plane[4],
372         float r_p1[3], float r_p2[3]);
373 bool clip_segment_v3_plane_n(
374         const float p1[3], const float p2[3], const float plane_array[][4], const int plane_tot,
375         float r_p1[3], float r_p2[3]);
376
377 bool point_in_slice_seg(float p[3], float l1[3], float l2[3]);
378
379 /****************************** Interpolation ********************************/
380 void interp_weights_tri_v3(float w[3], const float a[3], const float b[3], const float c[3], const float p[3]);
381 void interp_weights_quad_v3(float w[4], const float a[3], const float b[3], const float c[3], const float d[3], const float p[3]);
382 void interp_weights_poly_v3(float w[], float v[][3], const int n, const float co[3]);
383 void interp_weights_poly_v2(float w[], float v[][2], const int n, const float co[2]);
384
385 void interp_cubic_v3(float x[3], float v[3],
386                      const float x1[3], const float v1[3], const float x2[3], const float v2[3], const float t);
387
388 int interp_sparse_array(float *array, const int list_size, const float invalid);
389
390 void transform_point_by_tri_v3(
391         float pt_tar[3], float const pt_src[3],
392         const float tri_tar_p1[3], const float tri_tar_p2[3], const float tri_tar_p3[3],
393         const float tri_src_p1[3], const float tri_src_p2[3], const float tri_src_p3[3]);
394 void transform_point_by_seg_v3(
395         float p_dst[3], const float p_src[3],
396         const float l_dst_p1[3], const float l_dst_p2[3],
397         const float l_src_p1[3], const float l_src_p2[3]);
398
399 void barycentric_weights_v2(
400         const float v1[2], const float v2[2], const float v3[2],
401         const float co[2], float w[3]);
402 void barycentric_weights_v2_clamped(
403         const float v1[2], const float v2[2], const float v3[2],
404         const float co[2], float w[3]);
405 void barycentric_weights_v2_persp(
406         const float v1[4], const float v2[4], const float v3[4],
407         const float co[2], float w[3]);
408 void barycentric_weights_v2_quad(
409         const float v1[2], const float v2[2], const float v3[2], const float v4[2],
410         const float co[2], float w[4]);
411
412 bool barycentric_coords_v2(const float v1[2], const float v2[2], const float v3[2], const float co[2], float w[3]);
413 int barycentric_inside_triangle_v2(const float w[3]);
414
415 void resolve_tri_uv_v2(float r_uv[2], const float st[2], const float st0[2], const float st1[2], const float st2[2]);
416 void resolve_tri_uv_v3(float r_uv[2], const float st[3], const float st0[3], const float st1[3], const float st2[3]);
417 void resolve_quad_uv_v2(float r_uv[2], const float st[2], const float st0[2], const float st1[2], const float st2[2], const float st3[2]);
418 void resolve_quad_uv_v2_deriv(float r_uv[2], float r_deriv[2][2],
419                               const float st[2], const float st0[2], const float st1[2], const float st2[2], const float st3[2]);
420 float resolve_quad_u_v2(const float st[2], const float st0[2], const float st1[2], const float st2[2], const float st3[2]);
421
422 /* use to find the point of a UV on a face */
423 void interp_bilinear_quad_v3(float data[4][3], float u, float v, float res[3]);
424 void interp_barycentric_tri_v3(float data[3][3], float u, float v, float res[3]);
425
426 /***************************** View & Projection *****************************/
427
428 void lookat_m4(float mat[4][4], float vx, float vy,
429                float vz, float px, float py, float pz, float twist);
430 void polarview_m4(float mat[4][4], float dist, float azimuth,
431                   float incidence, float twist);
432
433 void perspective_m4(float mat[4][4], const float left, const float right,
434                     const float bottom, const float top, const float nearClip, const float farClip);
435 void orthographic_m4(float mat[4][4], const float left, const float right,
436                      const float bottom, const float top, const float nearClip, const float farClip);
437 void window_translate_m4(float winmat[4][4], float perspmat[4][4],
438                          const float x, const float y);
439
440 void planes_from_projmat(float mat[4][4], float left[4], float right[4], float top[4], float bottom[4],
441                          float front[4], float back[4]);
442
443 int box_clip_bounds_m4(float boundbox[2][3],
444                        const float bounds[4], float winmat[4][4]);
445 void box_minmax_bounds_m4(float min[3], float max[3],
446                           float boundbox[2][3], float mat[4][4]);
447
448 /********************************** Mapping **********************************/
449
450 void map_to_tube(float *r_u, float *r_v, const float x, const float y, const float z);
451 void map_to_sphere(float *r_u, float *r_v, const float x, const float y, const float z);
452 void map_to_plane_v2_v3v3(float r_co[2], const float co[3], const float no[3]);
453 void map_to_plane_axis_angle_v2_v3v3fl(float r_co[2], const float co[3], const float axis[3], const float angle);
454
455 /********************************** Normals **********************************/
456
457 void accumulate_vertex_normals_tri_v3(
458         float n1[3], float n2[3], float n3[3],
459         const float f_no[3],
460         const float co1[3], const float co2[3], const float co3[3]);
461
462 void accumulate_vertex_normals_v3(
463         float n1[3], float n2[3], float n3[3], float n4[3],
464         const float f_no[3],
465         const float co1[3], const float co2[3], const float co3[3], const float co4[3]);
466
467 void accumulate_vertex_normals_poly_v3(
468         float **vertnos, const float polyno[3],
469         const float **vertcos, float vdiffs[][3], const int nverts);
470
471 /********************************* Tangents **********************************/
472
473 void tangent_from_uv_v3(
474         const float uv1[2], const float uv2[2], const float uv3[2],
475         const float co1[3], const float co2[3], const float co3[3],
476         const float n[3],
477         float r_tang[3]);
478
479 /******************************** Vector Clouds ******************************/
480
481 void vcloud_estimate_transform_v3(
482         const int list_size, const float (*pos)[3], const float *weight, const float (*rpos)[3], const float *rweight,
483         float lloc[3], float rloc[3], float lrot[3][3], float lscale[3][3]);
484
485 /****************************** Spherical Harmonics *************************/
486
487 /* Uses 2nd order SH => 9 coefficients, stored in this order:
488  * 0 = (0, 0),
489  * 1 = (1, -1), 2 = (1, 0), 3 = (1, 1),
490  * 4 = (2, -2), 5 = (2, -1), 6 = (2, 0), 7 = (2, 1), 8 = (2, 2) */
491
492 MINLINE void zero_sh(float r[9]);
493 MINLINE void copy_sh_sh(float r[9], const float a[9]);
494 MINLINE void mul_sh_fl(float r[9], const float f);
495 MINLINE void add_sh_shsh(float r[9], const float a[9], const float b[9]);
496 MINLINE float dot_shsh(const float a[9], const float b[9]);
497
498 MINLINE float eval_shv3(float r[9], const float v[3]);
499 MINLINE float diffuse_shv3(float r[9], const float v[3]);
500 MINLINE void vec_fac_to_sh(float r[9], const float v[3], const float f);
501 MINLINE void madd_sh_shfl(float r[9], const float sh[3], const float f);
502
503 /********************************* Form Factor *******************************/
504
505 float form_factor_quad(const float p[3], const float n[3],
506                        const float q0[3], const float q1[3], const float q2[3], const float q3[3]);
507 bool form_factor_visible_quad(const float p[3], const float n[3],
508                               const float v0[3], const float v1[3], const float v2[3],
509                               float q0[3], float q1[3], float q2[3], float q3[3]);
510 float form_factor_hemi_poly(float p[3], float n[3],
511                             float v1[3], float v2[3], float v3[3], float v4[3]);
512
513 void axis_dominant_v3_to_m3_negate(float r_mat[3][3], const float normal[3]);
514 void axis_dominant_v3_to_m3(float r_mat[3][3], const float normal[3]);
515
516 MINLINE void  axis_dominant_v3(int *r_axis_a, int *r_axis_b, const float axis[3]);
517 MINLINE float axis_dominant_v3_max(int *r_axis_a, int *r_axis_b, const float axis[3]) ATTR_WARN_UNUSED_RESULT;
518 MINLINE int   axis_dominant_v3_single(const float vec[3]);
519 MINLINE int   axis_dominant_v3_ortho_single(const float vec[3]);
520
521 MINLINE int max_axis_v3(const float vec[3]);
522 MINLINE int min_axis_v3(const float vec[3]);
523
524 MINLINE int poly_to_tri_count(const int poly_count, const int corner_count);
525
526 MINLINE float shell_angle_to_dist(const float angle);
527 MINLINE float shell_v3v3_normalized_to_dist(const float a[3], const float b[3]);
528 MINLINE float shell_v2v2_normalized_to_dist(const float a[2], const float b[2]);
529 MINLINE float shell_v3v3_mid_normalized_to_dist(const float a[3], const float b[3]);
530 MINLINE float shell_v2v2_mid_normalized_to_dist(const float a[2], const float b[2]);
531
532 /********************************* Cubic (Bezier) *******************************/
533
534 float cubic_tangent_factor_circle_v3(const float tan_l[3], const float tan_r[3]);
535
536 /**************************** Inline Definitions ******************************/
537
538 #if BLI_MATH_DO_INLINE
539 #include "intern/math_geom_inline.c"
540 #endif
541
542 #ifdef BLI_MATH_GCC_WARN_PRAGMA
543 #  pragma GCC diagnostic pop
544 #endif
545
546 #ifdef __cplusplus
547 }
548 #endif
549
550 #endif /* __BLI_MATH_GEOM_H__ */