utils.hpp

Go to the documentation of this file.

#include <string>
#include <stdexcept>
#include "Triangle.hpp"
namespace Utils {
    int lineLength(const std::string &line);
 
    template<typename Stream>
    void loadFile(Stream &gStream, const std::string &fileName) {
        gStream.open(fileName.c_str(), std::ios::binary);
        if (not gStream.is_open()) {
            std::stringstream ss;
            ss << "Could not open file \"" << fileName << "\".";
            throw std::runtime_error(ss.str());
        }
    }
 
    // View documentation at https://mattnotmitt.github.io/DelaunayTriangulation/namespaceUtils.html#aea35446b52e730679709eefcdf6fc8e2
    template<typename T>
    double constantValueApprox(T func, const Triangle &triangle) {
        return triangle.area() * func(triangle.getCc().x, triangle.getCc().y);
    }
 
    // View documentation at https://mattnotmitt.github.io/DelaunayTriangulation/namespaceUtils.html#a5795be8719aeb33bdfe657d70ed61207
    template<typename T>
    double linearInterpolationApprox(T func, const Triangle &triangle) {
        std::vector<Vec> vecs = triangle.getOwner()->resolvePoints(triangle.getVertices());
        Eigen::Vector3d bar = triangle.barycentric(Eigen::Vector2d(triangle.getCc().x, triangle.getCc().y));
        double area = triangle.area();
        return area * ((bar(0)  * func(vecs[0].getX(), vecs[0].getY())) +
                (bar(1) * func(vecs[1].getX(), vecs[1].getY())) +
                (bar(2) * func(vecs[2].getX(), vecs[2].getY())));
    }
}