/Users/craigcornelius/Projects/SPRING Mac Release 0.2/tetra.h

Go to the documentation of this file.

// $Id: tetra.h,v 1.4 2006/05/24 16:52:40 sean Exp $
// $Copyright: (c)2001 National Biocomputation Center, Stanford University $

#ifndef Tetra_H
#define Tetra_H


#include "point3d.h"
#include "util.h"

#include "node.h"

#include "nodearray.h"
#include "edgearray.h"
#include "facearray.h"
#include "tetraarray.h"

class Node;
class Edge;
class Face;
class NodeArray;
class EdgeArray;
class FaceArray;
class TetraArray;
class BoundingSphereLeaf;

class Tetra {
private:

        Node *node[4];
        Edge *edge[6];
        Face *face[4];

  TetraArray* tetraarray_p;
        
  int marker; // 0 off, 1 on
    
  BoundingSphereLeaf* boundingsphere;
        
public:
  Tetra();
  ~Tetra();
  static const char* rcsid;
  static int debug;

  void init(NodeArray* nodearray_p_in, EdgeArray* edgearray_p_in, 
        FaceArray* facearray_p_in, TetraArray* tetraarray_p_in, 
                                int i0, int i1, int i2, int i3);
                                
  void init(NodeArray* nodearray_p_in, EdgeArray* edgearray_p_in, 
        FaceArray* facearray_p_in,  TetraArray* tetraarray_p_in, 
                                Node* n0, Node* n1, Node* n2, Node* n3);
                                
  void init();
                                        
  int getIndex(void);
  TetraArray* getTetraArray();
        
  BoundingSphereLeaf* getBoundingSphere();
  void        setBoundingSphere(BoundingSphereLeaf* sphere);

  // links to other arrays
  int hasNodeLink(Node* node);
  Node* getNodeLink(int which); // 0,1,2,3
        Node* getOtherNode(Node* n0, Node* n1, Node* n2);
        void deleteNodeLink(Node* n);
                
  int hasEdgeLink(Edge* edge);  
  Edge* getEdgeLink(int which); // 0,1,2,3,4,5
        void deleteEdgeLink(Edge* e);
        Edge* getEdge(int n0, int n1);
        int getEdgeIndex(Edge* e);

        
        int hasFaceLink(Face* face);
  Face* getFaceLink(int which); // 0,1,2,3
        void deleteFaceLink(Face* f);
        
  void unlink();  // unlinks Tetra from its nodes and edges
                
  int isAdjacent(Tetra *t);
  
  void draw(int do_texture);    
  void drawlabel(float offset);
        void drawnormals();
        
  void drawMarker();
  void setMarker(int value);
  int getMarker(void);
        
  Point3D getCenter();
  double getRadius();
        double getAveRadius();
        double getVolume();
  void computenormals();

        int isFacingInwards(Face* f);

        
        int intersectsFace(Face* f, Point3D* intPt);
        int intersectsEdge(Edge* e, Point3D* intPt);
                
  int SanityCheck(TetraArray* real_fa_p);       // returns 1 if ok
  friend ostream& operator<<(ostream& os, const Tetra& f);
        
};

inline BoundingSphereLeaf* Tetra::getBoundingSphere()
{ return(boundingsphere); }

inline void Tetra::setBoundingSphere(BoundingSphereLeaf* sphere)
{ boundingsphere = sphere; }

inline Node* Tetra::getNodeLink(int which)
{ if ((which < 0) || (which > 4)) return(NULL);
  return(node[which]);
}

inline void Tetra::deleteNodeLink(Node* n)
{ for (int i=0; i<4; i++) 
        if (node[i] == n)
                node[i] = NULL;
}

inline Edge* Tetra::getEdgeLink(int which)
{ if ((which < 0) || (which > 6)) return(NULL);
  return(edge[which]);
}

inline void Tetra::deleteEdgeLink(Edge* e)
{ for (int i=0; i<6; i++) 
        if (edge[i] == e)
                edge[i] = NULL;
}

inline int Tetra::getEdgeIndex(Edge* e)
{ for (int i=0; i<6; i++) 
        if (edge[i] == e)
                return i;
                return -1;
}

inline Face* Tetra::getFaceLink(int which)
{ if ((which < 0) || (which > 4)) return(NULL);
  return(face[which]);
}

inline void Tetra::deleteFaceLink(Face* f)
{ for (int i=0; i<4; i++) 
        if (face[i] == f)
                face[i] = NULL;
}

inline double Tetra::getRadius()
{ double radius = 0;
  Point3D cent = getCenter();
  for (int i=0; i<4; i++) {
    double currad = node[i]->p.Dist(cent);
    if (currad > radius)
      radius = currad;
  }
  return radius;
}

inline double Tetra::getAveRadius()
{ 
        double radius = 0;
        Point3D cent = getCenter();
        for(int i=0; i<4; i++) {
          radius += node[i]->p.Dist(cent);
        }
        radius /= 4;
        return radius;
}

inline int Tetra::isAdjacent(Tetra *t)
{ 
  if (node[0]->hasTetraLink(t) || node[1]->hasTetraLink(t) ||
      node[2]->hasTetraLink(t) || node[4]->hasTetraLink(t) )
    return(1);
  return(0);
}

inline TetraArray* Tetra::getTetraArray()
{ return(tetraarray_p); }

inline void Tetra::setMarker(int value)
{ marker = value; }

inline int Tetra::getMarker()
{ return(marker); }

inline Point3D Tetra::getCenter()
{ return ((node[0]->p + node[1]->p + node[2]->p + node[3]->p)/4.0); }

                

#endif

---