IJIntDerivation Derived Type

type, public, extends(Derivation) :: IJIntDerivation

Calculates the Shkarofsky I or J integrals for each passed harmonic.


Components

Type Visibility Attributes Name Initial
logical, public :: userDefined = .false.

True only if user explicitly sets it to true

type(SparseRowData), private :: ijMatBuffer

Sparse buffer used to calculate the integrals

integer(kind=ik), private :: numV

Type-Bound Procedures

procedure, public :: isDefined => isDefinedObject

  • interface

    private pure module function isDefinedObject(this) result(defined)

    Getter for userDefined

    Arguments

    Type IntentOptional Attributes Name
    class(Object), intent(in) :: this

    Return Value logical

procedure, public :: makeDefined => makeDefinedObject

  • interface

    private pure module subroutine makeDefinedObject(this)

    Set userDefined to .true.

    Arguments

    Type IntentOptional Attributes Name
    class(Object), intent(inout) :: this

procedure, public :: makeUndefined => makeUndefinedObject

  • interface

    private pure module subroutine makeUndefinedObject(this)

    Set userDefined to .false.

    Arguments

    Type IntentOptional Attributes Name
    class(Object), intent(inout) :: this

procedure, public :: init => initIJInt

  • interface

    private module subroutine initIJInt(this, vSpaceObj, ind, isJInt)

    Initialize Shkarofsky I/J integral derivation derivation

    Arguments

    Type IntentOptional Attributes Name
    class(IJIntDerivation), intent(inout) :: this
    type(VSpace), intent(in) :: vSpaceObj
    integer(kind=ik), intent(in) :: ind

    Index of indegral

    logical, intent(in), optional :: isJInt

    If true the lower triangular J integral is calculated instead of the I integral. Defaults to false.

procedure, public :: calculate => calculateIJInt

  • interface

    private module function calculateIJInt(this, inputArray, indices) result(output)

    Arguments

    Type IntentOptional Attributes Name
    class(IJIntDerivation), intent(inout) :: this
    type(RealArray), intent(in), dimension(:) :: inputArray
    integer(kind=ik), intent(in), dimension(:) :: indices

    Return Value real(kind=rk), allocatable, dimension(:)