Merge pull request #1 from Akari202/matrix

Merge matrix overhaul into main.

Author: Akari202

.github/workflows/rust.yml

@@ -17,6 +17,6 @@ jobs:
     steps:
     - uses: actions/checkout@v3
     - name: Build
-      run: cargo build --verbose
+      run: cargo build --verbose --all-features
     - name: Run tests
-      run: cargo test --verbose
+      run: cargo test --verbose --all-features

Cargo.lock

@@ -22,3 +22,5 @@ https://www.cs.utexas.edu/~ear/cs341/automatabook/AutomataTheoryBook.pdf
 https://dl.acm.org/doi/pdf/10.1145/360825.360855
 https://docs.google.com/document/d/1KkKC2-ozJkvbWQIAXeJ1MUGqxjn19c-Mmc7RtxFTA3c/edit?tab=t.0
 https://ntrs.nasa.gov/api/citations/19900013774/downloads/19900013774.pdf
+
+https://www.cs.unc.edu/techreports/96-043.pdf

README.md

@@ -1,4 +1,4 @@
-VERSION = "0.2.5"
+VERSION = "0.3.0"
 
 locations = [
     "./vec-utils/Cargo.toml",

update_version.py

@@ -1,6 +1,6 @@
 [package]
 name = "vec-utils-py"
-version = "0.2.5"
+version = "0.3.0"
 edition = "2024"
 authors = ["Akari Harada <akaritwo0two@gmail.com>"]
 license = "GPL-3"
@@ -17,4 +17,4 @@ pyo3 = { version = "0.27.2", features = [
   "generate-import-lib"
 ] }
 rayon = "1.11.0"
-vec-utils = { path = "../vec-utils" }
+vec-utils = { path = "../vec-utils", features = ["matrix"] }

vec-utils-py/Cargo.toml

@@ -1,6 +1,6 @@
 [project]
 name = "vec-utils-py"
-version = "0.2.5"
+version = "0.3.0"
 requires-python = ">=3.9"
 authors = [
   { name="Akari202", email="akaritwo0two@gmail.com" }

vec-utils-py/pyproject.toml

@@ -104,7 +104,7 @@ impl Quat {
     }
 
     fn to_rotation_matrix(&self) -> [[f64; 3]; 3] {
-        self.inner.to_rotation_matrix()
+        self.inner.to_rotation_matrix().to_nested_arr()
     }
 
     fn rotate(&self, v: &Vec3d) -> Vec3d {

vec-utils-py/src/quat.rs

@@ -1,6 +1,6 @@
 [package]
 name = "vec-utils"
-version = "0.2.5"
+version = "0.3.0"
 edition = "2024"
 authors = ["Akari Harada <akaritwo0two@gmail.com>"]
 license = "GPL-3"
@@ -10,6 +10,9 @@ name = "vec_utils"
 
 [dependencies]
 libm = { version = "0.2.15" }
+matrixmultiply = { version = "0.3.10", default-features = false, features = [
+  "cgemm"
+], optional = true }
 rkyv = { version = "0.8.12", default-features = false, features = ["bytecheck"], optional = true }
 serde = { version = "1.0.219", default-features = false, features = ["derive"], optional = true }
 
@@ -19,6 +22,7 @@ pretty_assertions = "1.4.1"
 
 [features]
 default = ["std"]
+matrix = ["dep:matrixmultiply"]
 rkyv = ["dep:rkyv"]
 serde = ["dep:serde"]
-std = ["rkyv?/std", "serde?/std"]
+std = ["matrixmultiply?/std", "rkyv?/std", "serde?/std"]

vec-utils/Cargo.toml

@@ -1,3 +1,4 @@
+use core::cmp::Ordering;
 use core::fmt;
 use core::ops::{Add, Div, Index, Mul, Sub};
 
@@ -7,6 +8,7 @@ use crate::{
 };
 
 /// A complex number
+#[repr(C)]
 #[cfg_attr(feature = "serde", derive(serde::Deserialize, serde::Serialize))]
 #[cfg_attr(
     feature = "rkyv",
@@ -241,6 +243,42 @@ impl_single_op_variants_other!(
     "Divide a real number by a complex number"
 );
 
+impl PartialOrd for Complex {
+    fn partial_cmp(&self, other: &Self) -> Option<Ordering> {
+        let self_mag_sq = self.real * self.real + self.imaginary * self.imaginary;
+        let other_mag_sq = other.real * other.real + other.imaginary * other.imaginary;
+
+        self_mag_sq.partial_cmp(&other_mag_sq)
+    }
+}
+
+impl PartialEq<f64> for Complex {
+    fn eq(&self, other: &f64) -> bool {
+        self.imaginary == 0.0 && self.real == *other
+    }
+}
+
+impl PartialOrd<f64> for Complex {
+    fn partial_cmp(&self, other: &f64) -> Option<Ordering> {
+        let self_mag_sq = self.real * self.real + self.imaginary * self.imaginary;
+        let other_mag_sq = other * other;
+
+        self_mag_sq.partial_cmp(&other_mag_sq)
+    }
+}
+
+impl PartialEq<Complex> for f64 {
+    fn eq(&self, other: &Complex) -> bool {
+        other.eq(self)
+    }
+}
+
+impl PartialOrd<Complex> for f64 {
+    fn partial_cmp(&self, other: &Complex) -> Option<Ordering> {
+        other.partial_cmp(self).map(Ordering::reverse)
+    }
+}
+
 impl Index<usize> for Complex {
     type Output = f64;
 

vec-utils/src/complex.rs

@@ -1,6 +1,3 @@
-//! Geometry module
-//! This module contains geometric shapes and operations on them.
-
 /// Circles
 pub mod circle;
 /// Intersections