Struct Matrix 

pub struct Matrix { /* private fields */ }

Sparse Matrix

Implementations

impl Matrix

pub fn new() -> Matrix

Create a new, initially empty Matrix

pub fn from_entries(entries: Vec<(usize, usize, Complex64)>) -> Matrix

Create a new Matrix from a vector of (row, col, val) entries.

pub fn identity(n: usize) -> Matrix

Create an n*n identity Matrix

pub fn add_element(&mut self, row: usize, col: usize, val: Complex64)

Add an element at location (row, col) with value val.

pub fn add_elements(&mut self, elements: Vec<(usize, usize, Complex64)>)

Add elements correspoding to each triplet (row, col, val) Rows and columns are usize, and vals are Complex64.

pub fn make(&mut self, row: usize, col: usize) -> Eindex

Create a zero-valued element at (row, col), or return existing Element index if present

pub fn reset(&mut self)

Reset all Elements to zero value.

pub fn update(&mut self, ei: Eindex, val: Complex64)

Update Element ei by val

pub fn vecmul(&self, x: &Vec<Complex64>) -> Result<Vec<Complex64>, &'static str>

Multiply by Vec

pub fn res( &self, x: &Vec<Complex64>, rhs: &Vec<Complex64>, ) -> Result<Vec<Complex64>, &'static str>

pub fn get_elem(&self, row: usize, col: usize) -> Option<Eindex>

Returns the Element-index at (row, col) if present, or None if not.

pub fn get(&self, row: usize, col: usize) -> Option<Complex64>

Returns the Element-value at (row, col) if present, or None if not.

pub fn solve( &mut self, rhs: &[Complex64], ) -> Result<Vec<Complex64>, &'static str>

Solve the system Ax=b, where:

• A is self
• b is argument rhs
• x is the return value.

Returns a Result containing the Vec<Complex64> representing x if successful. Returns an Err if unsuccessful.

Performs LU factorization, forward and backward substitution.

Trait Implementations

impl Debug for Matrix

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

Formats the value using the given formatter.

impl Index<Eindex> for Matrix

type Output = Element

The returned type after indexing.

fn index(&self, index: Eindex) -> &Self::Output

Performs the indexing (container[index]) operation.

impl IndexMut<Eindex> for Matrix

fn index_mut(&mut self, index: Eindex) -> &mut Self::Output

Performs the mutable indexing (container[index]) operation.