Base.m

Go to the documentation of this file.

namespace error{
namespace alpha{


/* (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 Base
  :public handle {
    protected:

        M1;

        
    public:

        Base(models.ReducedModel rm) {
            fm = rm.FullModel;
            if ~isempty(fm.Approx)
                /*  Get full d x N coeff matrix of approx function */
                Ma = fm.Approx.Expansion.Ma;
            else
                /*  Get full d x N coeff matrix of core function */
                Ma = fm.System.f.Expansion.Ma;
            end

            /*  Perform any offline computations/preparations
             * Only prepare matrices if projection is used */
            if ~isempty(rm.V)
                /*  Compute projection part matrices, without creating a
                 * d x d matrix (too big!) */
                M = Ma - rm.V*(rm.W^t*Ma);
                hlp = M^t*(rm.G*M);
                /*  Check if matrix needs to be made symmetric */
                if any(any(abs(hlp-hlp^t) > 1e-5))
                    hlp = (hlp + hlp^t)/2;
                    warning(" KerMor:errorest "," M1 matrix not sufficiently symmetric, updating (M+M "" )/2 ");
                end
                this.M1= hlp;

                this.inputOfflineComputations(rm, M);
                clear M;
            else
                /*  No projection means no projection error! */
                n = size(Ma,2);
                this.M1= zeros(n,n);
            end
        }
    public: /* ( Abstract ) */

        virtual function a = getAlpha(colvec phi,double ut,double t,colvec mu) = 0;
        virtual function  inputOfflineComputations(models.BaseFullModel model,matrix M) = 0;
};
}
}