Constant.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 Constant
  :public error.alpha.Base {
    private:

        M2 = "[]";

        M3 = "[]";

        
    public:


        Constant(models.BaseFullModel model) {
            this = this@error.alpha.Base(model);
        }


        function  inputOfflineComputations(models.ReducedModel rm,matrix<double> M) {
            fm = rm.FullModel;
            if ~isempty(fm.System.B)
                try
                    B = fm.System.B.evaluate([],[]);
                catch ME/* #ok */

                    B = fm.System.B.evaluate(0,rm.getRandomParam);
                    warning(" Some:Id "," Error estimator for current system will not work correctly! (B is not linear and mu-independent! ");
                end

                if ~isempty(rm.V)
                    B2 = B - rm.V*(rm.W^t*B);
                    this.M2= M^t*(rm.G*B2);
                    this.M3= B2^t*(rm.G*B2);
                    clear B2;
                else
                    n = size(this.M1,2);
                    b = size(B,2);
                    this.M2= zeros(n,b);
                    this.M3= zeros(b,b);
                end
            end
        }
        function a = getAlpha(colvec phi,double ut,double t,colvec mu) {
            a = phi*this.M1*phi^t;
            if ~isempty(ut) /*  An input function u is set */

                a = a + phi*this.M2*ut + ut^t*this.M3*ut;
            end
            a = sqrt(abs(a));
            /* a = sqrt(max(a,0)); */
        }
};
}
}