domain.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 <interface.h>
 
namespace OpenMEEG {

 
    class SimpleDomain {
    public:
 
        typedef enum { Inside, Outside } Side;
 
         SimpleDomain() { }
         SimpleDomain(Interface& i,const Side s): interf(i),side(s) { }
        ~SimpleDomain() { }
 
              Interface& interface()       { return interf;  }
        const Interface& interface() const { return interf;  }
 
        bool inside() const { return (side==Inside); }
 
        // omesh must belong to the interface.
 
        int mesh_orientation(const OrientedMesh& omesh) const {
            return (inside()) ? omesh.orientation() : -omesh.orientation();
        }
 
    private:
 
        Interface interf;
        Side      side;
    };

 
    class OPENMEEG_EXPORT Domain {
    public:
 
        typedef std::vector<SimpleDomain> Boundaries;
 
         Domain(const std::string& dname=""): domain_name(dname) { }
        ~Domain() { }

 
              Boundaries& boundaries()       { return bounds; }
        const Boundaries& boundaries() const { return bounds; }

 
              std::string& name()       { return domain_name; }
        const std::string& name() const { return domain_name; }

 
        void set_conductivity(const double c) { cond = c;    }
        bool has_conductivity() const { return cond != -1.0;    }
        const double& conductivity() const    { return cond; }

 
        void info(const bool outermost=false) const;
 
        bool contains(const Mesh& m) const { return mesh_orientation(m)!=0; }
 
        bool contains(const Vect3& point) const; 

 
        int mesh_orientation(const Mesh& m) const {
            for (const auto& boundary : boundaries())
                for (const auto& omesh : boundary.interface().oriented_meshes())
                    if (&omesh.mesh()==&m)
                        return boundary.mesh_orientation(omesh);
            return 0;
        }
 
    private:
 
        Boundaries  bounds;             
        std::string domain_name = "";   
        double      cond        = -1.0; 
    };

 
    typedef std::vector<Domain> Domains;
}