eigen_mod::solve Interface Reference

Solver for: More...

Collaboration diagram for eigen_mod::solve:

Private Member Functions
subroutine	eigen_solve (eig, A, B, Op)
	Eigen-problem resolution More...

Detailed Description

Solver for:

• standard eigenvalue problems: \(~~ A\,x = \lambda \,x \)
• geeralised eigenvalue problems: \(~~ A\,x = \lambda \, B \,x \)

call solve(eig, args) 

with args = arguments depending on the problem type and on the eigen-solver 'eig' mode:

• Standard problems.
  • Regular mode

    call solve(eig, A=A)  

  • Shift invert mode

    call solve(eig, Op=Op)  

    Op= \( (A-sId)^{-1} \),
    • \( s \) is the shift parameter
    • by default \( s=0 \) and Op= \( A^{-1} \)
• Geeralised problems.
  • inverse mode

    call solve(eig, A=A, B=B, Op=Op)  

    with Op= \( B^{-1} \)
  • Shift invert mode

    call solve(eig, B=B, Op=Op)  

    Op= \( (A-sB)^{-1} \),
    • \( s \) is the shift parameter
    • by default \( s=0 \) and Op= \( A^{-1} \) ...

Output:

• The computed eigen-values are in the 1D-array

  eig%lambda 

• The associated eigen-vectors are in the 2D-array

  eig%v 

Definition at line 239 of file eigen_mod.F90.

Member Function/Subroutine Documentation

eigen_solve()

subroutine eigen_mod::solve::eigen_solve ( type(eigen), intent(inout)  eig, procedure(rntorn), optional  A, procedure(rntorn), optional  B, procedure(rntornxl), optional  Op  )

private

Eigen-problem resolution

Parameters

[in,out]

eig

eigen-solver settings

See solve for the output description.

Optional arguments:

Parameters

[in]	A	eigen-problem matrix, procedural arguments with interface RnToRn
[in]	B	eigen-problem matrix (generalised problems), procedural arguments with interface RnToRn
[in]	Op	linear problem solver, procedural arguments with interface RnToRnxL.

See solve for the argument list description.

Definition at line 458 of file eigen_mod.F90.

---

The documentation for this interface was generated from the following file:

• /home/cpierre1/choral/SRC/eigen_mod.F90