PCDModel.m

Go to the documentation of this file.

namespace models{
namespace pcd{


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

        static const d1 = 1.8e-11;
        static const d2 = 1.86e-11;
        static const d3 = 1.89e-11;
        static const d4 = 2.27e-11;
        static const L = 1e-5;
        static const n = 2;
        static const Kc1_real = 0.08;
        static const Kc2_real = 0.08;
        static const Kd1_real = .0005;
        static const Kd2_real = .0005;
        static const Kd3_real = .0005;
        static const Kd4_real = .0005;
        static const Kp1_real = 0.0001;
        static const Kp2_real = 0.0001;
        static const xa0 = 1e-7;
        static const ya0 = 1e-7;
        static const xi0 = 1e-7;
        static const yi0 = 1e-7;
        static const SteadyStates = "[[0, 9.8153e-4, 0.1930]*models.pcd.PCDModel.xa0   \
                        [0, 3.0824e-5, 0.1713]*models.pcd.PCDModel.ya0   \
                        [.2, 0.1990, 0.0070]*models.pcd.PCDModel.xi0   \
                        [.2, 0.2, 0.0287]*models.pcd.PCDModel.yi0]";
    public: /* ( Dependent ) */

        Dimension;

        
    public: /* ( Dependent ) */

        PCDModel(dim) {

            if nargin == 0
                dim = 1;
            end

            this.T= 3600; /* [s] */

            this.dt= 5; /* [s]
             * time scaling */

            this.tau= this.L^2/this.d1;

            if dim > 10000
                this.Data.useFileTrajectoryData;
            end

/*              s = sampling.RandomSampler;
 *             s.Samples = 10; */
            s = sampling.GridSampler;
            this.Sampler= s;

/*              this.ODESolver = solvers.MLWrapper(@ode23);
 *             this.ODESolver = solvers.ExplEuler(this.dt);
 *             o = solvers.MLode15i;
 *             o.AbsTol = 1e-6;
 *             o.RelTol = 1e-5;
 *             o.MaxStep = [];
 *             this.ODESolver = o; */
            this.ODESolver= solvers.SemiImplicitEuler(this);

            switch dim
                case 3
                    error(" Fixme! Currently not working :-( ");
                    s = models.pcd.PCDSystem3D(this);
                    this.Name= " Programmed Cell Death 3D ";
                case 2 
                    s = models.pcd.PCDSystem2D(this);
                    this.Name= " Programmed Cell Death 2D ";
                otherwise
                    s = models.pcd.PCDSystem1D(this);
                    this.Name= " Programmed Cell Death 1D ";
            end
            s.MaxTimestep= this.dt;            
            this.System= s;

            /*  Space reduction setup */
            sr = spacereduction.PODGreedy;
            sr.Eps= 1e-5;
            this.SpaceReducer= sr;

            this.Approx= [];
        }
        function  plot(varargin) {
            this.System.plot(this, varargin[:]);
        }
        function  plotState(varargin) {
            this.System.plotState(this, varargin[:]);
        }
#if 0 //mtoc++: 'get.Dimension'
function value = Dimension() {
            value = this.System.Dims;
        }

#endif

};
}
}