fstr_eigen_setup.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
! Copyright (c) 2019 FrontISTR Commons
! This software is released under the MIT License, see LICENSE.txt
!-------------------------------------------------------------------------------
module m_fstr_eigen_setup
  use m_fstr
  use fstr_setup_util, only &
    : fstr_param_pack &
    , fstr_ctrl_err_stop
  use fstr_ctrl_eigen
  implicit none
 
contains
 
  subroutine fstr_eigen_init( fstrEIG )
    implicit none
    type(fstr_eigen) :: fstrEIG
 
    fstreig%nget        = 5
    fstreig%maxiter     = 60
    fstreig%iter        = 0
    !fstrEIG%sigma       = 0.0d0
    fstreig%sigma       = 0.01d0
    fstreig%tolerance   = 1.0d-6
    fstreig%totalmass   = 0.0d0
  end subroutine fstr_eigen_init
 
  subroutine fstr_eigen_finalize( fstrEIG )
    implicit none
    type(fstr_eigen) :: fstrEIG
    if(associated(fstreig%eigval)) deallocate(fstreig%eigval)
    if(associated(fstreig%eigvec)) deallocate(fstreig%eigvec)
    if(associated(fstreig%resid)) deallocate(fstreig%resid)
    if(associated(fstreig%filter)) deallocate(fstreig%filter)
    if(associated(fstreig%mass)) deallocate(fstreig%mass)
    if(associated(fstreig%effmass)) deallocate(fstreig%effmass)
    if(associated(fstreig%partfactor)) deallocate(fstreig%partfactor)
  end subroutine fstr_eigen_finalize
 
  !-----------------------------------------------------------------------------!
  !-----------------------------------------------------------------------------!
 
  subroutine fstr_setup_eigen( ctrl, P )
    implicit none
    integer(kind=kint) :: ctrl
    type(fstr_param_pack) :: P
 
    integer(kind=kint) :: rcode
 
    rcode = fstr_ctrl_get_eigen( ctrl, p%EIGEN%nget, p%EIGEN%tolerance, p%EIGEN%maxiter, p%EIGEN%sigma)
    if( rcode /= 0) call fstr_ctrl_err_stop
 
  end subroutine fstr_setup_eigen
 
end module m_fstr_eigen_setup