/*\file TimeAdapt.c *\brief: update time steps to respect CFL condition */ #include "./TimeAdapt.h" void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ /*input datasets: */ Elements* elements = NULL; Nodes *nodes = NULL; Vertices *vertices = NULL; Loads *loads = NULL; Materials *materials = NULL; Parameters *parameters = NULL; /*output*/ double dt; /*Boot module: */ MODULEBOOT(); /*checks on arguments on the matlab side: */ CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&TimeAdaptUsage); /*Input datasets: */ FetchData((DataSet**)&elements,ELEMENTSIN); FetchData((DataSet**)&nodes,NODESIN); FetchData((DataSet**)&vertices,VERTICESIN); FetchData((DataSet**)&loads,LOADSIN); FetchData((DataSet**)&materials,MATERIALSIN); FetchParams(¶meters,PARAMETERSIN); /*configure: */ elements-> Configure(elements,loads, nodes,vertices, materials,parameters); nodes-> Configure(elements,loads, nodes,vertices, materials,parameters); loads-> Configure(elements, loads, nodes,vertices, materials,parameters); /*call "x" code layer*/ TimeAdaptx(&dt,elements,nodes,vertices,loads, materials,parameters); /*write output datasets: */ WriteData(DT,dt); /*Free ressources: */ delete elements; delete nodes; delete vertices; delete loads; delete materials; delete parameters; /*end module: */ MODULEEND(); } void TimeAdaptUsage(void) { _printf_(true,"\n"); _printf_(true," usage: dt = %s(elements,nodes,vertices,loads,materials,parameters);\n",__FUNCT__); _printf_(true,"\n"); }