MapTriangle.hpp

Go to the documentation of this file.

// _LICENSE_HEADER_
//
// Copyright (C) 2019 - 2025.
// Terms register on the GPL-3.0 license.
//
// This file can be redistributed and/or modified under the license terms.
//
// See top level LICENSE file for more details.
//
// This file can be used citing references in CITATION.cff file.
 
#ifndef __MapTriangle_H
#define __MapTriangle_H
 
#include "Eigen/Eigen"
 
namespace Gedim
{
class MapTriangle
{
  public:
    struct MapTriangleData final
    {
        Eigen::Matrix3d BMatrix;
        Eigen::Matrix3d BMatrixInv;
        Eigen::Vector3d b;
        double DetBMatrix;
        double DetBMatrixInv;
    };
 
  private:
    static inline Eigen::Matrix3d BMatrix(const Eigen::Matrix3d &vertices)
    {
        Eigen::Matrix3d B;
        B.row(0) << vertices(0, 1) - vertices(0, 0), vertices(0, 2) - vertices(0, 0), 0.0;
        B.row(1) << vertices(1, 1) - vertices(1, 0), vertices(1, 2) - vertices(1, 0), 0.0;
        B.row(2) << 0.0, 0.0, 1.0;
        return B;
    }
 
    static inline Eigen::Vector3d b(const Eigen::Matrix3d &vertices)
    {
        return vertices.col(0);
    }
 
  public:
    MapTriangle()
    {
    }
    ~MapTriangle()
    {
    }
 
    MapTriangleData Compute(const Eigen::Matrix3d &vertices) const;
 
    static inline Eigen::MatrixXd F(const MapTriangleData &mapData, const Eigen::MatrixXd &x)
    {
        return (mapData.BMatrix * x).colwise() + mapData.b;
    }
    static inline Eigen::MatrixXd FInv(const MapTriangleData &mapData, const Eigen::MatrixXd &x)
    {
        return mapData.BMatrixInv * (x.colwise() - mapData.b);
    }
    static Eigen::MatrixXd J(const MapTriangleData &mapData, const Eigen::MatrixXd &x);
 
    static inline Eigen::VectorXd DetJ(const MapTriangleData &mapData, const Eigen::MatrixXd &x)
    {
        return Eigen::VectorXd::Constant(x.cols(), mapData.DetBMatrix);
    }
 
    static inline double DetJ(const MapTriangleData &mapData)
    {
        return mapData.DetBMatrix;
    };
};
} // namespace Gedim
 
#endif // __MapTriangle_H