Fixed impl for scalar / matrix and added matrix recip().

Author: bitshifter

src/f32/coresimd/mat2.rs

@@ -403,6 +403,13 @@ impl Mat2 {
         Self(self.0 / f32x4::splat(rhs))
     }
 
+    /// Returns a matrix containing the reciprocal `1.0/n` of each element of `self`.
+    #[inline]
+    #[must_use]
+    pub fn recip(&self) -> Self {
+        Self::from_cols(self.x_axis.recip(), self.y_axis.recip())
+    }
+
     /// Returns true if the absolute difference of all elements between `self` and `rhs`
     /// is less than or equal to `max_abs_diff`.
     ///
@@ -716,7 +723,7 @@ impl Div<Mat2> for f32 {
     type Output = Mat2;
     #[inline]
     fn div(self, rhs: Mat2) -> Self::Output {
-        rhs.div_scalar(self)
+        Mat2(f32x4::splat(self) / rhs.0)
     }
 }
 

src/f32/coresimd/mat3a.rs

@@ -826,6 +826,17 @@ impl Mat3A {
         )
     }
 
+    /// Returns a matrix containing the reciprocal `1.0/n` of each element of `self`.
+    #[inline]
+    #[must_use]
+    pub fn recip(&self) -> Self {
+        Self::from_cols(
+            self.x_axis.recip(),
+            self.y_axis.recip(),
+            self.z_axis.recip(),
+        )
+    }
+
     /// Returns true if the absolute difference of all elements between `self` and `rhs`
     /// is less than or equal to `max_abs_diff`.
     ///
@@ -1147,7 +1158,11 @@ impl Div<Mat3A> for f32 {
     type Output = Mat3A;
     #[inline]
     fn div(self, rhs: Mat3A) -> Self::Output {
-        rhs.div_scalar(self)
+        Mat3A::from_cols(
+            self.div(rhs.x_axis),
+            self.div(rhs.y_axis),
+            self.div(rhs.z_axis),
+        )
     }
 }
 

src/f32/coresimd/mat4.rs

@@ -1469,6 +1469,18 @@ impl Mat4 {
         )
     }
 
+    /// Returns a matrix containing the reciprocal `1.0/n` of each element of `self`.
+    #[inline]
+    #[must_use]
+    pub fn recip(&self) -> Self {
+        Self::from_cols(
+            self.x_axis.recip(),
+            self.y_axis.recip(),
+            self.z_axis.recip(),
+            self.w_axis.recip(),
+        )
+    }
+
     /// Returns true if the absolute difference of all elements between `self` and `rhs`
     /// is less than or equal to `max_abs_diff`.
     ///
@@ -1805,7 +1817,12 @@ impl Div<Mat4> for f32 {
     type Output = Mat4;
     #[inline]
     fn div(self, rhs: Mat4) -> Self::Output {
-        rhs.div_scalar(self)
+        Mat4::from_cols(
+            self.div(rhs.x_axis),
+            self.div(rhs.y_axis),
+            self.div(rhs.z_axis),
+            self.div(rhs.w_axis),
+        )
     }
 }
 

src/f32/mat3.rs

@@ -756,6 +756,17 @@ impl Mat3 {
         )
     }
 
+    /// Returns a matrix containing the reciprocal `1.0/n` of each element of `self`.
+    #[inline]
+    #[must_use]
+    pub fn recip(&self) -> Self {
+        Self::from_cols(
+            self.x_axis.recip(),
+            self.y_axis.recip(),
+            self.z_axis.recip(),
+        )
+    }
+
     /// Returns true if the absolute difference of all elements between `self` and `rhs`
     /// is less than or equal to `max_abs_diff`.
     ///
@@ -1077,7 +1088,11 @@ impl Div<Mat3> for f32 {
     type Output = Mat3;
     #[inline]
     fn div(self, rhs: Mat3) -> Self::Output {
-        rhs.div_scalar(self)
+        Mat3::from_cols(
+            self.div(rhs.x_axis),
+            self.div(rhs.y_axis),
+            self.div(rhs.z_axis),
+        )
     }
 }
 

src/f32/neon/mat2.rs

@@ -445,6 +445,13 @@ impl Mat2 {
         Self(unsafe { vdivq_f32(self.0, vld1q_dup_f32(&rhs)) })
     }
 
+    /// Returns a matrix containing the reciprocal `1.0/n` of each element of `self`.
+    #[inline]
+    #[must_use]
+    pub fn recip(&self) -> Self {
+        Self::from_cols(self.x_axis.recip(), self.y_axis.recip())
+    }
+
     /// Returns true if the absolute difference of all elements between `self` and `rhs`
     /// is less than or equal to `max_abs_diff`.
     ///
@@ -764,7 +771,7 @@ impl Div<Mat2> for f32 {
     type Output = Mat2;
     #[inline]
     fn div(self, rhs: Mat2) -> Self::Output {
-        rhs.div_scalar(self)
+        Mat2(unsafe { vdivq_f32(vld1q_dup_f32(&self), rhs.0) })
     }
 }
 

src/f32/neon/mat3a.rs

