ICIAMSystem.m

Go to the documentation of this file.

namespace models{
namespace iciam2011{


/* (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 ICIAMSystem
  :public models.BaseFirstOrderSystem {
    public: /* ( setObservable ) */

        svNum;
    public: /* ( setObservable ) */


        ICIAMSystem(models.BaseFullModel model) {
            
            this = this@models.BaseFirstOrderSystem(model);

            this.registerProps(" svNum ");

            /* % System settings */
            dims = model.dim;
            this.NumStateDofs= dims;
            this.MaxTimestep= model.dt;

            /*  Sample bases */
            this.svNum= 20;
            space = linspace(0,50,this.svNum);

            f = dscomponents.ParamTimeKernelCoreFun(this);
            fe = f.Expansion;
            fe.Centers.xi= repmat(space,dims,1);
            fe.Centers.ti= [];
            fe.Centers.mui= [];

            /* % BaseCompWiseKernelApprox settings
             * Choose the kernel with such that each kernel vanishes (<
             * kerneleps) after kernelsupport support vectors.  */
            kernelsupport = 2;
            kerneleps = 0.00001;

            d = sqrt(dims)*space(2);
            /* gamma = -((kernelsupport*d)^2)/log(kerneleps) preprint */
            gamma = -((kernelsupport*d))/log(kerneleps);
            fe.Kernel= kernels.GaussKernel(gamma);
            fe.Kernel.G= 1;
            /* this.SystemKernel = kernels.GaussKernel(1); */
            fe.TimeKernel= kernels.NoKernel;
            fe.ParamKernel= kernels.NoKernel;

            this.f= f;
            this.updateDimensions;
            this.updateSparsityPattern;
        }


};
}
}