ddfem 1.0.9__py3-none-any.whl → 1.0.11__py3-none-any.whl

ddfem/data/2.11.dev20250709/localfunction_dd5ab085a9d94f578c1f3d337ea360a3_a524c1196983e65de1c06d7d6afdeb44.cc

@@ -0,0 +1,161 @@
+#ifndef GUARD_dd5ab085a9d94f578c1f3d337ea360a3
+#define GUARD_dd5ab085a9d94f578c1f3d337ea360a3
+
+#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_dd5ab085a9d94f578c1f3d337ea360a3
+{
+
+  // UFLLocalFunction
+// ----------------
+
+template< class GridPart, class Coeffu_h >
+struct UFLLocalFunction
+  : public Dune::Fem::BindableGridFunctionWithSpace<GridPart,Dune::Dim<2>>
+{
+  typedef GridPart GridPartType;
+  typedef typename GridPartType::GridViewType GridView;
+  typedef typename GridView::ctype ctype;
+  typedef Dune::Fem::BindableGridFunctionWithSpace<GridPart,Dune::Dim<2>> BaseType;
+  typedef Dune::Fem::GridFunctionSpace<GridPartType,Dune::Dim<2>> 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, 2 > > > CoefficientFunctionSpaceTupleType;
+  typedef std::tuple< Coeffu_h > 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<Coeffu_h>>();
+  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 Coeffu_h &coeffu_h, const Dune::Fem::ParameterReader &parameter = Dune::Fem::Parameter::container() )
+    : BaseType(gridPart,name,order),
+      coefficients_( Dune::Fem::ConstLocalFunction< Coeffu_h >( coeffu_h ) )
+  {}
+
+  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 );
+    result[ 0 ] = (tmp0[ 0 ])[ 0 ];
+    result[ 1 ] = (tmp0[ 0 ])[ 1 ];
+  }
+
+  template< class Point >
+  void jacobian ( const Point &x, typename FunctionSpaceType::JacobianRangeType &result ) const
+  {
+    typename CoefficientFunctionSpaceType< 0 >::HessianRangeType tmp0 = hessianCoefficient< 0 >( x );
+    (result[ 0 ])[ 0 ] = ((tmp0[ 0 ])[ 0 ])[ 0 ];
+    (result[ 0 ])[ 1 ] = ((tmp0[ 0 ])[ 0 ])[ 1 ];
+    (result[ 1 ])[ 0 ] = ((tmp0[ 0 ])[ 1 ])[ 0 ];
+    (result[ 1 ])[ 1 ] = ((tmp0[ 0 ])[ 1 ])[ 1 ];
+  }
+
+  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< Coeffu_h > > coefficients_;
+};
+
+} // namespace UFLLocalFunctions_dd5ab085a9d94f578c1f3d337ea360a3
+
+PYBIND11_MODULE( localfunction_dd5ab085a9d94f578c1f3d337ea360a3_a524c1196983e65de1c06d7d6afdeb44, module )
+{
+  typedef UFLLocalFunctions_dd5ab085a9d94f578c1f3d337ea360a3::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, 2 > > > LocalFunctionType;
+  if constexpr( LocalFunctionType::gridPartValid )
+  {
+      auto cls = Dune::Python::insertClass<LocalFunctionType>(module,"UFLLocalFunction",Dune::Python::GenerateTypeName("UFLLocalFunctions_dd5ab085a9d94f578c1f3d337ea360a3::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, 2 > > >"), Dune::Python::IncludeFiles({"python/dune/generated/localfunction_dd5ab085a9d94f578c1f3d337ea360a3_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, 2 > > &coeffu_h ) {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, coeffu_h); } ), pybind11::keep_alive< 1, 2 >(), pybind11::keep_alive< 1, 3 >() );
+      cls.def_property_readonly( "virtualized", [] ( LocalFunctionType& ) -> bool { return true;});
+  }
+}
+
+#endif

ddfem/data/2.11.dev20250709/{localfunction_3e323e81c52891c0ecb656eac273e52f_a524c1196983e65de1c06d7d6afdeb44.cc → localfunction_e2b19c5987f9b6fb0cd5453bfa7f4f9f_a524c1196983e65de1c06d7d6afdeb44.cc}

@@ -1,5 +1,5 @@
-#ifndef GUARD_3e323e81c52891c0ecb656eac273e52f
-#define GUARD_3e323e81c52891c0ecb656eac273e52f
+#ifndef GUARD_e2b19c5987f9b6fb0cd5453bfa7f4f9f
+#define GUARD_e2b19c5987f9b6fb0cd5453bfa7f4f9f
 
 #define USING_DUNE_PYTHON 1
 #include <config.h>
@@ -19,7 +19,7 @@
 #include <dune/common/exceptions.hh>
 #include <dune/fempy/py/ufllocalfunction.hh>
 
-namespace UFLLocalFunctions_3e323e81c52891c0ecb656eac273e52f
+namespace UFLLocalFunctions_e2b19c5987f9b6fb0cd5453bfa7f4f9f
 {
 
   // UFLLocalFunction
@@ -27,22 +27,19 @@ namespace UFLLocalFunctions_3e323e81c52891c0ecb656eac273e52f
 
 template< class GridPart >
 struct UFLLocalFunction
-  : public Dune::Fem::BindableGridFunctionWithSpace<GridPart,Dune::Dim<1>>
+  : public Dune::Fem::BindableGridFunctionWithSpace<GridPart,Dune::Dim<2>>
 {
   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 Dune::Fem::BindableGridFunctionWithSpace<GridPart,Dune::Dim<2>> BaseType;
+  typedef Dune::Fem::GridFunctionSpace<GridPartType,Dune::Dim<2>> 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 double Conc5;
-  typedef std::tuple< std::shared_ptr< Conc5 > > ConstantTupleType;
-  template< std::size_t i >
-  using ConstantsRangeType = typename std::tuple_element_t< i, ConstantTupleType >::element_type;
+  typedef std::tuple<> ConstantTupleType;
   typedef std::tuple<> CoefficientTupleType;
   static constexpr bool gridPartValid = true;
   template< std::size_t i >
@@ -54,9 +51,7 @@ struct UFLLocalFunction
 
   UFLLocalFunction ( const GridPartType &gridPart, const std::string &name, int order, const Dune::Fem::ParameterReader &parameter = Dune::Fem::Parameter::container() )
     : BaseType(gridPart,name,order)
-  {
-    std::get< 0 >( constants_ ) = std::make_shared< Conc5 >( (Conc5(0)) );
-  }
+  {}
 
   void bind ( const IntersectionType &intersection, Side side )
   {
@@ -76,8 +71,8 @@ struct UFLLocalFunction
   template< class Point >
   void evaluate ( const Point &x, typename FunctionSpaceType::RangeType &result ) const
   {
-    double tmp0 = constant< 0 >();
-    result[ 0 ] = tmp0;
+    result[ 0 ] = 0.0;
+    result[ 1 ] = 0.0;
   }
 
   template< class Point >
@@ -103,31 +98,20 @@ struct UFLLocalFunction
   {
     return *std::get< i >( constants_ );
   }
-
-  const Conc5 &conc5 () const
-  {
-    return *std::get< 0 >( constants_ );
-  }
-
-  Conc5 &conc5 ()
-  {
-    return *std::get< 0 >( constants_ );
-  }
   ConstantTupleType constants_;
   std::tuple<  > coefficients_;
 };
 
-} // namespace UFLLocalFunctions_3e323e81c52891c0ecb656eac273e52f
+} // namespace UFLLocalFunctions_e2b19c5987f9b6fb0cd5453bfa7f4f9f
 
-PYBIND11_MODULE( localfunction_3e323e81c52891c0ecb656eac273e52f_a524c1196983e65de1c06d7d6afdeb44, module )
+PYBIND11_MODULE( localfunction_e2b19c5987f9b6fb0cd5453bfa7f4f9f_a524c1196983e65de1c06d7d6afdeb44, module )
 {
-  typedef UFLLocalFunctions_3e323e81c52891c0ecb656eac273e52f::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;
+  typedef UFLLocalFunctions_e2b19c5987f9b6fb0cd5453bfa7f4f9f::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_3e323e81c52891c0ecb656eac273e52f::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_3e323e81c52891c0ecb656eac273e52f_a524c1196983e65de1c06d7d6afdeb44.cc"})).first;
+      auto cls = Dune::Python::insertClass<LocalFunctionType>(module,"UFLLocalFunction",Dune::Python::GenerateTypeName("UFLLocalFunctions_e2b19c5987f9b6fb0cd5453bfa7f4f9f::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_e2b19c5987f9b6fb0cd5453bfa7f4f9f_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( "c5", [] ( LocalFunctionType &self ) -> UFLLocalFunctions_3e323e81c52891c0ecb656eac273e52f::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 > > >::Conc5 { return self.conc5(); }, [] ( LocalFunctionType &self, const UFLLocalFunctions_3e323e81c52891c0ecb656eac273e52f::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 > > >::Conc5 &v ) { self.conc5() = v; } );
       cls.def_property_readonly( "virtualized", [] ( LocalFunctionType& ) -> bool { return true;});
   }
 }

ddfem/data/2.11.dev20250709/{localfunction_e69118a51208cfab8a028711e966e76c_19659fe2f2f74bcfb5d532f1a0d38be9.cc → localfunction_e85b93fe9b8fcd40ee82782fc600b65b_a524c1196983e65de1c06d7d6afdeb44.cc}

@@ -1,10 +1,15 @@
-#ifndef GUARD_e69118a51208cfab8a028711e966e76c
-#define GUARD_e69118a51208cfab8a028711e966e76c
+#ifndef GUARD_e85b93fe9b8fcd40ee82782fc600b65b
+#define GUARD_e85b93fe9b8fcd40ee82782fc600b65b
 
 #define USING_DUNE_PYTHON 1
 #include <config.h>
-#include <dune/grid/io/file/dgfparser/dgfyasp.hh>
-#include <dune/grid/yaspgrid.hh>
+#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>
@@ -14,7 +19,7 @@
 #include <dune/common/exceptions.hh>
 #include <dune/fempy/py/ufllocalfunction.hh>
 
-namespace UFLLocalFunctions_e69118a51208cfab8a028711e966e76c
+namespace UFLLocalFunctions_e85b93fe9b8fcd40ee82782fc600b65b
 {
 
   // UFLLocalFunction
@@ -66,23 +71,23 @@ struct UFLLocalFunction
   template< class Point >
   void evaluate ( const Point &x, typename FunctionSpaceType::RangeType &result ) const
   {
-    using std::min;
+    using std::max;
     using std::sqrt;
     GlobalCoordinateType tmp0 = geometry().global( Dune::Fem::coordinate( x ) );
-    const auto tmp1 = 0.2 + tmp0[ 0 ];
+    const auto tmp1 = -0.8 + tmp0[ 1 ];
     const auto tmp2 = tmp1 * tmp1;
-    const auto tmp3 = tmp0[ 1 ] * tmp0[ 1 ];
+    const auto tmp3 = tmp0[ 0 ] * tmp0[ 0 ];
     const auto tmp4 = tmp3 + tmp2;
     const auto tmp5 = 1e-10 + tmp4;
     const auto tmp6 = std::sqrt( tmp5 );
     const auto tmp7 = -0.5 + tmp6;
-    const auto tmp8 = -0.2 + tmp0[ 0 ];
-    const auto tmp9 = tmp8 * tmp8;
+    const auto tmp8 = -1 * tmp7;
+    const auto tmp9 = tmp0[ 1 ] * tmp0[ 1 ];
     const auto tmp10 = tmp3 + tmp9;
     const auto tmp11 = 1e-10 + tmp10;
     const auto tmp12 = std::sqrt( tmp11 );
-    const auto tmp13 = -0.5 + tmp12;
-    const auto tmp14 = std::min( tmp13, tmp7 );
+    const auto tmp13 = -1 + tmp12;
+    const auto tmp14 = std::max( tmp13, tmp8 );
     result[ 0 ] = tmp14;
   }
 
@@ -91,38 +96,40 @@ struct UFLLocalFunction
   {
     using std::sqrt;
     GlobalCoordinateType tmp0 = geometry().global( Dune::Fem::coordinate( x ) );
-    const auto tmp1 = 0.2 + tmp0[ 0 ];
+    const auto tmp1 = -0.8 + tmp0[ 1 ];
     const auto tmp2 = tmp1 * tmp1;
-    const auto tmp3 = tmp0[ 1 ] * tmp0[ 1 ];
+    const auto tmp3 = tmp0[ 0 ] * tmp0[ 0 ];
     const auto tmp4 = tmp3 + tmp2;
     const auto tmp5 = 1e-10 + tmp4;
     const auto tmp6 = std::sqrt( tmp5 );
     const auto tmp7 = -0.5 + tmp6;
-    const auto tmp8 = -0.2 + tmp0[ 0 ];
-    const auto tmp9 = tmp8 * tmp8;
+    const auto tmp8 = -1 * tmp7;
+    const auto tmp9 = tmp0[ 1 ] * tmp0[ 1 ];
     const auto tmp10 = tmp3 + tmp9;
     const auto tmp11 = 1e-10 + tmp10;
     const auto tmp12 = std::sqrt( tmp11 );
-    const auto tmp13 = -0.5 + tmp12;
+    const auto tmp13 = -1 + tmp12;
     const auto tmp14 = 2 * tmp12;
-    const auto tmp15 = tmp8 + tmp8;
+    const auto tmp15 = tmp0[ 0 ] + tmp0[ 0 ];
     const auto tmp16 = tmp15 / tmp14;
-    const auto tmp17 = tmp16 * (tmp13 < tmp7 ? 1 : 0.0);
+    const auto tmp17 = tmp16 * (tmp13 > tmp8 ? 1 : 0.0);
     const auto tmp18 = 2 * tmp6;
-    const auto tmp19 = tmp1 + tmp1;
-    const auto tmp20 = tmp19 / tmp18;
-    const auto tmp21 = -1 * (tmp13 < tmp7 ? 1 : 0.0);
+    const auto tmp19 = tmp15 / tmp18;
+    const auto tmp20 = -1 * tmp19;
+    const auto tmp21 = -1 * (tmp13 > tmp8 ? 1 : 0.0);
     const auto tmp22 = 1.0 + tmp21;
     const auto tmp23 = tmp22 * tmp20;
     const auto tmp24 = tmp23 + tmp17;
     const auto tmp25 = tmp0[ 1 ] + tmp0[ 1 ];
     const auto tmp26 = tmp25 / tmp14;
-    const auto tmp27 = tmp26 * (tmp13 < tmp7 ? 1 : 0.0);
-    const auto tmp28 = tmp25 / tmp18;
-    const auto tmp29 = tmp22 * tmp28;
-    const auto tmp30 = tmp29 + tmp27;
+    const auto tmp27 = tmp26 * (tmp13 > tmp8 ? 1 : 0.0);
+    const auto tmp28 = tmp1 + tmp1;
+    const auto tmp29 = tmp28 / tmp18;
+    const auto tmp30 = -1 * tmp29;
+    const auto tmp31 = tmp22 * tmp30;
+    const auto tmp32 = tmp31 + tmp27;
     (result[ 0 ])[ 0 ] = tmp24;
-    (result[ 0 ])[ 1 ] = tmp30;
+    (result[ 0 ])[ 1 ] = tmp32;
   }
 
   template< class Point >
@@ -130,37 +137,37 @@ struct UFLLocalFunction
   {
     using std::sqrt;
     GlobalCoordinateType tmp0 = geometry().global( Dune::Fem::coordinate( x ) );
-    const auto tmp1 = 0.2 + tmp0[ 0 ];
+    const auto tmp1 = -0.8 + tmp0[ 1 ];
     const auto tmp2 = tmp1 * tmp1;
-    const auto tmp3 = tmp0[ 1 ] * tmp0[ 1 ];
+    const auto tmp3 = tmp0[ 0 ] * tmp0[ 0 ];
     const auto tmp4 = tmp3 + tmp2;
     const auto tmp5 = 1e-10 + tmp4;
     const auto tmp6 = std::sqrt( tmp5 );
     const auto tmp7 = -0.5 + tmp6;
-    const auto tmp8 = -0.2 + tmp0[ 0 ];
-    const auto tmp9 = tmp8 * tmp8;
+    const auto tmp8 = -1 * tmp7;
+    const auto tmp9 = tmp0[ 1 ] * tmp0[ 1 ];
     const auto tmp10 = tmp3 + tmp9;
     const auto tmp11 = 1e-10 + tmp10;
     const auto tmp12 = std::sqrt( tmp11 );
-    const auto tmp13 = -0.5 + tmp12;
+    const auto tmp13 = -1 + tmp12;
     const auto tmp14 = 2 * tmp12;
-    const auto tmp15 = tmp8 + tmp8;
+    const auto tmp15 = tmp0[ 0 ] + tmp0[ 0 ];
     const auto tmp16 = tmp15 / tmp14;
     const auto tmp17 = 2 * tmp16;
     const auto tmp18 = tmp17 * tmp16;
     const auto tmp19 = -1 * tmp18;
     const auto tmp20 = 2 + tmp19;
     const auto tmp21 = tmp20 / tmp14;
-    const auto tmp22 = tmp21 * (tmp13 < tmp7 ? 1 : 0.0);
+    const auto tmp22 = tmp21 * (tmp13 > tmp8 ? 1 : 0.0);
     const auto tmp23 = 2 * tmp6;
-    const auto tmp24 = tmp1 + tmp1;
-    const auto tmp25 = tmp24 / tmp23;
-    const auto tmp26 = 2 * tmp25;
-    const auto tmp27 = tmp26 * tmp25;
-    const auto tmp28 = -1 * tmp27;
-    const auto tmp29 = 2 + tmp28;
-    const auto tmp30 = tmp29 / tmp23;
-    const auto tmp31 = -1 * (tmp13 < tmp7 ? 1 : 0.0);
+    const auto tmp24 = tmp15 / tmp23;
+    const auto tmp25 = 2 * tmp24;
+    const auto tmp26 = tmp25 * tmp24;
+    const auto tmp27 = -1 * tmp26;
+    const auto tmp28 = 2 + tmp27;
+    const auto tmp29 = tmp28 / tmp23;
+    const auto tmp30 = -1 * tmp29;
+    const auto tmp31 = -1 * (tmp13 > tmp8 ? 1 : 0.0);
     const auto tmp32 = 1.0 + tmp31;
     const auto tmp33 = tmp32 * tmp30;
     const auto tmp34 = tmp33 + tmp22;
@@ -170,38 +177,42 @@ struct UFLLocalFunction
     const auto tmp38 = tmp37 * tmp16;
     const auto tmp39 = -1 * tmp38;
     const auto tmp40 = tmp39 / tmp14;
-    const auto tmp41 = tmp40 * (tmp13 < tmp7 ? 1 : 0.0);
-    const auto tmp42 = tmp35 / tmp23;
-    const auto tmp43 = 2 * tmp42;
-    const auto tmp44 = tmp43 * tmp25;
-    const auto tmp45 = -1 * tmp44;
-    const auto tmp46 = tmp45 / tmp23;
-    const auto tmp47 = tmp32 * tmp46;
-    const auto tmp48 = tmp47 + tmp41;
-    const auto tmp49 = tmp17 * tmp36;
-    const auto tmp50 = -1 * tmp49;
-    const auto tmp51 = tmp50 / tmp14;
-    const auto tmp52 = tmp51 * (tmp13 < tmp7 ? 1 : 0.0);
-    const auto tmp53 = tmp26 * tmp42;
-    const auto tmp54 = -1 * tmp53;
-    const auto tmp55 = tmp54 / tmp23;
-    const auto tmp56 = tmp32 * tmp55;
-    const auto tmp57 = tmp56 + tmp52;
-    const auto tmp58 = tmp37 * tmp36;
-    const auto tmp59 = -1 * tmp58;
-    const auto tmp60 = 2 + tmp59;
-    const auto tmp61 = tmp60 / tmp14;
-    const auto tmp62 = tmp61 * (tmp13 < tmp7 ? 1 : 0.0);
-    const auto tmp63 = tmp43 * tmp42;
-    const auto tmp64 = -1 * tmp63;
-    const auto tmp65 = 2 + tmp64;
-    const auto tmp66 = tmp65 / tmp23;
-    const auto tmp67 = tmp32 * tmp66;
-    const auto tmp68 = tmp67 + tmp62;
+    const auto tmp41 = tmp40 * (tmp13 > tmp8 ? 1 : 0.0);
+    const auto tmp42 = tmp1 + tmp1;
+    const auto tmp43 = tmp42 / tmp23;
+    const auto tmp44 = 2 * tmp43;
+    const auto tmp45 = tmp44 * tmp24;
+    const auto tmp46 = -1 * tmp45;
+    const auto tmp47 = tmp46 / tmp23;
+    const auto tmp48 = -1 * tmp47;
+    const auto tmp49 = tmp32 * tmp48;
+    const auto tmp50 = tmp49 + tmp41;
+    const auto tmp51 = tmp17 * tmp36;
+    const auto tmp52 = -1 * tmp51;
+    const auto tmp53 = tmp52 / tmp14;
+    const auto tmp54 = tmp53 * (tmp13 > tmp8 ? 1 : 0.0);
+    const auto tmp55 = tmp25 * tmp43;
+    const auto tmp56 = -1 * tmp55;
+    const auto tmp57 = tmp56 / tmp23;
+    const auto tmp58 = -1 * tmp57;
+    const auto tmp59 = tmp32 * tmp58;
+    const auto tmp60 = tmp59 + tmp54;
+    const auto tmp61 = tmp37 * tmp36;
+    const auto tmp62 = -1 * tmp61;
+    const auto tmp63 = 2 + tmp62;
+    const auto tmp64 = tmp63 / tmp14;
+    const auto tmp65 = tmp64 * (tmp13 > tmp8 ? 1 : 0.0);
+    const auto tmp66 = tmp44 * tmp43;
+    const auto tmp67 = -1 * tmp66;
+    const auto tmp68 = 2 + tmp67;
+    const auto tmp69 = tmp68 / tmp23;
+    const auto tmp70 = -1 * tmp69;
+    const auto tmp71 = tmp32 * tmp70;
+    const auto tmp72 = tmp71 + tmp65;
     ((result[ 0 ])[ 0 ])[ 0 ] = tmp34;
-    ((result[ 0 ])[ 0 ])[ 1 ] = tmp48;
-    ((result[ 0 ])[ 1 ])[ 0 ] = tmp57;
-    ((result[ 0 ])[ 1 ])[ 1 ] = tmp68;
+    ((result[ 0 ])[ 0 ])[ 1 ] = tmp50;
+    ((result[ 0 ])[ 1 ])[ 0 ] = tmp60;
+    ((result[ 0 ])[ 1 ])[ 1 ] = tmp72;
   }
 
   template< std::size_t i >
@@ -219,16 +230,16 @@ struct UFLLocalFunction
   std::tuple<  > coefficients_;
 };
 
-} // namespace UFLLocalFunctions_e69118a51208cfab8a028711e966e76c
+} // namespace UFLLocalFunctions_e85b93fe9b8fcd40ee82782fc600b65b
 
-PYBIND11_MODULE( localfunction_e69118a51208cfab8a028711e966e76c_19659fe2f2f74bcfb5d532f1a0d38be9, module )
+PYBIND11_MODULE( localfunction_e85b93fe9b8fcd40ee82782fc600b65b_a524c1196983e65de1c06d7d6afdeb44, module )
 {
-  typedef UFLLocalFunctions_e69118a51208cfab8a028711e966e76c::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::YaspGrid< 2, Dune::EquidistantOffsetCoordinates< double, 2 > >::LeafGridView > > LocalFunctionType;
+  typedef UFLLocalFunctions_e85b93fe9b8fcd40ee82782fc600b65b::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_e69118a51208cfab8a028711e966e76c::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::YaspGrid< 2, Dune::EquidistantOffsetCoordinates< double, 2 > >::LeafGridView > >"), Dune::Python::IncludeFiles({"python/dune/generated/localfunction_e69118a51208cfab8a028711e966e76c_19659fe2f2f74bcfb5d532f1a0d38be9.cc"})).first;
+      auto cls = Dune::Python::insertClass<LocalFunctionType>(module,"UFLLocalFunction",Dune::Python::GenerateTypeName("UFLLocalFunctions_e85b93fe9b8fcd40ee82782fc600b65b::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_e85b93fe9b8fcd40ee82782fc600b65b_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::YaspGrid< 2, Dune::EquidistantOffsetCoordinates< double, 2 > >::LeafGridView>(gridView),name,order); } ), pybind11::keep_alive< 1, 2 >() );
+      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;});
   }
 }

ddfem/data/extra.modules

@@ -0,0 +1,13 @@
+femspace_11b3436cb033df683615c50d5d4239bd_f14dd0b3b474c0cd816334330f55ad72.cc
+localfunction_3d957c5b2628e5120636b8e3f461bc0d_a524c1196983e65de1c06d7d6afdeb44.cc
+localfunction_6363af791dde11d699a9a3d876cc5d9a_a524c1196983e65de1c06d7d6afdeb44.cc
+localfunction_762aa49bc5a47ad3d8d8bbde36b42ee5_a524c1196983e65de1c06d7d6afdeb44.cc
+localfunction_e2b19c5987f9b6fb0cd5453bfa7f4f9f_a524c1196983e65de1c06d7d6afdeb44.cc
+femscheme_765ef00d995a939c34bdc0d0cc4847ca.cc
+localfunction_22e797a8f3e8c322f33c41e07bd15165_a524c1196983e65de1c06d7d6afdeb44.cc
+localfunction_32613626342183eafa360dc1f0f18924_a524c1196983e65de1c06d7d6afdeb44.cc
+localfunction_4976bbef3267846b59f3a26497ffe9db_a524c1196983e65de1c06d7d6afdeb44.cc
+localfunction_5591a2c803b9658084d8f11d44bf3ae1_a524c1196983e65de1c06d7d6afdeb44.cc
+localfunction_6363af791dde11d699a9a3d876cc5d9a_af122c1df944c95cd395ec0f91d0f970.cc
+localfunction_d4b9ca85d48735ee01f395a42eaa0cee_af122c1df944c95cd395ec0f91d0f970.cc
+localfunction_dd1fb4bd60f2d8b3a062dbbdaedb9fee_a524c1196983e65de1c06d7d6afdeb44.cc