Math Lib: rotate matrix cleanup
authorCampbell Barton <ideasman42@gmail.com>
Fri, 25 Nov 2016 05:20:30 +0000 (16:20 +1100)
committerCampbell Barton <ideasman42@gmail.com>
Fri, 25 Nov 2016 05:39:34 +0000 (16:39 +1100)
- Remove 'rotate_m2', unlike 'rotate_m4' it created a new matrix
  duplicating 'angle_to_mat2' - now used instead.
  (better avoid matching functions having different behavior).

- Add 'axis_angle_to_mat4_single',
  convenience wrapper for 'axis_angle_to_mat3_single'.

- Replace 'unit_m4(), rotate_m4()' with a single call to 'axis_angle_to_mat4_single'.

15 files changed:
source/blender/alembic/intern/abc_transform.cc
source/blender/alembic/intern/abc_util.cc
source/blender/blenkernel/intern/object_dupli.c
source/blender/blenkernel/intern/tracking_stabilize.c
source/blender/blenlib/BLI_math_matrix.h
source/blender/blenlib/BLI_math_rotation.h
source/blender/blenlib/intern/math_geom.c
source/blender/blenlib/intern/math_matrix.c
source/blender/blenlib/intern/math_rotation.c
source/blender/collada/collada_internal.cpp
source/blender/editors/mesh/editmesh_knife.c
source/blender/editors/object/object_warp.c
source/blender/editors/sculpt_paint/sculpt.c
source/blender/editors/uvedit/uvedit_smart_stitch.c
source/blender/modifiers/intern/MOD_screw.c

index 7f8984f9970d303178685d2e4a22c257e1187927..e2fc7674c4e728dfb30c7384280587fa5ff446eb 100644 (file)
@@ -92,8 +92,7 @@ void AbcTransformWriter::do_write()
        /* Only apply rotation to root camera, parenting will propagate it. */
        if (m_object->type == OB_CAMERA && !has_parent_camera(m_object)) {
                float rot_mat[4][4];
-               unit_m4(rot_mat);
-               rotate_m4(rot_mat, 'X', -M_PI_2);
+               axis_angle_to_mat4_single(rot_mat, 'X', -M_PI_2);
                mul_m4_m4m4(mat, mat, rot_mat);
        }
 
index f87d18605d4a228874decf4526cdb61cee503bf8..979a72fcf4e37f8ccd92ed1312d4180b5ffc0a08 100644 (file)
@@ -221,8 +221,7 @@ void convert_matrix(const Imath::M44d &xform, Object *ob,
 
        if (ob->type == OB_CAMERA) {
                float cam_to_yup[4][4];
-               unit_m4(cam_to_yup);
-               rotate_m4(cam_to_yup, 'X', M_PI_2);
+               axis_angle_to_mat4_single(cam_to_yup, 'X', M_PI_2);
                mul_m4_m4m4(r_mat, r_mat, cam_to_yup);
        }
 
index 14cc5ec084901be61a5b169846159157ee2de371..e3b801b319336a749ce4c53f05ab1211331c5a1d 100644 (file)
@@ -642,8 +642,7 @@ static void make_duplis_font(const DupliContext *ctx)
                                float rmat[4][4];
 
                                zero_v3(obmat[3]);
-                               unit_m4(rmat);
-                               rotate_m4(rmat, 'Z', -ct->rot);
+                               axis_angle_to_mat4_single(rmat, 'Z', -ct->rot);
                                mul_m4_m4m4(obmat, obmat, rmat);
                        }
 
