Utils.h

Go to the documentation of this file.

 
#pragma once
 
#include <Eigen/Core>
 
#include "MeshLib/Elements/Element.h"
#include "MeshLib/Node.h"
 
namespace ProcessLib
{
namespace LIE
{
template <typename Derived>
Eigen::Vector3d computePhysicalCoordinates(
    MeshLib::Element const& e, Eigen::MatrixBase<Derived> const& shape)
{
    Eigen::Vector3d pt = Eigen::Vector3d::Zero();
    for (unsigned i = 0; i < e.getNumberOfNodes(); i++)
    {
        MeshLib::Node const& node = *e.getNode(i);
        for (unsigned j = 0; j < 3; j++)
        {
            pt[j] += shape[i] * node[j];
        }
    }
    return pt;
}
 
}  // namespace LIE
}  // namespace ProcessLib