impl Outermorphism for Matrix4x4

ega/src/lib.rs

@@ -5,11 +5,43 @@ use ::libm::Libm;
 
 mod operators;
 mod optional_features;
+mod outermorphisms;
 mod values;
 
 pub use operators::*;
+pub use outermorphisms::*;
 pub use values::*;
 
+/// return `Empty`
+#[inline]
+fn return_empty<Lhs, Rhs>(_: Lhs, _: Rhs) -> Empty {
+ Empty
+}
+
+/// return the left-hand-side
+#[inline]
+fn return_lhs<Lhs, Rhs>(lhs: Lhs, _: Rhs) -> Lhs {
+ lhs
+}
+
+/// return the right-hand-side
+#[inline]
+fn return_rhs<Lhs, Rhs>(_: Lhs, rhs: Rhs) -> Rhs {
+ rhs
+}
+
+/// negate and return the right-hand-side
+#[inline]
+fn return_neg_rhs<Lhs, Rhs: Neg<Output = Rhs>>(_: Lhs, rhs: Rhs) -> Rhs {
+ -rhs
+}
+
+/// return zero
+#[inline]
+fn return_zero<Lhs, Rhs, Output: Zero>(_: Lhs, _: Rhs) -> Output {
+ Output::zero()
+}
+
 #[rustfmt::skip]
 #[cfg(any(test, doctest))]
 mod test_values {
@@ -83,4 +115,11 @@ mod test_values {
  pub const PSEUDOSCALAR_A: Pseudoscalar = Pseudoscalar { e0123: 397. };
  pub const PSEUDOSCALAR_B: Pseudoscalar = Pseudoscalar { e0123: 401. };
  pub const PSEUDOSCALAR_C: Pseudoscalar = Pseudoscalar { e0123: -409. };
+
+ pub const MATRIX_4X4_A: Matrix4x4 = Matrix4x4 {
+ m00: 1.0, m01: 0.5, m02: 0.5, m03: 0.5,
+ m10: 0.5, m11: 0.9, m12: -0.5, m13: -1.5,
+ m20: 0.5, m21: -1.0, m22: 0.5, m23: -2.0,
+ m30: 0.5, m31: -0.9, m32: 2.0, m33: -1.0,
+ };
 }

ega/src/operators/add.rs

@@ -1,4 +1,3 @@
-use super::{return_lhs, return_rhs};
 use crate::*;
 
 pub use ::core::ops::Add;

ega/src/operators/dot.rs

@@ -1,4 +1,3 @@
-use super::return_empty;
 use crate::*;
 
 /// The inner product

ega/src/operators/geometric_product.rs

@@ -1,4 +1,3 @@
-use super::return_empty;
 use crate::*;
 
 /// The geometric product

ega/src/operators/join.rs

@@ -1,4 +1,3 @@
-use super::return_empty;
 use crate::*;
 
 /// The regressive product

ega/src/operators/meet.rs

@@ -1,4 +1,3 @@
-use super::return_empty;
 use crate::*;
 
 /// The outer product

ega/src/operators/mod.rs

@@ -38,8 +38,6 @@ pub use reverse::Reverse;
 pub use scalar_product::ScalarProduct;
 pub use sub::Sub;
 
-use crate::{values::Empty, Zero};
-
 /// Grade Involution
 pub trait Involution {
  type Output;
@@ -55,33 +53,3 @@ pub trait Exponent<Rhs> {
  /// Exponentiation
  fn exp(self, rhs: Rhs) -> Self::Output;
 }
-
-/// return `Empty`
-#[inline]
-fn return_empty<Lhs, Rhs>(_: Lhs, _: Rhs) -> Empty {
- Empty
-}
-
-/// return the left-hand-side
-#[inline]
-fn return_lhs<Lhs, Rhs>(lhs: Lhs, _: Rhs) -> Lhs {
- lhs
-}
-
-/// return the right-hand-side
-#[inline]
-fn return_rhs<Lhs, Rhs>(_: Lhs, rhs: Rhs) -> Rhs {
- rhs
-}
-
-/// negate and return the right-hand-side
-#[inline]
-fn return_neg_rhs<Lhs, Rhs: Neg<Output = Rhs>>(_: Lhs, rhs: Rhs) -> Rhs {
- -rhs
-}
-
-/// return zero
-#[inline]
-fn return_zero<Lhs, Rhs, Output: Zero>(_: Lhs, _: Rhs) -> Output {
- Output::zero()
-}

ega/src/operators/scalar_product.rs

@@ -1,4 +1,3 @@
-use super::return_zero;
 use crate::*;
 
 /// Scalar Product

ega/src/operators/sub.rs

@@ -1,4 +1,3 @@
-use super::{return_lhs, return_neg_rhs};
 use crate::*;
 
 pub use core::ops::Sub;