ddfem 1.0.3__tar.gz → 1.0.6__tar.gz

{ddfem-1.0.3 → ddfem-1.0.6}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ddfem
-Version: 1.0.3
+Version: 1.0.6
 Summary: Diffuse domain finite element solver
 Author-email: Luke Benfield <luke.benfield@warwick.ac.uk>, Andreas Dedner <a.s.dedner@warwick.ac.uk>
 License-Expression: MIT
@@ -15,8 +15,17 @@ Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: fenics-ufl==2022.2.0
 Requires-Dist: numpy
+Requires-Dist: setuptools==79
 Provides-Extra: dune
 Requires-Dist: dune-fem; extra == "dune"
+Requires-Dist: pygmsh; extra == "dune"
+Provides-Extra: doc
+Requires-Dist: dune-fem; extra == "doc"
+Requires-Dist: pygmsh; extra == "doc"
+Requires-Dist: jupyter; extra == "doc"
+Requires-Dist: jupytext; extra == "doc"
+Requires-Dist: nbsphinx; extra == "doc"
+Requires-Dist: sphinx_rtd_theme; extra == "doc"
 Dynamic: license-file
 
 # Diffuse Domain Finite Element Methods

ddfem-1.0.6/ddfem/__init__.py

@@ -0,0 +1,77 @@
+from . import geometry
+from . import transformers
+from . import boundary
+from .model2ufl import model2ufl
+from .plot import plotSolution
+
+# a utility method to construct a suitable grid
+def GridView(domain, omega=None, *,
+                filterTolerance=0,
+                factor=0, refinements=0,
+                adaptLevels=0):
+    import dune.ufl
+    from dune.fem.function import gridFunction
+    from dune.fem import adapt, mark, markNeighbors
+    from dune.fem.view import adaptiveLeafGridView, filteredGridView
+    from dune.grid import cartesianDomain, Marker
+    from dune.alugrid import aluConformGrid as GridView
+    import ufl
+    if hasattr(domain,"indexSet"):
+        return domain
+    if type(domain) == list:
+        domain = cartesianDomain(*domain)
+    gridView = adaptiveLeafGridView( GridView(domain) )
+    dims = (gridView.dimGrid, gridView.dimWorld)
+    x = ufl.SpatialCoordinate(dune.ufl.domain(dims))
+
+    if omega and factor>0 and refinements>0:
+        startEps = omega.epsilon.value * factor
+        region = gridFunction(
+            omega.phi(x) * (1 - omega.phi(x)), gridView=gridView
+        )  # interface
+        for j in range(1, refinements + 1):
+            omega.epsilon.value = startEps * j / refinements
+            mark(region, 0.00247262315663, maxLevel=refinements)  # 1 epsilon
+            adapt(gridView.hierarchicalGrid)
+        omega.epsilon.value = startEps/factor * 0.5**(refinements/dims[0])
+    if adaptLevels>0:
+        sd = gridFunction(omega(x), gridView=gridView)
+        bary = dims[0] * [1.0 / (dims[0]+1)]
+        lowerTol=-0.1
+        upperTol=0.1
+        for _ in range(adaptLevels):
+            def markStrategy(e):
+                v = sd(e, bary)
+                return (
+                    Marker.refine
+                    if v > lowerTol and v < upperTol
+                    else Marker.keep
+                )
+            gridView.hierarchicalGrid.adapt(markStrategy)
+            lowerTol /= 2
+            upperTol /= 2
+
+    if filterTolerance<0:
+        return gridView
+    if omega and filterTolerance>0:
+        sd = gridFunction(omega(x), gridView=gridView)
+        bary = dims[0] * [1.0 / (dims[0]+1)]
+        filter = lambda e: 1 if sd(e, bary) < filterTolerance else 2
+    else:
+        filter = lambda e: 1
+    return filteredGridView(
+               gridView, filter, domainId=1, useFilteredIndexSet=True)
+
+def solve(uflModel, domain):
+    # we use a Lagrange space
+    from dune.fem.space import lagrange
+    # and the build-in solver
+    from dune.fem.scheme import galerkin
+
+    gridView = GridView(domain)
+    space = lagrange(gridView, dimRange=1)
+    scheme = galerkin(uflModel, space=space,
+                      solver=("suitesparse","umfpack"))
+    uh = space.interpolate(0,name="ddm")
+    scheme.solve(target=uh)
+    return uh

