BLI: Add support for non-square matrix multiplication.

remove the square matrix operator and correct a few naming
mistake.

@@ -384,1 +382,4 @@
 
+  /** Multiply two compatible matrices using matrix multiplication. */
+  template<int OtherNumRow>
+  MatBase<T, OtherNumRow, NumRow> operator*(const MatBase<T, OtherNumRow, NumCol> &b) const

@@ -458,6 +473,13 @@ struct alignas(Alignment) MatBase : public vec_struct_base<VecBase<T, NumRow>, N
     return result;
   }
 
+  static MatBase diagonal(const VecBase<T, min_dim> value)

@@ -371,6 +371,30 @@ TEST(math_matrix, MatrixMultiply)
   }
 }
 
+TEST(math_matrix, RectangularMatrixMultiply)

@@ -374,0 +375,4 @@
+{
+
+  {
+    const float3x4 matrix_a = transpose(float4x3(float3(1.0, 0.0, 1.0),

@@ -374,0 +376,4 @@
+
+  {
+    const float3x4 matrix_a = transpose(float4x3(float3(1.0, 0.0, 1.0),
+                                                 float3(2.0, 1.0, 1.0),

@@ -374,0 +387,4 @@
+                                                 float3(6.0f, 5.0f, 3.0f),
+                                                 float3(11.0f, 9.0f, 6.0f)));
+
+    const float3x4 result = matrix_a * matrix_b;