/*\file InterpFromMeshToMesh2d.c
 *\brief: bamg module.
 */
#include "./InterpFromMeshToMesh2d.h"

void InterpFromMeshToMesh2dUsage(void){/*{{{*/
	_pprintLine_("INTERFROMMESHTOMESH2D - interpolation from a 2d triangular mesh onto a list of point");
	_pprintLine_("");
	_pprintLine_("   This function is a multi-threaded mex file that interpolates a field");
	_pprintLine_("   defined on a Delaunay triangulation onto a list of point");
	_pprintLine_("");
	_pprintLine_("   Usage:");
	_pprintLine_("         data_interp=InterpFromMeshToMesh2d(index,x,y,data,x_interp,y_interp);");
	_pprintLine_("      or data_interp=InterpFromMeshToMesh2d(index,x,y,data,x_interp,y_interp,OPTIONS);");
	_pprintLine_("");
	_pprintLine_("      index             : index of the mesh where data is defined");
	_pprintLine_("      x,y               : coordinates of the nodes where data is defined");
	_pprintLine_("      data              : matrix holding the data to be interpolated onto the mesh. (one column per field)");
	_pprintLine_("      x_interp,y_interp : coordinates of the points onto which we interpolate.");
	_pprintLine_("      data_interp       : vector of mesh interpolated data.");
	_pprintLine_("      Available options :");
	_pprintLine_("         - 'default' : default value if point is outsite of triangulation (instead of linear interolation)");
	_pprintLine_("");
	_pprintLine_("   Example:");
	_pprintLine_("      load('temperature.mat');");
	_pprintLine_("      md.initialization.temperature=InterpFromMeshToMesh2d(index,x,y,temperature,md.mesh.x,md.mesh.y);");
	_pprintLine_("      md.initialization.temperature=InterpFromMeshToMesh2d(index,x,y,temperature,md.mesh.x,md.mesh.y,'default',253);");
	_pprintLine_("");
}/*}}}*/
WRAPPER(InterpFromMeshToMesh2d){

	/*Intermediaties*/
	int     *index              = NULL;
	double  *x_data             = NULL;
	double  *y_data             = NULL;
	double  *data               = NULL;
	int      nods_data,nels_data;
	int      M_data,N_data;
	double  *x_interp           = NULL;
	double  *y_interp           = NULL;
	int      N_interp;
	Options *options   = NULL;
	double  *data_interp = NULL;
	int      test1,test2,test;

	/*Boot module: */
	MODULEBOOT();

	/*checks on output arguments on the matlab side: */
	#ifdef _HAVE_MATLAB_MODULES_
	if(nlhs!=NLHS){
		InterpFromMeshToMesh2dUsage();
		_error_("InterpFromMeshToMesh2dUsage usage error");
	}
	#endif
	/*check on input arguments: */
	if(nrhs<NRHS){
		InterpFromMeshToMesh2dUsage();
		_error_("InterpFromMeshToMesh2dUsage usage error");
	}

	/*Fetch inputs: */
	FetchData(&index,&nels_data,&test,INDEX); if(test!=3) _error_("index should have 3 columns");
	FetchData(&x_data,&nods_data,X);          if(nods_data<3) _error_("there should be at least three points");
	FetchData(&y_data,&test,Y);               if(test!=nods_data) _error_("vectors x and y should have the same length");
	FetchData(&data,&M_data,&N_data,DATA);    if(M_data*N_data<1) _error_("data is empty");
	FetchData(&x_interp,&N_interp,XINTERP);   if(N_interp<1) _error_("no interpolation requested");
	FetchData(&y_interp,&test,YINTERP);       if(test!=N_interp) _error_("vectors x_interp and y_interp should have the same length");
	FetchData(&options,NRHS,nrhs,ARGUMENTS);

	/*Run core computations*/
	InterpFromMeshToMesh2dx(&data_interp,index,x_data,y_data,nods_data,nels_data,data,M_data,N_data,x_interp,y_interp,N_interp,options);

	/*Write data: */
	WriteData(DATAINTERP,data_interp,N_interp,N_data);

	/*end module: */
	MODULEEND();
}