index b8949f9a0def29afc6976c7ba69e27a5be241d1d..36b24fbb2dc3819586ade8bd7c1e27251b581f4c 100644 (file)
@@ -587,7 +587,7 @@ static void compensate_rotation_center(const int size, float aspect,
 
        copy_v2_v2(intended_pivot, pivot);
        copy_v2_v2(rotated_pivot, pivot);
-       rotate_m2(rotation_mat, +angle);
+       angle_to_mat2(rotation_mat, +angle);
        sub_v2_v2(rotated_pivot, origin);
        mul_m2v2(rotation_mat, rotated_pivot);
        mul_v2_fl(rotated_pivot, scale);
@@ -967,7 +967,7 @@ static void initialize_track_for_stabilization(StabContext *ctx,
 
        pos[0] *= aspect;
        angle = average_angle - atan2f(pos[1],pos[0]);
-       rotate_m2(local_data->stabilization_rotation_base, angle);
+       angle_to_mat2(local_data->stabilization_rotation_base, angle);
 
        /* Per track baseline value for zoom. */
        len = len_v2(pos) + SCALE_ERROR_LIMIT_BIAS;
index 8124e07dd4796c0013ce1c241ca1a8288cc8143d..d0dfad2a02ffd4e32dae5ff3a4d35e5e424198a6 100644 (file)
@@ -219,7 +219,6 @@ void mat4_to_size(float r[3], float M[4][4]);
 
 void translate_m4(float mat[4][4], float tx, float ty, float tz);
 void rotate_m4(float mat[4][4], const char axis, const float angle);
-void rotate_m2(float mat[2][2], const float angle);
 void transform_pivot_set_m4(float mat[4][4], const float pivot[3]);
 
 void mat3_to_rot_size(float rot[3][3], float size[3], float mat3[3][3]);
index 24c20ee7b5009535f73d9ee08436a9312ddb7483..d60be30e10dbcb78a277a311299aa9d4f5d43b19 100644 (file)
@@ -122,8 +122,9 @@ void mat3_to_axis_angle(float axis[3], float *angle, float M[3][3]);
 void mat4_to_axis_angle(float axis[3], float *angle, float M[4][4]);
 void quat_to_axis_angle(float axis[3], float *angle, const float q[4]);
 
-void axis_angle_to_mat3_single(float R[3][3], const char axis, const float angle);
 void      angle_to_mat2(float R[2][2], const float angle);
+void axis_angle_to_mat3_single(float R[3][3], const char axis, const float angle);
+void axis_angle_to_mat4_single(float R[4][4], const char axis, const float angle);
 
 void axis_angle_to_quat_single(float q[4], const char axis, const float angle);
 
index f31d0935b773c23149353e22e934e6e6d4e2bd72..0790d65345d3bcc703e0999d5401afd3164f7228 100644 (file)
@@ -3917,10 +3917,9 @@ void lookat_m4(float mat[4][4], float vx, float vy, float vz, float px, float py
        float sine, cosine, hyp, hyp1, dx, dy, dz;
        float mat1[4][4];
 
-       unit_m4(mat);
        unit_m4(mat1);
 
-       rotate_m4(mat, 'Z', -twist);
+       axis_angle_to_mat3_single(mat, 'Z', -twist);
 
        dx = px - vx;
        dy = py - vy;
index c9c61d5c8787681d4874916401db6e67d44180e3..7176686df54c3e62a2927bcf8c49937fbf06f2d4 100644 (file)
@@ -1625,6 +1625,13 @@ void translate_m4(float mat[4][4], float Tx, float Ty, float Tz)
        mat[3][2] += (Tx * mat[0][2] + Ty * mat[1][2] + Tz * mat[2][2]);
 }
 
+/**
+ * Rotate a matrix in-place.
+ *
+ * \note To create a new rotation matrix see:
+ * #axis_angle_to_mat4_single, #axis_angle_to_mat3_single, #angle_to_mat2
+ * (axis & angle args are compatible).
+ */
 void rotate_m4(float mat[4][4], const char axis, const float angle)
 {
        int col;
@@ -1665,13 +1672,6 @@ void rotate_m4(float mat[4][4], const char axis, const float angle)
        }
 }
 
-void rotate_m2(float mat[2][2], const float angle)
-{
-       mat[0][0] = mat[1][1] = cosf(angle);
-       mat[0][1] = sinf(angle);
-       mat[1][0] = -mat[0][1];
-}
-
 /**
  * Scale or rotate around a pivot point,
  * a convenience function to avoid having to do inline.
index b285a74b8acc0b12ca4113f4dd165b3b54192112..9b5dad3f7670d600b2db6b4bf6b4fc816ad31603 100644 (file)
@@ -1009,6 +1009,13 @@ void mat4_to_axis_angle(float axis[3], float *angle, float mat[4][4])
        quat_to_axis_angle(axis, angle, q);
 }
 
+void axis_angle_to_mat4_single(float mat[4][4], const char axis, const float angle)
+{
+       float mat3[3][3];
+       axis_angle_to_mat3_single(mat3, axis, angle);
+       copy_m4_m3(mat, mat3);
+}
+
 /* rotation matrix from a single axis */
 void axis_angle_to_mat3_single(float mat[3][3], const char axis, const float angle)
 {
index 38855013ee1bba01e7c170d0b1de2bec1c6bafce..e1a13559b082ce33fa726001f7c6329e54f8f4ba 100644 (file)
 
 UnitConverter::UnitConverter() : unit(), up_axis(COLLADAFW::FileInfo::Z_UP)
 {
-       unit_m4(x_up_mat4);
-       rotate_m4(x_up_mat4, 'Y', -0.5 * M_PI);
-
-       unit_m4(y_up_mat4);
-       rotate_m4(y_up_mat4, 'X', 0.5 * M_PI);
+       axis_angle_to_mat4_single(x_up_mat4, 'Y', -0.5 * M_PI);
+       axis_angle_to_mat4_single(y_up_mat4, 'X', 0.5 * M_PI);
 
        unit_m4(z_up_mat4);
        unit_m4(scale_mat4);
index a84b8d9dcc82f21644b5c82fb2793698aa754a2c..44453d03ade70a1a2dc5399d8a8c04b020a4f2da 100644 (file)
@@ -2140,7 +2140,7 @@ static float snap_v2_angle(float r[2], const float v[2], const float v_ref[2], f
        normalize_v2_v2(v_unit, v);
        angle = angle_signed_v2v2(v_unit, v_ref);
        angle_delta = (roundf(angle / angle_snap) * angle_snap) - angle;
-       rotate_m2(m2, angle_delta);
+       angle_to_mat2(m2, angle_delta);
 
        mul_v2_m2v2(r, m2, v);
        return angle + angle_delta;
index 9f4da87903d7e4884ace16fdfe7c4e73dc0f655e..92b82e2a31b85ca6ff72d5d97e5be167928514a0 100644 (file)
@@ -53,8 +53,7 @@ static void object_warp_calc_view_matrix(float r_mat_view[4][4], float r_center_
        float viewmat_roll[4][4];
 
        /* apply the rotation offset by rolling the view */
-       unit_m4(mat_offset);
-       rotate_m4(mat_offset, 'Z', offset_angle);
+       axis_angle_to_mat4_single(mat_offset, 'Z', offset_angle);
        mul_m4_m4m4(viewmat_roll, mat_offset, viewmat);
 
        /* apply the view and the object matrix */
index 53434b18d063cb6b531ee2edaad77a487a31749a..ef2f2d36ab7a38c5a09eb82e8f0caa37dbff4a31 100644 (file)
@@ -821,10 +821,9 @@ static float calc_overlap(StrokeCache *cache, const char symm, const char axis,
        flip_v3_v3(mirror, cache->true_location, symm);
 
        if (axis != 0) {
-               float mat[4][4];
-               unit_m4(mat);
-               rotate_m4(mat, axis, angle);
-               mul_m4_v3(mat, mirror);
+               float mat[3][3];
+               axis_angle_to_mat3_single(mat, axis, angle);
+               mul_m3_v3(mat, mirror);
        }
 
        /* distsq = len_squared_v3v3(mirror, cache->traced_location); */
index 59442e8978747f1b8d26cb6dbf7598eb652fc5cf..50aec737c8ef0ebe05be7434b4009cd10be8e66c 100644 (file)
@@ -459,7 +459,7 @@ static void stitch_calculate_island_snapping(
                                            island_stitch_data[i].num_rot_elements_neg) / totelem;
                        }
 
-                       rotate_m2(rotation_mat, rotation);
+                       angle_to_mat2(rotation_mat, rotation);
                        numOfIslandUVs = getNumOfIslandUvs(state->element_map, i);
                        element = &state->element_map->buf[state->element_map->islandIndices[i]];
                        for (j = 0; j < numOfIslandUVs; j++, element++) {
index df94975e2748333fe78f055e7148c9a4cb6d3f09..290e19736bb8d4c34d19f330320d22e958c9bd1c 100644 (file)
@@ -798,13 +798,11 @@ static DerivedMesh *applyModifier(ModifierData *md, Object *ob,
 
                if (ltmd->ob_axis) {
                        axis_angle_normalized_to_mat3(mat3, axis_vec, step_angle);
-                       copy_m4_m3(mat, mat3);
                }
                else {
-                       unit_m4(mat);
-                       rotate_m4(mat, axis_char, step_angle);
-                       copy_m3_m4(mat3, mat);
+                       axis_angle_to_mat3_single(mat3, axis_char, step_angle);
                }
+               copy_m4_m3(mat, mat3);
 
                if (screw_ofs)
                        madd_v3_v3fl(mat[3], axis_vec, screw_ofs * ((float)step / (float)(step_tot - 1)));