ddfem 1.0.4__py3-none-any.whl → 1.0.6__py3-none-any.whl

ddfem/data/localfunction_a8c726f728cd35d137188b33301aeef1_a524c1196983e65de1c06d7d6afdeb44.cc

@@ -0,0 +1,176 @@
+#ifndef GUARD_a8c726f728cd35d137188b33301aeef1
+#define GUARD_a8c726f728cd35d137188b33301aeef1
+
+#define USING_DUNE_PYTHON 1
+#include <config.h>
+#include <dune/alugrid/dgf.hh>
+#include <dune/alugrid/grid.hh>
+#include <dune/fem/gridpart/adaptiveleafgridpart.hh>
+#include <dune/fem/gridpart/filter/simple.hh>
+#include <dune/fem/gridpart/filteredgridpart.hh>
+#include <dune/fempy/py/gridview.hh>
+#include <dune/python/grid/gridview.hh>
+#include <dune/python/grid/hierarchical.hh>
+#include <dune/fem/function/localfunction/bindable.hh>
+#include <dune/fem/common/intersectionside.hh>
+#include <dune/python/pybind11/pybind11.h>
+#include <dune/python/pybind11/extensions.h>
+#include <dune/fempy/py/grid/gridpart.hh>
+#include <dune/common/exceptions.hh>
+#include <dune/fempy/function/virtualizedgridfunction.hh>
+#include <dune/fempy/py/ufllocalfunction.hh>
+
+namespace UFLLocalFunctions_a8c726f728cd35d137188b33301aeef1
+{
+
+  // UFLLocalFunction
+// ----------------
+
+template< class GridPart, class Coeffddm >
+struct UFLLocalFunction
+  : public Dune::Fem::BindableGridFunctionWithSpace<GridPart,Dune::Dim<1>>
+{
+  typedef GridPart GridPartType;
+  typedef typename GridPartType::GridViewType GridView;
+  typedef typename GridView::ctype ctype;
+  typedef Dune::Fem::BindableGridFunctionWithSpace<GridPart,Dune::Dim<1>> BaseType;
+  typedef Dune::Fem::GridFunctionSpace<GridPartType,Dune::Dim<1>> FunctionSpaceType;
+  typedef typename GridPartType::template Codim< 0 >::EntityType EntityType;
+  typedef typename GridPartType::IntersectionType IntersectionType;
+  typedef typename EntityType::Geometry Geometry;
+  typedef typename Geometry::GlobalCoordinate GlobalCoordinateType;
+  typedef Dune::Fem::IntersectionSide Side;
+  typedef std::tuple<> ConstantTupleType;
+  typedef std::tuple< Dune::Fem::GridFunctionSpace< GridPartType, Dune::FieldVector< double, 1 > > > CoefficientFunctionSpaceTupleType;
+  typedef std::tuple< Coeffddm > CoefficientTupleType;
+  template< std::size_t i >
+  using CoefficientFunctionSpaceType = std::tuple_element_t< i, CoefficientFunctionSpaceTupleType >;
+  template< std::size_t i >
+  using CoefficientRangeType = typename CoefficientFunctionSpaceType< i >::RangeType;
+  template< std::size_t i >
+  using CoefficientJacobianRangeType = typename CoefficientFunctionSpaceType< i >::JacobianRangeType;
+  static constexpr bool gridPartValid = Dune::Fem::checkGridPartValid<GridPartType,Dune::Fem::ConstLocalFunction<Coeffddm>>();
+  template< std::size_t i >
+  using CoefficientType = std::tuple_element_t< i, CoefficientTupleType >;
+  template< std::size_t i >
+  using ConstantType = typename std::tuple_element_t< i, ConstantTupleType >::element_type;
+  using BaseType::entity;
+  using BaseType::geometry;
+
+  UFLLocalFunction ( const GridPartType &gridPart, const std::string &name, int order, const Coeffddm &coeffddm, const Dune::Fem::ParameterReader &parameter = Dune::Fem::Parameter::container() )
+    : BaseType(gridPart,name,order),
+      coefficients_( Dune::Fem::ConstLocalFunction< Coeffddm >( coeffddm ) )
+  {}
+
+  void bind ( const IntersectionType &intersection, Side side )
+  {
+    BaseType::bind(intersection, side);
+    std::get< 0 >( coefficients_ ).bind( this->entity() );
+  }
+
+  void bind ( const EntityType &entity )
+  {
+    BaseType::bind(entity);
+    std::get< 0 >( coefficients_ ).bind( this->entity() );
+  }
+
+  void unbind ()
+  {
+    BaseType::unbind();
+    std::get< 0 >( coefficients_ ).unbind();
+  }
+
+  template< class Point >
+  void evaluate ( const Point &x, typename FunctionSpaceType::RangeType &result ) const
+  {
+    typename CoefficientFunctionSpaceType< 0 >::JacobianRangeType tmp0 = jacobianCoefficient< 0 >( x );
+    const auto tmp1 = (tmp0[ 0 ])[ 1 ] * (tmp0[ 0 ])[ 1 ];
+    const auto tmp2 = (tmp0[ 0 ])[ 0 ] * (tmp0[ 0 ])[ 0 ];
+    const auto tmp3 = tmp2 + tmp1;
+    result[ 0 ] = tmp3;
+  }
+
+  template< class Point >
+  void jacobian ( const Point &x, typename FunctionSpaceType::JacobianRangeType &result ) const
+  {
+    typename CoefficientFunctionSpaceType< 0 >::HessianRangeType tmp0 = hessianCoefficient< 0 >( x );
+    typename CoefficientFunctionSpaceType< 0 >::JacobianRangeType tmp1 = jacobianCoefficient< 0 >( x );
+    const auto tmp2 = (tmp1[ 0 ])[ 1 ] * ((tmp0[ 0 ])[ 1 ])[ 0 ];
+    const auto tmp3 = ((tmp0[ 0 ])[ 1 ])[ 0 ] * (tmp1[ 0 ])[ 1 ];
+    const auto tmp4 = tmp3 + tmp2;
+    const auto tmp5 = (tmp1[ 0 ])[ 0 ] * ((tmp0[ 0 ])[ 0 ])[ 0 ];
+    const auto tmp6 = ((tmp0[ 0 ])[ 0 ])[ 0 ] * (tmp1[ 0 ])[ 0 ];
+    const auto tmp7 = tmp6 + tmp5;
+    const auto tmp8 = tmp7 + tmp4;
+    const auto tmp9 = (tmp1[ 0 ])[ 1 ] * ((tmp0[ 0 ])[ 1 ])[ 1 ];
+    const auto tmp10 = ((tmp0[ 0 ])[ 1 ])[ 1 ] * (tmp1[ 0 ])[ 1 ];
+    const auto tmp11 = tmp10 + tmp9;
+    const auto tmp12 = (tmp1[ 0 ])[ 0 ] * ((tmp0[ 0 ])[ 0 ])[ 1 ];
+    const auto tmp13 = ((tmp0[ 0 ])[ 0 ])[ 1 ] * (tmp1[ 0 ])[ 0 ];
+    const auto tmp14 = tmp13 + tmp12;
+    const auto tmp15 = tmp14 + tmp11;
+    (result[ 0 ])[ 0 ] = tmp8;
+    (result[ 0 ])[ 1 ] = tmp15;
+  }
+
+  template< class Point >
+  void hessian ( const Point &x, typename FunctionSpaceType::HessianRangeType &result ) const
+  {
+    DUNE_THROW(Dune::NotImplemented,"hessian method could not be generated for local function (TooHighDerivative('CodeGenerator does not allow for third order derivatives, yet.'))");
+    result=typename FunctionSpaceType::HessianRangeType(0);
+  }
+
+  template< std::size_t i >
+  const ConstantType< i > &constant () const
+  {
+    return *std::get< i >( constants_ );
+  }
+
+  template< std::size_t i >
+  ConstantType< i > &constant ()
+  {
+    return *std::get< i >( constants_ );
+  }
+
+  template< std::size_t i, class Point >
+  typename CoefficientFunctionSpaceType< i >::RangeType evaluateCoefficient ( const Point &x ) const
+  {
+    typename CoefficientFunctionSpaceType< i >::RangeType result;
+    std::get< i >( coefficients_ ).evaluate( x, result );;
+    return result;
+  }
+
+  template< std::size_t i, class Point >
+  typename CoefficientFunctionSpaceType< i >::JacobianRangeType jacobianCoefficient ( const Point &x ) const
+  {
+    typename CoefficientFunctionSpaceType< i >::JacobianRangeType result;
+    std::get< i >( coefficients_ ).jacobian( x, result );;
+    return result;
+  }
+
+  template< std::size_t i, class Point >
+  typename CoefficientFunctionSpaceType< i >::HessianRangeType hessianCoefficient ( const Point &x ) const
+  {
+    typename CoefficientFunctionSpaceType< i >::HessianRangeType result;
+    std::get< i >( coefficients_ ).hessian( x, result );;
+    return result;
+  }
+  ConstantTupleType constants_;
+  std::tuple< Dune::Fem::ConstLocalFunction< Coeffddm > > coefficients_;
+};
+
+} // namespace UFLLocalFunctions_a8c726f728cd35d137188b33301aeef1
+
+PYBIND11_MODULE( localfunction_a8c726f728cd35d137188b33301aeef1_a524c1196983e65de1c06d7d6afdeb44, module )
+{
+  typedef UFLLocalFunctions_a8c726f728cd35d137188b33301aeef1::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > >, Dune::FemPy::VirtualizedGridFunction< Dune::FemPy::GridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>, Dune::FieldVector< double, 1 > > > LocalFunctionType;
+  if constexpr( LocalFunctionType::gridPartValid )
+  {
+      auto cls = Dune::Python::insertClass<LocalFunctionType>(module,"UFLLocalFunction",Dune::Python::GenerateTypeName("UFLLocalFunctions_a8c726f728cd35d137188b33301aeef1::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > >, Dune::FemPy::VirtualizedGridFunction< Dune::FemPy::GridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>, Dune::FieldVector< double, 1 > > >"), Dune::Python::IncludeFiles({"python/dune/generated/localfunction_a8c726f728cd35d137188b33301aeef1_a524c1196983e65de1c06d7d6afdeb44.cc"})).first;
+      Dune::FemPy::registerUFLLocalFunction( module, cls );
+      cls.def( pybind11::init( [] ( pybind11::object gridView, const std::string &name, int order, const Dune::FemPy::VirtualizedGridFunction< Dune::FemPy::GridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>, Dune::FieldVector< double, 1 > > &coeffddm ) {return new LocalFunctionType( Dune::FemPy::gridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>(gridView),name,order, coeffddm); } ), pybind11::keep_alive< 1, 2 >(), pybind11::keep_alive< 1, 3 >() );
+      cls.def_property_readonly( "virtualized", [] ( LocalFunctionType& ) -> bool { return true;});
+  }
+}
+
+#endif

