mne_proj_data.cpp

Go to the documentation of this file.

//=============================================================================================================
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "mne_proj_data.h"
#include "mne_surface_old.h"
#include "mne_triangle.h"
 
#define MALLOC_46(x,t) (t *)malloc((x)*sizeof(t))
 
#ifndef TRUE
#define TRUE 1
#endif
 
#ifndef FALSE
#define FALSE 0
#endif
 
#define X_46 0
#define Y_46 1
#define Z_46 2
 
#define VEC_DOT_46(x,y) ((x)[X_46]*(y)[X_46] + (x)[Y_46]*(y)[Y_46] + (x)[Z_46]*(y)[Z_46])
 
#define FREE_46(x) if ((char *)(x) != NULL) free((char *)(x))
 
//=============================================================================================================
// USED NAMESPACES
//=============================================================================================================
 
using namespace MNELIB;
 
//=============================================================================================================
// DEFINE MEMBER METHODS
//=============================================================================================================
 
MneProjData::MneProjData(MNELIB::MneSurfaceOld* s)
{
    int k;
    MneTriangle* tri = Q_NULLPTR;
 
    a   = MALLOC_46(s->ntri,float);
    b   = MALLOC_46(s->ntri,float);
    c   = MALLOC_46(s->ntri,float);
    act = MALLOC_46(s->ntri,int);
 
    for (k = 0, tri = s->tris; k < s->ntri; k++, tri++) {
      a[k] =  VEC_DOT_46(tri->r12,tri->r12);
      b[k] =  VEC_DOT_46(tri->r13,tri->r13);
      c[k] =  VEC_DOT_46(tri->r12,tri->r13);
 
      act[k] = TRUE;
    }
    nactive = s->ntri;
}
 
//=============================================================================================================
 
MneProjData::~MneProjData()
{
    FREE_46(a);
    FREE_46(b);
    FREE_46(c);
}