scale_gmres Module Reference

module common / GMRES More...

Data Types
type	gmres

Functions/Subroutines
subroutine	gmres_init (this, n, m, eps, eps0)

Detailed Description

module common / GMRES

Description

A module to provide a iterative solver for system of linear equations using generalized minimal residual method (GMRES)

Reference

• Y. Saad and M.H. Schultz, 1986: GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Stat. Comput., 7:856–869

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

gmres_init()

subroutine scale_gmres::gmres_init	(	class(gmres), intent(inout)	this,
		integer, intent(in)	n,
		integer, intent(in)	m,
		real(rp), intent(in)	eps,
		real(rp), intent(in)	eps0 )

Definition at line 78 of file scale_gmres.F90.

 
  implicit none
  class(GMRES), intent(inout) :: this
  integer, intent(in) :: N
  integer, intent(in) :: m
  real(RP), intent(in) :: eps
  real(RP), intent(in) :: eps0
 
  !--------------------------------------
 
  this%N = n
  this%m = m
  this%EPS = eps
  this%EPS0 = eps0
 
  allocate( this%v(n,m+1) )
  allocate( this%hj(m+1) )
  allocate( this%g(m+1) )
  allocate( this%r(m+1,m) )
  allocate( this%co(m), this%si(m) )
  allocate( this%y(m) )
 
  return