Mesh.hpp

Go to the documentation of this file.

//
// Created by matt on 17/11/2019.
//
 
#include <map>
#include <memory>
#include <ostream>
#include "Vertex.hpp"
#include "Triangle.hpp"
#include "IMesh.hpp"
 
#pragma once
// Whole program could (?) be sped up using DCEL https://en.wikipedia.org/wiki/Doubly_connected_edge_list
class Mesh : public IMesh {
private:
    int vertexAttributes; 
    int triangleAttributes; 
    int dimensions; 
    std::map<triInd, Triangle> triangles; 
    std::map<vertInd, Vertex> vertices; 

    std::map<edge, std::vector<triInd> > edges;
    std::map<vertInd, std::vector<triInd> > vertTri;
public:
    // Constructors
    Mesh() : vertexAttributes(0), triangleAttributes(0), dimensions(3) {};
    ~Mesh() {};
    // Functionality
    std::vector<Vec> resolvePoints(std::vector<vertInd> vertInds);
    bool isDelaunay();
    std::vector<edge> newEdges(triInd triInd, const std::vector<vertInd> &vert);
    void removeEdges(triInd triInd, const std::vector<edge> &rEdge);
    int containingTriangle(double x, double y);
    std::vector<int> adjacentTriangles(triInd triInd);
    void recalcCircum(vertInd vertInd);
    void updateVertTri(triInd triInd, std::vector<vertInd> vertInds);
    void removeVertTri(triInd triInd, std::vector<vertInd> rVertInds);
    template <typename T>
    double integrate(T func, double (*interp)(T, const Triangle&)) {
        double res = 0;
        for (int i = 0; i < triangles.size(); i++) {
            res += std::abs(interp(func, triangles[i]));
        }
        return res;
    };
    // Getters & Setters
    const std::map<int, Vertex> &getVertices() const;
    void setVertices(const std::map<int, Vertex> &vertices);
    const std::map<int, Triangle> &getTriangles() const;
    void setTriangles(const std::map<int, Triangle> &triangles);
    int getVertexAttributes() const;
    void setVertexAttributes(int verAttr);
    int getTriangleAttributes() const;
    void setTriangleAttributes(int triAttr);
    int getDimensions() const;
    const std::map<std::pair<int, int>, std::vector<int> > &getEdges() const;
    // Equality Operators

    bool operator==(const Mesh &rhs) const;
    bool operator!=(const Mesh &rhs) const;
    // Stream Operators
    friend std::ifstream &operator>>(std::ifstream &ifs, Mesh &mesh);
    friend std::ofstream &operator<<(std::ofstream &ofs, Mesh &mesh);
};