smod_qr_invert Submodule

Submodule containing the implementation of the QR-invert algorithm. The original problem is written as a standard eigenvalue problem through $$\mathcal{B}^{-1}\mathcal{A}\textbf{X} = \omega\textbf{X}\ $$ . This is done using a LU decomposition via LAPACKS's zgbsv. Eventually a call to LAPACK's zgeev routine is done to obtain all eigenvalues and eigenvectors.

Contents

Module Procedures

qr_invert

Module Procedures

module procedure qr_invert module module subroutine qr_invert(matrix_A, matrix_B, settings, omega, vr) Interface →

Solves the eigenvalue problem by rewriting it to a standard form through inversion of the B-matrix.

Read more…

Arguments

Type	Intent	Optional	Attributes		Name
type(matrix_t),	intent(in)			::	matrix_A	matrix A
type(matrix_t),	intent(in)			::	matrix_B	matrix B
type(settings_t),	intent(in)			::	settings	settings object
complex(kind=dp),	intent(out)			::	omega(:)	array with eigenvalues
complex(kind=dp),	intent(out)			::	vr(:,:)	array with right eigenvectors