MapTetrahedron.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 __MapTetrahedron_H
#define __MapTetrahedron_H
 
#include "Eigen/Eigen"
#include "GeometryUtilities.hpp"
 
namespace Gedim
{
class MapTetrahedron
{
  public:
    struct MapTetrahedronData final
    {
        Eigen::Matrix3d Q;
        Eigen::Vector3d b;
        Eigen::Matrix3d QInv;
        double DetQ;
        double DetQInv;
    };
 
  private:
    const Gedim::GeometryUtilities &geometryUtilities;
 
    bool TestMapConfiguration(const Eigen::MatrixXd &vertices,
                              const Eigen::MatrixXd &referencePoints,
                              const unsigned int &secondVertexIndex,
                              const unsigned int &thirdVertexIndex,
                              const unsigned int &fourthVertexIndex,
                              MapTetrahedron::MapTetrahedronData &result) const;
 
  public:
    MapTetrahedron(const Gedim::GeometryUtilities &geometryUtilities);
 
    ~MapTetrahedron()
    {
    }
 
    static inline Eigen::Matrix3d Q(const Eigen::Vector3d &firstVertex,
                                    const Eigen::Vector3d &secondVertex,
                                    const Eigen::Vector3d &thirdVertex,
                                    const Eigen::Vector3d &fourthVertex)
    {
        Eigen::Matrix3d Q;
        Q.row(0) << secondVertex.x() - firstVertex.x(), thirdVertex.x() - firstVertex.x(),
            fourthVertex.x() - firstVertex.x();
        Q.row(1) << secondVertex.y() - firstVertex.y(), thirdVertex.y() - firstVertex.y(),
            fourthVertex.y() - firstVertex.y();
        Q.row(2) << secondVertex.z() - firstVertex.z(), thirdVertex.z() - firstVertex.z(),
            fourthVertex.z() - firstVertex.z();
 
        return Q;
    }
 
    static inline Eigen::Vector3d b(const Eigen::Vector3d &firstVertex)
    {
        return firstVertex;
    }
 
    MapTetrahedronData Compute(const Eigen::MatrixXd &vertices) const;
 
    static inline Eigen::MatrixXd F(const MapTetrahedronData &mapData, const Eigen::MatrixXd &x)
    {
        return (mapData.Q * x).colwise() + mapData.b;
    }
 
    static inline Eigen::MatrixXd FInv(const MapTetrahedronData &mapData, const Eigen::MatrixXd &x)
    {
        return mapData.QInv * (x.colwise() - mapData.b);
    }
 
    static Eigen::MatrixXd J(const MapTetrahedronData &mapData, const Eigen::MatrixXd &x);
    static Eigen::VectorXd DetJ(const MapTetrahedronData &mapData, const Eigen::MatrixXd &x);
 
    static inline double DetJ(const MapTetrahedronData &mapData)
    {
        return mapData.DetQ;
    };
};
} // namespace Gedim
 
#endif // __MapTetrahedron_H