ddfem-1.0.6/ddfem/__main__.py

@@ -0,0 +1,24 @@
+import sys, os
+try:
+    import dune
+except ImportError:
+    dune = None
+import ddfem
+
+def main():
+    if dune:
+        from dune.fem.utility import FemThreadPoolExecutor
+        from dune.generator import builder
+        builder.initialize()
+
+        dataPath = os.path.join(ddfem.__path__[0],"data")
+        with FemThreadPoolExecutor(max_workers=8) as executor:
+            for file in os.listdir(dataPath):
+                filename = os.path.join(dataPath,os.fsdecode(file))
+                if filename.endswith(".cc"):
+                    with open(filename,"r") as f:
+                        executor.submit( builder.load, file.replace(".cc",""), f.read(), None )
+
+if __name__ == "__main__":
+    main()
+

ddfem-1.0.6/ddfem/data/adapt_010e1455a389aa01992af42c0a3f4b14.cc

@@ -0,0 +1,27 @@
+#ifndef Guard_adapt_010e1455a389aa01992af42c0a3f4b14
+#define Guard_adapt_010e1455a389aa01992af42c0a3f4b14
+
+#include <config.h>
+
+#define USING_DUNE_PYTHON 1
+
+#include <dune/alugrid/dgf.hh>
+#include <dune/alugrid/grid.hh>
+#include <dune/fempy/py/grid/adaptation.hh>
+#include <dune/python/grid/hierarchical.hh>
+
+#include <dune/python/common/typeregistry.hh>
+#include <dune/python/pybind11/pybind11.h>
+#include <dune/python/pybind11/stl.h>
+
+PYBIND11_MODULE( adapt_010e1455a389aa01992af42c0a3f4b14, module )
+{
+  using pybind11::operator""_a;
+  pybind11::module cls0 = module;
+  {
+    using DuneType = Dune::FemPy::GridAdaptation< Dune::ALUGrid< 2, 2, Dune::simplex > >;
+    auto cls = Dune::Python::insertClass< DuneType >( cls0, "GridAdaptation", Dune::Python::GenerateTypeName("Dune::FemPy::GridAdaptation< Dune::ALUGrid< 2, 2, Dune::simplex > >"), Dune::Python::IncludeFiles{"dune/alugrid/dgf.hh","dune/alugrid/grid.hh","dune/fempy/py/grid/adaptation.hh","dune/python/grid/hierarchical.hh"}).first;
+    Dune::FemPy::registerGridAdaptation( cls0, cls );
+  }
+}
+#endif

ddfem-1.0.6/ddfem/data/femscheme_f2a3cb32a864b08c0236366e732fb908.cc

