Revision 2a26041b8cfe3af78c295ee4d7f0d4ef17c7a807 authored by Dmitry Yu. Naumov on 16 August 2021, 06:32:22 UTC, committed by Dmitry Yu. Naumov on 16 August 2021, 06:32:22 UTC
Staggered: check all processes for convergence in coupling iterations See merge request ogs/ogs!3740
ElementCoordinatesMappingLocal.h
/**
* \file
* \copyright
* Copyright (c) 2012-2021, 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
Computing file changes ...