ddfem/data/localfunction_a9bb07439605d560b583608f20441342_a524c1196983e65de1c06d7d6afdeb44.cc

@@ -0,0 +1,173 @@
+#ifndef GUARD_a9bb07439605d560b583608f20441342
+#define GUARD_a9bb07439605d560b583608f20441342
+
+#define USING_DUNE_PYTHON 1
+#include <config.h>
+#include <dune/alugrid/dgf.hh>
+#include <dune/alugrid/grid.hh>
+#include <dune/fem/gridpart/adaptiveleafgridpart.hh>
+#include <dune/fem/gridpart/filter/simple.hh>
+#include <dune/fem/gridpart/filteredgridpart.hh>
+#include <dune/fempy/py/gridview.hh>
+#include <dune/python/grid/gridview.hh>
+#include <dune/python/grid/hierarchical.hh>
+#include <dune/fem/function/localfunction/bindable.hh>
+#include <dune/fem/common/intersectionside.hh>
+#include <dune/python/pybind11/pybind11.h>
+#include <dune/python/pybind11/extensions.h>
+#include <dune/fempy/py/grid/gridpart.hh>
+#include <dune/common/exceptions.hh>
+#include <dune/fempy/py/ufllocalfunction.hh>
+
+namespace UFLLocalFunctions_a9bb07439605d560b583608f20441342
+{
+
+  // UFLLocalFunction
+// ----------------
+
+template< class GridPart >
+struct UFLLocalFunction
+  : public Dune::Fem::BindableGridFunctionWithSpace<GridPart,Dune::Dim<1>>
+{
+  typedef GridPart GridPartType;
+  typedef typename GridPartType::GridViewType GridView;
+  typedef typename GridView::ctype ctype;
+  typedef Dune::Fem::BindableGridFunctionWithSpace<GridPart,Dune::Dim<1>> BaseType;
+  typedef Dune::Fem::GridFunctionSpace<GridPartType,Dune::Dim<1>> FunctionSpaceType;
+  typedef typename GridPartType::template Codim< 0 >::EntityType EntityType;
+  typedef typename GridPartType::IntersectionType IntersectionType;
+  typedef typename EntityType::Geometry Geometry;
+  typedef typename Geometry::GlobalCoordinate GlobalCoordinateType;
+  typedef Dune::Fem::IntersectionSide Side;
+  typedef std::tuple<> ConstantTupleType;
+  typedef std::tuple<> CoefficientTupleType;
+  static constexpr bool gridPartValid = true;
+  template< std::size_t i >
+  using CoefficientType = std::tuple_element_t< i, CoefficientTupleType >;
+  template< std::size_t i >
+  using ConstantType = typename std::tuple_element_t< i, ConstantTupleType >::element_type;
+  using BaseType::entity;
+  using BaseType::geometry;
+
+  UFLLocalFunction ( const GridPartType &gridPart, const std::string &name, int order, const Dune::Fem::ParameterReader &parameter = Dune::Fem::Parameter::container() )
+    : BaseType(gridPart,name,order)
+  {}
+
+  void bind ( const IntersectionType &intersection, Side side )
+  {
+    BaseType::bind(intersection, side);
+  }
+
+  void bind ( const EntityType &entity )
+  {
+    BaseType::bind(entity);
+  }
+
+  void unbind ()
+  {
+    BaseType::unbind();
+  }
+
+  template< class Point >
+  void evaluate ( const Point &x, typename FunctionSpaceType::RangeType &result ) const
+  {
+    using std::sqrt;
+    GlobalCoordinateType tmp0 = geometry().global( Dune::Fem::coordinate( x ) );
+    const auto tmp1 = -1 + tmp0[ 0 ];
+    const auto tmp2 = tmp1 * tmp1;
+    const auto tmp3 = tmp0[ 1 ] * tmp0[ 1 ];
+    const auto tmp4 = tmp3 + tmp2;
+    const auto tmp5 = 1e-10 + tmp4;
+    const auto tmp6 = std::sqrt( tmp5 );
+    const auto tmp7 = -0.5 + tmp6;
+    result[ 0 ] = tmp7;
+  }
+
+  template< class Point >
+  void jacobian ( const Point &x, typename FunctionSpaceType::JacobianRangeType &result ) const
+  {
+    using std::sqrt;
+    GlobalCoordinateType tmp0 = geometry().global( Dune::Fem::coordinate( x ) );
+    const auto tmp1 = -1 + tmp0[ 0 ];
+    const auto tmp2 = tmp1 * tmp1;
+    const auto tmp3 = tmp0[ 1 ] * tmp0[ 1 ];
+    const auto tmp4 = tmp3 + tmp2;
+    const auto tmp5 = 1e-10 + tmp4;
+    const auto tmp6 = std::sqrt( tmp5 );
+    const auto tmp7 = 2 * tmp6;
+    const auto tmp8 = tmp1 + tmp1;
+    const auto tmp9 = tmp8 / tmp7;
+    const auto tmp10 = tmp0[ 1 ] + tmp0[ 1 ];
+    const auto tmp11 = tmp10 / tmp7;
+    (result[ 0 ])[ 0 ] = tmp9;
+    (result[ 0 ])[ 1 ] = tmp11;
+  }
+
+  template< class Point >
+  void hessian ( const Point &x, typename FunctionSpaceType::HessianRangeType &result ) const
+  {
+    using std::sqrt;
+    GlobalCoordinateType tmp0 = geometry().global( Dune::Fem::coordinate( x ) );
+    const auto tmp1 = -1 + tmp0[ 0 ];
+    const auto tmp2 = tmp1 * tmp1;
+    const auto tmp3 = tmp0[ 1 ] * tmp0[ 1 ];
+    const auto tmp4 = tmp3 + tmp2;
+    const auto tmp5 = 1e-10 + tmp4;
+    const auto tmp6 = std::sqrt( tmp5 );
+    const auto tmp7 = 2 * tmp6;
+    const auto tmp8 = tmp1 + tmp1;
+    const auto tmp9 = tmp8 / tmp7;
+    const auto tmp10 = 2 * tmp9;
+    const auto tmp11 = tmp10 * tmp9;
+    const auto tmp12 = -1 * tmp11;
+    const auto tmp13 = 2 + tmp12;
+    const auto tmp14 = tmp13 / tmp7;
+    const auto tmp15 = tmp0[ 1 ] + tmp0[ 1 ];
+    const auto tmp16 = tmp15 / tmp7;
+    const auto tmp17 = 2 * tmp16;
+    const auto tmp18 = tmp17 * tmp9;
+    const auto tmp19 = -1 * tmp18;
+    const auto tmp20 = tmp19 / tmp7;
+    const auto tmp21 = tmp10 * tmp16;
+    const auto tmp22 = -1 * tmp21;
+    const auto tmp23 = tmp22 / tmp7;
+    const auto tmp24 = tmp17 * tmp16;
+    const auto tmp25 = -1 * tmp24;
+    const auto tmp26 = 2 + tmp25;
+    const auto tmp27 = tmp26 / tmp7;
+    ((result[ 0 ])[ 0 ])[ 0 ] = tmp14;
+    ((result[ 0 ])[ 0 ])[ 1 ] = tmp20;
+    ((result[ 0 ])[ 1 ])[ 0 ] = tmp23;
+    ((result[ 0 ])[ 1 ])[ 1 ] = tmp27;
+  }
+
+  template< std::size_t i >
+  const ConstantType< i > &constant () const
+  {
+    return *std::get< i >( constants_ );
+  }
+
+  template< std::size_t i >
+  ConstantType< i > &constant ()
+  {
+    return *std::get< i >( constants_ );
+  }
+  ConstantTupleType constants_;
+  std::tuple<  > coefficients_;
+};
+
+} // namespace UFLLocalFunctions_a9bb07439605d560b583608f20441342
+
+PYBIND11_MODULE( localfunction_a9bb07439605d560b583608f20441342_a524c1196983e65de1c06d7d6afdeb44, module )
+{
+  typedef UFLLocalFunctions_a9bb07439605d560b583608f20441342::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > > > LocalFunctionType;
+  if constexpr( LocalFunctionType::gridPartValid )
+  {
+      auto cls = Dune::Python::insertClass<LocalFunctionType>(module,"UFLLocalFunction",Dune::Python::GenerateTypeName("UFLLocalFunctions_a9bb07439605d560b583608f20441342::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > > >"), Dune::Python::IncludeFiles({"python/dune/generated/localfunction_a9bb07439605d560b583608f20441342_a524c1196983e65de1c06d7d6afdeb44.cc"})).first;
+      Dune::FemPy::registerUFLLocalFunction( module, cls );
+      cls.def( pybind11::init( [] ( pybind11::object gridView, const std::string &name, int order ) {return new LocalFunctionType( Dune::FemPy::gridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>(gridView),name,order); } ), pybind11::keep_alive< 1, 2 >() );
+      cls.def_property_readonly( "virtualized", [] ( LocalFunctionType& ) -> bool { return true;});
+  }
+}
+
+#endif