RotationDecorator.m

Go to the documentation of this file.

namespace spacereduction{


/* (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 RotationDecorator
  :public spacereduction.BaseSpaceReducer {
    public:

        Degree = .05;
        Dims = 5;
    private:

        sp;
    public:

        RotationDecorator(s) {
            if isempty(s)
                error(" Rotation decorator must be given an underlying space reduction class. ");
            end
            if ~strcmp(s.TargetDimensions," : ")
                error(" Rotation decorator not implemented for partial reduction (ReducableDims != :) ");
            end
            this.sp= s;
            this.registerProps(" Degree "," Dims ");
        }
    protected:


        function [V , W ] = generateReducedSpaceImpl(models.BaseFullModel model,subset) {

            /*  Call subclass reduction */
            [V, W] = this.sp.generateReducedSpaceImpl(model, subset);

            /* rnd = RandStream('mt19937ar','Seed',2564); */

            n = size(V,1); /* #ok<*PROP> */

            if this.Dims >= n/2
                warning(" KerMor:spacereduction:RotationDecorator "," Too large Dims property (%d) for full space dimension %d, setting to %d ",this.Dims,n,floor(n/2));
                this.Dims= floor(n/2);                
            end
            R = spdiags(ones(n,1),0,n,n);
            for idx = 1:this.Dims
                Q = spdiags(ones(n,1),0,n,n);
                /* idx1 = rnd.randi(n); */
                idx1 = 2*(idx-1)+1;
                /* idx2 = rnd.randi(n); */
                idx2 = idx1+1;

                Q(idx1,idx1) = cos(this.Degree);
                Q(idx1,idx2) = -sin(this.Degree);
                Q(idx2,idx1) = sin(this.Degree);
                Q(idx2,idx2) = cos(this.Degree);
                R = R*Q;
            end
            V = R*V;
            W = R*W;
        }
};
}