mne_triangle.cpp

Go to the documentation of this file.

//=============================================================================================================
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "mne_triangle.h"
 
#include <Eigen/Geometry>
 
//=============================================================================================================
// USED NAMESPACES
//=============================================================================================================
 
using namespace Eigen;
using namespace MNELIB;
 
//=============================================================================================================
// DEFINE MEMBER METHODS
//=============================================================================================================
 
MNETriangle::MNETriangle() = default;
 
//=============================================================================================================
 
void MNETriangle::compute_data()
{
    r12 = Eigen::Map<const Eigen::Vector3f>(r2) - Eigen::Map<const Eigen::Vector3f>(r1);
    r13 = Eigen::Map<const Eigen::Vector3f>(r3) - Eigen::Map<const Eigen::Vector3f>(r1);
 
    nn = r12.cross(r13);
    float size = nn.norm();
    if (size > 0)
        nn /= size;
    area = size / 2.0f;
 
    float sizey = r13.norm();
    if (sizey <= 0)
        sizey = 1.0f;
    ey = r13 / sizey;
 
    cent = (Eigen::Map<const Eigen::Vector3f>(r1)
          + Eigen::Map<const Eigen::Vector3f>(r2)
          + Eigen::Map<const Eigen::Vector3f>(r3)) / 3.0f;
 
    ex = ey.cross(nn);
}