@@ -0,0 +1,43 @@
+#ifndef Guard_femscheme_f2a3cb32a864b08c0236366e732fb908
+#define Guard_femscheme_f2a3cb32a864b08c0236366e732fb908
+
+#include <config.h>
+
+#define USING_DUNE_PYTHON 1
+
+#include <dune/alugrid/dgf.hh>
+#include <dune/alugrid/grid.hh>
+#include <dune/fem/function/adaptivefunction.hh>
+#include <dune/fem/gridpart/adaptiveleafgridpart.hh>
+#include <dune/fem/gridpart/filter/simple.hh>
+#include <dune/fem/gridpart/filteredgridpart.hh>
+#include <dune/fem/operator/linear/spoperator.hh>
+#include <dune/fem/schemes/dirichletwrapper.hh>
+#include <dune/fem/schemes/molgalerkin.hh>
+#include <dune/fem/solver/ldlsolver.hh>
+#include <dune/fem/solver/spqrsolver.hh>
+#include <dune/fem/solver/umfpacksolver.hh>
+#include <dune/fem/space/lagrange.hh>
+#include <dune/fempy/parameter.hh>
+#include <dune/fempy/py/discretefunction.hh>
+#include <dune/fempy/py/gridview.hh>
+#include <dune/fempy/py/scheme.hh>
+#include <dune/fempy/py/space.hh>
+#include <dune/python/grid/gridview.hh>
+#include <dune/python/grid/hierarchical.hh>
+
+#include <dune/python/common/typeregistry.hh>
+#include <dune/python/pybind11/pybind11.h>
+#include <dune/python/pybind11/stl.h>
+
+PYBIND11_MODULE( femscheme_f2a3cb32a864b08c0236366e732fb908, module )
+{
+  using pybind11::operator""_a;
+  pybind11::module cls0 = module;
+  {
+    using DuneType = Dune::Fem::GalerkinScheme< Dune::Fem::VirtualizedIntegrands< typename Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage>::GridPartType, std::tuple< Dune::FieldVector< double, 1 >, Dune::FieldMatrix< double, 1, 2 > >, std::tuple< Dune::FieldVector< double, 1 >, Dune::FieldMatrix< double, 1, 2 > > >, Dune::Fem::SparseRowLinearOperator< Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >,Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >,Dune::Fem::SparseRowMatrix<double,int>>, Dune::Fem::UMFPACKInverseOperator< Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >, typename Dune::Fem::SparseRowLinearOperator< Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >,Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >,Dune::Fem::SparseRowMatrix<double,int>>::MatrixType >, true >;
+    auto cls = Dune::Python::insertClass< DuneType >( cls0, "Scheme",pybind11::dynamic_attr(), Dune::Python::GenerateTypeName("Dune::Fem::GalerkinScheme< Dune::Fem::VirtualizedIntegrands< typename Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage>::GridPartType, std::tuple< Dune::FieldVector< double, 1 >, Dune::FieldMatrix< double, 1, 2 > >, std::tuple< Dune::FieldVector< double, 1 >, Dune::FieldMatrix< double, 1, 2 > > >, Dune::Fem::SparseRowLinearOperator< Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >,Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >,Dune::Fem::SparseRowMatrix<double,int>>, Dune::Fem::UMFPACKInverseOperator< Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >, typename Dune::Fem::SparseRowLinearOperator< Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >,Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >,Dune::Fem::SparseRowMatrix<double,int>>::MatrixType >, true >"), Dune::Python::IncludeFiles{"dune/alugrid/dgf.hh","dune/alugrid/grid.hh","dune/fem/function/adaptivefunction.hh","dune/fem/gridpart/adaptiveleafgridpart.hh","dune/fem/gridpart/filter/simple.hh","dune/fem/gridpart/filteredgridpart.hh","dune/fem/operator/linear/spoperator.hh","dune/fem/schemes/dirichletwrapper.hh","dune/fem/schemes/molgalerkin.hh","dune/fem/solver/ldlsolver.hh","dune/fem/solver/spqrsolver.hh","dune/fem/solver/umfpacksolver.hh","dune/fem/space/lagrange.hh","dune/fempy/parameter.hh","dune/fempy/py/discretefunction.hh","dune/fempy/py/gridview.hh","dune/fempy/py/scheme.hh","dune/fempy/py/space.hh","dune/python/grid/gridview.hh","dune/python/grid/hierarchical.hh"}).first;
+    Dune::FemPy::registerScheme( cls0, cls );
+  }
+}
+#endif

ddfem-1.0.6/ddfem/data/femspace_4ff9049b6f297245e5d33296d14cd684_d5b136dbe3c5077b69c99b8c322eb563.cc

