Tetrahedron.hpp

Go to the documentation of this file.

//  This file is part of ff3d - http://www.freefem.org/ff3d
//  Copyright (C) 2001, 2002, 2003 Stéphane Del Pino

//  This program is free software; you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation; either version 2, or (at your option)
//  any later version.

//  This program is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.

//  You should have received a copy of the GNU General Public License
//  along with this program; if not, write to the Free Software Foundation,
//  Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.  

//  $Id: Tetrahedron.hpp,v 1.9 2007/06/09 10:37:06 delpinux Exp $

#ifndef TETRAHEDRON_HPP
#define TETRAHEDRON_HPP

#include <Cell.hpp>
#include <Triangle.hpp>

#include <TinyMatrix.hpp>

class Tetrahedron
  : public Cell
{
public:
  enum {
    NumberOfVertices   = 4,     
    NumberOfFaces      = 4,     
    NumberOfEdges      = 6      
  };

  typedef Triangle FaceType;

  static const size_t faces[NumberOfFaces][FaceType::NumberOfVertices];
  static const size_t edges[NumberOfEdges][Edge::NumberOfVertices];

  virtual Cell::Type type() const
  {
    return Cell::tetrahedron;
  }

  inline void getBarycentricCoordinates(const TinyVector<3, real_t>& X,
                                        TinyVector<4, real_t>& lambda) const
  {
    TinyVector<4, real_t> X2;
    for (size_t i=0; i<3; ++i) X2[i] = X[i];
    X2[3] = 1;

    TinyMatrix<4, 4, real_t> M = 0;
    for (size_t i=0; i<4; ++i) {
      const Vertex& V = (*this)(i);
      for (size_t j=0; j<3; ++j) {
        M(j, i) = V[j];
      }
    }

    for (size_t i=0; i<4; ++i) {
      M(3,i) = 1;
    }

    lambda = X2/M;
  }

  size_t numberOfVertices() const
  {
    return NumberOfVertices;
  }

  size_t numberOfEdges() const
  {
    return NumberOfEdges;
  }

  inline const Tetrahedron& operator=(const Tetrahedron& T)
  {
    Cell::operator=(T);
    return *this;
  }

  Tetrahedron()
    : Cell(Tetrahedron::NumberOfVertices)
  {
    ;
  }

  Tetrahedron(Vertex& x0,
              Vertex& x1,
              Vertex& x2,
              Vertex& x3,
              const size_t& ref = 0);

  ~Tetrahedron()
  {
    ;
  }

};

#endif // TETRAHEDRON_HPP

---