[6130] | 1 | /*\file TimeAdapt.c
|
---|
| 2 | *\brief: update time steps to respect CFL condition
|
---|
| 3 | */
|
---|
| 4 |
|
---|
| 5 | #include "./TimeAdapt.h"
|
---|
| 6 |
|
---|
| 7 | void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){
|
---|
| 8 |
|
---|
| 9 | /*input datasets: */
|
---|
| 10 | Elements* elements = NULL;
|
---|
| 11 | Nodes *nodes = NULL;
|
---|
| 12 | Vertices *vertices = NULL;
|
---|
| 13 | Loads *loads = NULL;
|
---|
| 14 | Materials *materials = NULL;
|
---|
| 15 | Parameters *parameters = NULL;
|
---|
| 16 |
|
---|
| 17 | /*output*/
|
---|
| 18 | double dt;
|
---|
| 19 |
|
---|
| 20 | /*Boot module: */
|
---|
| 21 | MODULEBOOT();
|
---|
| 22 |
|
---|
| 23 | /*checks on arguments on the matlab side: */
|
---|
| 24 | CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&TimeAdaptUsage);
|
---|
| 25 |
|
---|
| 26 | /*Input datasets: */
|
---|
| 27 | FetchData((DataSet**)&elements,ELEMENTSIN);
|
---|
| 28 | FetchData((DataSet**)&nodes,NODESIN);
|
---|
| 29 | FetchData((DataSet**)&vertices,VERTICESIN);
|
---|
| 30 | FetchData((DataSet**)&loads,LOADSIN);
|
---|
| 31 | FetchData((DataSet**)&materials,MATERIALSIN);
|
---|
| 32 | FetchParams(¶meters,PARAMETERSIN);
|
---|
| 33 |
|
---|
| 34 | /*configure: */
|
---|
| 35 | elements-> Configure(elements,loads, nodes,vertices, materials,parameters);
|
---|
| 36 | nodes-> Configure(elements,loads, nodes,vertices, materials,parameters);
|
---|
| 37 | loads-> Configure(elements, loads, nodes,vertices, materials,parameters);
|
---|
| 38 |
|
---|
| 39 | /*call "x" code layer*/
|
---|
| 40 | TimeAdaptx(&dt,elements,nodes,vertices,loads, materials,parameters);
|
---|
| 41 |
|
---|
| 42 | /*write output datasets: */
|
---|
| 43 | WriteData(DT,dt);
|
---|
| 44 |
|
---|
| 45 | /*Free ressources: */
|
---|
| 46 | delete elements;
|
---|
| 47 | delete nodes;
|
---|
| 48 | delete vertices;
|
---|
| 49 | delete loads;
|
---|
| 50 | delete materials;
|
---|
| 51 | delete parameters;
|
---|
| 52 |
|
---|
| 53 | /*end module: */
|
---|
| 54 | MODULEEND();
|
---|
| 55 | }
|
---|
| 56 |
|
---|
| 57 | void TimeAdaptUsage(void)
|
---|
| 58 | {
|
---|
| 59 | _printf_("\n");
|
---|
| 60 | _printf_(" usage: dt = %s(elements,nodes,vertices,loads,materials,parameters);\n",__FUNCT__);
|
---|
| 61 | _printf_("\n");
|
---|
| 62 | }
|
---|