ngh_interface.hh

Go to the documentation of this file.

/*
 * ngh_interface.hh
 *
 *  Created on: Apr 18, 2014
 *      Author: jb
 */

#ifndef NGH_INTERFACE_HH_
#define NGH_INTERFACE_HH_

#include "mesh/elements.h"
#include "mesh/ngh/include/abscissa.h"
#include "mesh/ngh/include/point.h"
#include "mesh/ngh/include/triangle.h"
#include "mesh/ngh/include/tetrahedron.h"

namespace ngh {
    /**
     * Create tetrahedron from element
     */
    inline void set_tetrahedron_from_element( TTetrahedron &te, Element *ele) {
        ASSERT(( ele->dim() == 3 ), "Dimension of element must be 3!\n");

         te.SetPoints(TPoint(ele->node[0]->point()(0), ele->node[0]->point()(1), ele->node[0]->point()(2)),
                     TPoint(ele->node[1]->point()(0), ele->node[1]->point()(1), ele->node[1]->point()(2)),
                     TPoint(ele->node[2]->point()(0), ele->node[2]->point()(1), ele->node[2]->point()(2)),
                     TPoint(ele->node[3]->point()(0), ele->node[3]->point()(1), ele->node[3]->point()(2)) );
    }

    /**
     * Create triangle from element
     */
    inline void set_triangle_from_element(TTriangle &tr, const Element *ele) {
        ASSERT(( ele->dim() == 2 ), "Dimension of element must be 2!\n");

        tr.SetPoints(TPoint(ele->node[0]->point()(0), ele->node[0]->point()(1), ele->node[0]->point()(2)),
                     TPoint(ele->node[1]->point()(0), ele->node[1]->point()(1), ele->node[1]->point()(2)),
                     TPoint(ele->node[2]->point()(0), ele->node[2]->point()(1), ele->node[2]->point()(2)) );
    }

    /**
     * Create abscissa from element
     */
    inline void set_abscissa_from_element(TAbscissa &ab, const Element *ele) {
        ASSERT(( ele->dim() == 1 ), "Dimension of element must be 1!\n");

        ab.SetPoints(TPoint(ele->node[0]->point()(0), ele->node[0]->point()(1), ele->node[0]->point()(2)),
                     TPoint(ele->node[1]->point()(0), ele->node[1]->point()(1), ele->node[1]->point()(2)) );
    }

    /**
     * Create point from element
     */
    inline void set_point_from_element(TPoint &p, const Element *ele) {
        ASSERT(( ele->dim() == 0 ), "Dimension of element must be 0!\n");

        p.SetCoord( ele->node[0]->point()(0), ele->node[0]->point()(1), ele->node[0]->point()(2) );
    }

} // namespace ngh

#endif /* NGH_INTERFACE_HH_ */