unitcube.cpp

Go to the documentation of this file.

#include "unitcube.h"

//UnitCube *UnitCube::_singleton=NULL;

 UnitCube::UnitCube() 
{
  //Set the unit vectors
  e1.reinit(3);
  e2.reinit(3);
  e3.reinit(3);

  e1=0;
  e2=0;
  e3=0;

  e1(0)=1;
  e2(1)=1;
  e3(2)=1;

  //Set the vertices
  _vertices.reinit(8,3);
  _vertices(0,0) = 0;  _vertices(0,1) = 0;  _vertices(0,2) = 0;
  _vertices(1,0) = 1;  _vertices(1,1) = 0;  _vertices(1,2) = 0;
  _vertices(2,0) = 0;  _vertices(2,1) = 1;  _vertices(2,2) = 0;
  _vertices(3,0) = 1;  _vertices(3,1) = 1;  _vertices(3,2) = 0;
  _vertices(4,0) = 0;  _vertices(4,1) = 0;  _vertices(4,2) = 1;
  _vertices(5,0) = 1;  _vertices(5,1) = 0;  _vertices(5,2) = 1;
  _vertices(6,0) = 0;  _vertices(6,1) = 1;  _vertices(6,2) = 1;
  _vertices(7,0) = 1;  _vertices(7,1) = 1;  _vertices(7,2) = 1;


  


}

 UnitCube::~UnitCube() 
{

}


unsigned proj_info[6][4]={
  {1,2,0,0},
  {1,2,0,1},
  {0,2,1,0},
  {0,2,1,1},
  {0,1,2,0},
  {0,1,2,1}};
  
  


void UnitCube::project_2d_points_into_faces(const VecDouble &pt2D,VecDouble *ppt3D,FaceDirection3D dir) 
{
  assert(pt2D.size() == 2);
  assert(ppt3D->size() == 3);
  VecDouble &pt3D = *ppt3D;
  pt3D(proj_info[dir][0])=pt2D(0);
  pt3D(proj_info[dir][1])=pt2D(1);
  pt3D(proj_info[dir][2])=static_cast<double>(proj_info[dir][3]);
  
  
}


void  UnitCube::face_base_vectors(Index faceId,VecDouble *v1,VecDouble *v2) 
{
  if (faceId < 2)
  {
    (*v1)=e2;
    (*v2)=e3;
    return;
  }
  if (faceId < 4)
  {
    (*v1)=e1;
    (*v2)=e3;
    return;
  }
  if (faceId < 7)
  {
    (*v1)=e1;
    (*v2)=e2;
    return;
  }
  assert(0); //Invalid faceId
  return;
}


// UnitCube& UnitCube::getUnitCube() 
// {
//   if (_singleton)
//     return _singleton;
//   else
//   {
//     _singleton = new UnitCube();
//     return _singleton;
//   }
// }

---