#include "mex.h"
#include <stdio.h>
/*#define NUMBER_OF_STRUCTS (sizeof(regularity)/sizeof(struct reg))
#define NUMBER_OF_FIELDS (sizeof(field_reg)/sizeof(*field_reg))*/
#define NUMBER_OF_STRUCTS_POS (sizeof(position)/sizeof(struct pos))
#define NUMBER_OF_FIELDS_POS (sizeof(field_pos)/sizeof(*field_pos))
/*struct reg
{
double *time,*amp2,*amp3,*amp4,*alphalinha,*alpha1,*alpha2,*alpha3,*alpha4;
double *thalpha1,*thalpha2,*thalpha3,*th2,*th3,*th4,*thlinha,*escolhidos;
};*/
struct pos
{
double Pon, *P, *Poff, *QRSon, *Q, *R, *Fiducial, *qrs, *Rprima, *S, *QRSoff;
double *Ton, *T, *Tprima, *Toff, *Ttipo;
};
void mexFunction(
int nlhs,
mxArray *plhs[],
int nrhs,
const mxArray *prhs[]
)
{
/* Input variables declaration */
char *sigdir, *headir, *ecgnr, *matdir, *anot;
double ft, lead, qrs_flag;
double *time;
/* Struct array "regularity" declaration */
/*const char *field_reg[] = {"time", "amp2", "amp3", "amp4", "alphalinha",
"alpha1", "alpha2", "alpha3", "alpha4", "thalpha1", "thalpha2", "thalpha3",
"th2", "th3", "th4", "thlinha", "escolhidos"};
struct reg regularity[] = {{NULL},{NULL},{NULL},{NULL},{NULL},{NULL},
{NULL},{NULL},{NULL},{NULL},{NULL},{NULL},{NULL},{NULL},{NULL},{NULL},
{NULL}};
int dims[2] = {1, NUMBER_OF_STRUCTS};
double *time_field, *amp2_field, *amp3_field, *amp4_field, *alphalinha_field;
double *alpha1_field, *alpha2_field, *alpha3_field, *alpha4_field, *thalpha1_field, *thalpha2_field, *thalpha3_field;
double *th2_field, *th3_field, *th4_field, *thlinha_field, *escolhidos_field;*/
/* Struct array "position" output declaration */
const char *field_pos[] = {"Pon", "P", "Poff", "QRSon", "Q", "R", "Fiducial", "qrs", "Rprima",
"S", "QRSoff", "Ton", "T", "Tprima", "Toff", "Ttipo"};
struct pos position[] = {{NULL},{NULL},{NULL},{NULL},{NULL},{NULL},{NULL},{NULL},
{NULL},{NULL},{NULL},{NULL},{NULL},{NULL},{NULL},{NULL}};
int dims_pos[2] = {1, NUMBER_OF_STRUCTS_POS};
double *Pon_field, *P_field, *Poff_field, *QRSon_field, *Q_field, *R_field;
double *Fiducial_field, *qrs_field, *Rprima_field, *S_field, *QRSoff_field;
double *Ton_field, *T_field, *Tprima_field, *Toff_field, *Ttipo_field;
/* Checking inputs */
if (nrhs < 9)
{
mexErrMsgTxt ("Wavedet: nine input arguments required.");
}
if (!mxIsChar(prhs[0]))
mexErrMsgTxt("Signal directory must be a string.");
else if (mxGetM(prhs[0]) != 1)
mexErrMsgTxt("Signal directory must be a string in a row vector.");
if (!mxIsChar(prhs[1]))
mexErrMsgTxt("Header signal directory must be a string.");
else if (mxGetM(prhs[1]) != 1)
mexErrMsgTxt("Header signal directory must be a string in a row vector.");
if (!mxIsChar(prhs[2]))
mexErrMsgTxt("Output directory must be a string.");
else if (mxGetM(prhs[2]) != 1)
mexErrMsgTxt("Output directory must be a string in a row vector.");
if (!mxIsChar(prhs[3]))
mexErrMsgTxt("Name of the signal must be a string.");
else if (mxGetM(prhs[3]) != 1)
mexErrMsgTxt("Name of the signal must be a string in a row vector.");
if (mxIsChar(prhs[4]))
mexErrMsgTxt("Format file must be an integer: 0 (MIT), 1 (LUND) or 2 (mat files).");
else if ((mxGetM(prhs[4]) != 1)||(mxGetN(prhs[4]) != 1))
mexErrMsgTxt("Format file must be an integer!!");
else
{
ft = mxGetScalar(prhs[4]);
if ((ft != 1)&&(ft != 0)&&(ft != 2))
mexErrMsgTxt("Format file must be an integer: 0, 1 or 2.");
}
if (!mxIsChar(prhs[5]))
mexErrMsgTxt("Annotation output file must be a string.");
else if (mxGetM(prhs[5]) != 1)
mexErrMsgTxt("Annotation output file must be a string in a row vector.");
if (mxIsChar(prhs[6]))
mexErrMsgTxt("Number of lead to analyze must be an integer.");
else if ((mxGetM(prhs[6]) != 1)||(mxGetN(prhs[6]) != 1))
mexErrMsgTxt("Number of lead to analyze must be an integer.");
lead = mxGetScalar(prhs[6]);
if (mxIsChar(prhs[7]))
mexErrMsgTxt("Time vector must be a 1x2 matrix.");
else if ((mxGetM(prhs[7]) != 1)||(mxGetN(prhs[7]) != 2))
mexErrMsgTxt("Time vector must be a 1x2 matrix.");
time = mxGetPr(prhs[7]);
if (time[1]<=time[0])
mexErrMsgTxt("Bad time vector, first element greater or equal than second one");
if (mxIsChar(prhs[8]))
mexErrMsgTxt("Qrs_flag must be an integer: 1 (External), 0 (Internal) QRS Detector.");
else
{
qrs_flag = mxGetScalar(prhs[8]);
if ((qrs_flag != 1)&&(qrs_flag != 0))
mexErrMsgTxt("Qrs_flag must be a number: 0 or 1.");
}
/* End of checking inputs */
sigdir = mxArrayToString(prhs[0]);
headir = mxArrayToString(prhs[1]);
matdir = mxArrayToString(prhs[2]);
ecgnr = mxArrayToString(prhs[3]);
anot = mxArrayToString(prhs[5]);
// plhs[0] = mxCreateString(anot);
//plhs[0] = mxCreateDoubleMatrix(1,1,mxREAL);
//lead = mxGetPr(plhs[0]);
//lead = mxGetPr(prhs[6]);
//mexPrintf("%f & %f",time[0],time[1]);
/* Create an structure output array */
plhs[0] = mxCreateStructArray(2, dims_pos, NUMBER_OF_FIELDS_POS, NULL);
wavedet(sigdir,headir,matdir,ecgnr,anot,ft,lead,qrs_flag,time);
mexPrintf("%f & %f",time[0],time[1]);
}
void wavedet(char *sigdir, char *headir, char *matdir, char *ecgnr, char *anot,
double ft, double lead, double qrs_flag, double *time)
{
}