diff options
author | Campbell Barton <ideasman42@gmail.com> | 2019-04-17 07:17:24 +0300 |
---|---|---|
committer | Campbell Barton <ideasman42@gmail.com> | 2019-04-17 07:21:24 +0300 |
commit | e12c08e8d170b7ca40f204a5b0423c23a9fbc2c1 (patch) | |
tree | 8cf3453d12edb177a218ef8009357518ec6cab6a /source/blender/blenlib/BLI_math_geom.h | |
parent | b3dabc200a4b0399ec6b81f2ff2730d07b44fcaa (diff) |
ClangFormat: apply to source, most of intern
Apply clang format as proposed in T53211.
For details on usage and instructions for migrating branches
without conflicts, see:
https://wiki.blender.org/wiki/Tools/ClangFormat
Diffstat (limited to 'source/blender/blenlib/BLI_math_geom.h')
-rw-r--r-- | source/blender/blenlib/BLI_math_geom.h | 850 |
1 files changed, 532 insertions, 318 deletions
diff --git a/source/blender/blenlib/BLI_math_geom.h b/source/blender/blenlib/BLI_math_geom.h index 2f81fbf37f6..484d5af194d 100644 --- a/source/blender/blenlib/BLI_math_geom.h +++ b/source/blender/blenlib/BLI_math_geom.h @@ -42,7 +42,8 @@ extern "C" { /********************************** Polygons *********************************/ float normal_tri_v3(float r[3], const float a[3], const float b[3], const float c[3]); -float normal_quad_v3(float r[3], const float a[3], const float b[3], const float c[3], const float d[3]); +float normal_quad_v3( + float r[3], const float a[3], const float b[3], const float c[3], const float d[3]); float normal_poly_v3(float r[3], const float verts[][3], unsigned int nr); MINLINE float area_tri_v2(const float a[2], const float b[2], const float c[2]); @@ -50,7 +51,10 @@ MINLINE float area_squared_tri_v2(const float a[2], const float b[2], const floa MINLINE float area_tri_signed_v2(const float v1[2], const float v2[2], const float v3[2]); float area_tri_v3(const float a[3], const float b[3], const float c[3]); float area_squared_tri_v3(const float a[3], const float b[3], const float c[3]); -float area_tri_signed_v3(const float v1[3], const float v2[3], const float v3[3], const float normal[3]); +float area_tri_signed_v3(const float v1[3], + const float v2[3], + const float v3[3], + const float normal[3]); float area_quad_v3(const float a[3], const float b[3], const float c[3], const float d[3]); float area_squared_quad_v3(const float a[3], const float b[3], const float c[3], const float d[3]); float area_poly_v3(const float verts[][3], unsigned int nr); @@ -60,136 +64,174 @@ float area_squared_poly_v2(const float verts[][2], unsigned int nr); float area_poly_signed_v2(const float verts[][2], unsigned int nr); float cotangent_tri_weight_v3(const float v1[3], const float v2[3], const float v3[3]); -void cross_tri_v3(float n[3], const float v1[3], const float v2[3], const float v3[3]); +void cross_tri_v3(float n[3], const float v1[3], const float v2[3], const float v3[3]); MINLINE float cross_tri_v2(const float v1[2], const float v2[2], const float v3[2]); void cross_poly_v3(float n[3], const float verts[][3], unsigned int nr); float cross_poly_v2(const float verts[][2], unsigned int nr); /********************************* Planes **********************************/ -void plane_from_point_normal_v3(float r_plane[4], const float plane_co[3], const float plane_no[3]); -void plane_to_point_vector_v3(const float plane[4], float r_plane_co[3], float r_plane_no[3]); -void plane_to_point_vector_v3_normalized(const float plane[4], float r_plane_co[3], float r_plane_no[3]); +void plane_from_point_normal_v3(float r_plane[4], + const float plane_co[3], + const float plane_no[3]); +void plane_to_point_vector_v3(const float plane[4], float r_plane_co[3], float r_plane_no[3]); +void plane_to_point_vector_v3_normalized(const float plane[4], + float r_plane_co[3], + float r_plane_no[3]); MINLINE float plane_point_side_v3(const float plane[4], const float co[3]); /********************************* Volume **********************************/ -float volume_tetrahedron_v3(const float v1[3], const float v2[3], const float v3[3], const float v4[3]); -float volume_tetrahedron_signed_v3(const float v1[3], const float v2[3], const float v3[3], const float v4[3]); +float volume_tetrahedron_v3(const float v1[3], + const float v2[3], + const float v3[3], + const float v4[3]); +float volume_tetrahedron_signed_v3(const float v1[3], + const float v2[3], + const float v3[3], + const float v4[3]); bool is_quad_convex_v3(const float v1[3], const float v2[3], const float v3[3], const float v4[3]); bool is_quad_convex_v2(const float v1[2], const float v2[2], const float v3[2], const float v4[2]); bool is_poly_convex_v2(const float verts[][2], unsigned int nr); -int is_quad_flip_v3(const float v1[3], const float v2[3], const float v3[3], const float v4[3]); -bool is_quad_flip_v3_first_third_fast(const float v0[3], const float v1[3], const float v2[3], const float v3[3]); +int is_quad_flip_v3(const float v1[3], const float v2[3], const float v3[3], const float v4[3]); +bool is_quad_flip_v3_first_third_fast(const float v0[3], + const float v1[3], + const float v2[3], + const float v3[3]); /********************************* Distance **********************************/ float dist_squared_to_line_v2(const float p[2], const float l1[2], const float l2[2]); -float dist_to_line_v2(const float p[2], const float l1[2], const float l2[2]); +float dist_to_line_v2(const float p[2], const float l1[2], const float l2[2]); float dist_squared_to_line_segment_v2(const float p[2], const float l1[2], const float l2[2]); -float dist_to_line_segment_v2(const float p[2], const float l1[2], const float l2[2]); +float dist_to_line_segment_v2(const float p[2], const float l1[2], const float l2[2]); float dist_signed_squared_to_plane_v3(const float p[3], const float plane[4]); -float dist_squared_to_plane_v3(const float p[3], const float plane[4]); +float dist_squared_to_plane_v3(const float p[3], const float plane[4]); float dist_signed_to_plane_v3(const float p[3], const float plane[4]); -float dist_to_plane_v3(const float p[3], const float plane[4]); +float dist_to_plane_v3(const float p[3], const float plane[4]); /* plane3 versions */ float dist_signed_squared_to_plane3_v3(const float p[3], const float plane[4]); -float dist_squared_to_plane3_v3(const float p[3], const float plane[4]); +float dist_squared_to_plane3_v3(const float p[3], const float plane[4]); float dist_signed_to_plane3_v3(const float p[3], const float plane[4]); -float dist_to_plane3_v3(const float p[3], const float plane[4]); +float dist_to_plane3_v3(const float p[3], const float plane[4]); float dist_squared_to_line_segment_v3(const float p[3], const float l1[3], const float l2[3]); -float dist_to_line_segment_v3(const float p[3], const float l1[3], const float l2[3]); +float dist_to_line_segment_v3(const float p[3], const float l1[3], const float l2[3]); float dist_squared_to_line_v3(const float p[3], const float l1[3], const float l2[3]); -float dist_to_line_v3(const float p[3], const float l1[3], const float l2[3]); -float dist_signed_squared_to_corner_v3v3v3( - const float p[3], - const float v1[3], const float v2[3], const float v3[3], - const float axis_ref[3]); -float dist_squared_to_ray_v3_normalized( - const float ray_origin[3], const float ray_direction[3], - const float co[3]); -float dist_squared_ray_to_seg_v3( - const float ray_origin[3], const float ray_direction[3], - const float v0[3], const float v1[3], - float r_point[3], float *r_depth); - -void aabb_get_near_far_from_plane( - const float plane_no[3], const float bbmin[3], const float bbmax[3], - float bb_near[3], float bb_afar[3]); +float dist_to_line_v3(const float p[3], const float l1[3], const float l2[3]); +float dist_signed_squared_to_corner_v3v3v3(const float p[3], + const float v1[3], + const float v2[3], + const float v3[3], + const float axis_ref[3]); +float dist_squared_to_ray_v3_normalized(const float ray_origin[3], + const float ray_direction[3], + const float co[3]); +float dist_squared_ray_to_seg_v3(const float ray_origin[3], + const float ray_direction[3], + const float v0[3], + const float v1[3], + float r_point[3], + float *r_depth); + +void aabb_get_near_far_from_plane(const float plane_no[3], + const float bbmin[3], + const float bbmax[3], + float bb_near[3], + float bb_afar[3]); struct DistRayAABB_Precalc { - float ray_origin[3]; - float ray_direction[3]; - float ray_inv_dir[3]; + float ray_origin[3]; + float ray_direction[3]; + float ray_inv_dir[3]; }; -void dist_squared_ray_to_aabb_v3_precalc( - struct DistRayAABB_Precalc *neasrest_precalc, - const float ray_origin[3], const float ray_direction[3]); -float dist_squared_ray_to_aabb_v3( - const struct DistRayAABB_Precalc *data, - const float bb_min[3], const float bb_max[3], - float r_point[3], float *r_depth); +void dist_squared_ray_to_aabb_v3_precalc(struct DistRayAABB_Precalc *neasrest_precalc, + const float ray_origin[3], + const float ray_direction[3]); +float dist_squared_ray_to_aabb_v3(const struct DistRayAABB_Precalc *data, + const float bb_min[3], + const float bb_max[3], + float r_point[3], + float *r_depth); /* when there is no advantage to precalc. */ -float dist_squared_ray_to_aabb_v3_simple( - const float ray_origin[3], const float ray_direction[3], - const float bb_min[3], const float bb_max[3], - float r_point[3], float *r_depth); +float dist_squared_ray_to_aabb_v3_simple(const float ray_origin[3], + const float ray_direction[3], + const float bb_min[3], + const float bb_max[3], + float r_point[3], + float *r_depth); struct DistProjectedAABBPrecalc { - float ray_origin[3]; - float ray_direction[3]; - float ray_inv_dir[3]; - float pmat[4][4]; - float mval[2]; + float ray_origin[3]; + float ray_direction[3]; + float ray_inv_dir[3]; + float pmat[4][4]; + float mval[2]; }; -void dist_squared_to_projected_aabb_precalc( - struct DistProjectedAABBPrecalc *precalc, - const float projmat[4][4], const float winsize[2], const float mval[2]); -float dist_squared_to_projected_aabb( - struct DistProjectedAABBPrecalc *data, - const float bbmin[3], const float bbmax[3], - bool r_axis_closest[3]); -float dist_squared_to_projected_aabb_simple( - const float projmat[4][4], const float winsize[2], const float mval[2], - const float bbmin[3], const float bbmax[3]); +void dist_squared_to_projected_aabb_precalc(struct DistProjectedAABBPrecalc *precalc, + const float projmat[4][4], + const float winsize[2], + const float mval[2]); +float dist_squared_to_projected_aabb(struct DistProjectedAABBPrecalc *data, + const float bbmin[3], + const float bbmax[3], + bool r_axis_closest[3]); +float dist_squared_to_projected_aabb_simple(const float projmat[4][4], + const float winsize[2], + const float mval[2], + const float bbmin[3], + const float bbmax[3]); float closest_to_line_v2(float r_close[2], const float p[2], const float l1[2], const float l2[2]); float closest_to_line_v3(float r_close[3], const float p[3], const float l1[3], const float l2[3]); -void closest_to_line_segment_v2(float r_close[2], const float p[2], const float l1[2], const float l2[2]); -void closest_to_line_segment_v3(float r_close[3], const float p[3], const float l1[3], const float l2[3]); +void closest_to_line_segment_v2(float r_close[2], + const float p[2], + const float l1[2], + const float l2[2]); +void closest_to_line_segment_v3(float r_close[3], + const float p[3], + const float l1[3], + const float l2[3]); void closest_to_plane_normalized_v3(float r_close[3], const float plane[4], const float pt[3]); void closest_to_plane_v3(float r_close[3], const float plane[4], const float pt[3]); void closest_to_plane3_normalized_v3(float r_close[3], const float plane[3], const float pt[3]); void closest_to_plane3_v3(float r_close[3], const float plane[3], const float pt[3]); /* Set 'r' to the point in triangle (t1, t2, t3) closest to point 'p' */ -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]); - -float ray_point_factor_v3_ex( - const float p[3], const float ray_origin[3], const float ray_direction[3], - const float epsilon, const float fallback); -float ray_point_factor_v3( - const float p[3], const float ray_origin[3], const float ray_direction[3]); - -float line_point_factor_v3_ex( - const float p[3], const float l1[3], const float l2[3], - const float epsilon, const float fallback); -float line_point_factor_v3( - const float p[3], const float l1[3], const float l2[3]); - -float line_point_factor_v2_ex( - const float p[2], const float l1[2], const float l2[2], - const float epsilon, const float fallback); -float line_point_factor_v2( - const float p[2], const float l1[2], const float l2[2]); - -float line_plane_factor_v3(const float plane_co[3], const float plane_no[3], - const float l1[3], const float l2[3]); +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]); + +float ray_point_factor_v3_ex(const float p[3], + const float ray_origin[3], + const float ray_direction[3], + const float epsilon, + const float fallback); +float ray_point_factor_v3(const float p[3], + const float ray_origin[3], + const float ray_direction[3]); + +float line_point_factor_v3_ex(const float p[3], + const float l1[3], + const float l2[3], + const float epsilon, + const float fallback); +float line_point_factor_v3(const float p[3], const float l1[3], const float l2[3]); + +float line_point_factor_v2_ex(const float p[2], + const float l1[2], + const float l2[2], + const float epsilon, + const float fallback); +float line_point_factor_v2(const float p[2], const float l1[2], const float l2[2]); + +float line_plane_factor_v3(const float plane_co[3], + const float plane_no[3], + const float l1[3], + const float l2[3]); void limit_dist_v3(float v1[3], float v2[3], const float dist); @@ -198,227 +240,351 @@ void limit_dist_v3(float v1[3], float v2[3], const float dist); /* TODO int return value consistency */ /* line-line */ -#define ISECT_LINE_LINE_COLINEAR -1 -#define ISECT_LINE_LINE_NONE 0 -#define ISECT_LINE_LINE_EXACT 1 -#define ISECT_LINE_LINE_CROSS 2 - -int isect_seg_seg_v2(const float a1[2], const float a2[2], const float b1[2], const float b2[2]); -void isect_seg_seg_v3( - const float a0[3], const float a1[3], - const float b0[3], const float b1[3], - float r_a[3], float r_b[3]); - -int isect_seg_seg_v2_int(const int a1[2], const int a2[2], const int b1[2], const int b2[2]); -int isect_seg_seg_v2_point_ex( - const float v0[2], const float v1[2], const float v2[2], const float v3[2], const float endpoint_bias, - float vi[2]); -int isect_seg_seg_v2_point( - const float v0[2], const float v1[2], const float v2[2], const float v3[2], - float vi[2]); -bool isect_seg_seg_v2_simple( - const float v1[2], const float v2[2], const float v3[2], const float v4[2]); - -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]); -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]); +#define ISECT_LINE_LINE_COLINEAR -1 +#define ISECT_LINE_LINE_NONE 0 +#define ISECT_LINE_LINE_EXACT 1 +#define ISECT_LINE_LINE_CROSS 2 + +int isect_seg_seg_v2(const float a1[2], const float a2[2], const float b1[2], const float b2[2]); +void isect_seg_seg_v3(const float a0[3], + const float a1[3], + const float b0[3], + const float b1[3], + float r_a[3], + float r_b[3]); + +int isect_seg_seg_v2_int(const int a1[2], const int a2[2], const int b1[2], const int b2[2]); +int isect_seg_seg_v2_point_ex(const float v0[2], + const float v1[2], + const float v2[2], + const float v3[2], + const float endpoint_bias, + float vi[2]); +int isect_seg_seg_v2_point( + const float v0[2], const float v1[2], const float v2[2], const float v3[2], float vi[2]); +bool isect_seg_seg_v2_simple(const float v1[2], + const float v2[2], + const float v3[2], + const float v4[2]); + +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]); +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]); int isect_line_line_v2_point( - const float v0[2], const float v1[2], - const float v2[2], const float v3[2], - float r_vi[2]); -int isect_line_line_epsilon_v3( - const float v1[3], const float v2[3], - const float v3[3], const float v4[3], - float i1[3], float i2[3], - const float epsilon); -int isect_line_line_v3( - const float v1[3], const float v2[3], - const float v3[3], const float v4[3], - float r_i1[3], float r_i2[3]); -bool isect_line_line_strict_v3( - const float v1[3], const float v2[3], - const float v3[3], const float v4[3], - float vi[3], float *r_lambda); - -bool isect_ray_plane_v3( - const float ray_origin[3], const float ray_direction[3], - const float plane[4], - float *r_lambda, const bool clip); + const float v0[2], const float v1[2], const float v2[2], const float v3[2], float r_vi[2]); +int isect_line_line_epsilon_v3(const float v1[3], + const float v2[3], + const float v3[3], + const float v4[3], + float i1[3], + float i2[3], + const float epsilon); +int isect_line_line_v3(const float v1[3], + const float v2[3], + const float v3[3], + const float v4[3], + float r_i1[3], + float r_i2[3]); +bool isect_line_line_strict_v3(const float v1[3], + const float v2[3], + const float v3[3], + const float v4[3], + float vi[3], + float *r_lambda); + +bool isect_ray_plane_v3(const float ray_origin[3], + const float ray_direction[3], + const float plane[4], + float *r_lambda, + const bool clip); bool isect_point_planes_v3(float (*planes)[4], int totplane, const float p[3]); -bool isect_point_planes_v3_negated( - const float (*planes)[4], const int totplane, const float p[3]); - -bool isect_line_plane_v3( - float r_isect_co[3], const float l1[3], const float l2[3], - const float plane_co[3], const float plane_no[3]) ATTR_WARN_UNUSED_RESULT; - -bool isect_plane_plane_plane_v3( - const float plane_a[4], const float plane_b[4], const float plane_c[4], - float r_isect_co[3]) ATTR_WARN_UNUSED_RESULT; -bool isect_plane_plane_v3( - const float plane_a[4], const float plane_b[4], - float r_isect_co[3], float r_isect_no[3]) ATTR_WARN_UNUSED_RESULT; +bool isect_point_planes_v3_negated(const float (*planes)[4], const int totplane, const float p[3]); + +bool isect_line_plane_v3(float r_isect_co[3], + const float l1[3], + const float l2[3], + const float plane_co[3], + const float plane_no[3]) ATTR_WARN_UNUSED_RESULT; + +bool isect_plane_plane_plane_v3(const float plane_a[4], + const float plane_b[4], + const float plane_c[4], + float r_isect_co[3]) ATTR_WARN_UNUSED_RESULT; +bool isect_plane_plane_v3(const float plane_a[4], + const float plane_b[4], + float r_isect_co[3], + float r_isect_no[3]) ATTR_WARN_UNUSED_RESULT; /* line/ray triangle */ -bool isect_line_segment_tri_v3( - const float p1[3], const float p2[3], - const float v0[3], const float v1[3], const float v2[3], - float *r_lambda, float r_uv[2]); -bool isect_line_segment_tri_epsilon_v3( - const float p1[3], const float p2[3], - const float v0[3], const float v1[3], const float v2[3], - float *r_lambda, float r_uv[2], const float epsilon); -bool isect_axial_line_segment_tri_v3( - const int axis, const float p1[3], const float p2[3], - const float v0[3], const float v1[3], const float v2[3], float *r_lambda); - -bool isect_ray_tri_v3( - const float ray_origin[3], const float ray_direction[3], - const float v0[3], const float v1[3], const float v2[3], - float *r_lambda, float r_uv[2]); -bool isect_ray_tri_threshold_v3( - const float ray_origin[3], const float ray_direction[3], - const float v0[3], const float v1[3], const float v2[3], - float *r_lambda, float r_uv[2], const float threshold); -bool isect_ray_tri_epsilon_v3( - const float ray_origin[3], const float ray_direction[3], - const float v0[3], const float v1[3], const float v2[3], - float *r_lambda, float r_uv[2], const float epsilon); -bool isect_tri_tri_epsilon_v3( - const float t_a0[3], const float t_a1[3], const float t_a2[3], - const float t_b0[3], const float t_b1[3], const float t_b2[3], - float r_i1[3], float r_i2[3], - const float epsilon); +bool isect_line_segment_tri_v3(const float p1[3], + const float p2[3], + const float v0[3], + const float v1[3], + const float v2[3], + float *r_lambda, + float r_uv[2]); +bool isect_line_segment_tri_epsilon_v3(const float p1[3], + const float p2[3], + const float v0[3], + const float v1[3], + const float v2[3], + float *r_lambda, + float r_uv[2], + const float epsilon); +bool isect_axial_line_segment_tri_v3(const int axis, + const float p1[3], + const float p2[3], + const float v0[3], + const float v1[3], + const float v2[3], + float *r_lambda); + +bool isect_ray_tri_v3(const float ray_origin[3], + const float ray_direction[3], + const float v0[3], + const float v1[3], + const float v2[3], + float *r_lambda, + float r_uv[2]); +bool isect_ray_tri_threshold_v3(const float ray_origin[3], + const float ray_direction[3], + const float v0[3], + const float v1[3], + const float v2[3], + float *r_lambda, + float r_uv[2], + const float threshold); +bool isect_ray_tri_epsilon_v3(const float ray_origin[3], + const float ray_direction[3], + const float v0[3], + const float v1[3], + const float v2[3], + float *r_lambda, + float r_uv[2], + const float epsilon); +bool isect_tri_tri_epsilon_v3(const float t_a0[3], + const float t_a1[3], + const float t_a2[3], + const float t_b0[3], + const float t_b1[3], + const float t_b2[3], + float r_i1[3], + float r_i2[3], + const float epsilon); /* water-tight raycast (requires pre-calculation) */ struct IsectRayPrecalc { - /* Maximal dimension kz, and orthogonal dimensions. */ - int kx, ky, kz; + /* Maximal dimension kz, and orthogonal dimensions. */ + int kx, ky, kz; - /* Shear constants. */ - float sx, sy, sz; + /* Shear constants. */ + float sx, sy, sz; }; -void isect_ray_tri_watertight_v3_precalc( - struct IsectRayPrecalc *isect_precalc, const float ray_direction[3]); -bool isect_ray_tri_watertight_v3( - const float ray_origin[3], const struct IsectRayPrecalc *isect_precalc, - const float v0[3], const float v1[3], const float v2[3], - float *r_dist, float r_uv[2]); +void isect_ray_tri_watertight_v3_precalc(struct IsectRayPrecalc *isect_precalc, + const float ray_direction[3]); +bool isect_ray_tri_watertight_v3(const float ray_origin[3], + const struct IsectRayPrecalc *isect_precalc, + const float v0[3], + const float v1[3], + const float v2[3], + float *r_dist, + float r_uv[2]); /* slower version which calculates IsectRayPrecalc each time */ -bool isect_ray_tri_watertight_v3_simple( - const float ray_origin[3], const float ray_direction[3], - const float v0[3], const float v1[3], const float v2[3], - float *r_lambda, float r_uv[2]); - -bool isect_ray_seg_v2( - const float ray_origin[2], const float ray_direction[2], - const float v0[2], const float v1[2], - float *r_lambda, float *r_u); - -bool isect_ray_seg_v3( - const float ray_origin[3], const float ray_direction[3], - const float v0[3], const float v1[3], - float *r_lambda); +bool isect_ray_tri_watertight_v3_simple(const float ray_origin[3], + const float ray_direction[3], + const float v0[3], + const float v1[3], + const float v2[3], + float *r_lambda, + float r_uv[2]); + +bool isect_ray_seg_v2(const float ray_origin[2], + const float ray_direction[2], + const float v0[2], + const float v1[2], + float *r_lambda, + float *r_u); + +bool isect_ray_seg_v3(const float ray_origin[3], + const float ray_direction[3], + const float v0[3], + const float v1[3], + float *r_lambda); /* point in polygon */ -bool isect_point_poly_v2(const float pt[2], const float verts[][2], const unsigned int nr, const bool use_holes); -bool isect_point_poly_v2_int(const int pt[2], const int verts[][2], const unsigned int nr, const bool use_holes); - -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]); - -int isect_point_tri_v2(const float pt[2], const float v1[2], const float v2[2], const float v3[2]); -bool isect_point_tri_v2_cw(const float pt[2], const float v1[2], const float v2[2], const float v3[2]); -int isect_point_tri_v2_int(const int x1, const int y1, const int x2, const int y2, const int a, const int b); -bool isect_point_tri_prism_v3(const float p[3], const float v1[3], const float v2[3], const float v3[3]); -bool isect_point_tri_v3( - const float p[3], const float v1[3], const float v2[3], const float v3[3], - float r_isect_co[3]); +bool isect_point_poly_v2(const float pt[2], + const float verts[][2], + const unsigned int nr, + const bool use_holes); +bool isect_point_poly_v2_int(const int pt[2], + const int verts[][2], + const unsigned int nr, + const bool use_holes); + +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]); + +int isect_point_tri_v2(const float pt[2], const float v1[2], const float v2[2], const float v3[2]); +bool isect_point_tri_v2_cw(const float pt[2], + const float v1[2], + const float v2[2], + const float v3[2]); +int isect_point_tri_v2_int( + const int x1, const int y1, const int x2, const int y2, const int a, const int b); +bool isect_point_tri_prism_v3(const float p[3], + const float v1[3], + const float v2[3], + const float v3[3]); +bool isect_point_tri_v3(const float p[3], + const float v1[3], + const float v2[3], + const float v3[3], + float r_isect_co[3]); /* axis-aligned bounding box */ -bool isect_aabb_aabb_v3(const float min1[3], const float max1[3], const float min2[3], const float max2[3]); +bool isect_aabb_aabb_v3(const float min1[3], + const float max1[3], + const float min2[3], + const float max2[3]); struct IsectRayAABB_Precalc { - float ray_origin[3]; - float ray_inv_dir[3]; - int sign[3]; + float ray_origin[3]; + float ray_inv_dir[3]; + int sign[3]; }; -void isect_ray_aabb_v3_precalc( - struct IsectRayAABB_Precalc *data, - const float ray_origin[3], const float ray_direction[3]); -bool isect_ray_aabb_v3( - const struct IsectRayAABB_Precalc *data, - const float bb_min[3], const float bb_max[3], float *tmin); -bool isect_ray_aabb_v3_simple( - const float orig[3], const float dir[3], - const float bb_min[3], const float bb_max[3], - float *tmin, float *tmax); +void isect_ray_aabb_v3_precalc(struct IsectRayAABB_Precalc *data, + const float ray_origin[3], + const float ray_direction[3]); +bool isect_ray_aabb_v3(const struct IsectRayAABB_Precalc *data, + const float bb_min[3], + const float bb_max[3], + float *tmin); +bool isect_ray_aabb_v3_simple(const float orig[3], + const float dir[3], + const float bb_min[3], + const float bb_max[3], + float *tmin, + float *tmax); /* other */ -#define ISECT_AABB_PLANE_BEHIND_ANY 0 -#define ISECT_AABB_PLANE_CROSS_ANY 1 +#define ISECT_AABB_PLANE_BEHIND_ANY 0 +#define ISECT_AABB_PLANE_CROSS_ANY 1 #define ISECT_AABB_PLANE_IN_FRONT_ALL 2 -int isect_aabb_planes_v3( - const float (*planes)[4], const int totplane, - const float bbmin[3], const float bbmax[3]); +int isect_aabb_planes_v3(const float (*planes)[4], + const int totplane, + const float bbmin[3], + const float bbmax[3]); -bool isect_sweeping_sphere_tri_v3(const float p1[3], const float p2[3], const float radius, - const float v0[3], const float v1[3], const float v2[3], float *r_lambda, float ipoint[3]); +bool isect_sweeping_sphere_tri_v3(const float p1[3], + const float p2[3], + const float radius, + const float v0[3], + const float v1[3], + const float v2[3], + float *r_lambda, + float ipoint[3]); bool clip_segment_v3_plane( - const float p1[3], const float p2[3], const float plane[4], - float r_p1[3], float r_p2[3]); -bool clip_segment_v3_plane_n( - const float p1[3], const float p2[3], const float plane_array[][4], const int plane_tot, - float r_p1[3], float r_p2[3]); + const float p1[3], const float p2[3], const float plane[4], float r_p1[3], float r_p2[3]); +bool clip_segment_v3_plane_n(const float p1[3], + const float p2[3], + const float plane_array[][4], + const int plane_tot, + float r_p1[3], + float r_p2[3]); bool point_in_slice_seg(float p[3], float l1[3], float l2[3]); /****************************** Interpolation ********************************/ -void interp_weights_tri_v3(float w[3], const float a[3], const float b[3], const float c[3], const float p[3]); -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]); +void interp_weights_tri_v3( + float w[3], const float a[3], const float b[3], const float c[3], const float p[3]); +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]); void interp_weights_poly_v3(float w[], float v[][3], const int n, const float co[3]); void interp_weights_poly_v2(float w[], float v[][2], const int n, const float co[2]); -void interp_cubic_v3(float x[3], float v[3], - const float x1[3], const float v1[3], const float x2[3], const float v2[3], const float t); +void interp_cubic_v3(float x[3], + float v[3], + const float x1[3], + const float v1[3], + const float x2[3], + const float v2[3], + const float t); int interp_sparse_array(float *array, const int list_size, const float invalid); -void transform_point_by_tri_v3( - float pt_tar[3], float const pt_src[3], - const float tri_tar_p1[3], const float tri_tar_p2[3], const float tri_tar_p3[3], - const float tri_src_p1[3], const float tri_src_p2[3], const float tri_src_p3[3]); -void transform_point_by_seg_v3( - float p_dst[3], const float p_src[3], - const float l_dst_p1[3], const float l_dst_p2[3], - const float l_src_p1[3], const float l_src_p2[3]); +void transform_point_by_tri_v3(float pt_tar[3], + float const pt_src[3], + const float tri_tar_p1[3], + const float tri_tar_p2[3], + const float tri_tar_p3[3], + const float tri_src_p1[3], + const float tri_src_p2[3], + const float tri_src_p3[3]); +void transform_point_by_seg_v3(float p_dst[3], + const float p_src[3], + const float l_dst_p1[3], + const float l_dst_p2[3], + const float l_src_p1[3], + const float l_src_p2[3]); void barycentric_weights_v2( - const float v1[2], const float v2[2], const float v3[2], - const float co[2], float w[3]); + const float v1[2], const float v2[2], const float v3[2], const float co[2], float w[3]); void barycentric_weights_v2_clamped( - const float v1[2], const float v2[2], const float v3[2], - const float co[2], float w[3]); + const float v1[2], const float v2[2], const float v3[2], const float co[2], float w[3]); void barycentric_weights_v2_persp( - const float v1[4], const float v2[4], const float v3[4], - const float co[2], float w[3]); -void barycentric_weights_v2_quad( - const float v1[2], const float v2[2], const float v3[2], const float v4[2], - const float co[2], float w[4]); - -bool barycentric_coords_v2(const float v1[2], const float v2[2], const float v3[2], const float co[2], float w[3]); + const float v1[4], const float v2[4], const float v3[4], const float co[2], float w[3]); +void barycentric_weights_v2_quad(const float v1[2], + const float v2[2], + const float v3[2], + const float v4[2], + const float co[2], + float w[4]); + +bool barycentric_coords_v2( + const float v1[2], const float v2[2], const float v3[2], const float co[2], float w[3]); int barycentric_inside_triangle_v2(const float w[3]); -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]); -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]); -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]); -void resolve_quad_uv_v2_deriv(float r_uv[2], float r_deriv[2][2], - const float st[2], const float st0[2], const float st1[2], const float st2[2], const float st3[2]); -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]); +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]); +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]); +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]); +void resolve_quad_uv_v2_deriv(float r_uv[2], + float r_deriv[2][2], + const float st[2], + const float st0[2], + const float st1[2], + const float st2[2], + const float st3[2]); +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]); /* use to find the point of a UV on a face */ void interp_bilinear_quad_v3(float data[4][3], float u, float v, float res[3]); @@ -426,67 +592,103 @@ void interp_barycentric_tri_v3(float data[3][3], float u, float v, float res[3]) /***************************** View & Projection *****************************/ -void lookat_m4(float mat[4][4], float vx, float vy, - float vz, float px, float py, float pz, float twist); -void polarview_m4(float mat[4][4], float dist, float azimuth, - float incidence, float twist); - -void perspective_m4(float mat[4][4], const float left, const float right, - const float bottom, const float top, const float nearClip, const float farClip); -void orthographic_m4(float mat[4][4], const float left, const float right, - const float bottom, const float top, const float nearClip, const float farClip); -void window_translate_m4(float winmat[4][4], float perspmat[4][4], - const float x, const float y); - -void planes_from_projmat(float mat[4][4], float left[4], float right[4], float top[4], float bottom[4], - float front[4], float back[4]); +void lookat_m4( + float mat[4][4], float vx, float vy, float vz, float px, float py, float pz, float twist); +void polarview_m4(float mat[4][4], float dist, float azimuth, float incidence, float twist); + +void perspective_m4(float mat[4][4], + const float left, + const float right, + const float bottom, + const float top, + const float nearClip, + const float farClip); +void orthographic_m4(float mat[4][4], + const float left, + const float right, + const float bottom, + const float top, + const float nearClip, + const float farClip); +void window_translate_m4(float winmat[4][4], float perspmat[4][4], const float x, const float y); + +void planes_from_projmat(float mat[4][4], + float left[4], + float right[4], + float top[4], + float bottom[4], + float front[4], + float back[4]); void projmat_dimensions(const float projmat[4][4], - float *r_left, float *r_right, - float *r_bottom, float *r_top, - float *r_near, float *r_far); + float *r_left, + float *r_right, + float *r_bottom, + float *r_top, + float *r_near, + float *r_far); -int box_clip_bounds_m4(float boundbox[2][3], - const float bounds[4], float winmat[4][4]); -void box_minmax_bounds_m4(float min[3], float max[3], - float boundbox[2][3], float mat[4][4]); +int box_clip_bounds_m4(float boundbox[2][3], const float bounds[4], float winmat[4][4]); +void box_minmax_bounds_m4(float min[3], float max[3], float boundbox[2][3], float mat[4][4]); /********************************** Mapping **********************************/ void map_to_tube(float *r_u, float *r_v, const float x, const float y, const float z); void map_to_sphere(float *r_u, float *r_v, const float x, const float y, const float z); void map_to_plane_v2_v3v3(float r_co[2], const float co[3], const float no[3]); -void map_to_plane_axis_angle_v2_v3v3fl(float r_co[2], const float co[3], const float axis[3], const float angle); +void map_to_plane_axis_angle_v2_v3v3fl(float r_co[2], + const float co[3], + const float axis[3], + const float angle); /********************************** Normals **********************************/ -void accumulate_vertex_normals_tri_v3( - float n1[3], float n2[3], float n3[3], - const float f_no[3], - const float co1[3], const float co2[3], const float co3[3]); - -void accumulate_vertex_normals_v3( - float n1[3], float n2[3], float n3[3], float n4[3], - const float f_no[3], - const float co1[3], const float co2[3], const float co3[3], const float co4[3]); - -void accumulate_vertex_normals_poly_v3( - float **vertnos, const float polyno[3], - const float **vertcos, float vdiffs[][3], const int nverts); +void accumulate_vertex_normals_tri_v3(float n1[3], + float n2[3], + float n3[3], + const float f_no[3], + const float co1[3], + const float co2[3], + const float co3[3]); + +void accumulate_vertex_normals_v3(float n1[3], + float n2[3], + float n3[3], + float n4[3], + const float f_no[3], + const float co1[3], + const float co2[3], + const float co3[3], + const float co4[3]); + +void accumulate_vertex_normals_poly_v3(float **vertnos, + const float polyno[3], + const float **vertcos, + float vdiffs[][3], + const int nverts); /********************************* Tangents **********************************/ -void tangent_from_uv_v3( - const float uv1[2], const float uv2[2], const float uv3[2], - const float co1[3], const float co2[3], const float co3[3], - const float n[3], - float r_tang[3]); +void tangent_from_uv_v3(const float uv1[2], + const float uv2[2], + const float uv3[2], + const float co1[3], + const float co2[3], + const float co3[3], + const float n[3], + float r_tang[3]); /******************************** Vector Clouds ******************************/ -void vcloud_estimate_transform_v3( - const int list_size, const float (*pos)[3], const float *weight, const float (*rpos)[3], const float *rweight, - float lloc[3], float rloc[3], float lrot[3][3], float lscale[3][3]); +void vcloud_estimate_transform_v3(const int list_size, + const float (*pos)[3], + const float *weight, + const float (*rpos)[3], + const float *rweight, + float lloc[3], + float rloc[3], + float lrot[3][3], + float lscale[3][3]); /****************************** Spherical Harmonics *************************/ @@ -508,21 +710,33 @@ MINLINE void madd_sh_shfl(float r[9], const float sh[3], const float f); /********************************* Form Factor *******************************/ -float form_factor_quad(const float p[3], const float n[3], - const float q0[3], const float q1[3], const float q2[3], const float q3[3]); -bool form_factor_visible_quad(const float p[3], const float n[3], - const float v0[3], const float v1[3], const float v2[3], - float q0[3], float q1[3], float q2[3], float q3[3]); -float form_factor_hemi_poly(float p[3], float n[3], - float v1[3], float v2[3], float v3[3], float v4[3]); +float form_factor_quad(const float p[3], + const float n[3], + const float q0[3], + const float q1[3], + const float q2[3], + const float q3[3]); +bool form_factor_visible_quad(const float p[3], + const float n[3], + const float v0[3], + const float v1[3], + const float v2[3], + float q0[3], + float q1[3], + float q2[3], + float q3[3]); +float form_factor_hemi_poly( + float p[3], float n[3], float v1[3], float v2[3], float v3[3], float v4[3]); void axis_dominant_v3_to_m3_negate(float r_mat[3][3], const float normal[3]); void axis_dominant_v3_to_m3(float r_mat[3][3], const float normal[3]); -MINLINE void axis_dominant_v3(int *r_axis_a, int *r_axis_b, const float axis[3]); -MINLINE float axis_dominant_v3_max(int *r_axis_a, int *r_axis_b, const float axis[3]) ATTR_WARN_UNUSED_RESULT; -MINLINE int axis_dominant_v3_single(const float vec[3]); -MINLINE int axis_dominant_v3_ortho_single(const float vec[3]); +MINLINE void axis_dominant_v3(int *r_axis_a, int *r_axis_b, const float axis[3]); +MINLINE float axis_dominant_v3_max(int *r_axis_a, + int *r_axis_b, + const float axis[3]) ATTR_WARN_UNUSED_RESULT; +MINLINE int axis_dominant_v3_single(const float vec[3]); +MINLINE int axis_dominant_v3_ortho_single(const float vec[3]); MINLINE int max_axis_v3(const float vec[3]); MINLINE int min_axis_v3(const float vec[3]); @@ -542,7 +756,7 @@ float cubic_tangent_factor_circle_v3(const float tan_l[3], const float tan_r[3]) /**************************** Inline Definitions ******************************/ #if BLI_MATH_DO_INLINE -#include "intern/math_geom_inline.c" +# include "intern/math_geom_inline.c" #endif #ifdef BLI_MATH_GCC_WARN_PRAGMA |