D:/Basile/VisualCppProjects/hydroMEX3/hydroMEX3/hydroMEX3.cpp

Go to the documentation of this file.

//----------------Includes for graphics------------------------------------------
#include "oGL_graphics/graphics_relative.h"
#include "oGL_graphics/solid_objects.h"
//---------------------------------------------------------------------------------

//----------------Includes for MATLAB MEX------------------------------------------
#include "math.h"
#include "mex.h"   
//---------------------------------------------------------------------------------

#include "string.h"

//----------------Includes for hydropter-------------------------------------------
#include "hydro_source/hydro_includes.h"
//---------------------------------------------------------------------------------


//-------------Includes for SUNDIALS CVODE solver:---------------------------------
#include "CVODEincludes.h"  // ODE solver includes
//---------------------------------------------------------------------------------

//-------------Include for timing functions----------------------------------------
#include "timing/timing.h"
//---------------------------------------------------------------------------------

#include "f.h"


#define MAXFLAGLEN 100

#ifndef PI_DEF
#define pi 3.141592653589793
#define PI 3.141592653589793
#endif





void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
        parameters Par; 
        parameters *P = &Par; 
        joy_parameters jPar;  
        joy_parameters *jP = &jPar;  
        datas dat;  
        datas *d = &dat; 
        ODE_data f_dats;  
        ODE_data *f_dat = &f_dats;  


        int kk;

        double *state;
        char* flag;
        char def_flag[]  = "";
        char output_flag[]  = "output";
        long tLength;

        float Time;
        bool go=1;

        state = get_state(nlhs, plhs, nrhs, prhs);

        f_dat->d = d;
        f_dat->P = P;




        f_dat->state = state;


        initParameters(P);
        get_Parameters(nlhs, plhs, nrhs, prhs, P);

        get_time(nlhs, plhs, nrhs, prhs, f_dat);
        
        get_joystick_parameters(nlhs, plhs, nrhs, prhs, jP);
        pass_joy_parameters(jP);

        //flag = get_flag(nlhs, plhs, nrhs, prhs); //no more flag argument...


        data_init(d);
        data_update_state(d, state);
        data_update_param(d, P);



        solver_init(f_dat);

        kk = draw_graphics_init(d);//, gd);


        wave_init();

        Boat_create(d);


        tic();
        while (go)
        {
                Time = toc();
                solve_ODE(f_dat, Time);
                go = draw_graphics(d);
                if ((f_dat->time_param != 0.0) && (d->t > f_dat->time_param)) go=0;
        }

        free_wave_variables();

        kk = draw_graphics_kill();

        solver_free(f_dat);

        return_args(nlhs, plhs, nrhs, prhs, f_dat);


        delete_state(f_dat->state);
        Boat_delete();
        

        return;
}

---