AliasSampler_1D_t Derived Type

type, public :: AliasSampler_1D_t


Contents

Variables

Type-Bound Procedures

Source Code

AliasSampler_1D_t

Components

Type Visibility Attributes Name Initial
type(AliasSampler_3D_t), private :: alias_sampler

Type-Bound Procedures

procedure, public, :: shared_alloc => setupSamplerArray_1D

  • private subroutine setupSamplerArray_1D(this, nEntries, entrySize, name)

    Setup an array of samplers using a single shared resource (split into parts associated with one of them each). This only does the allocation. @param[in] nEntries number of samplers to initialise @param[in] entrySize number of values per sampler

    Arguments

    Type IntentOptional Attributes Name
    class(AliasSampler_1D_t) :: this
    integer, intent(in) :: nEntries
    integer, intent(in) :: entrySize
    character(len=*), intent(in) :: name

procedure, public, :: setup_entry => setupEntry_1D

  • private subroutine setupEntry_1D(this, iEntry, rank_with_info, arr)

    Initialise one sampler of an array @param[in] iEntry index of the entry to initialize @param[in] arr data to be loaded by that entry

    Arguments

    Type IntentOptional Attributes Name
    class(AliasSampler_1D_t), intent(inout) :: this
    integer, intent(in) :: iEntry
    integer, intent(in) :: rank_with_info
    real(kind=dp), intent(in) :: arr(:)

procedure, public, :: finalize => samplerArrayDestructor_1D

procedure, public, :: sample => aSample_1D

  • private subroutine aSample_1D(this, iEntry, tgt, prob)

    Draw a random element from 1:size(this%probs) with the probabilities listed in prob

    Arguments

    Type IntentOptional Attributes Name
    class(AliasSampler_1D_t), intent(in) :: this
    integer, intent(in) :: iEntry

    The index of the sampler.
    integer, intent(out) :: tgt

    The sampled value tgt.
    real(kind=dp), intent(out) :: prob

    The probability of sampling tgt.

