EpsSelector.m

Go to the documentation of this file.

namespace data{
namespace selection{


/* (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 EpsSelector
  :public data.selection.ASelector {
    public: /* ( setObservable ) */

        EpsRad = 5;
        SubspaceProject = true;
    public: /* ( setObservable ) */

        EpsSelector() {
            this = this@data.selection.ASelector;
            this.registerProps(" EpsRad "," SubspaceProject ");
            this.EpsRad= 1;
        }


        function copy = clone() {
            copy = data.selection.EpsSelector;
            /* copy = clone@data.selection.ASelector(this, copy); */
            copy.EpsRad= this.EpsRad;
            copy.SubspaceProject= this.SubspaceProject;
        }


        
#if 0 //mtoc++: 'set.EpsRad'
function  EpsRad(value) {
            if ~isposrealscalar(value)
                error(" EpsRad must be a positive real scalar. ");
            end
            this.EpsRad= value;
        }

#endif

    
    protected: /* ( setObservable ) */

        function [matrixxi , rowvecti , matrixmui , fxi ] = select(models.BaseFullModel model) {
            td = model.Data.TrajectoryData;
            if td.getNumTrajectories > 1
                warning(" KerMor:approx:selection:EpsSelector "," The epsilon-selector can only be used with one trajectory at the moment. Ignoring all but the first. ");
            end
            [x, mu] = td.getTrajectoryNr(1); 
            if this.SubspaceProject && ~isempty(model.SpaceReducer)
                [V, W] = model.assembleProjectionMatrices;
                x = V*(W^t*x);
            end
            /* [m,M] = Utils.getBoundingBox(x);
             *d = norm(M-m) / 10;
             *this.EpsRad = d; */
            selidx = 1;
            idx = 1; 
            cur = x(:,idx);
            while idx < size(x,2)
                while(idx < size(x,2) && norm(x(:,idx)-cur) < this.EpsRad)
                    idx = idx+1;
                end
                selidx(end+1) = idx;/* #ok */

                cur = x(:,idx);
                idx = idx+1;
            end
            xi = x(:,selidx);
            ti = model.Times(selidx);
            mui = repmat(mu,1,length(selidx));
            fxi = [];
        }
};
}
}