read_hsx.cpp

Go to the documentation of this file.

#include "mex.h"
#include <unistd.h>
#include <string.h>
#include "EQuUs_MATLAB.h"

#ifdef __cplusplus
extern "C" {
#endif

#ifdef __INTEL_COMPILER
void hsx_equus_mp_read_hsx_( char* fname, long *fname_len, EQstruct* eq_struct);
#else
void __hsx_equus_MOD_read_hsx( char* fname, long *fname_len, EQstruct* eq_struct);
#endif  


#ifdef __cplusplus
 }
#endif




void mexFunction(int nlhs, mxArray *plhs[],
                 int nrhs, const mxArray *prhs[])
{
      /* input filename */
      char *filename;
      /* length of the filename */
      long buflen;
      /* field names of the output structure */
      char **fieldnames;  



/*-----------------------------------------------------------------------*/
/*     Check for proper number of arguments. */
      if(nrhs != 1) {
         mexErrMsgIdAndTxt ("siesta_EQuUs:hsx_m:nInput", "One input required.");
      }
      else if (nlhs > 1) {
         mexErrMsgIdAndTxt ("siesta_EQuUs:hsx_m:nOutput", "Too many output arguments.");
      }

/*     Validate inputs */
      
/*     Check that the input is a double.*/
      if ((mxIsChar(prhs[0]) != 1)) {
         mexErrMsgIdAndTxt ("siesta_EQuUs:hsx_m:NonString", "Input must be a string.");
      }

      /* get the length of the input string */
      buflen = (mxGetM(prhs[0]) * mxGetN(prhs[0])) + 1;


      /* copy the string data from prhs[0] into a C string.    */
      filename = mxArrayToString(prhs[0]);
#ifdef DEBUG
      mexPrintf( "%s\n", filename );
#endif

      if (access(filename, F_OK) == -1) {
        mexErrMsgIdAndTxt ("siesta_EQuUs:hsx_m:FileNotExist", "The given file %s does not exist.", filename);
      }

      EQelement eq_field, *eq_fields;
      EQstruct eq_hsx;


      /* Retriving the HSX structure */
#ifdef __INTEL_COMPILER
      hsx_equus_mp_read_hsx_( filename, &buflen, &eq_hsx);  
#else
       __hsx_equus_MOD_read_hsx( filename, &buflen, &eq_hsx);
#endif  

      /* obtaining the fieldnames */
      fieldnames =  getFieldnames( &eq_hsx );

      /* creating the MATLAB version of the output structure */
      plhs[0] = mxCreateStructMatrix(1, 1, eq_hsx.numfields, (const char**)fieldnames);

#ifdef DEBUG
      mexPrintf( "\nSetting the fields\n*********************************\n" );
#endif
      int idx;
      for( idx=0; idx<eq_hsx.numfields; idx++ ) {
#ifdef DEBUG
        mexPrintf( "setting field: %s:\n", fieldnames[idx] );
#endif
        setField( &eq_hsx.fields[idx], plhs[0] );
    }

}