helperfunctions.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2021-2025 The Ikarus Developers mueller@ibb.uni-stuttgart.de
// SPDX-License-Identifier: LGPL-3.0-or-later
 
#pragma once
 
#include <dune/localfefunctions/eigenDuneTransformations.hh>
 
#include <ikarus/utils/tensorutils.hh>
 
namespace Ikarus::EAS::Impl {
template <typename GEO>
auto transformationMatrixAtCenterWithDetJ(const GEO& geometry) {
  const auto& referenceElement = Dune::ReferenceElements<double, GEO::mydimension>::general(geometry.type());
  const auto quadPos0          = referenceElement.position(0, 0);
  const auto detJ0             = geometry.integrationElement(quadPos0);
 
  return (transformationMatrix(geometry, quadPos0) * detJ0).eval();
}
 
template <typename GEO>
void transformDisplacementGradient(const GEO& geometry, Eigen::Matrix<double, GEO::mydimension, GEO::mydimension>& H,
                                   const Dune::FieldVector<double, GEO::mydimension>& quadPos) {
  const auto& referenceElement = Dune::ReferenceElements<double, GEO::mydimension>::general(geometry.type());
  const auto quadPos0          = referenceElement.position(0, 0);
  const auto detJ0             = geometry.integrationElement(quadPos0);
  const auto detJ              = geometry.integrationElement(quadPos);
  const auto JInv              = Dune::toEigen(geometry.jacobianInverseTransposed(quadPos0)).eval(); // pre-transposed
  H                            = ((detJ0 / detJ) * JInv * H * JInv.transpose()).eval();
}
} // namespace Ikarus::EAS::Impl