Struct totsu_f64lapack::F64LAPACK

pub struct F64LAPACK;

f64-specific LinAlgEx implementation using cblas-sys and lapacke-sys

You need a BLAS/LAPACK source to link.

Trait Implementations

impl Clone for F64LAPACK

fn clone(&self) -> F64LAPACK

Returns a copy of the value.

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

Performs copy-assignment from source.

impl LinAlg for F64LAPACK

type F = f64

Floating point data type used as scalars.

type Sl = [f64]

Data type of slice of F used as vectors.

fn norm(x: &[f64]) -> f64

Calculates 2-norm (or euclidean norm) \(\|x\|_2=\sqrt{\sum_i x_i^2}\).

fn copy(x: &[f64], y: &mut [f64])

Copies from a vector to another vector.

fn scale(alpha: f64, x: &mut [f64])

Calculates \(\alpha x\).

fn add(alpha: f64, x: &[f64], y: &mut [f64])

Calculates \(\alpha x + y\).

fn adds(s: f64, y: &mut [f64])

Calculates \(s\mathbb{1} + y\).

fn abssum(x: &[f64], incx: usize) -> f64

Calculates 1-norm (or sum of absolute values) \(\|x\|_1=\sum_i |x_i|\).

fn transform_di(alpha: f64, mat: &[f64], x: &[f64], beta: f64, y: &mut [f64])

Calculates \(\alpha D x + \beta y\), where \(D={\bf diag}(d)\) is a diagonal matrix.

impl LinAlgEx for F64LAPACK

fn transform_ge(
    transpose: bool,
    n_row: usize,
    n_col: usize,
    alpha: f64,
    mat: &[f64],
    x: &[f64],
    beta: f64,
    y: &mut [f64]
)

Calculates \(\alpha G x + \beta y\).

fn transform_sp(
    n: usize,
    alpha: f64,
    mat: &[f64],
    x: &[f64],
    beta: f64,
    y: &mut [f64]
)

Calculates \(\alpha S x + \beta y\), where \(S\) is a symmetric matrix, supplied in packed form.

fn map_eig_worklen(n: usize) -> usize

Query of a length of work slice that LinAlgEx::map_eig requires.

fn map_eig<M>(
    mat: &mut [f64],
    scale_diag: Option<f64>,
    eps_zero: f64,
    work: &mut [f64],
    map: M
)where
    M: Fn(f64) -> Option<f64>,

Applies a map to eigenvalues of a symmetric matrix \(S\) supplied in packed form.