@@ -0,0 +1,45 @@
+#ifndef Guard_femspace_4ff9049b6f297245e5d33296d14cd684_d5b136dbe3c5077b69c99b8c322eb563
+#define Guard_femspace_4ff9049b6f297245e5d33296d14cd684_d5b136dbe3c5077b69c99b8c322eb563
+
+#include <config.h>
+
+#define USING_DUNE_PYTHON 1
+
+#define BACKENDNAME "as_numpy"
+#include <dune/alugrid/dgf.hh>
+#include <dune/alugrid/grid.hh>
+#include <dune/fem/function/adaptivefunction.hh>
+#include <dune/fem/gridpart/adaptiveleafgridpart.hh>
+#include <dune/fem/gridpart/filter/simple.hh>
+#include <dune/fem/gridpart/filteredgridpart.hh>
+#include <dune/fem/space/lagrange.hh>
+#include <dune/fempy/py/discretefunction.hh>
+#include <dune/fempy/py/gridview.hh>
+#include <dune/fempy/py/space.hh>
+#include <dune/python/grid/gridview.hh>
+#include <dune/python/grid/hierarchical.hh>
+
+#include <dune/python/common/typeregistry.hh>
+#include <dune/python/pybind11/pybind11.h>
+#include <dune/python/pybind11/stl.h>
+
+PYBIND11_MODULE( femspace_4ff9049b6f297245e5d33296d14cd684_d5b136dbe3c5077b69c99b8c322eb563, module )
+{
+  using pybind11::operator""_a;
+  pybind11::module cls0 = module;
+  {
+    using DuneType = Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage>;
+    auto cls = Dune::Python::insertClass< DuneType, std::shared_ptr<DuneType> >( cls0, "Space",pybind11::dynamic_attr(), Dune::Python::GenerateTypeName("Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage>"), Dune::Python::IncludeFiles{"dune/alugrid/dgf.hh","dune/alugrid/grid.hh","dune/fem/function/adaptivefunction.hh","dune/fem/gridpart/adaptiveleafgridpart.hh","dune/fem/gridpart/filter/simple.hh","dune/fem/gridpart/filteredgridpart.hh","dune/fem/space/lagrange.hh","dune/fempy/py/discretefunction.hh","dune/fempy/py/gridview.hh","dune/fempy/py/space.hh","dune/python/grid/gridview.hh","dune/python/grid/hierarchical.hh"}).first;
+    Dune::FemPy::registerSpace( cls0, cls );
+  }
+  using pybind11::operator""_a;
+  {
+    using DuneType = Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >;
+    auto cls = Dune::Python::insertClass< DuneType >( cls0, "DiscreteFunction",pybind11::buffer_protocol(),pybind11::dynamic_attr(), Dune::Python::GenerateTypeName("Dune::Fem::AdaptiveDiscreteFunction< Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage> >"), Dune::Python::IncludeFiles{"dune/alugrid/dgf.hh","dune/alugrid/grid.hh","dune/fem/function/adaptivefunction.hh","dune/fem/gridpart/adaptiveleafgridpart.hh","dune/fem/gridpart/filter/simple.hh","dune/fem/gridpart/filteredgridpart.hh","dune/fem/space/lagrange.hh","dune/fempy/py/discretefunction.hh","dune/fempy/py/gridview.hh","dune/fempy/py/space.hh","dune/python/grid/gridview.hh","dune/python/grid/hierarchical.hh"}).first;
+    Dune::FemPy::registerDiscreteFunction( cls0, cls );
+    cls.def( pybind11::init( [] ( const std::string &name, const Dune::Fem::DynamicLagrangeDiscreteFunctionSpace< Dune::Fem::FunctionSpace< double, double, 2, 1 >, 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::Fem::CodegenStorage>&space, pybind11::array_t<double> dofVector ) {double *dof = static_cast< double* >( dofVector.request(false).ptr );
+        return new DuneType(name,space,dof);
+      } ), "name"_a, "space"_a, "dofVector"_a, pybind11::keep_alive< 1, 3 >(), pybind11::keep_alive< 1, 4 >() );
+  }
+}
+#endif

ddfem-1.0.6/ddfem/data/generate.py

