Struct Matrix 

pub struct Matrix { /* private fields */ }

Represents a square matrix over C of maximum size MAX_SIZE.

Each element is an instance of Complex, and we store the elements internally in an array of size MAX_SIZE^2 * sizeof(Complex).

In practice, this means each matrix occupies around 16KiB.

Implementations

impl Matrix

pub fn new(size: usize) -> Matrix

Construct a new zero-initialized matrix of given size.

Panics

We panic if the given size exceeds MAX_SIZE.

pub fn new_from_elements(size: usize, elements: Vec<Complex>) -> Matrix

Construct a new matrix of given size from elements.

Panics

We panic if the given size exceeds MAX_SIZE.

pub fn identity(size: usize) -> Matrix

Construct a new identity matrix of given size.

Panics

We panic if the given size exceeds MAX_SIZE.

pub fn embed(&mut self, other: &Matrix, i: usize, j: usize)

Embed another matrix into this one, overrising elements.

Embed with top-left position at (i, j).

Panics

We panic if this matrix isn’t large enough.

pub fn permute_rows(&self, permutation: Vec<usize>) -> Matrix

Permute the rows to generate a new matrix.

Row i goes to row perutation[i].

Panics

We panic if set(permutation) != {0, …, self.size - 1}.

pub fn permute_columns(&self, permutation: Vec<usize>) -> Matrix

Permute the columns to generate a new matrix.

Column i goes to column perutation[i].

Panics

We panic if set(permutation) != {0, …, self.size - 1}.

pub fn size(&self) -> usize

Size of the matrix.

pub fn get(&self, i: usize, j: usize) -> Complex

Get the element in position (i, j).

pub fn set(&mut self, i: usize, j: usize, value: Complex)

Set the element in position (i, j) to value.

pub fn approx_eq(&self, other: &Matrix) -> bool

Approximately equal test.

Trait Implementations

impl<'a> Add for &'a Matrix

type Output = Matrix

The resulting type after applying the + operator.

fn add(self, rhs: &'a Matrix) -> Matrix

Performs the + operation.

impl Debug for Matrix

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

Formats the value using the given formatter.

impl<'a> Mul for &'a Matrix

type Output = Matrix

The resulting type after applying the * operator.

fn mul(self, rhs: &Matrix) -> Matrix

Performs the * operation.

impl<'a> Mul<&'a [Complex; 32]> for &'a Matrix

Implements standard matrix vector multiplication.

Panics

We panic if the vector contains non-zero elements in positions self.size or beyond.

type Output = [Complex; 32]

The resulting type after applying the * operator.

fn mul(self, rhs: &Vector) -> Vector

Performs the * operation.

impl PartialEq for Matrix

fn eq(&self, other: &Matrix) -> bool

Equality operator ==.

fn ne(&self, other: &Rhs) -> bool

Inequality operator !=.