Burgers.m

Go to the documentation of this file.

namespace models{
namespace burgers{


/* (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 Burgers
  :public models.BaseFullModel {
    public: /* ( Constant ) */

        static const Omega = "[0 1]";

    
    public: /* ( Dependent ) */

        Dimension;

        
    private:

        fDim;

        
    public:

        PlotAzEl = "[]";
    public:

        Burgers(integer dim,integer version) {
            this = this@models.BaseFullModel;
            if nargin < 2
                version = 1;
                if nargin < 1
                    dim = 2000;
                end
            end
            this.T= 1;
            this.dt= .01;

            if version == 1
                this.System= models.burgers.BurgersSys(this);
                this.ODESolver= solvers.MLode15i;
                this.SaveTag= sprintf(" burgers_1D_d%d_combRHS ",dim);
                this.Name= sprintf(" 1D-%dd unsteady Burgers equation (combined RHS) ",dim);
            elseif version == 2
                this.System= models.burgers.BurgersSys_A(this);
                this.System.MaxTimestep= this.dt;
                this.ODESolver= solvers.SemiImplicitEuler(this);
                this.System.MaxTimestep= this.dt;
                this.SaveTag= sprintf(" burgers_1D_d%d ",dim);
                this.Name= sprintf(" 1D-%dd unsteady Burgers equation ",dim);
            end
            this.Dimension= dim;

            this.SpaceReducer= spacereduction.PODGreedy;
            this.SpaceReducer.Eps= 1e-9;

            a = approx.DEIM(this.System);
            a.MaxOrder= 80;
            this.Approx= a;

            this.ErrorEstimator= error.DEIMEstimator;
        }
        function  plot(double t,matrix<double> y,pm_ax) {
            if nargin < 4
                pm_ax = PlotManager;
                pm_ax.LeaveOpen= true;
            end
            nt = length(t);
            y  = [zeros(nt,1) y^t zeros(nt,1)]; /*  add boundaries */

            xx = linspace(this.Omega(1), this.Omega(2), this.fDim+2);

            if ~ishandle(pm_ax)
                pm_ax = pm_ax.nextPlot(" burgers ",sprintf(" %s: dim=%d ",this.Name,this.fDim)," x "," t ");
            end
            surfc(pm_ax,xx,t,y);
            shading interp;
            zlabel(" y ");
            rotate3d on;
            if ~isempty(this.PlotAzEl)
                view(this.PlotAzEl);
            end
        }


        
#if 0 //mtoc++: 'set.Dimension'
function  Dimension(value) {
            this.fDim= value;
            this.System.newDim;
            this.Data.SimCache.clearTrajectories;
        }

#endif


        
#if 0 //mtoc++: 'get.Dimension'
function dim = Dimension() {
            dim = this.fDim;
        }

#endif

};
}
}