AliasSampler_1D_t Derived Type

type, public :: 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.

    Arguments

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

    number of samplers to initialise
    integer, intent(in) :: entrySize

    number of samplers to initialise number of values per sampler
    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

    Arguments

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

    index of the entry to initialize
    integer, intent(in) :: rank_with_info

    index of the entry to initialize
    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

    Arguments

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

    index of the sampler to use
    integer, intent(in) :: tgt

    the number for which we request the probability of sampling

    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

    Arguments

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

    Index of the sampler to use
    integer, intent(in) :: contain(:)
    real(kind=dp), intent(in) :: renorm
    integer, intent(in) :: tgt

    the number for which we request the probability of sampling

    Return Value real(kind=dp)

procedure, public :: get_memory_demand => get_memory_demand_AliasSampler_1D_t

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
        procedure, public :: get_memory_demand => get_memory_demand_AliasSampler_1D_t
    end type AliasSampler_1D_t