Blender V5.0
BLI_math_matrix.h
Go to the documentation of this file.
1/* SPDX-FileCopyrightText: 2001-2002 NaN Holding BV. All rights reserved.
2 *
3 * SPDX-License-Identifier: GPL-2.0-or-later */
4
5#pragma once
6
10
11#include "BLI_compiler_attrs.h"
12
13/* -------------------------------------------------------------------- */
16
17void zero_m3(float m[3][3]);
18void zero_m4(float m[4][4]);
19
20void unit_m2(float m[2][2]);
21void unit_m3(float m[3][3]);
22void unit_m4(float m[4][4]);
23void unit_m4_db(double m[4][4]);
24
25void copy_m2_m2(float m1[2][2], const float m2[2][2]);
26void copy_m3_m3(float m1[3][3], const float m2[3][3]);
27void copy_m4_m4(float m1[4][4], const float m2[4][4]);
28void copy_m3_m4(float m1[3][3], const float m2[4][4]);
29void copy_m4_m3(float m1[4][4], const float m2[3][3]);
30
31void copy_m4_m4_db(double m1[4][4], const double m2[4][4]);
32
33/* double->float */
34
35void copy_m3_m3d(float m1[3][3], const double m2[3][3]);
36
37/* float->double */
38
39void copy_m3d_m3(double m1[3][3], const float m2[3][3]);
40void copy_m4d_m4(double m1[4][4], const float m2[4][4]);
41
42void swap_m4m4(float m1[4][4], float m2[4][4]);
43
45void shuffle_m4(float R[4][4], const int index[4]);
46
48
49/* -------------------------------------------------------------------- */
52
53void add_m3_m3m3(float R[3][3], const float A[3][3], const float B[3][3]);
54void add_m4_m4m4(float R[4][4], const float A[4][4], const float B[4][4]);
55
56void madd_m3_m3m3fl(float R[3][3], const float A[3][3], const float B[3][3], float f);
57void madd_m4_m4m4fl(float R[4][4], const float A[4][4], const float B[4][4], float f);
58
59void sub_m3_m3m3(float R[3][3], const float A[3][3], const float B[3][3]);
60
61void mul_m3_m3m3(float R[3][3], const float A[3][3], const float B[3][3]);
62void mul_m4_m3m4(float R[4][4], const float A[3][3], const float B[4][4]);
63void mul_m4_m4m3(float R[4][4], const float A[4][4], const float B[3][3]);
64void mul_m4_m4m4(float R[4][4], const float A[4][4], const float B[4][4]);
68void mul_m3_m3m4(float R[3][3], const float A[3][3], const float B[4][4]);
72void mul_m3_m4m3(float R[3][3], const float A[4][4], const float B[3][3]);
73void mul_m3_m4m4(float R[3][3], const float A[4][4], const float B[4][4]);
74
80void mul_m3_m3_pre(float R[3][3], const float A[3][3]);
81void mul_m3_m3_post(float R[3][3], const float B[3][3]);
82void mul_m4db_m4db_m4fl(double R[4][4], const double A[4][4], const float B[4][4]);
83void mul_m4_m4_pre(float R[4][4], const float A[4][4]);
84void mul_m4_m4_post(float R[4][4], const float B[4][4]);
85
86/* Implement #mul_m3_series macro. */
87
88void _va_mul_m3_series_3(float r[3][3], const float m1[3][3], const float m2[3][3]) ATTR_NONNULL();
89void _va_mul_m3_series_4(float r[3][3],
90 const float m1[3][3],
91 const float m2[3][3],
92 const float m3[3][3]) ATTR_NONNULL();
93void _va_mul_m3_series_5(float r[3][3],
94 const float m1[3][3],
95 const float m2[3][3],
96 const float m3[3][3],
97 const float m4[3][3]) ATTR_NONNULL();
98void _va_mul_m3_series_6(float r[3][3],
99 const float m1[3][3],
100 const float m2[3][3],
101 const float m3[3][3],
102 const float m4[3][3],
103 const float m5[3][3]) ATTR_NONNULL();
104void _va_mul_m3_series_7(float r[3][3],
105 const float m1[3][3],
106 const float m2[3][3],
107 const float m3[3][3],
108 const float m4[3][3],
109 const float m5[3][3],
110 const float m6[3][3]) ATTR_NONNULL();
111void _va_mul_m3_series_8(float r[3][3],
112 const float m1[3][3],
113 const float m2[3][3],
114 const float m3[3][3],
115 const float m4[3][3],
116 const float m5[3][3],
117 const float m6[3][3],
118 const float m7[3][3]) ATTR_NONNULL();
119void _va_mul_m3_series_9(float r[3][3],
120 const float m1[3][3],
121 const float m2[3][3],
122 const float m3[3][3],
123 const float m4[3][3],
124 const float m5[3][3],
125 const float m6[3][3],
126 const float m7[3][3],
127 const float m8[3][3]) ATTR_NONNULL();
128
129/* Implement #mul_m4_series macro. */
130
131void _va_mul_m4_series_3(float r[4][4], const float m1[4][4], const float m2[4][4]) ATTR_NONNULL();
132void _va_mul_m4_series_4(float r[4][4],
133 const float m1[4][4],
134 const float m2[4][4],
135 const float m3[4][4]) ATTR_NONNULL();
136void _va_mul_m4_series_5(float r[4][4],
137 const float m1[4][4],
138 const float m2[4][4],
139 const float m3[4][4],
140 const float m4[4][4]) ATTR_NONNULL();
141void _va_mul_m4_series_6(float r[4][4],
142 const float m1[4][4],
143 const float m2[4][4],
144 const float m3[4][4],
145 const float m4[4][4],
146 const float m5[4][4]) ATTR_NONNULL();
147void _va_mul_m4_series_7(float r[4][4],
148 const float m1[4][4],
149 const float m2[4][4],
150 const float m3[4][4],
151 const float m4[4][4],
152 const float m5[4][4],
153 const float m6[4][4]) ATTR_NONNULL();
154void _va_mul_m4_series_8(float r[4][4],
155 const float m1[4][4],
156 const float m2[4][4],
157 const float m3[4][4],
158 const float m4[4][4],
159 const float m5[4][4],
160 const float m6[4][4],
161 const float m7[4][4]) ATTR_NONNULL();
162void _va_mul_m4_series_9(float r[4][4],
163 const float m1[4][4],
164 const float m2[4][4],
165 const float m3[4][4],
166 const float m4[4][4],
167 const float m5[4][4],
168 const float m6[4][4],
169 const float m7[4][4],
170 const float m8[4][4]) ATTR_NONNULL();
171
172#define mul_m3_series(...) VA_NARGS_CALL_OVERLOAD(_va_mul_m3_series_, __VA_ARGS__)
173#define mul_m4_series(...) VA_NARGS_CALL_OVERLOAD(_va_mul_m4_series_, __VA_ARGS__)
174
175void mul_m4_v3(const float M[4][4], float r[3]);
176void mul_v3_m4v3(float r[3], const float mat[4][4], const float vec[3]);
177void mul_v3_m4v3_db(double r[3], const double mat[4][4], const double vec[3]);
178void mul_v4_m4v3_db(double r[4], const double mat[4][4], const double vec[3]);
179void mul_v2_m4v3(float r[2], const float mat[4][4], const float vec[3]);
180void mul_v2_m2v2(float r[2], const float mat[2][2], const float vec[2]);
181void mul_m2_v2(const float mat[2][2], float vec[2]);
183void mul_mat3_m4_v3(const float mat[4][4], float r[3]);
184void mul_v3_mat3_m4v3(float r[3], const float mat[4][4], const float vec[3]);
185void mul_v3_mat3_m4v3_db(double r[3], const double mat[4][4], const double vec[3]);
186void mul_m4_v4(const float mat[4][4], float r[4]);
187void mul_v4_m4v4(float r[4], const float mat[4][4], const float v[4]);
188void mul_v4_m4v3(float r[4], const float M[4][4], const float v[3]); /* v has implicit w = 1.0f */
189void mul_project_m4_v3(const float mat[4][4], float vec[3]);
190void mul_v3_project_m4_v3(float r[3], const float mat[4][4], const float vec[3]);
191void mul_v2_project_m4_v3(float r[2], const float mat[4][4], const float vec[3]);
192
193void mul_m3_v2(const float m[3][3], float r[2]);
194void mul_v2_m3v2(float r[2], const float m[3][3], const float v[2]);
195void mul_m3_v3(const float M[3][3], float r[3]);
196void mul_v3_m3v3(float r[3], const float M[3][3], const float a[3]);
197void mul_v2_m3v3(float r[2], const float M[3][3], const float a[3]);
198void mul_transposed_m3_v3(const float M[3][3], float r[3]);
199void mul_transposed_mat3_m4_v3(const float M[4][4], float r[3]);
200
209void mul_m4_m4m4_aligned_scale(float R[4][4], const float A[4][4], const float B[4][4]);
213void mul_m4_m4m4_split_channels(float R[4][4], const float A[4][4], const float B[4][4]);
214
215void mul_m3_fl(float R[3][3], float f);
216void mul_m4_fl(float R[4][4], float f);
217void mul_mat3_m4_fl(float R[4][4], float f);
218
219void negate_m3(float R[3][3]);
220void negate_mat3_m4(float R[4][4]);
221void negate_m4(float R[4][4]);
222
223bool invert_m3(float mat[3][3]);
224bool invert_m2_m2(float inverse[2][2], const float mat[2][2]);
225bool invert_m3_m3(float inverse[3][3], const float mat[3][3]);
226bool invert_m4(float mat[4][4]);
227bool invert_m4_m4(float inverse[4][4], const float mat[4][4]);
238bool invert_m4_m4_fallback(float inverse[4][4], const float mat[4][4]);
239
240/* Double matrix functions (no mixing types). */
241
242void mul_v3_m3v3_db(double r[3], const double M[3][3], const double a[3]);
243void mul_m3_v3_db(const double M[3][3], double r[3]);
244
246
247/* -------------------------------------------------------------------- */
250
251void transpose_m3(float R[3][3]);
252void transpose_m3_m3(float R[3][3], const float M[3][3]);
256void transpose_m3_m4(float R[3][3], const float M[4][4]);
257void transpose_m4(float R[4][4]);
258void transpose_m4_m4(float R[4][4], const float M[4][4]);
259
260bool compare_m4m4(const float mat1[4][4], const float mat2[4][4], float limit);
261
262void normalize_m2_m2(float R[2][2], const float M[2][2]) ATTR_NONNULL();
263void normalize_m3(float R[3][3]) ATTR_NONNULL();
264void normalize_m3_m3(float R[3][3], const float M[3][3]) ATTR_NONNULL();
265void normalize_m4_ex(float R[4][4], float r_scale[3]) ATTR_NONNULL();
266void normalize_m4(float R[4][4]) ATTR_NONNULL();
267void normalize_m4_m4(float rmat[4][4], const float mat[4][4]) ATTR_NONNULL();
268
274void orthogonalize_m3(float R[3][3], int axis);
280void orthogonalize_m4(float R[4][4], int axis);
281
290void orthogonalize_m4_stable(float R[4][4], int axis, bool normalize);
291
292bool orthogonalize_m3_zero_axes(float m[3][3], float unit_length);
293bool orthogonalize_m4_zero_axes(float m[4][4], float unit_length);
294
295bool is_orthogonal_m3(const float m[3][3]);
296bool is_orthogonal_m4(const float m[4][4]);
297bool is_orthonormal_m3(const float m[3][3]);
298bool is_orthonormal_m4(const float m[4][4]);
299
300bool is_identity_m4(const float m[4][4]);
301bool is_uniform_scaled_m3(const float m[3][3]);
302bool is_uniform_scaled_m4(const float m[4][4]);
303
304/* NOTE: 'adjoint' here means the adjugate (adjunct, "classical adjoint") matrix!
305 * Nowadays 'adjoint' usually refers to the conjugate transpose,
306 * which for real-valued matrices is simply the transpose. */
307
308void adjoint_m2_m2(float R[2][2], const float M[2][2]);
309void adjoint_m3_m3(float R[3][3], const float M[3][3]);
310void adjoint_m4_m4(float R[4][4], const float M[4][4]);
311
312float determinant_m2(float a, float b, float c, float d);
313float determinant_m3(
314 float a1, float a2, float a3, float b1, float b2, float b3, float c1, float c2, float c3);
315float determinant_m3_array(const float m[3][3]);
316float determinant_m4_mat3_array(const float m[4][4]);
317float determinant_m4(const float m[4][4]);
318
319#define PSEUDOINVERSE_EPSILON 1e-8f
320
328void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4]);
329void pseudoinverse_m4_m4(float inverse[4][4], const float mat[4][4], float epsilon);
330void pseudoinverse_m3_m3(float inverse[3][3], const float mat[3][3], float epsilon);
331
332bool has_zero_axis_m4(const float matrix[4][4]);
333
334void invert_m4_m4_safe(float inverse[4][4], const float mat[4][4]);
335
336void invert_m3_m3_safe_ortho(float inverse[3][3], const float mat[3][3]);
343void invert_m4_m4_safe_ortho(float inverse[4][4], const float mat[4][4]);
344
346
347/* -------------------------------------------------------------------- */
350
351void scale_m3_fl(float R[3][3], float scale);
352void scale_m4_fl(float R[4][4], float scale);
353
359float mat4_to_volume_scale(const float mat[4][4]);
360
366float mat3_to_scale(const float mat[3][3]);
367float mat4_to_scale(const float mat[4][4]);
368
369void size_to_mat3(float R[3][3], const float size[3]);
370void size_to_mat4(float R[4][4], const float size[3]);
371
372void mat3_to_size(float size[3], const float M[3][3]);
373void mat4_to_size(float size[3], const float M[4][4]);
374
384float mat4_to_size_max_axis(const float M[4][4]);
385
390void mat4_to_size_fix_shear(float size[3], const float M[4][4]);
391
392void translate_m4(float mat[4][4], float Tx, float Ty, float Tz);
400void rotate_m4(float mat[4][4], char axis, float angle);
402void rescale_m4(float mat[4][4], const float scale[3]);
411void transform_pivot_set_m4(float mat[4][4], const float pivot[3]);
412
417void mat3_to_rot_size(float rot[3][3], float size[3], const float mat3[3][3]);
422void mat4_to_loc_rot_size(float loc[3], float rot[3][3], float size[3], const float wmat[4][4]);
423void mat4_to_loc_quat(float loc[3], float quat[4], const float wmat[4][4]);
424void mat4_decompose(float loc[3], float quat[4], float size[3], const float wmat[4][4]);
425
426void mat3_polar_decompose(const float mat3[3][3], float r_U[3][3], float r_P[3][3]);
427
432void loc_rot_size_to_mat4(float R[4][4],
433 const float loc[3],
434 const float rot[3][3],
435 const float size[3]);
442void loc_eul_size_to_mat4(float R[4][4],
443 const float loc[3],
444 const float eul[3],
445 const float size[3]);
451 float R[4][4], const float loc[3], const float eul[3], const float size[3], short order);
456void loc_quat_size_to_mat4(float R[4][4],
457 const float loc[3],
458 const float quat[4],
459 const float size[3]);
460
461void blend_m3_m3m3(float out[3][3], const float dst[3][3], const float src[3][3], float srcweight);
462void blend_m4_m4m4(float out[4][4], const float dst[4][4], const float src[4][4], float srcweight);
463
480void interp_m3_m3m3(float R[3][3], const float A[3][3], const float B[3][3], float t);
490void interp_m4_m4m4(float R[4][4], const float A[4][4], const float B[4][4], float t);
491
502bool is_negative_m3(const float mat[3][3]);
504bool is_negative_m4(const float mat[4][4]);
505
506bool is_zero_m4(const float mat[4][4]);
507
508bool equals_m3m3(const float mat1[3][3], const float mat2[3][3]);
509bool equals_m4m4(const float mat1[4][4], const float mat2[4][4]);
510
529typedef struct SpaceTransform {
530 float local2target[4][4];
531 float target2local[4][4];
532
534
548 const float local[4][4],
549 const float target[4][4]);
563 const float local[4][4],
564 const float target[4][4]);
565void BLI_space_transform_apply(const struct SpaceTransform *data, float co[3]);
566void BLI_space_transform_invert(const struct SpaceTransform *data, float co[3]);
567void BLI_space_transform_apply_normal(const struct SpaceTransform *data, float no[3]);
569
570#define BLI_SPACE_TRANSFORM_SETUP(data, local, target) \
571 BLI_space_transform_from_matrices( \
572 (data), (local)->object_to_world().ptr(), (target)->object_to_world().ptr())
573
575
576/* -------------------------------------------------------------------- */
579
580void print_m3(const char *str, const float m[3][3]);
581void print_m4(const char *str, const float m[4][4]);
582
583#define print_m3_id(M) print_m3(STRINGIFY(M), M)
584#define print_m4_id(M) print_m4(STRINGIFY(M), M)
585
#define ATTR_NONNULL(...)
void mul_v4_m4v3(float r[4], const float M[4][4], const float v[3])
bool invert_m2_m2(float inverse[2][2], const float mat[2][2])
void mul_v2_m4v3(float r[2], const float mat[4][4], const float vec[3])
float mat3_to_scale(const float mat[3][3])
float mat4_to_volume_scale(const float mat[4][4])
bool is_negative_m3(const float mat[3][3])
void mul_v2_project_m4_v3(float r[2], const float mat[4][4], const float vec[3])
void BLI_space_transform_from_matrices(struct SpaceTransform *data, const float local[4][4], const float target[4][4])
bool is_orthogonal_m3(const float m[3][3])
void sub_m3_m3m3(float R[3][3], const float A[3][3], const float B[3][3])
void unit_m2(float m[2][2])
void mul_v4_m4v3_db(double r[4], const double mat[4][4], const double vec[3])
void negate_m3(float R[3][3])
void interp_m4_m4m4(float R[4][4], const float A[4][4], const float B[4][4], float t)
float mat4_to_scale(const float mat[4][4])
void orthogonalize_m4(float R[4][4], int axis)
void mul_v3_project_m4_v3(float r[3], const float mat[4][4], const float vec[3])
void mul_m3_v3(const float M[3][3], float r[3])
void madd_m4_m4m4fl(float R[4][4], const float A[4][4], const float B[4][4], float f)
void zero_m4(float m[4][4])
void BLI_space_transform_apply_normal(const struct SpaceTransform *data, float no[3])
void mul_m4_fl(float R[4][4], float f)
void mul_m4_m4m4(float R[4][4], const float A[4][4], const float B[4][4])
void mul_mat3_m4_fl(float R[4][4], float f)
void mul_v4_m4v4(float r[4], const float mat[4][4], const float v[4])
void mul_m3_m3_pre(float R[3][3], const float A[3][3])
void pseudoinverse_m3_m3(float inverse[3][3], const float mat[3][3], float epsilon)
void BLI_space_transform_apply(const struct SpaceTransform *data, float co[3])
void size_to_mat3(float R[3][3], const float size[3])
void _va_mul_m4_series_7(float r[4][4], const float m1[4][4], const float m2[4][4], const float m3[4][4], const float m4[4][4], const float m5[4][4], const float m6[4][4]) ATTR_NONNULL()
void normalize_m3_m3(float R[3][3], const float M[3][3]) ATTR_NONNULL()
void mul_v3_m3v3_db(double r[3], const double M[3][3], const double a[3])
void normalize_m2_m2(float R[2][2], const float M[2][2]) ATTR_NONNULL()
void copy_m3_m3(float m1[3][3], const float m2[3][3])
void mat4_decompose(float loc[3], float quat[4], float size[3], const float wmat[4][4])
bool is_orthogonal_m4(const float m[4][4])
void adjoint_m3_m3(float R[3][3], const float M[3][3])
void _va_mul_m3_series_6(float r[3][3], const float m1[3][3], const float m2[3][3], const float m3[3][3], const float m4[3][3], const float m5[3][3]) ATTR_NONNULL()
void unit_m3(float m[3][3])
void mul_m3_v2(const float m[3][3], float r[2])
void mat3_to_rot_size(float rot[3][3], float size[3], const float mat3[3][3])
void _va_mul_m3_series_3(float r[3][3], const float m1[3][3], const float m2[3][3]) ATTR_NONNULL()
void blend_m4_m4m4(float out[4][4], const float dst[4][4], const float src[4][4], float srcweight)
void add_m4_m4m4(float R[4][4], const float A[4][4], const float B[4][4])
void scale_m3_fl(float R[3][3], float scale)
void copy_m3_m4(float m1[3][3], const float m2[4][4])
float determinant_m4_mat3_array(const float m[4][4])
void mul_m4_m4m3(float R[4][4], const float A[4][4], const float B[3][3])
void normalize_m4_m4(float rmat[4][4], const float mat[4][4]) ATTR_NONNULL()
void mul_v3_m4v3_db(double r[3], const double mat[4][4], const double vec[3])
void transpose_m4_m4(float R[4][4], const float M[4][4])
void mul_m4_m4_pre(float R[4][4], const float A[4][4])
void mul_m3_fl(float R[3][3], float f)
void mul_v2_m3v3(float r[2], const float M[3][3], const float a[3])
void mul_m4db_m4db_m4fl(double R[4][4], const double A[4][4], const float B[4][4])
void orthogonalize_m3(float R[3][3], int axis)
void copy_m3_m3d(float m1[3][3], const double m2[3][3])
bool invert_m3_m3(float inverse[3][3], const float mat[3][3])
void mul_m3_m3_post(float R[3][3], const float B[3][3])
void add_m3_m3m3(float R[3][3], const float A[3][3], const float B[3][3])
void copy_m4_m3(float m1[4][4], const float m2[3][3])
bool is_uniform_scaled_m3(const float m[3][3])
void loc_rot_size_to_mat4(float R[4][4], const float loc[3], const float rot[3][3], const float size[3])
void mat4_to_size_fix_shear(float size[3], const float M[4][4])
bool is_uniform_scaled_m4(const float m[4][4])
void mat4_to_loc_rot_size(float loc[3], float rot[3][3], float size[3], const float wmat[4][4])
void svd_m4(float U[4][4], float s[4], float V[4][4], float A_[4][4])
bool is_orthonormal_m3(const float m[3][3])
void zero_m3(float m[3][3])
void translate_m4(float mat[4][4], float Tx, float Ty, float Tz)
bool has_zero_axis_m4(const float matrix[4][4])
void normalize_m3(float R[3][3]) ATTR_NONNULL()
void transform_pivot_set_m4(float mat[4][4], const float pivot[3])
float determinant_m2(float a, float b, float c, float d)
void _va_mul_m4_series_6(float r[4][4], const float m1[4][4], const float m2[4][4], const float m3[4][4], const float m4[4][4], const float m5[4][4]) ATTR_NONNULL()
bool equals_m3m3(const float mat1[3][3], const float mat2[3][3])
void _va_mul_m4_series_5(float r[4][4], const float m1[4][4], const float m2[4][4], const float m3[4][4], const float m4[4][4]) ATTR_NONNULL()
void rescale_m4(float mat[4][4], const float scale[3])
void loc_eulO_size_to_mat4(float R[4][4], const float loc[3], const float eul[3], const float size[3], short order)
bool is_identity_m4(const float m[4][4])
void size_to_mat4(float R[4][4], const float size[3])
void mul_project_m4_v3(const float mat[4][4], float vec[3])
void unit_m4_db(double m[4][4])
void _va_mul_m3_series_4(float r[3][3], const float m1[3][3], const float m2[3][3], const float m3[3][3]) ATTR_NONNULL()
void copy_m4d_m4(double m1[4][4], const float m2[4][4])
void BLI_space_transform_global_from_matrices(struct SpaceTransform *data, const float local[4][4], const float target[4][4])
void mul_m4_v3(const float M[4][4], float r[3])
void invert_m4_m4_safe(float inverse[4][4], const float mat[4][4])
void orthogonalize_m4_stable(float R[4][4], int axis, bool normalize)
void scale_m4_fl(float R[4][4], float scale)
void adjoint_m4_m4(float R[4][4], const float M[4][4])
bool equals_m4m4(const float mat1[4][4], const float mat2[4][4])
void invert_m3_m3_safe_ortho(float inverse[3][3], const float mat[3][3])
void mul_m4_m4m4_split_channels(float R[4][4], const float A[4][4], const float B[4][4])
float determinant_m4(const float m[4][4])
void swap_m4m4(float m1[4][4], float m2[4][4])
void transpose_m3_m4(float R[3][3], const float M[4][4])
void _va_mul_m3_series_8(float r[3][3], const float m1[3][3], const float m2[3][3], const float m3[3][3], const float m4[3][3], const float m5[3][3], const float m6[3][3], const float m7[3][3]) ATTR_NONNULL()
void copy_m2_m2(float m1[2][2], const float m2[2][2])
void mat3_polar_decompose(const float mat3[3][3], float r_U[3][3], float r_P[3][3])
void mul_v2_m2v2(float r[2], const float mat[2][2], const float vec[2])
void _va_mul_m4_series_4(float r[4][4], const float m1[4][4], const float m2[4][4], const float m3[4][4]) ATTR_NONNULL()
void _va_mul_m4_series_9(float r[4][4], const float m1[4][4], const float m2[4][4], const float m3[4][4], const float m4[4][4], const float m5[4][4], const float m6[4][4], const float m7[4][4], const float m8[4][4]) ATTR_NONNULL()
void normalize_m4_ex(float R[4][4], float r_scale[3]) ATTR_NONNULL()
void _va_mul_m3_series_9(float r[3][3], const float m1[3][3], const float m2[3][3], const float m3[3][3], const float m4[3][3], const float m5[3][3], const float m6[3][3], const float m7[3][3], const float m8[3][3]) ATTR_NONNULL()
void _va_mul_m3_series_7(float r[3][3], const float m1[3][3], const float m2[3][3], const float m3[3][3], const float m4[3][3], const float m5[3][3], const float m6[3][3]) ATTR_NONNULL()
void copy_m4_m4(float m1[4][4], const float m2[4][4])
void mat4_to_loc_quat(float loc[3], float quat[4], const float wmat[4][4])
void interp_m3_m3m3(float R[3][3], const float A[3][3], const float B[3][3], float t)
void mul_transposed_mat3_m4_v3(const float M[4][4], float r[3])
void mul_m3_m3m4(float R[3][3], const float A[3][3], const float B[4][4])
void _va_mul_m4_series_3(float r[4][4], const float m1[4][4], const float m2[4][4]) ATTR_NONNULL()
bool is_orthonormal_m4(const float m[4][4])
void mul_m3_v3_db(const double M[3][3], double r[3])
void madd_m3_m3m3fl(float R[3][3], const float A[3][3], const float B[3][3], float f)
float mat4_to_size_max_axis(const float M[4][4])
void mul_v3_m4v3(float r[3], const float mat[4][4], const float vec[3])
float determinant_m3_array(const float m[3][3])
void _va_mul_m4_series_8(float r[4][4], const float m1[4][4], const float m2[4][4], const float m3[4][4], const float m4[4][4], const float m5[4][4], const float m6[4][4], const float m7[4][4]) ATTR_NONNULL()
void shuffle_m4(float R[4][4], const int index[4])
bool is_negative_m4(const float mat[4][4])
void negate_mat3_m4(float R[4][4])
bool invert_m4_m4(float inverse[4][4], const float mat[4][4])
void mul_v2_m3v2(float r[2], const float m[3][3], const float v[2])
void mul_m4_m4m4_aligned_scale(float R[4][4], const float A[4][4], const float B[4][4])
void print_m3(const char *str, const float m[3][3])
void mul_m4_v4(const float mat[4][4], float r[4])
void rotate_m4(float mat[4][4], char axis, float angle)
void loc_quat_size_to_mat4(float R[4][4], const float loc[3], const float quat[4], const float size[3])
void mul_v3_m3v3(float r[3], const float M[3][3], const float a[3])
void loc_eul_size_to_mat4(float R[4][4], const float loc[3], const float eul[3], const float size[3])
bool invert_m4(float mat[4][4])
void BLI_space_transform_invert_normal(const struct SpaceTransform *data, float no[3])
void mat4_to_size(float size[3], const float M[4][4])
void transpose_m3(float R[3][3])
void BLI_space_transform_invert(const struct SpaceTransform *data, float co[3])
float determinant_m3(float a1, float a2, float a3, float b1, float b2, float b3, float c1, float c2, float c3)
void mul_m2_v2(const float mat[2][2], float vec[2])
void invert_m4_m4_safe_ortho(float inverse[4][4], const float mat[4][4])
void transpose_m3_m3(float R[3][3], const float M[3][3])
void blend_m3_m3m3(float out[3][3], const float dst[3][3], const float src[3][3], float srcweight)
bool orthogonalize_m3_zero_axes(float m[3][3], float unit_length)
void negate_m4(float R[4][4])
void mul_transposed_m3_v3(const float M[3][3], float r[3])
bool compare_m4m4(const float mat1[4][4], const float mat2[4][4], float limit)
void copy_m3d_m3(double m1[3][3], const float m2[3][3])
void pseudoinverse_m4_m4(float inverse[4][4], const float mat[4][4], float epsilon)
void mul_v3_mat3_m4v3_db(double r[3], const double mat[4][4], const double vec[3])
void _va_mul_m3_series_5(float r[3][3], const float m1[3][3], const float m2[3][3], const float m3[3][3], const float m4[3][3]) ATTR_NONNULL()
void mul_m3_m3m3(float R[3][3], const float A[3][3], const float B[3][3])
void print_m4(const char *str, const float m[4][4])
bool orthogonalize_m4_zero_axes(float m[4][4], float unit_length)
bool is_zero_m4(const float mat[4][4])
void mat3_to_size(float size[3], const float M[3][3])
void mul_m4_m3m4(float R[4][4], const float A[3][3], const float B[4][4])
void transpose_m4(float R[4][4])
void mul_m4_m4_post(float R[4][4], const float B[4][4])
void mul_v3_mat3_m4v3(float r[3], const float mat[4][4], const float vec[3])
bool invert_m3(float mat[3][3])
void mul_mat3_m4_v3(const float mat[4][4], float r[3])
void mul_m3_m4m4(float R[3][3], const float A[4][4], const float B[4][4])
void copy_m4_m4_db(double m1[4][4], const double m2[4][4])
bool invert_m4_m4_fallback(float inverse[4][4], const float mat[4][4])
void mul_m3_m4m3(float R[3][3], const float A[4][4], const float B[3][3])
void normalize_m4(float R[4][4]) ATTR_NONNULL()
void adjoint_m2_m2(float R[2][2], const float M[2][2])
void unit_m4(float m[4][4])
static double angle(const Eigen::Vector3d &v1, const Eigen::Vector3d &v2)
Definition IK_Math.h:117
#define U
BMesh const char void * data
ATTR_WARN_UNUSED_RESULT const BMVert * v
#define A
static DBVT_INLINE btScalar size(const btDbvtVolume &a)
Definition btDbvt.cpp:52
#define rot(x, k)
#define str(s)
#define out
MatBase< C, R > inverse(MatBase< C, R >) RET
VecBase< float, D > normalize(VecOp< float, D >) RET
#define M
#define B
#define R
float local2target[4][4]
float target2local[4][4]
CCL_NAMESPACE_BEGIN struct Window V