https://gitlab.opengeosys.org/ogs/ogs.git
Tip revision: 3bc70067ad8267c465bb475948f4d8b4f03e7bb6 authored by Wenqing Wang on 20 November 2020, 10:38:43 UTC
[TRM] Update for the changes in RM for using drho_l/dp
[TRM] Update for the changes in RM for using drho_l/dp
Tip revision: 3bc7006
ElementCoordinatesMappingLocal.h
/**
* \file
* \copyright
* Copyright (c) 2012-2020, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*/
#pragma once
#include <vector>
#include <Eigen/Eigen>
#include "MathLib/Point3d.h"
namespace MeshLib
{
class Element;
}
namespace MeshLib
{
using RotationMatrix = Eigen::Matrix<double, 3u, 3u, Eigen::RowMajor>;
/**
* This class maps node coordinates on intrinsic coordinates of the given element.
*/
class ElementCoordinatesMappingLocal final
{
public:
/**
* Constructor
* \param e Mesh element whose node coordinates are mapped
* \param global_dim Global dimension
*/
ElementCoordinatesMappingLocal(const Element &e, const unsigned global_dim);
/// return the global dimension
unsigned getGlobalDimension() const { return _global_dim; }
/// return mapped coordinates of the node
MathLib::Point3d const& getMappedCoordinates(std::size_t node_id) const
{
return _points[node_id];
}
/// return a rotation matrix converting to global coordinates
const RotationMatrix& getRotationMatrixToGlobal() const {return _matR2global;}
private:
const unsigned _global_dim;
std::vector<MathLib::Point3d> _points;
RotationMatrix _matR2global;
};
} // namespace MeshLib