KernelLS.m

Go to the documentation of this file.

namespace general{
namespace regression{


/* (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 KernelLS
  :public KerMorObject,
   public IKernelCoeffComp {
        /*  @change{0,4,sa,2011-05-07} Implemented Setter for the properties G and Epsilon */
    
    public: /* ( setObservable ) */

        K;
        double lambda = 1;
        CGMaxIt = "[]";
        CGTol = 1e-6;
        MaxStraightInvDim = 2000;
    public: /* ( setObservable ) */


        KernelLS() {
            this = this@KerMorObject;
            this.registerProps(" K "," lambda "," CGMaxIt "," CGTol "," MaxStraightInvDim ");
        }


        function copy = clone() {
            copy = general.regression.KernelLS;
            copy = clone@IKernelCoeffComp(this, copy);
            copy.MaxStraightInvDim= this.MaxStraightInvDim;
            copy.CGTol= this.CGTol;
            copy.CGMaxIt= this.CGMaxIt;
            copy.lambda= this.lambda;
            copy.K= this.K;
        }


        function [a , sf ] = regress(fx,ainit) {
            n = size(this.K,1);

            y = (fx*this.K)^t;
            M = this.K^t*this.K + this.lambda*speye(n, n);

            if length(y) <= this.MaxStraightInvDim
                a = M\y;
                sf = StopFlag.SUCCESS;
            else
                [a, flag] = bicg(M,y,this.CGTol, this.CGMaxIt);/* #ok
                 *| @todo return correct stop flag */

                sf = StopFlag.SUCCESS;
            end
        }


        function  init(kernels.KernelExpansion kexp) {
            this.K= kexp.getKernelMatrix;
        }
        function [ai , svidx , sf ] = computeKernelCoefficients(yi,initialai) {
            [ai, sf] = this.regress(yi);
            svidx = 1:size(yi,2);
        }


        
#if 0 //mtoc++: 'set.K'
function  K(value) {
            if (isa(value," handle ") && ~isa(value, " data.FileMatrix ")) || (~ismatrix(value) || ~isa(value," double "))
                error(" Value must be a data.FileMatrix or a double matrix. ");
            end
            this.K= value;
        }

#endif

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

#endif


        
#if 0 //mtoc++: 'set.CGMaxIt'
function  CGMaxIt(value) {
            if ~isempty(value) && ~isposintscalar(value)
                error(" value must be a positive integer scalar ");                
            end
            this.CGMaxIt= value;
        }

#endif


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

#endif


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

#endif

};
}
}