Skip to main content

SpectralRadiusComputer

Struct SpectralRadiusComputer 

Source
pub struct SpectralRadiusComputer {
    pub max_iters: u32,
    pub tol: f64,
}
Expand description

Computes the spectral radius of a square matrix via the power method.

The power method applied to ||A^k||^{1/k} converges monotonically down to r(A) by the Gelfand formula.

Fields§

§max_iters: u32

Maximum number of iterations for the power iteration.

§tol: f64

Convergence tolerance: stop when successive estimates differ by less than tol.

Implementations§

Source§

impl SpectralRadiusComputer

Source

pub fn new(max_iters: u32, tol: f64) -> Self

Create a new computer with given iteration count and tolerance.

Source

pub fn default() -> Self

Default: 30 iterations, tolerance 1e-8.

Source

pub fn compute(&self, mat: &SquareMatrix) -> f64

Compute the spectral radius of mat using the Gelfand formula r(A) = inf_k ||A^k||^{1/k}.

Returns the best estimate found within max_iters steps.

Source

pub fn power_vector_method(&self, mat: &SquareMatrix, init: &[f64]) -> f64

Use the power-vector method: iterate v_{k+1} = A v_k / ||A v_k|| and track ||A v_k|| to estimate the dominant eigenvalue magnitude.

Trait Implementations§

Source§

impl Clone for SpectralRadiusComputer

Source§

fn clone(&self) -> SpectralRadiusComputer

Returns a duplicate of the value. Read more
1.0.0 · Source§

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

Performs copy-assignment from source. Read more
Source§

impl Debug for SpectralRadiusComputer

Source§

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

Formats the value using the given formatter. Read more

Auto Trait Implementations§

Blanket Implementations§

Source§

impl<T> Any for T
where T: 'static + ?Sized,

Source§

fn type_id(&self) -> TypeId

Gets the TypeId of self. Read more
Source§

impl<T> Borrow<T> for T
where T: ?Sized,

Source§

fn borrow(&self) -> &T

Immutably borrows from an owned value. Read more
Source§

impl<T> BorrowMut<T> for T
where T: ?Sized,

Source§

fn borrow_mut(&mut self) -> &mut T

Mutably borrows from an owned value. Read more
Source§

impl<T> CloneToUninit for T
where T: Clone,

Source§

unsafe fn clone_to_uninit(&self, dest: *mut u8)

🔬This is a nightly-only experimental API. (clone_to_uninit)
Performs copy-assignment from self to dest. Read more
Source§

impl<T> From<T> for T

Source§

fn from(t: T) -> T

Returns the argument unchanged.

Source§

impl<T, U> Into<U> for T
where U: From<T>,

Source§

fn into(self) -> U

Calls U::from(self).

That is, this conversion is whatever the implementation of From<T> for U chooses to do.

Source§

impl<T> ToOwned for T
where T: Clone,

Source§

type Owned = T

The resulting type after obtaining ownership.
Source§

fn to_owned(&self) -> T

Creates owned data from borrowed data, usually by cloning. Read more
Source§

fn clone_into(&self, target: &mut T)

Uses borrowed data to replace owned data, usually by cloning. Read more
Source§

impl<T, U> TryFrom<U> for T
where U: Into<T>,

Source§

type Error = Infallible

The type returned in the event of a conversion error.
Source§

fn try_from(value: U) -> Result<T, <T as TryFrom<U>>::Error>

Performs the conversion.
Source§

impl<T, U> TryInto<U> for T
where U: TryFrom<T>,

Source§

type Error = <U as TryFrom<T>>::Error

The type returned in the event of a conversion error.
Source§

fn try_into(self) -> Result<U, <U as TryFrom<T>>::Error>

Performs the conversion.