interface.h

Go to the documentation of this file.

// Project Name: OpenMEEG (http://openmeeg.github.io)
// © INRIA and ENPC under the French open source license CeCILL-B.
// See full copyright notice in the file LICENSE.txt
// If you make a copy of this file, you must either:
// - provide also LICENSE.txt and modify this header to refer to it.
// - replace this header by the LICENSE.txt content.
 
#pragma once
 
#include <string>
#include <sstream>
#include <vector>
#include <cstdlib>
#include <ctime>
#include <om_common.h>
#include <mesh.h>
 
namespace OpenMEEG {
 
    class Mesh;

 
    class OrientedMesh {
    public:
 
        typedef enum { Normal=1, Opposite=-1 } Orientation ;
 
        OrientedMesh() {}
 
        OrientedMesh(Mesh& m,const Orientation o): meshptr(&m),orient(o) { }
 
              Mesh&  mesh()       { return *meshptr; } 
        const Mesh&  mesh() const { return *meshptr; } 
 
        int orientation() const { return orient; } 
 
        void change_orientation() { orient = -orient; }
 
    private:
 
        Mesh* meshptr;
        int   orient;
    };

 
    class OPENMEEG_EXPORT Interface {
    public:
 
        //        using OrientedMesh::Orientation;
 
        typedef std::vector<OrientedMesh> OrientedMeshes;

 
        Interface() { }

 
        Interface(const std::string& interfname): interface_name(interfname) { }
 
        const std::string& name() const { return interface_name; } 
 
              OrientedMeshes& oriented_meshes()       { return orientedmeshes; }
        const OrientedMeshes& oriented_meshes() const { return orientedmeshes; }
 
        bool outermost() const { return outermost_interface; } 
        void set_to_outermost(); 
 
        bool contains(const Vect3& p) const; 
 
        bool is_mesh_orientations_coherent(const bool doublechecked=false); 

 
        size_t nb_vertices() const {
            size_t nb = 0;
            for (const auto& omesh : oriented_meshes())
                nb += omesh.mesh().vertices().size();
            return nb;
        }

 
        size_t nb_triangles() const {
            size_t nb = 0;
            for (const auto& omesh : oriented_meshes())
                nb += omesh.mesh().triangles().size();
            return nb;
        }

 
        TrianglesRefs adjacent_triangles(const Triangle& t) const {
            TrianglesRefs triangles;
            for (const auto& omesh : oriented_meshes()) {
                const TrianglesRefs& tri = omesh.mesh().adjacent_triangles(t);
                triangles.insert(triangles.end(),tri.begin(),tri.end());
            }
            return triangles;
        }
 
    private:
 
        double solid_angle(const Vect3& p) const; 
 
        std::string    interface_name      = "";    
        bool           outermost_interface = false; 
        OrientedMeshes orientedmeshes;
    };

 
    typedef std::vector<Interface> Interfaces;
}