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_creatematrix_and_dampingforce
    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
    use hecmw_ebc_defer
 
    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
    type(hecmwst_ebc) :: hecEBC
    real(kind=kreal), pointer :: mass_orig(:)
    real(kind=kreal) :: t1, t2
 
    t1 = hecmw_wtime()
 
    call hecmw_mpc_mat_init(hecmesh, hecmat, hecmeshmpc, hecmatmpc)
 
    fstrsolid%dunode = 0.0d0
    call fstr_creatematrix_and_dampingforce(hecmesh, hecmat, fstrsolid, 0.0d0, 0.0d0)
    call fstr_addspring(1, hecmesh, hecmat, fstrsolid, fstrparam)
 
    call hecmw_ebc_init(hecmat, hecebc)
    call fstr_addbc(1,  hecmesh, hecmat, fstrsolid, fstrparam, heclagmat, 2, hecebc)
    call hecmw_mpc_mat_ass(hecmesh, hecmat, hecmeshmpc, hecmatmpc)
    call hecmw_mpc_trans_rhs(hecmesh, hecmat, hecmatmpc)
    call hecmw_ebc_apply(hecmeshmpc, hecmatmpc, hecebc)
 
    call setmass(fstrsolid, hecmesh, hecmat, fstreig)
    mass_orig => fstreig%mass
    allocate(fstreig%mass(size(mass_orig)))
    fstreig%mass(:) = mass_orig(:)
    call hecmw_mpc_trans_mass(hecmesh, hecmat, hecmatmpc, fstreig%mass)
 
    call fstr_solve_lanczos(hecmeshmpc, hecmatmpc, fstrsolid, fstreig, hecebc)
    call hecmw_ebc_finalize(hecebc)
    call fstr_eigen_residual(hecmeshmpc, hecmatmpc, fstreig)
 
    call hecmw_mpc_tback_eigvec(hecmesh, hecmat, fstreig%iter, fstreig%eigvec)
    ! the eigenvectors are back in the original space, so is the mass used by the postprocessing
    deallocate(fstreig%mass)
    fstreig%mass => mass_orig
 
    call fstr_eigen_output(hecmesh, hecmat, fstreig)
 
    if( iresult.eq.1 ) then
      call fstr_eigen_output_result(hecmesh, fstreig)
    end if
    if( ivisual.eq.1 ) then
      call fstr_eigen_output_visual(hecmesh, fstreig, fstrresult)
    end if
 
    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