mne_proj_data.cpp

Go to the documentation of this file.

//=============================================================================================================
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "mne_proj_data.h"
#include "mne_surface.h"
#include "mne_triangle.h"
 
#ifndef TRUE
#define TRUE 1
#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])
 
//=============================================================================================================
// USED NAMESPACES
//=============================================================================================================
 
using namespace MNELIB;
 
//=============================================================================================================
// DEFINE MEMBER METHODS
//=============================================================================================================
 
MNEProjData::MNEProjData(const MNELIB::MNESurface* s)
{
    a.resize(s->ntri);
    b.resize(s->ntri);
    c.resize(s->ntri);
    act.resize(s->ntri);
 
    const MNETriangle* tri = s->tris.data();
    for (int k = 0; 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;
}