src/tetmesh/face.cpp

Go to the documentation of this file.

/***************************************************************************
                          face.cpp  -  description
                             -------------------
    begin                : Fri Dec 12 2003
    copyright            : (C) 2003 by Roman Geus
    email                : roman.geus@psi.ch
***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   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 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/

#include "matrix.h"
#include "linalg.h"
#include "vector3.h"
#include "point.h"
#include "face.h"
#include "tetmesh.h"

namespace mesh {

using namespace colarray;

Face::Face(id_t id, id_t node0, id_t node1, id_t node2, id_t tet) :
    Entity(id)
{ 
    _tets[0] = tet;
    _tets[1] = ID_NONE;
}

Point* Face::get_corner(int i) const {
    return get_mesh()->get_point(get_corner_id(i)); 
}

Edge* Face::get_edge(int i) const {
        return get_mesh()->get_edge(get_edge_id(i));    
}

id_t Face::get_corner_id(int i) const {
    Tet* tet = get_mesh()->get_tet(get_tet_id());
    
    // Query info of parent Tet 
    int face_lid;
    for (face_lid = 0; face_lid < 4; ++ face_lid)
        if (tet->get_face_id(face_lid) == get_id())
            break;
    assert(face_lid < 4);
    int node_lids[3];
    tet->get_face_point_lids(face_lid, node_lids[0], node_lids[1], node_lids[2]);
    return tet->get_corner_id(node_lids[i]);
}

id_t Face::get_edge_id(int i) const {
    Tet* tet = get_mesh()->get_tet(get_tet_id());
    
    // Query info of parent Tet 
    int face_lid;
    for (face_lid = 0; face_lid < 4; ++ face_lid)
        if (tet->get_face_id(face_lid) == get_id())
            break;
    assert(face_lid < 4);
    int edge_lids[3];
    tet->get_face_edge_lids(face_lid, edge_lids[0], edge_lids[1], edge_lids[2]);
    return tet->get_edge_id(edge_lids[i]);
}

TetMesh* Face::get_mesh() const {
    return TetMesh::get_instance();
}

} // namespace mesh

---