/*\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(&parameters,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");
}