pf_misdcQ_t Derived Type

  1. pf_misdcQ.f90
  2. pf_mod_misdcQ
  3. pf_misdcQ_t

type, public, abstract, extends(pf_sweeper_t) :: pf_misdcQ_t

Multi-implicit SDC sweeper type, extends abstract sweeper


Inherits

type~~pf_misdcq_t~~InheritsGraph type~pf_misdcq_t pf_misdcQ_t type~pf_sweeper_t pf_sweeper_t type~pf_misdcq_t->type~pf_sweeper_t type~pf_encap_t pf_encap_t type~pf_misdcq_t->type~pf_encap_t I3, rhs

Contents

Variables

Type-Bound Procedures

Source Code

pf_misdcQ_t

Components

TypeVisibility AttributesNameInitial
integer, public :: npieces
logical, public :: use_LUq
real(kind=pfdp), public, allocatable:: QdiffE(:,:)
real(kind=pfdp), public, allocatable:: QdiffI(:,:)
real(kind=pfdp), public, allocatable:: QtilE(:,:)
real(kind=pfdp), public, allocatable:: QtilI(:,:)
real(kind=pfdp), public, allocatable:: dtsdc(:)
class(pf_encap_t), public, allocatable:: I3(:)
class(pf_encap_t), public, allocatable:: rhs

Type-Bound Procedures

procedure(pf_f_eval_p), public, deferred :: f_eval

  • subroutine pf_f_eval_p(this, y, t, level_index, f, piece) Prototype

    This is the interface for the routine to compute the RHS function values Evaluate f_piece(y), where piece is one or two Evaluate f_piece(y), where piece is one or two

    Arguments

    Type IntentOptional AttributesName
    class(pf_misdcQ_t), intent(inout) :: this
    class(pf_encap_t), intent(in) :: y

    Argument for evaluation
    real(kind=pfdp), intent(in) :: t

    Time at evaluation
    integer, intent(in) :: level_index

    Level index
    class(pf_encap_t), intent(inout) :: f

    RHS function value
    integer, intent(in) :: piece

    Which piece to evaluate

procedure(pf_f_comp_p), public, deferred :: f_comp

  • subroutine pf_f_comp_p(this, y, t, dtq, rhs, level_index, f, piece) Prototype

    Solve the equation where n is given by the argument piece

    Arguments

    Type IntentOptional AttributesName
    class(pf_misdcQ_t), intent(inout) :: this
    class(pf_encap_t), intent(inout) :: y

    Solution of implicit solve

    real(kind=pfdp), intent(in) :: t

    Time of solve
    real(kind=pfdp), intent(in) :: dtq

    dt*quadrature weight
    class(pf_encap_t), intent(in) :: rhs

    RHS for solve
    integer, intent(in) :: level_index

    Level index
    class(pf_encap_t), intent(inout) :: f

    f_n of solution y
    integer, intent(in) :: piece

    Which piece to evaluate

procedure, public :: sweep => misdcQ_sweep

  • public subroutine misdcQ_sweep(this, pf, level_index, t0, dt, nsweeps, flags)

    Assign level pointer

    Read more…

    Arguments

    Type IntentOptional AttributesName
    class(pf_misdcQ_t), intent(inout) :: this
    type(pf_pfasst_t), intent(inout), target:: pf

    PFASST structure
    integer, intent(in) :: level_index

    which level to sweep on
    real(kind=pfdp), intent(in) :: t0

    Time at beginning of time step
    real(kind=pfdp), intent(in) :: dt

    time step size
    integer, intent(in) :: nsweeps

    number of sweeps to do
    integer, intent(in), optional :: flags

procedure, public :: initialize => misdcQ_initialize

procedure, public :: integrate => misdcQ_integrate

  • public subroutine misdcQ_integrate(this, lev, qSDC, fSDC, dt, fintSDC, flags)

    Arguments

    Type IntentOptional AttributesName
    class(pf_misdcQ_t), intent(inout) :: this
    class(pf_level_t), intent(in) :: lev
    class(pf_encap_t), intent(in) :: qSDC(:)
    class(pf_encap_t), intent(in) :: fSDC(:,:)
    real(kind=pfdp), intent(in) :: dt
    class(pf_encap_t), intent(inout) :: fintSDC(:)
    integer, intent(in), optional :: flags

procedure, public :: residual => misdcQ_residual

  • public subroutine misdcQ_residual(this, lev, dt, flags)

    Subroutine to compute Residual

    Arguments

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

    Current level
    real(kind=pfdp), intent(in) :: dt

    Time step
    integer, intent(in), optional :: flags

procedure, public :: spreadq0 => misdcQ_spreadq0

  • public subroutine misdcQ_spreadq0(this, lev, t0, flags, step)

    Arguments

    Type IntentOptional AttributesName
    class(pf_misdcQ_t), intent(inout) :: this
    class(pf_level_t), intent(inout) :: lev
    real(kind=pfdp), intent(in) :: t0
    integer, intent(in), optional :: flags
    integer, intent(in), optional :: step

procedure, public :: evaluate_all => misdcQ_evaluate_all

  • public subroutine misdcQ_evaluate_all(this, lev, t, flags, step)

    Subroutine to evaluate the function values at all nodes

    Arguments

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

    Current level
    real(kind=pfdp), intent(in) :: t(:)

    Array of times at each node
    integer, intent(in), optional :: flags
    integer, intent(in), optional :: step

procedure, public :: evaluate => misdcQ_evaluate

  • public subroutine misdcQ_evaluate(this, lev, t, m, flags, step)

    Subroutine to evaluate function value at node m

    Arguments

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

    Current level
    real(kind=pfdp), intent(in) :: t

    Time at which to evaluate
    integer, intent(in) :: m

    Node at which to evaluate
    integer, intent(in), optional :: flags
    integer, intent(in), optional :: step

procedure, public :: destroy => misdcQ_destroy

procedure, public :: misdcQ_destroy

Source Code

  type, extends(pf_sweeper_t), abstract :: pf_misdcQ_t
     real(pfdp), allocatable :: QdiffE(:,:)
     real(pfdp), allocatable :: QdiffI(:,:)
     real(pfdp), allocatable :: QtilE(:,:)
     real(pfdp), allocatable :: QtilI(:,:)
     real(pfdp), allocatable :: dtsdc(:)

     class(pf_encap_t), allocatable :: I3(:)
     class(pf_encap_t), allocatable :: rhs
   contains 
     procedure(pf_f_eval_p), deferred :: f_eval
     procedure(pf_f_comp_p), deferred :: f_comp
     procedure :: sweep        => misdcQ_sweep
     procedure :: initialize   => misdcQ_initialize
     procedure :: integrate    => misdcQ_integrate
     procedure :: residual   => misdcQ_residual
     procedure :: spreadq0   => misdcQ_spreadq0
     procedure :: evaluate_all => misdcQ_evaluate_all
     procedure :: evaluate   => misdcQ_evaluate
     procedure :: destroy      => misdcQ_destroy
     procedure :: misdcQ_destroy
  end type pf_misdcQ_t