gateway.c

著名的Jade算法

#include <math.h>
#include "mex.h"
  
#include <stdio.h>
#include <stdlib.h>


/* #include "Matutil.c" repository of once useful matrix-related functions */

#include "JnS.c"  /* The real thing is maintained in an other file */


/* -- The gateway ------------------------------------------------ */


/* I/O Arguments */
#define	Xin	prhs[0]
#define	Bout	plhs[0]


void mexFunction(
                 int nlhs,       mxArray *plhs[],
                 int nrhs, const mxArray *prhs[]
		 )
{
  double	*X;  /* the input data matrix */
  double	*B;  /* the output separating matrix */

  unsigned int	datadim ;    /* number of sensors */
  unsigned int	samplesize ; /* number of data points */
  
  /* Check for proper number of arguments */
  if (nrhs != 1) mexErrMsgTxt("Requires a unique input argument.");
  if (nlhs >  1) mexErrMsgTxt("Returns a single matrix.");
  
  
  /* Check the dimensions of X  */
  datadim    = mxGetM(Xin);
  samplesize = mxGetN(Xin);

  /*   if ( */
  /*       !mxIsNumeric(Y_IN)  */
  /*       || mxIsComplex(Y_IN) */
  /*       || mxIsSparse(Y_IN) */
  /*       || !mxIsDouble(Y_IN)  */
  /*       )  */
  /*     mexErrMsgTxt("Ragnagna..."); */
  
  
  /* Create a matrix for the return argument */
  Bout = mxCreateDoubleMatrix(datadim, datadim, mxREAL);
  
  
  /* Get pointers to the I/O parameters */
  B = mxGetPr(Bout);
  X = mxGetPr(Xin);
  
  /* Do it */
  ICA_ALGORITHM (B, X, (int) datadim, (int) samplesize);

  /* Note: this is a stupid trick to avoid duplicating the gateway code, */

  return;

}