Function power_iteration

pub fn power_iteration(
    a: ArrayView2<'_, f64>,
    x: &mut Array1<f64>,
    options: &Options,
) -> (f64, usize)

Computes the dominant eigenvector of the given matrix a using the power iteration algorithm.

The power iteration algorithm is an iterative method for finding the dominant eigenvector of a matrix. This function takes the matrix a, an initial guess for the eigenvector x, and a set of options controlling the algorithm’s behavior.

The function returns the dominant eigenvalue and the number of iterations performed.