@@ -0,0 +1,30 @@
+# Step 1: in the doc folder run
+#         make clean
+#         DUNE_LOGMODULES=8 make -ij8
+# Step 2: in this folder run
+#         python generate.py
+
+import os, shutil
+from dune.common.module import getDunePyDir
+def main():
+    dunepy = getDunePyDir()
+    docPath = os.path.join("..","..","doc")
+    names = set()
+    for file in os.listdir(docPath):
+        filename = os.path.join(docPath,os.fsdecode(file))
+        if filename.endswith(".modules"):
+            with open(filename,"r") as f:
+                names.update(f.readlines())
+    for name in names:
+        name = name.strip()
+        print(name)
+        src = os.path.join(dunepy,"python","dune","generated",name+".cc")
+        try:
+            shutil.copy(src, ".")
+        except FileNotFoundError:
+            print(f"Error: can't copy {src}",flush=True)
+            pass
+
+if __name__ == "__main__":
+    main()
+

ddfem-1.0.6/ddfem/data/hierarchicalgrid_d827e80243feae36e8b9006bb24207e7.cc

@@ -0,0 +1,26 @@
+#ifndef Guard_hierarchicalgrid_d827e80243feae36e8b9006bb24207e7
+#define Guard_hierarchicalgrid_d827e80243feae36e8b9006bb24207e7
+
+#include <config.h>
+
+#define USING_DUNE_PYTHON 1
+
+#include <dune/alugrid/dgf.hh>
+#include <dune/alugrid/grid.hh>
+#include <dune/python/grid/hierarchical.hh>
+
+#include <dune/python/common/typeregistry.hh>
+#include <dune/python/pybind11/pybind11.h>
+#include <dune/python/pybind11/stl.h>
+
+PYBIND11_MODULE( hierarchicalgrid_d827e80243feae36e8b9006bb24207e7, module )
+{
+  using pybind11::operator""_a;
+  pybind11::module cls0 = module;
+  {
+    using DuneType = Dune::ALUGrid< 2, 2, Dune::simplex >;
+    auto cls = Dune::Python::insertClass< DuneType, std::shared_ptr<DuneType> >( cls0, "HierarchicalGrid",pybind11::dynamic_attr(), Dune::Python::GenerateTypeName("Dune::ALUGrid< 2, 2, Dune::simplex >"), Dune::Python::IncludeFiles{"dune/alugrid/dgf.hh","dune/alugrid/grid.hh","dune/python/grid/hierarchical.hh"}).first;
+    Dune::Python::registerHierarchicalGrid( cls0, cls );
+  }
+}
+#endif

ddfem-1.0.6/ddfem/data/indexset_e9c0cdc96a4fc702a0969c29910954e6.cc

@@ -0,0 +1,32 @@
+#ifndef Guard_indexset_e9c0cdc96a4fc702a0969c29910954e6
+#define Guard_indexset_e9c0cdc96a4fc702a0969c29910954e6
+
+#include <config.h>
+
+#define USING_DUNE_PYTHON 1
+
+#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/python/grid/indexset.hh>
+
+#include <dune/python/common/typeregistry.hh>
+#include <dune/python/pybind11/pybind11.h>
+#include <dune/python/pybind11/stl.h>
+
+PYBIND11_MODULE( indexset_e9c0cdc96a4fc702a0969c29910954e6, module )
+{
+  using pybind11::operator""_a;
+  pybind11::module cls0 = module;
+  {
+    using DuneType = 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 >::IndexSet;
+    auto cls = Dune::Python::insertClass< DuneType >( cls0, "GridViewIndexSet", Dune::Python::GenerateTypeName("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 >::IndexSet"), Dune::Python::IncludeFiles{"dune/alugrid/dgf.hh","dune/alugrid/grid.hh","dune/fem/gridpart/adaptiveleafgridpart.hh","dune/fem/gridpart/filter/simple.hh","dune/fem/gridpart/filteredgridpart.hh","dune/fempy/py/gridview.hh","dune/python/grid/gridview.hh","dune/python/grid/hierarchical.hh","dune/python/grid/indexset.hh"}).first;
+    Dune::Python::registerGridViewIndexSet( cls0, cls );
+  }
+}
+#endif