#ifndef OERSTED_FINITEELEMENTMESH_H
#define OERSTED_FINITEELEMENTMESH_H

#include <vector>

#include "Eigen"

#include "Node.h"
#include "Triangle.h"

template<size_t M, size_t Q>
class QuadratureDerivatives {
private:
    SmallVector<M> Value[Q];
};

template<size_t Q>
class QuadratureDerivatives<2, Q> {
public:
    double dx(size_t i) { return d[i].x(); };

    double dy(size_t i) { return d[i].y(); };

private:
    XY d[Q];
};

template<size_t Q>
class QuadratureDerivatives<3, Q> {
public:
    double dx(size_t i) { return d[i].x(); };

    double dy(size_t i) { return d[i].y(); };

    double dz(size_t i) { return d[i].z(); };

private:
    XYZ d[Q];
};

template<size_t M, size_t N>
class FiniteElementMesh {
};

template<size_t P>
class FiniteElementMesh<2, P> {
public:
    FiniteElementMesh(std::vector<XY> nodes, std::vector<Triangle<P>> tris) : Nodes(nodes), Triangles(tris) {};

    std::vector<XY> const nodes() { return Nodes; };

    std::vector<Triangle<P>> const triangles() { return Triangles; };

protected:
    std::vector<XY> Nodes;
    std::vector<Triangle<P>> Triangles;
};

#endif //OERSTED_FINITEELEMENTMESH_H