fstr_solve_eigen.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_solve_eigen
contains
 
  subroutine fstr_solve_eigen( hecMESH, hecMAT, fstrEIG, fstrSOLID, &
      & fstrRESULT, fstrPARAM, hecLagMAT)
    use hecmw_util
    use m_fstr
    use m_fstr_stiffmatrix
    use m_fstr_addbc
    use m_fstr_eig_setmass
    use m_fstr_eig_lanczos
    use m_fstr_eig_output
    use m_static_lib
    use m_hecmw2fstr_mesh_conv
    use m_fstr_spring
 
    implicit none
 
    type(hecmwst_local_mesh)  :: hecmesh
    type(hecmwst_matrix)      :: hecMAT
    type(fstr_solid)          :: fstrSOLID
    type(hecmwst_result_data) :: fstrRESULT
    type(fstr_param)          :: fstrPARAM
    type(fstr_eigen)          :: fstrEIG
    type(hecmwst_matrix_lagrange) :: hecLagMAT
 
    type(hecmwst_local_mesh), pointer :: hecMESHmpc
    type(hecmwst_matrix), pointer :: hecMATmpc
    real(kind=kreal) :: t1, t2
 
    t1 = hecmw_wtime()
 
    call hecmw_mpc_mat_init(hecmesh, hecmat, hecmeshmpc, hecmatmpc)
 
    fstrsolid%dunode = 0.0d0
    call fstr_stiffmatrix(hecmesh, hecmat, fstrsolid, 0.0d0, 0.0d0)
    call fstr_addspring(1, hecmesh, hecmat, fstrsolid, fstrparam)
 
    call fstr_addbc(1,  hecmesh, hecmat, fstrsolid, fstrparam, heclagmat, 2)
    call hecmw_mpc_mat_ass(hecmesh, hecmat, hecmeshmpc, hecmatmpc)
    call hecmw_mpc_trans_rhs(hecmesh, hecmat, hecmatmpc)
 
    call setmass(fstrsolid, hecmesh, hecmat, fstreig)
    call hecmw_mpc_trans_mass(hecmesh, hecmat, fstreig%mass)
 
    call fstr_solve_lanczos(hecmeshmpc, hecmatmpc, fstrsolid, fstreig)
 
    call hecmw_mpc_tback_eigvec(hecmesh, hecmat, fstreig%iter, fstreig%eigvec)
 
    call fstr_eigen_output(hecmesh, hecmat, fstreig)
 
    call fstr_eigen_make_result(hecmesh, hecmat, fstreig, fstrresult)
 
    call hecmw_mpc_mat_finalize(hecmesh, hecmat, hecmeshmpc, hecmatmpc)
 
    t2 = hecmw_wtime()
 
    if(myrank == 0)then
      write(imsg,'("### FSTR_SOLVE_EIGEN FINISHED!")')
      write(*,'("### FSTR_SOLVE_EIGEN FINISHED!")')
    endif
 
  end subroutine fstr_solve_eigen
end module m_fstr_solve_eigen