u3d_triangle.hpp

/* u3d_triangle.hpp */
/* $Id: u3d_triangle.hpp,v 1.2 2008/05/13 12:07:38 sad Exp $ */
#ifndef __u3d_triangle_hpp
#define __u3d_triangle_hpp

#include "u3d_point.hpp"
#include <cmath>  // for sqrt()

#define COS_TRIANGLE_EPSILON 1.7e-7

//Temporaly:
#include <iostream>
using std::cout;
using std::endl;


struct U3dTriangle
{
  U3dTriangle(){};
  
  inline U3dTriangle(const int& iind1, const int& iind2, const int& iind3);


#if 0
  bool intersectWithSegmentOld(const U3dPoint& PTr0, const U3dPoint& PTr1, const U3dPoint& PTr2,
                const U3dPoint& pA, const U3dPoint& pB, 
                const double& A, const double& B, const double& C,
                double& outX, double& outY, double& outZ) const;




  bool intersectWithSegmentNew(const U3dPoint& PTr0, const U3dPoint& PTr1, const U3dPoint& PTr2,
                const U3dPoint& pA, const U3dPoint& pB, 
                const double& A, const double& B, const double& C,
                double& outX, double& outY, double& outZ) const;
#endif 


  bool intersectWithSegment(const U3dPoint& PTr0, const U3dPoint& PTr1, const U3dPoint& PTr2,
                const U3dPoint& pA, const U3dPoint& pB, 
                const double& A, const double& B, const double& C,
                double& outX, double& outY, double& outZ) const;



  int ind1, ind2, ind3; 
  double a, b, c, d;
};


inline U3dTriangle::U3dTriangle(const int& iind1, const int& iind2, const int& iind3) : 
  ind1(iind1), ind2(iind2), ind3(iind3)
{};



inline double vectorLength(const double& a1, const double& a2, const double& a3)
{
  return sqrt(a1*a1 + a2*a2 + a3*a3);
}


void calculateTrigABC(U3dTriangle& trig, const U3dPoint&, const U3dPoint&, const U3dPoint&);



typedef std::vector<U3dTriangle*> U3dTrianglesAr;



#endif /* u3d_triangle.hpp */