ARBFKernel.m

Go to the documentation of this file.

namespace kernels{


/* (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 ARBFKernel
  :public kernels.BaseKernel {
    public: /* ( setObservable ) */

        double Gamma = 1;
    public: /* ( setObservable ) */

        ARBFKernel() {
            this = this@kernels.BaseKernel;
            this.IsRBF= true;
        }


        function K = evaluate(matrix<double> x,matrix<double> y) {
            K = this.evaluateScalar(sqrt(this.getSqDiffNorm(x, y)));
        }
        function bool = eq(B) {
            bool = eq@kernels.BaseKernel(A, B) && A.Gamma == B.Gamma;
        }


        function copy = clone(copy) {
            copy = clone@kernels.BaseKernel(this, copy);
            copy.Gamma= this.Gamma;
        }

    
    public: /* ( Sealed ) */ /* ( setObservable ) */

        function r = getSqDiffNorm(matrix<double> x,matrix<double> y) {
            if ~isempty(this.fP)
                x = x(this.fP,:);
            end
            sx = this.fG*x;
            n1sq = sum(x.*sx,1);
            n1 = size(x,2);
            if isempty(y)
                n2sq = n1sq;
                n2 = n1;
                y = x;
            else
                if ~isempty(this.fP)
                    y = y(this.fP,:);
                end
                n2sq = sum(y.*(this.fG*y),1);
                n2 = size(y,2);
            end;
            r = ((ones(n2,1)*n1sq)"  + ones(n1,1)*n2sq - 2*sx "*y);
        }
        /* % Getter & setter */
    public: /* ( setObservable ) */

        
#if 0 //mtoc++: 'set.Gamma'
function  Gamma(value) {
            if ~isreal(value) || ~isscalar(value) || value <= 0
                error(" Only positive scalar values allowed for Gamma. ");
            end
            this.Gamma= value;
        }

#endif

    public: /* ( Abstract ) */ /* ( setObservable ) */

        virtual function phir = evaluateScalar(matrix<double> r) = 0;
};
}