Struct MatrixOperator 

pub struct MatrixOperator<const P: usize, const Q: usize, const R: usize> { /* private fields */ }

A linear operator represented as a matrix.

The matrix is stored in row-major order.

Implementations

impl<const P: usize, const Q: usize, const R: usize> MatrixOperator<P, Q, R>

pub const DIM: usize

The dimension of the Clifford algebra.

pub fn new(entries: Vec<f64>, rows: usize, cols: usize) -> Result<Self>

Create a new matrix operator from entries.

Entries should be in row-major order.

pub fn identity() -> Self

Create an identity matrix.

pub fn zeros() -> Self

Create a zero matrix.

pub fn diagonal(diag: &[f64]) -> Result<Self>

Create a diagonal matrix from diagonal entries.

pub fn get(&self, row: usize, col: usize) -> f64

Get a matrix entry.

pub fn set(&mut self, row: usize, col: usize, value: f64)

Set a matrix entry.

pub fn rows(&self) -> usize

Get the number of rows.

pub fn cols(&self) -> usize

Get the number of columns.

pub fn transpose(&self) -> Self

Compute the transpose of the matrix.

pub fn trace(&self) -> f64

Compute the trace of the matrix.

pub fn multiply(&self, other: &Self) -> Result<Self>

Multiply this matrix by another.

pub fn add(&self, other: &Self) -> Result<Self>

Add another matrix.

pub fn scale(&self, scalar: f64) -> Self

Scale the matrix by a scalar.

pub fn is_symmetric(&self, tolerance: f64) -> bool

Check if the matrix is symmetric.

Trait Implementations

impl<const P: usize, const Q: usize, const R: usize> AdjointableOperator<Multivector<P, Q, R>> for MatrixOperator<P, Q, R>

type Adjoint = MatrixOperator<P, Q, R>

The type of the adjoint operator.

fn adjoint(&self) -> Self::Adjoint

Compute the adjoint operator.

fn is_self_adjoint(&self) -> bool

Check if this operator is self-adjoint (T = T*).

fn is_normal(&self) -> bool

Check if this operator is normal (TT* = T*T).

impl<const P: usize, const Q: usize, const R: usize> BoundedOperator<Multivector<P, Q, R>> for MatrixOperator<P, Q, R>

fn operator_norm(&self) -> f64

Compute or estimate the operator norm ||T||.

fn is_bounded_by(&self, bound: f64) -> bool

Check if the operator is bounded by a given constant.

impl<const P: usize, const Q: usize, const R: usize> Clone for MatrixOperator<P, Q, R>

fn clone(&self) -> MatrixOperator<P, Q, R>

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl<const P: usize, const Q: usize, const R: usize> Debug for MatrixOperator<P, Q, R>

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl<const P: usize, const Q: usize, const R: usize> LinearOperator<Multivector<P, Q, R>> for MatrixOperator<P, Q, R>

fn apply(&self, x: &Multivector<P, Q, R>) -> Result<Multivector<P, Q, R>>

Apply the operator to an element.

fn domain_dimension(&self) -> Option<usize>

Get the domain dimension (None if infinite).

fn codomain_dimension(&self) -> Option<usize>

Get the codomain dimension (None if infinite).

impl<const P: usize, const Q: usize, const R: usize> OperatorNorm for MatrixOperator<P, Q, R>

fn norm(&self) -> f64

Compute the operator norm.

fn frobenius_norm(&self) -> Option<f64>

Compute the Frobenius norm (for matrix operators).