propagator_exp_gauss2.F90

Go to the documentation of this file.

!! Copyright (C) 2024 Sebastian Ohlmann
!!
!! This program is free software; you can redistribute it and/or modify
!! it under the terms of the GNU General Public License as published by
!! the Free Software Foundation; either version 2, or (at your option)
!! any later version.
!!
!! This program is distributed in the hope that it will be useful,
!! but WITHOUT ANY WARRANTY; without even the implied warranty of
!! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!! GNU General Public License for more details.
!!
!! You should have received a copy of the GNU General Public License
!! along with this program; if not, write to the Free Software
!! Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
!! 02110-1301, USA.
!!
 
#include "global.h"
 
module propagator_exp_gauss2_oct_m
  use algorithm_oct_m
  use debug_oct_m
  use global_oct_m
  use propagator_oct_m
 
  implicit none
 
  private
  public ::                            &
    propagator_exp_gauss2_t
 
  type, extends(propagator_t) :: propagator_exp_gauss2_t
    private
  end type propagator_exp_gauss2_t
 
  interface propagator_exp_gauss2_t
    procedure propagator_exp_gauss2_constructor
  end interface propagator_exp_gauss2_t
 
  !# doc_start exp_gauss2_propagation_operations
  ! Specific exponential mid-point propagation operations identifiers
  character(len=ALGO_LABEL_LEN), public, parameter :: &
    EXP_GAUSS2_START        = 'EXP_GAUSS2_START',             &
    exp_gauss2_finish       = 'EXP_GAUSS2_FINISH',            &
    exp_gauss2_extrapolate  = 'EXP_GAUSS2_EXTRAPOLATE',      &
    exp_gauss2_propagate    = 'EXP_GAUSS2_PROPAGATE'
 
  ! Specific exponential mid-point propagation operations
  type(algorithmic_operation_t), public, parameter :: &
    OP_EXP_GAUSS2_START        = algorithmic_operation_t(exp_gauss2_start,        'Starting exponential with Gauss order 1'),  &
    op_exp_gauss2_finish       = algorithmic_operation_t(exp_gauss2_finish,       'Finishing exponential with Gauss order 1'), &
    op_exp_gauss2_extrapolate  = algorithmic_operation_t(exp_gauss2_extrapolate,  'Extrapolation step'), &
    op_exp_gauss2_propagate    = algorithmic_operation_t(exp_gauss2_propagate,    'Propagation step for exponential midpoint')
  !# doc_end
 
contains
 
  ! ---------------------------------------------------------
  function propagator_exp_gauss2_constructor(dt) result(this)
    real(real64),           intent(in) :: dt
    type(propagator_exp_gauss2_t), pointer   :: this
 
    push_sub(propagator_exp_gauss2_constructor)
 
    allocate(this)
 
    this%predictor_corrector = .false.
    this%start_operation = op_exp_gauss2_start
    this%final_operation = op_exp_gauss2_finish
 
    call this%add_operation(op_exp_gauss2_extrapolate)
    call this%add_operation(op_exp_gauss2_propagate)
    call this%add_operation(op_update_couplings)
    call this%add_operation(op_update_interactions)
    call this%add_operation(op_iteration_done)
    call this%add_operation(op_rewind_algorithm)
 
    this%algo_steps = 1
    this%dt = dt
 
    pop_sub(propagator_exp_gauss2_constructor)
  end function propagator_exp_gauss2_constructor
 
end module propagator_exp_gauss2_oct_m
 
 
!! Local Variables:
!! mode: f90
!! coding: utf-8
!! End: