Class EigenvaluesProblems::QEPOperator

class QEPOperator

Operator for structured quadratic eigenvalue problem (QEP) linearization Implements the matrix-vector product for the linearized system C1^{-1}C0 where the original QEP is: (A + λB + λ²C)v = 0

Public Types

using Scalar = NT

Scalar type required by Spectra library.

Public Functions

inline QEPOperator(const MT &A, const MT &B, const MT &C)

Construct QEP operator from coefficient matrices

Parameters:

  • A – [in] Constant term matrix

  • B – [in] Linear term matrix

  • C – [in] Quadratic term matrix (should be positive definite)

inline int rows() const

Number of rows in linearized system (twice original dimension)

inline int cols() const

Number of columns in linearized system (twice original dimension)

inline void perform_op(const NT *x_in, NT *y_out) const

Perform matrix-vector product y = Op * x for Spectra eigenvalue solver Implements the linearized QEP system operator using efficient factorizations

Parameters:

  • x_in – [in] Input vector of size 2*n

  • y_out – [out] Output vector of size 2*n