generic, public, :: constrained_sample => constrained_sample_1D_nI, constrained_sample_1D_fast

  • private subroutine constrained_sample_1D_nI(this, i, contain, renorm, pos, tgt, prob)

    Draw a random element from 1:size(this%probs) with the probabilities listed in prob while adherring to constraints

    Arguments

    Type IntentOptional Attributes Name
    class(AliasSampler_1D_t), intent(in) :: this
    integer, intent(in) :: i

    The index of the sampler.
    integer, intent(in) :: contain(:)

    The constraint in nI format.
    real(kind=dp), intent(in) :: renorm

    The renormalization. (i.e. sum(this%get_prob(… contain…))
    integer, intent(out) :: pos

    The sampled value tgt and its position pos in `contain.
    integer, intent(out) :: tgt

    The sampled value tgt and its position pos in `contain.
    real(kind=dp), intent(out) :: prob

    The probability of sampling tgt from contain

  • private subroutine constrained_sample_1D_fast(this, i, contain, contain_ilut, renormalization, pos, tgt, prob)

    Draw a random element from 1:size(this%probs) with the probabilities listed in prob while adherring to constraints

    Arguments

    Type IntentOptional Attributes Name
    class(AliasSampler_1D_t), intent(in) :: this
    integer, intent(in) :: i

    The index of the sampler.
    integer, intent(in) :: contain(:)

    The constraint in nI format.
    integer(kind=n_int), intent(in) :: contain_ilut(0:)

    The constraint in ilut (bitmask) format
    real(kind=dp), intent(in) :: renormalization

    The renormalization. (i.e. sum(this%get_prob(… contain…))
    integer, intent(out) :: pos

    The sampled value tgt and its position pos in `contain.
    integer, intent(out) :: tgt

    The sampled value tgt and its position pos in `contain.
    real(kind=dp), intent(out) :: prob

    The probability of sampling tgt from contain

procedure, private :: constrained_sample_1D_fast

  • private subroutine constrained_sample_1D_fast(this, i, contain, contain_ilut, renormalization, pos, tgt, prob)

    Draw a random element from 1:size(this%probs) with the probabilities listed in prob while adherring to constraints

    Arguments

    Type IntentOptional Attributes Name
    class(AliasSampler_1D_t), intent(in) :: this
    integer, intent(in) :: i

    The index of the sampler.
    integer, intent(in) :: contain(:)

    The constraint in nI format.
    integer(kind=n_int), intent(in) :: contain_ilut(0:)

    The constraint in ilut (bitmask) format
    real(kind=dp), intent(in) :: renormalization

    The renormalization. (i.e. sum(this%get_prob(… contain…))
    integer, intent(out) :: pos

    The sampled value tgt and its position pos in `contain.
    integer, intent(out) :: tgt

    The sampled value tgt and its position pos in `contain.
    real(kind=dp), intent(out) :: prob

    The probability of sampling tgt from contain

procedure, private :: constrained_sample_1D_nI

  • private subroutine constrained_sample_1D_nI(this, i, contain, renorm, pos, tgt, prob)

    Draw a random element from 1:size(this%probs) with the probabilities listed in prob while adherring to constraints

    Arguments

    Type IntentOptional Attributes Name
    class(AliasSampler_1D_t), intent(in) :: this
    integer, intent(in) :: i

    The index of the sampler.
    integer, intent(in) :: contain(:)

    The constraint in nI format.
    real(kind=dp), intent(in) :: renorm

    The renormalization. (i.e. sum(this%get_prob(… contain…))
    integer, intent(out) :: pos

    The sampled value tgt and its position pos in `contain.
    integer, intent(out) :: tgt

    The sampled value tgt and its position pos in `contain.
    real(kind=dp), intent(out) :: prob

    The probability of sampling tgt from contain

procedure, public, :: get_prob => aGetProb_1D

  • private elemental function aGetProb_1D(this, iEntry, tgt) result(prob)

    Returns the probability to draw tgt from the sampler with index iEntry @param[in] iEntry index of the sampler to use @param[in] tgt the number for which we request the probability of sampling @return prob the probability of drawing tgt with the sample routine

    Arguments

    Type IntentOptional Attributes Name
    class(AliasSampler_1D_t), intent(in) :: this
    integer, intent(in) :: iEntry
    integer, intent(in) :: tgt

    Return Value real(kind=dp)

procedure, public, :: constrained_getProb => constrained_get_prob_1D

  • private pure function constrained_get_prob_1D(this, i, contain, renorm, tgt) result(prob)

    Returns the probability to draw tgt from the sampler with index iEntry @param[in] i Index of the sampler to use @param[in] constraint pick only elements from constraint @param[in] tgt the number for which we request the probability of sampling @return prob the probability of drawing tgt with the sample routine from constraint

    Arguments

    Type IntentOptional Attributes Name
    class(AliasSampler_1D_t), intent(in) :: this
    integer, intent(in) :: i
    integer, intent(in) :: contain(:)
    real(kind=dp), intent(in) :: renorm
    integer, intent(in) :: tgt

    Return Value real(kind=dp)

Source Code

    type AliasSampler_1D_t
        private
        type(AliasSampler_3D_t) :: alias_sampler
    contains
        ! constructor
        procedure :: shared_alloc => setupSamplerArray_1D
        procedure :: setup_entry => setupEntry_1D
        ! destructor
        procedure :: finalize => samplerArrayDestructor_1D
        ! get a random element and the generation probability from one of the samplers
        procedure :: sample => aSample_1D
        generic :: constrained_sample => constrained_sample_1D_nI, constrained_sample_1D_fast
        procedure, private :: constrained_sample_1D_nI, constrained_sample_1D_fast
        procedure :: get_prob => aGetProb_1D
        procedure :: constrained_getProb => constrained_get_prob_1D
    end type AliasSampler_1D_t