System.m

Go to the documentation of this file.

namespace models{
namespace golf{


/* (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 System
  :public models.BaseSecondOrderSystem {
    public:

        maxvforhole;

        
    public:

        System(models.BaseFullModel model) {
            this = this@models.BaseSecondOrderSystem(model);

            /*  Initial position */
            this.addParam(" x_0 ",1);
            this.addParam(" y_0 ",1);

            /*  Initial velocity */
            this.addParam(" vx_0 ",1);
            this.addParam(" vy_0 ",0);

            /*  Hole position */
            this.addParam(" hole_x ",0);
            this.addParam(" hole_y ",0);

            /*  #7: Gravity */
            this.addParam(" Gravity ",9.81);
            /*  #8: Friction on grass */
            this.addParam(" Grass friction ",0.03);

            this.NumStateDofs= 2;
            this.NumDerivativeDofs= 2;
            this.updateDimensions;

            x0 = dscomponents.AffineInitialValue;
            x0.addMatrix(" mu(1) ",[1;0]);
            x0.addMatrix(" mu(2) ",[0;1]);
            this.x0= x0;

            this.f= models.golf.Force(this);

            this.updateSparsityPattern;
        }


        function  prepareSimulation(colvec<double> mu,integer inputidx) {
            prepareSimulation@models.BaseSecondOrderSystem(this, mu, inputidx);
            m = this.Model;

            maxv = sqrt((1/(2*m.rad_ball)) * mu(7) * (2*m.rad_hole - m.rad_ball)^2 ); /*  m/s */

            this.maxvforhole= maxv;

            /*  Evtl. den Zeitschritt verkleinern, so dass der Ball nicht pro Schritt
             * die Lochweite überspringen kann ODER
             * Mindestens eine Auflösung von 24 Bildern pro sekunde, damit beim
             * Aviexport ein flüssiges Video entsteht! */
            dt = m.dt;
            if (dt > 2*m.rad_hole/maxv) || (dt > 1/24)
                this.MaxTimestep= min(0.9*2*m.rad_hole/maxv, 1/24);
            else
                this.MaxTimestep= dt;
            end
        }

    
    protected:

        function dv = getDerivativeDirichletValues() {
            dv = [];
        }


        function val = getDerivativeInitialValues(colvec<double> mu) {
            val = mu(3:4);
        }

    

};
}
}