BaseSampler.m

Go to the documentation of this file.

namespace sampling{


/* (Autoinserted by mtoc++)
 * This source code has been filtered by the mtoc++ executable,
 * which generates code that can be processed by the doxygen documentation tool.
 *
 * On the other hand, it can neither be interpreted by MATLAB, nor can it be compiled with a C++ compiler.
 * Except for the comments, the function bodies of your M-file functions are untouched.
 * Consequently, the FILTER_SOURCE_FILES doxygen switch (default in our Doxyfile.template) will produce
 * attached source files that are highly readable by humans.
 *
 * Additionally, links in the doxygen generated documentation to the source code of functions and class members refer to
 * the correct locations in the source code browser.
 * However, the line numbers most likely do not correspond to the line numbers in the original MATLAB source files.
 */

class BaseSampler
  :public KerMorObject {
    public:

        sampling.Domain Domain = "[]";
    public:

        function samples = generateSamples(models.BaseFullModel model) {
            sys = model.System;
            if sys.ParamCount == 0
                samples = [];
            else
                params = 1:sys.ParamCount;
                if ~isempty(model.TrainingParams)
                    params = params(model.TrainingParams);
                end
                samples = this.performSampling(sys.Params(params));
                if ~isempty(model.TrainingParams)
                    full = repmat(model.DefaultMu,1,size(samples,2));
                    full(model.TrainingParams,:) = samples;
                    samples = full;
                end
            end
        }


        
#if 0 //mtoc++: 'set.Domain'
function  Domain(value) {
            if ~isempty(value) && ~isa(value," sampling.Domain ")
                error(" The Domain property must be a sampling.Domain instance. ");
            end
            this.Domain= value;
        }

#endif

    
    public: /* ( Abstract ) */

        virtual function samples = performSampling() = 0;
    protected: /* ( Static ) */

        static function obj = loadobj(obj,from) {
            if ~isa(obj," sampling.BaseSampler ")
                if nargin < 2
                    error(" Invalid call to loadobject method. Must pass the struct to instantiate the class from. ");
                end
                if isfield(from," Domain ")
                    obj.Domain= from.Domain;
                else
                    obj.Domain= [];
                end
                obj = loadobj@KerMorObject(obj, from);
            else
                obj = loadobj@KerMorObject(obj);
            end
        }

    
    public: /* ( Static ) */

        static function res = test_SubsetSampling() {
            m = models.BaseFullModel;
            s = models.BaseFirstOrderSystem(m);
            m.System=s;
            s.addParam(" param_a ",.5);
            s.addParam(" param_b ",1);
            s.addParam(" param_c ",0," Range ",[0,1]);
            s.addParam(" param_d ",2);

            m.TrainingParams= [1 3];
            m.DefaultMu= rand(4,1);

            m.Sampler= sampling.RandomSampler;
            m.Sampler.Samples= 500;
            samples = m.Sampler.generateSamples(m);

            res = size(samples,1) == 4;
            res = res & all(samples(2,:) == m.DefaultMu(2));
            res = res & all(samples(4,:) == m.DefaultMu(4));

            m.Sampler= sampling.GridSampler;
            samples = m.Sampler.generateSamples(m);
            res = res & size(samples,1) == 4;
            res = res & all(samples(2,:) == m.DefaultMu(2));
            res = res & all(samples(4,:) == m.DefaultMu(4));
        }
};
}