eigen_driver.f

! Subroutine for computing time evolution through time slices
module QMDynamicSlice_driver_m

    use type_m
    use constants_m
    use parameters_m        , only : survival , driver , n_part
    use Data_Output         , only : Dump_stuff 
    use Schroedinger_m      , only : DeAllocate_QDyn
    use TDSE_adiabatic_m    , only : TDSE_adiabatic
    use CSDM_adiabatic_m    , only : CSDM_adiabatic
    use Chebyshev_driver_m  , only : Chebyshev_driver

    public :: QMDynamicSlice_driver

    private

contains
!
!
!================================
 subroutine QMDynamicSlice_driver
!================================
implicit none

! local variables ...
integer                :: it
real*8  , allocatable  :: QDyn_temp(:,:,:)
type(f_time)           :: QDyn

If( .NOT. survival ) then
    Print*, ">>> Survival = .FALSE. <<<"
    Stop
end If

select case ( DRIVER )

    case( "slice_AO" , "slice_FSSH" )

        CALL TDSE_adiabatic ( QDyn , it )

    case( "slice_CSDM" )

        CALL CSDM_adiabatic ( QDyn , it )

    case( "slice_Cheb" )

        CALL Chebyshev_driver ( QDyn , it )

    case default
        Print*, " >>> Check your QMDynamicSlice_driver options <<< :" , driver

end select

! prepare data for survival probability ...
allocate ( QDyn_temp( it , 0:size(QDyn%fragments)+1 , n_part ) , source=QDyn%dyn( 1:it , 0:size(QDyn%fragments)+1 , : ) )
CALL move_alloc( from=QDyn_temp , to=QDyn%dyn )

CALL Dump_stuff( QDyn=QDyn ) 

! final procedures ...
CALL DeAllocate_QDyn( QDyn , flag="dealloc" )

include 'formats.h'

end subroutine QMDynamicSlice_driver
!
!
!
end module QMDynamicSlice_driver_m