@@ -844,6 +844,17 @@ impl Mat3A {
         )
     }
 
+    /// Returns a matrix containing the reciprocal `1.0/n` of each element of `self`.
+    #[inline]
+    #[must_use]
+    pub fn recip(&self) -> Self {
+        Self::from_cols(
+            self.x_axis.recip(),
+            self.y_axis.recip(),
+            self.z_axis.recip(),
+        )
+    }
+
     /// Returns true if the absolute difference of all elements between `self` and `rhs`
     /// is less than or equal to `max_abs_diff`.
     ///
@@ -1165,7 +1176,11 @@ impl Div<Mat3A> for f32 {
     type Output = Mat3A;
     #[inline]
     fn div(self, rhs: Mat3A) -> Self::Output {
-        rhs.div_scalar(self)
+        Mat3A::from_cols(
+            self.div(rhs.x_axis),
+            self.div(rhs.y_axis),
+            self.div(rhs.z_axis),
+        )
     }
 }
 

src/f32/neon/mat4.rs

@@ -1481,6 +1481,18 @@ impl Mat4 {
         )
     }
 
+    /// Returns a matrix containing the reciprocal `1.0/n` of each element of `self`.
+    #[inline]
+    #[must_use]
+    pub fn recip(&self) -> Self {
+        Self::from_cols(
+            self.x_axis.recip(),
+            self.y_axis.recip(),
+            self.z_axis.recip(),
+            self.w_axis.recip(),
+        )
+    }
+
     /// Returns true if the absolute difference of all elements between `self` and `rhs`
     /// is less than or equal to `max_abs_diff`.
     ///
@@ -1817,7 +1829,12 @@ impl Div<Mat4> for f32 {
     type Output = Mat4;
     #[inline]
     fn div(self, rhs: Mat4) -> Self::Output {
-        rhs.div_scalar(self)
+        Mat4::from_cols(
+            self.div(rhs.x_axis),
+            self.div(rhs.y_axis),
+            self.div(rhs.z_axis),
+            self.div(rhs.w_axis),
+        )
     }
 }
 

src/f32/scalar/mat2.rs

@@ -412,6 +412,13 @@ impl Mat2 {
         Self::from_cols(self.x_axis.div(rhs), self.y_axis.div(rhs))
     }
 
+    /// Returns a matrix containing the reciprocal `1.0/n` of each element of `self`.
+    #[inline]
+    #[must_use]
+    pub fn recip(&self) -> Self {
+        Self::from_cols(self.x_axis.recip(), self.y_axis.recip())
+    }
+
     /// Returns true if the absolute difference of all elements between `self` and `rhs`
     /// is less than or equal to `max_abs_diff`.
     ///
@@ -716,7 +723,7 @@ impl Div<Mat2> for f32 {
     type Output = Mat2;
     #[inline]
     fn div(self, rhs: Mat2) -> Self::Output {
-        rhs.div_scalar(self)
+        Mat2::from_cols(self.div(rhs.x_axis), self.div(rhs.y_axis))
     }
 }
 

src/f32/scalar/mat3a.rs

@@ -828,6 +828,17 @@ impl Mat3A {
         )
     }
 
+    /// Returns a matrix containing the reciprocal `1.0/n` of each element of `self`.
+    #[inline]
+    #[must_use]
+    pub fn recip(&self) -> Self {
+        Self::from_cols(
+            self.x_axis.recip(),
+            self.y_axis.recip(),
+            self.z_axis.recip(),
+        )
+    }
+
     /// Returns true if the absolute difference of all elements between `self` and `rhs`
     /// is less than or equal to `max_abs_diff`.
     ///
@@ -1149,7 +1160,11 @@ impl Div<Mat3A> for f32 {
     type Output = Mat3A;
     #[inline]
     fn div(self, rhs: Mat3A) -> Self::Output {
-        rhs.div_scalar(self)
+        Mat3A::from_cols(
+            self.div(rhs.x_axis),
+            self.div(rhs.y_axis),
+            self.div(rhs.z_axis),
+        )
     }
 }
 

src/f32/scalar/mat4.rs

@@ -1374,6 +1374,18 @@ impl Mat4 {
         )
     }
 
+    /// Returns a matrix containing the reciprocal `1.0/n` of each element of `self`.
+    #[inline]
+    #[must_use]
+    pub fn recip(&self) -> Self {
+        Self::from_cols(
+            self.x_axis.recip(),
+            self.y_axis.recip(),
+            self.z_axis.recip(),
+            self.w_axis.recip(),
+        )
+    }
+
     /// Returns true if the absolute difference of all elements between `self` and `rhs`
     /// is less than or equal to `max_abs_diff`.
     ///
@@ -1710,7 +1722,12 @@ impl Div<Mat4> for f32 {
     type Output = Mat4;
     #[inline]
     fn div(self, rhs: Mat4) -> Self::Output {
-        rhs.div_scalar(self)
+        Mat4::from_cols(
+            self.div(rhs.x_axis),
+            self.div(rhs.y_axis),
+            self.div(rhs.z_axis),
+            self.div(rhs.w_axis),
+        )
     }
 }