pf_generic_spreadq0 Subroutine

  1. pf_utils.f90
  2. pf_mod_utils
  3. pf_generic_spreadq0

public subroutine pf_generic_spreadq0(this, lev, t0)

Generic routine to spread initial conditions Each sweeper can define its own spreadq0 or use this generic one

Arguments

Type IntentOptional AttributesName
class(pf_sweeper_t), intent(in) :: this
class(pf_level_t), intent(inout) :: lev

Level on which to spread
real(kind=pfdp), intent(in) :: t0

time at beginning of interval

Called by

proc~~pf_generic_spreadq0~~CalledByGraph proc~pf_generic_spreadq0 pf_generic_spreadq0 proc~imexq_spreadq0 imexQ_spreadq0 proc~imexq_spreadq0->proc~pf_generic_spreadq0 proc~imex_spreadq0 imex_spreadq0 proc~imex_spreadq0->proc~pf_generic_spreadq0 proc~misdcq_spreadq0 misdcQ_spreadq0 proc~misdcq_spreadq0->proc~pf_generic_spreadq0 proc~magpicard_spreadq0 magpicard_spreadq0 proc~magpicard_spreadq0->proc~pf_generic_spreadq0 proc~verlet_spreadq0 verlet_spreadq0 proc~verlet_spreadq0->proc~pf_generic_spreadq0

Contents

Source Code

pf_generic_spreadq0

Source Code

  subroutine pf_generic_spreadq0(this,lev, t0)
    class(pf_sweeper_t), intent(in)  :: this
    class(pf_level_t), intent(inout) :: lev  !!  Level on which to spread
    real(pfdp),       intent(in)    :: t0    !!  time at beginning of interval

    integer :: m, p

    !  Stick initial condition into first node slot
    call lev%Q(1)%copy(lev%q0)

    !  Evaluate F at first spot
    call lev%ulevel%sweeper%evaluate(lev, t0, 1)

    ! Spread F and solution to all nodes
    do m = 2, lev%nnodes
       call lev%Q(m)%copy(lev%Q(1))
       do p = 1, lev%ulevel%sweeper%npieces
         call lev%F(m,p)%copy(lev%F(1,p))
       end do
    end do
  end subroutine pf_generic_spreadq0