registerpreelement.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/common/fvector.hh>
#include <dune/python/pybind11/pybind11.h>
 
#include <Eigen/Core>
 
namespace Ikarus::Python {
 
template <class NonLinearElasticPre, class... options>
void registerNonLinearElasticPre(pybind11::handle scope, pybind11::class_<NonLinearElasticPre, options...> cls) {
  using Material = typename NonLinearElasticPre::Material;
  cls.def(pybind11::init([](const Material& mat) { return new NonLinearElasticPre(mat); }));
}
 
template <class LinearElasticPre, class... options>
void registerLinearElasticPre(pybind11::handle scope, pybind11::class_<LinearElasticPre, options...> cls) {
  using Material = typename LinearElasticPre::Material;
  cls.def(pybind11::init([](const Material& mat) { return new LinearElasticPre(mat); }));
}
 
template <class TrussPre, class... options>
void registerTrussPre(pybind11::handle scope, pybind11::class_<TrussPre, options...> cls) {
  cls.def(pybind11::init([](double emod, double area) { return new TrussPre({emod, area}); }));
}
 
template <class KirchhoffLoveShellPre, class... options>
void registerKirchhoffLoveShellPre(pybind11::handle scope, pybind11::class_<KirchhoffLoveShellPre, options...> cls) {
  cls.def(pybind11::init(
      [](const double& E, const double& nu, const double& h) { return new KirchhoffLoveShellPre({E, nu}, h); }));
}
 
template <class EASPre, class... options>
void registerEnhancedAssumedStrainsPre(pybind11::handle scope, pybind11::class_<EASPre, options...> cls) {
  cls.def(pybind11::init([](int numberOfParameter) { return new EASPre(numberOfParameter); }));
}
 
template <class NeumannBoundaryLoadPre, class... options>
void registerNeumannBoundaryLoadPre(pybind11::handle scope, pybind11::class_<NeumannBoundaryLoadPre, options...> cls) {
  using BoundaryPatchType = typename NeumannBoundaryLoadPre::BoundaryPatchType;
  using GridView          = typename NeumannBoundaryLoadPre::GridView;
 
  using LoadFunction = std::function<Eigen::Vector<double, NeumannBoundaryLoadPre::worldDim>(
      Dune::FieldVector<double, NeumannBoundaryLoadPre::worldDim>, const double&)>;
  cls.def(pybind11::init([](const BoundaryPatchType& patch, LoadFunction volumeLoad) {
            return new NeumannBoundaryLoadPre(&patch, volumeLoad);
          }),
          pybind11::keep_alive<1, 2>());
}
 
template <class VolumeLoadPre, class... options>
void registerVolumeLoadPre(pybind11::handle scope, pybind11::class_<VolumeLoadPre, options...> cls) {
  using LoadFunction = std::function<Eigen::Vector<double, VolumeLoadPre::worldDim>(
      Dune::FieldVector<double, VolumeLoadPre::worldDim>, const double&)>;
  cls.def(pybind11::init([](LoadFunction volumeLoad) { return new VolumeLoadPre(volumeLoad); }));
}
 
} // namespace Ikarus::Python