SampleRandomizedSequence.h

Go to the documentation of this file.

//
// Class SampleRandomizedSequence
//   This class creates randomized sequences, i.e., it randomly picks values from
//   a sequence.
//
// Copyright (c) 2019, Matthias Frey, Paul Scherrer Institut, Villigen PSI, Switzerland
// All rights reserved
//
// Implemented as part of the PhD thesis
// "Precise Simulations of Multibunches in High Intensity Cyclotrons"
//
// This file is part of OPAL.
//
// OPAL is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// You should have received a copy of the GNU General Public License
// along with OPAL. If not, see <https://www.gnu.org/licenses/>.
//
#ifndef OPAL_SAMPLE_WEIGHTED_SEQUENCE_H
#define OPAL_SAMPLE_WEIGHTED_SEQUENCE_H
 
#include "Sample/Uniform.h"
 
template <typename T>
class SampleRandomizedSequence : public SamplingMethod
{
 
public:
 
    SampleRandomizedSequence(T lower, T upper, double step)
        : unif_m(0, size_t((upper - lower) / step))
        , lower_m(lower)
        , step_m(step)
    { }
 
    SampleRandomizedSequence(T lower, T upper, double step, size_t seed)
        : unif_m(0, size_t((upper - lower) / step), seed)
        , lower_m(lower)
        , step_m(step)
    { }
 
    void create(std::shared_ptr<SampleIndividual>& ind, size_t i) {
        size_t idx = unif_m.getNext();
        ind->genes[i] = static_cast<T>(lower_m + idx * step_m);
    }
 
    void allocate(const CmdArguments_t& args, const Comm::Bundle_t& comm) {
        unif_m.allocate(args, comm);
    }
 
private:
    std::vector<T> points_m;
    Uniform<size_t> unif_m;
    T lower_m;
    double step_m;
};
 
#endif