ngsolve 6.2.2506.post60.dev0__cp312-cp312-macosx_10_15_universal2.whl → 6.2.2506.post74.dev0__cp312-cp312-macosx_10_15_universal2.whl

netgen/include/kernels.hpp

@@ -94,14 +94,14 @@ namespace ngsbem
     
     Array<KernelTerm> terms;
     
-    auto CreateMultipoleExpansion (Vec<3> c, double r) const
+    auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<SingularMLExpansion<mp_type>> (c, r, 1e-16);
+      return make_shared<SingularMLExpansion<mp_type>> (c, r, 1e-16, fmm_params);
     }
 
-    auto CreateLocalExpansion (Vec<3> c, double r) const
+    auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<RegularMLExpansion<mp_type>> (c, r, 1e-16);
+      return make_shared<RegularMLExpansion<mp_type>> (c, r, 1e-16, fmm_params);
     }
 
     void AddSource (SingularMLExpansion<mp_type> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<double> val) const
@@ -157,14 +157,14 @@ namespace ngsbem
 
     Array<KernelTerm> terms;
 
-    auto CreateMultipoleExpansion (Vec<3> c, double r) const
+    auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<SingularMLExpansion<mp_type>> (c, r, 1e-16);
+      return make_shared<SingularMLExpansion<mp_type>> (c, r, 1e-16, fmm_params);
     }
 
-    auto CreateLocalExpansion (Vec<3> c, double r) const
+    auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<RegularMLExpansion<mp_type>> (c, r, 1e-16);
+      return make_shared<RegularMLExpansion<mp_type>> (c, r, 1e-16, fmm_params);
     }
 
     void AddSource (SingularMLExpansion<mp_type> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<double> val) const
@@ -231,14 +231,14 @@ namespace ngsbem
     double GetKappa() const { return kappa; }
     Array<KernelTerm> terms = { KernelTerm{1.0, 0, 0, 0}, };    
 
-    auto CreateMultipoleExpansion (Vec<3> c, double r) const
+    auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<SingularMLExpansion<Complex>> (c, r, kappa);
+      return make_shared<SingularMLExpansion<Complex>> (c, r, kappa, fmm_params);
     }
 
-    auto CreateLocalExpansion (Vec<3> c, double r) const
+    auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<RegularMLExpansion<Complex>> (c, r, kappa);
+      return make_shared<RegularMLExpansion<Complex>> (c, r, kappa, fmm_params);
     }
 
     void AddSource (SingularMLExpansion<Complex> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
@@ -285,14 +285,14 @@ namespace ngsbem
     double GetKappa() const { return kappa; }
     Array<KernelTerm> terms = { KernelTerm{1.0, 0, 0, 0}, };    
 
-    auto CreateMultipoleExpansion (Vec<3> c, double r) const
+    auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<SingularMLExpansion<Complex>> (c, r, kappa);
+      return make_shared<SingularMLExpansion<Complex>> (c, r, kappa, fmm_params);
     }
 
-    auto CreateLocalExpansion (Vec<3> c, double r) const
+    auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<RegularMLExpansion<Complex>> (c, r, kappa);
+      return make_shared<RegularMLExpansion<Complex>> (c, r, kappa, fmm_params);
     }
 
     void AddSource (SingularMLExpansion<Complex> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
@@ -345,14 +345,14 @@ namespace ngsbem
         KernelTerm{1.0, 0, 2, 2},
       };
 
-    auto CreateMultipoleExpansion (Vec<3> c, double r) const
+    auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<SingularMLExpansion<Vec<3,Complex>>> (c, r, kappa);
+      return make_shared<SingularMLExpansion<Vec<3,Complex>>> (c, r, kappa, fmm_params);
     }
 
-    auto CreateLocalExpansion (Vec<3> c, double r) const
+    auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<RegularMLExpansion<Vec<3,Complex>>> (c, r, kappa);
+      return make_shared<RegularMLExpansion<Vec<3,Complex>>> (c, r, kappa, fmm_params);
     }
 
     void AddSource (SingularMLExpansion<Vec<3,Complex>> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
@@ -398,14 +398,14 @@ namespace ngsbem
 	KernelTerm{1.0, 1, 3, 3},
       };
 
-    auto CreateMultipoleExpansion (Vec<3> c, double r) const
+    auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<SingularMLExpansion<Vec<6,Complex>>> (c, r, kappa);
+      return make_shared<SingularMLExpansion<Vec<6,Complex>>> (c, r, kappa, fmm_params);
     }
  
-    auto CreateLocalExpansion (Vec<3> c, double r) const
+    auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<RegularMLExpansion<Vec<6,Complex>>> (c, r, kappa);
+      return make_shared<RegularMLExpansion<Vec<6,Complex>>> (c, r, kappa, fmm_params);
     }
 
     void AddSource (SingularMLExpansion<Vec<6,Complex>> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
@@ -461,14 +461,14 @@ namespace ngsbem
     double GetKappa() const { return kappa; }
     Array<KernelTerm> terms = { KernelTerm{1.0, 0, 0, 0}, };    
 
-    auto CreateMultipoleExpansion (Vec<3> c, double r) const
+    auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<SingularMLExpansion<Complex>> (c, r, kappa);
+      return make_shared<SingularMLExpansion<Complex>> (c, r, kappa, fmm_params);
     }
 
-    auto CreateLocalExpansion (Vec<3> c, double r) const
+    auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<RegularMLExpansion<Complex>> (c, r, kappa);
+      return make_shared<RegularMLExpansion<Complex>> (c, r, kappa, fmm_params);
     }
 
     void AddSource (SingularMLExpansion<Complex> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
@@ -519,14 +519,14 @@ namespace ngsbem
     double GetKappa() const { return kappa; }
     Array<KernelTerm> terms;
 
-    auto CreateMultipoleExpansion (Vec<3> c, double r) const
+    auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<SingularMLExpansion<Vec<4,Complex>>> (c, r, kappa);
+      return make_shared<SingularMLExpansion<Vec<4,Complex>>> (c, r, kappa, fmm_params);
     }
 
-    auto CreateLocalExpansion (Vec<3> c, double r) const
+    auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<RegularMLExpansion<Vec<4,Complex>>> (c, r, kappa);
+      return make_shared<RegularMLExpansion<Vec<4,Complex>>> (c, r, kappa, fmm_params);
     }
 
     void AddSource (SingularMLExpansion<Vec<4,Complex>> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
@@ -586,14 +586,14 @@ namespace ngsbem
         KernelTerm{-1.0, 2, 1, 0},
       };    
 
-    auto CreateMultipoleExpansion (Vec<3> c, double r) const
+    auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<SingularMLExpansion<Vec<3,Complex>>> (c, r, kappa);
+      return make_shared<SingularMLExpansion<Vec<3,Complex>>> (c, r, kappa, fmm_params);
     }
 
-    auto CreateLocalExpansion (Vec<3> c, double r) const
+    auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<RegularMLExpansion<Vec<3,Complex>>> (c, r, kappa);
+      return make_shared<RegularMLExpansion<Vec<3,Complex>>> (c, r, kappa, fmm_params);
     }
 
     void AddSource (SingularMLExpansion<Vec<3,Complex>> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
@@ -700,14 +700,14 @@ namespace ngsbem
     }
 
     
-    auto CreateMultipoleExpansion (Vec<3> c, double r) const
+    auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<SingularMLExpansion<Vec<6,Complex>>> (c, r, 1e-16);
+      return make_shared<SingularMLExpansion<Vec<6,Complex>>> (c, r, 1e-16, fmm_params);
     }
 
-    auto CreateLocalExpansion (Vec<3> c, double r) const
+    auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<RegularMLExpansion<Vec<6,Complex>>> (c, r, 1e-16);
+      return make_shared<RegularMLExpansion<Vec<6,Complex>>> (c, r, 1e-16, fmm_params);
     }
 
     

netgen/include/mp_coefficient.hpp

@@ -119,9 +119,9 @@ namespace ngsbem
     shared_ptr<RegularMLExpansion<entry_type>> mlmp;
   public:
     RegularMLExpansionCF (shared_ptr<SingularMLExpansionCF<entry_type>> asingmp, Vec<3> center, double r)
-      : CoefficientFunction(sizeof(entry_type)/sizeof(Complex), true), mlmp{make_shared<RegularMLExpansion<entry_type>>(asingmp->MLExpansion(), center, r)} { } 
+      : CoefficientFunction(sizeof(entry_type)/sizeof(Complex), true), mlmp{make_shared<RegularMLExpansion<entry_type>>(asingmp->MLExpansion(), center, r, FMM_Parameters())} { } 
     RegularMLExpansionCF (shared_ptr<SingularMLExpansion<entry_type>> asingmp, Vec<3> center, double r)
-      : CoefficientFunction(sizeof(entry_type)/sizeof(Complex), true), mlmp{make_shared<RegularMLExpansion<entry_type>>(asingmp, center, r)} { } 
+      : CoefficientFunction(sizeof(entry_type)/sizeof(Complex), true), mlmp{make_shared<RegularMLExpansion<entry_type>>(asingmp, center, r, FMM_Parameters())} { } 
     
     virtual double Evaluate (const BaseMappedIntegrationPoint & ip) const override
     { throw Exception("real eval not available"); }

netgen/include/mptools.hpp

@@ -58,6 +58,15 @@ namespace ngsbem
   
 
 
+  class FMM_Parameters
+  {
+  public:
+    int maxdirect = 100;
+    int minorder = 20;    // order = minorder + 2 kappa r 
+  };
+
+
+  
   
   inline std::tuple<double, double, double> SphericalCoordinates(Vec<3> dist){
     double len, theta, phi;
@@ -516,12 +525,14 @@ namespace ngsbem
 
   // ***************** parameters ****************
 
+  /*
   static constexpr int MPOrder (double rho_kappa)
   {
     // return max (20, int(2*rho_kappa));
     return 20+int(2*rho_kappa);
   }
   static constexpr int maxdirect = 100;
+  */
 
 
   template <typename SCAL, auto S>
@@ -665,11 +676,13 @@ namespace ngsbem
                   Vec<3,SIMD<double,FMM_SW>>, simd_entry_type>> simd_chargedipoles;
       
       int total_sources;
+      const FMM_Parameters & fmm_params;
       std::mutex node_mutex;
       atomic<bool> have_childs{false};
       
-      Node (Vec<3> acenter, double ar, int alevel, double akappa)
-        : center(acenter), r(ar), level(alevel), mp(MPOrder(ar*akappa), akappa, ar) // min(1.0, ar*akappa))
+      Node (Vec<3> acenter, double ar, int alevel, double akappa, const FMM_Parameters & afmm_params)
+      // : center(acenter), r(ar), level(alevel), mp(MPOrder(ar*akappa), akappa, ar), fmm_params(afmm_params)
+        : center(acenter), r(ar), level(alevel), mp(afmm_params.minorder+2*ar*akappa, akappa, ar), fmm_params(afmm_params)
       {
         if (level < nodes_on_level.Size())
           nodes_on_level[level]++;
@@ -693,7 +706,7 @@ namespace ngsbem
             cc(0) += (i&1) ? r/2 : -r/2;
             cc(1) += (i&2) ? r/2 : -r/2;
             cc(2) += (i&4) ? r/2 : -r/2;
-            childs[i] = make_unique<Node> (cc, r/2, level+1, mp.Kappa());
+            childs[i] = make_unique<Node> (cc, r/2, level+1, mp.Kappa(), fmm_params);
           }
         have_childs = true;
       }
@@ -742,7 +755,7 @@ namespace ngsbem
 
         // if (r*mp.Kappa() < 1e-8) return;
         if (level > 20) return;
-        if (charges.Size() < maxdirect && r*mp.Kappa() < 5)
+        if (charges.Size() < fmm_params.maxdirect && r*mp.Kappa() < 5)
           return;
         
         SendSourcesToChilds();
@@ -772,7 +785,7 @@ namespace ngsbem
         dipoles.Append (tuple{x,d,c});
         
         if (level > 20) return;
-        if (dipoles.Size() < maxdirect)
+        if (dipoles.Size() < fmm_params.maxdirect)
           return;
 
         SendSourcesToChilds();
@@ -801,7 +814,7 @@ namespace ngsbem
         
         chargedipoles.Append (tuple{x,c,dir,c2});
 
-        if (chargedipoles.Size() < maxdirect || r < 1e-8)
+        if (chargedipoles.Size() < fmm_params.maxdirect || r < 1e-8)
           return;
 
         SendSourcesToChilds();
@@ -1269,12 +1282,13 @@ namespace ngsbem
       }
     };
     
-    Node root;
+    FMM_Parameters fmm_params;
+    Node root;    
     bool havemp = false;
     
   public:
-    SingularMLExpansion (Vec<3> center, double r, double kappa)
-      : root(center, r, 0, kappa)
+    SingularMLExpansion (Vec<3> center, double r, double kappa, FMM_Parameters _params = FMM_Parameters())
+      : fmm_params(_params), root(center, r, 0, kappa, fmm_params)
     {
       nodes_on_level = 0;
       nodes_on_level[0] = 1;
@@ -1661,11 +1675,13 @@ namespace ngsbem
       atomic<bool> have_childs{false};
 
       Array<const typename SingularMLExpansion<elem_type>::Node*> singnodes;
+      const FMM_Parameters & params;
 
-      Node (Vec<3> acenter, double ar, int alevel, double kappa)
+      
+      Node (Vec<3> acenter, double ar, int alevel, double kappa, const FMM_Parameters & _params)
         : center(acenter), r(ar), level(alevel),
           // mp(MPOrder(ar*kappa), kappa, ar) // 1.0/min(1.0, 0.25*r*kappa))
-          mp(-1, kappa, ar)
+          mp(-1, kappa, ar), params(_params)
           // : center(acenter), r(ar), level(alevel), mp(MPOrder(ar*kappa), kappa, 1.0)
       {
         if (level < nodes_on_level.Size())
@@ -1674,7 +1690,8 @@ namespace ngsbem
 
       void Allocate()
       {
-        mp = SphericalExpansion<Regular,elem_type>(MPOrder(r*mp.Kappa()), mp.Kappa(), r);
+        // mp = SphericalExpansion<Regular,elem_type>(MPOrder(r*mp.Kappa()), mp.Kappa(), r);
+        mp = SphericalExpansion<Regular,elem_type>(params.minorder+2*r*mp.Kappa(), mp.Kappa(), r);        
       }
       
       
@@ -1688,7 +1705,7 @@ namespace ngsbem
             cc(0) += (i&1) ? r/2 : -r/2;
             cc(1) += (i&2) ? r/2 : -r/2;
             cc(2) += (i&4) ? r/2 : -r/2;
-            childs[i] = make_unique<Node> (cc, r/2, level+1, mp.Kappa());
+            childs[i] = make_unique<Node> (cc, r/2, level+1, mp.Kappa(), params);
             if (allocate)
               childs[i] -> Allocate();
           }
@@ -1909,7 +1926,7 @@ namespace ngsbem
 
         // if (r*mp.Kappa() < 1e-8) return;
         if (level > 20) return;        
-        if (targets.Size() < maxdirect && r*mp.Kappa() < 5)
+        if (targets.Size() < params.maxdirect && r*mp.Kappa() < 5)
           return;
 
         CreateChilds();
@@ -1949,7 +1966,7 @@ namespace ngsbem
         vol_targets.Append (tuple(x,tr));
 
         if (level > 20) return;
-        if (vol_targets.Size() < maxdirect && (r*mp.Kappa() < 5))
+        if (vol_targets.Size() < params.maxdirect && (r*mp.Kappa() < 5))
           return;
 
         CreateChilds();
@@ -2018,14 +2035,16 @@ namespace ngsbem
       }
 
     };
-    
+
+    FMM_Parameters fmm_params;
     Node root;
     shared_ptr<SingularMLExpansion<elem_type>> singmp;
     
   public:
-  RegularMLExpansion (shared_ptr<SingularMLExpansion<elem_type>> asingmp, Vec<3> center, double r)
-      : root(center, r, 0, asingmp->Kappa()), singmp(asingmp)
-    {
+  RegularMLExpansion (shared_ptr<SingularMLExpansion<elem_type>> asingmp, Vec<3> center, double r,
+                      const FMM_Parameters & _params)
+  : fmm_params(_params), root(center, r, 0, asingmp->Kappa(), fmm_params), singmp(asingmp)
+  {
       if (!singmp->havemp) throw Exception("first call Calc for singular MP");
       root.Allocate();
       
@@ -2054,13 +2073,13 @@ namespace ngsbem
       }
     }
 
-    RegularMLExpansion (Vec<3> center, double r, double kappa)
-      : root(center, r, 0, kappa)
-    {
-      nodes_on_level = 0;
-      nodes_on_level[0] = 1;
-    }
-    
+  RegularMLExpansion (Vec<3> center, double r, double kappa, const FMM_Parameters & _params)
+  : fmm_params(_params), root(center, r, 0, kappa, fmm_params)
+  {
+    nodes_on_level = 0;
+    nodes_on_level[0] = 1;
+  }
+  
     void AddTarget (Vec<3> t)
     {
       root.AddTarget (t);

netgen/include/preconditioner.hpp

@@ -414,7 +414,7 @@ namespace ngcomp
     ///
     // MGPreconditioner * low_order_preconditioner;
     ///
-    shared_ptr<Preconditioner> coarse_pre;
+    shared_ptr<BaseMatrix> coarse_pre;
     ///
     int finesmoothingsteps;
     ///
@@ -462,7 +462,7 @@ namespace ngcomp
     void MgTest () const;
 
     void SetDirectSolverCluster(shared_ptr<Array<int>> cluster);
-    void SetCoarsePreconditioner(shared_ptr<Preconditioner> prec);
+    void SetCoarsePreconditioner(shared_ptr<BaseMatrix> prec);
   };
 
   class CommutingAMGPreconditioner : public Preconditioner

ngsolve/__init__.pyi

@@ -16,18 +16,15 @@ from builtins import sum as builtin_sum
 import netgen as netgen
 from netgen import Redraw
 from netgen import TimeFunction
-from ngsolve.bla import InnerProduct
-from ngsolve.bla import Matrix
-from ngsolve.bla import Norm
-from ngsolve.bla import Vector
+from ngsolve.bla.pybind11_detail_function_record_v1_system_libcpp_abi1 import InnerProduct
+from ngsolve.bla.pybind11_detail_function_record_v1_system_libcpp_abi1 import Matrix
+from ngsolve.bla.pybind11_detail_function_record_v1_system_libcpp_abi1 import Norm
+from ngsolve.bla.pybind11_detail_function_record_v1_system_libcpp_abi1 import Vector
 from ngsolve.comp import APhiHCurlAMG
 from ngsolve.comp import BilinearForm
-from ngsolve.comp import BoundaryFromVolumeCF
 from ngsolve.comp import COUPLING_TYPE
 from ngsolve.comp import Compress
-from ngsolve.comp import CompressCompound
 from ngsolve.comp import ContactBoundary
-from ngsolve.comp import ConvertOperator
 from ngsolve.comp import Discontinuous
 from ngsolve.comp import ElementId
 from ngsolve.comp import FESpace
@@ -44,10 +41,8 @@ from ngsolve.comp import HDiv
 from ngsolve.comp import HDivDiv
 from ngsolve.comp import HDivDivSurface
 from ngsolve.comp import HDivSurface
-from ngsolve.comp import Hidden as PrivateSpace
 from ngsolve.comp import Hidden
-from ngsolve.comp import Integrate
-from ngsolve.comp import Interpolate
+from ngsolve.comp import Hidden as PrivateSpace
 from ngsolve.comp import L2
 from ngsolve.comp import LinearForm
 from ngsolve.comp import MatrixValued
@@ -59,18 +54,12 @@ from ngsolve.comp import NormalFacetFESpace
 from ngsolve.comp import NormalFacetSurface
 from ngsolve.comp import NumberSpace
 from ngsolve.comp import ORDER_POLICY
-from ngsolve.comp import PatchwiseSolve
 from ngsolve.comp import Periodic
 from ngsolve.comp import PlateauFESpace
 from ngsolve.comp import Preconditioner
 from ngsolve.comp import ProductSpace
 from ngsolve.comp import Region
-from ngsolve.comp import SetHeapSize
-from ngsolve.comp import SetTestoutFile
 from ngsolve.comp import SurfaceL2
-from ngsolve.comp import SymbolicBFI
-from ngsolve.comp import SymbolicEnergy
-from ngsolve.comp import SymbolicLFI
 from ngsolve.comp import TangentialFacetFESpace
 from ngsolve.comp import TangentialSurfaceL2
 from ngsolve.comp import VTKOutput
@@ -83,70 +72,81 @@ from ngsolve.comp import VectorNodalFESpace
 from ngsolve.comp import VectorSurfaceL2
 from ngsolve.comp import VectorValued
 from ngsolve.comp import pml
+from ngsolve.comp.pybind11_detail_function_record_v1_system_libcpp_abi1 import BoundaryFromVolumeCF
+from ngsolve.comp.pybind11_detail_function_record_v1_system_libcpp_abi1 import CompressCompound
+from ngsolve.comp.pybind11_detail_function_record_v1_system_libcpp_abi1 import ConvertOperator
+from ngsolve.comp.pybind11_detail_function_record_v1_system_libcpp_abi1 import Integrate
+from ngsolve.comp.pybind11_detail_function_record_v1_system_libcpp_abi1 import Interpolate
+from ngsolve.comp.pybind11_detail_function_record_v1_system_libcpp_abi1 import PatchwiseSolve
+from ngsolve.comp.pybind11_detail_function_record_v1_system_libcpp_abi1 import SetHeapSize
+from ngsolve.comp.pybind11_detail_function_record_v1_system_libcpp_abi1 import SetTestoutFile
+from ngsolve.comp.pybind11_detail_function_record_v1_system_libcpp_abi1 import SymbolicBFI
+from ngsolve.comp.pybind11_detail_function_record_v1_system_libcpp_abi1 import SymbolicEnergy
+from ngsolve.comp.pybind11_detail_function_record_v1_system_libcpp_abi1 import SymbolicLFI
 from ngsolve.fem import BFI
 from ngsolve.fem import BSpline
-from ngsolve.fem import BlockBFI
-from ngsolve.fem import BlockLFI
-from ngsolve.fem import CacheCF
 from ngsolve.fem import CoefficientFunction
 from ngsolve.fem import CoefficientFunction as CF
-from ngsolve.fem import Cof
-from ngsolve.fem import CompoundBFI
-from ngsolve.fem import CompoundLFI
-from ngsolve.fem import Conj
-from ngsolve.fem import Cross
-from ngsolve.fem import Det
 from ngsolve.fem import ET
-from ngsolve.fem import Id
-from ngsolve.fem import IfPos
 from ngsolve.fem import IntegrationRule
-from ngsolve.fem import Inv
 from ngsolve.fem import LFI
 from ngsolve.fem import Parameter
 from ngsolve.fem import ParameterC
 from ngsolve.fem import PlaceholderCF
-from ngsolve.fem import Skew
-from ngsolve.fem import Sym
-from ngsolve.fem import Trace
-from ngsolve.fem import VoxelCoefficient
-from ngsolve.fem import Zero as ZeroCF
-from ngsolve.fem import acos
-from ngsolve.fem import asin
-from ngsolve.fem import atan
-from ngsolve.fem import atan2
-from ngsolve.fem import ceil
-from ngsolve.fem import cos
-from ngsolve.fem import cosh
-from ngsolve.fem import erf
-from ngsolve.fem import exp
-from ngsolve.fem import floor
-from ngsolve.fem import log
-from ngsolve.fem import pow
-from ngsolve.fem import sin
-from ngsolve.fem import sinh
-from ngsolve.fem import sqrt
-from ngsolve.fem import tan
-from ngsolve.la import ArnoldiSolver
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import BlockBFI
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import BlockLFI
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import CacheCF
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import Cof
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import CompoundBFI
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import CompoundLFI
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import Conj
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import Cross
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import Det
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import Id
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import IfPos
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import Inv
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import Skew
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import Sym
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import Trace
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import VoxelCoefficient
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import Zero as ZeroCF
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import acos
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import asin
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import atan
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import atan2
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import ceil
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import cos
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import cosh
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import erf
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import exp
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import floor
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import log
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import pow
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import sin
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import sinh
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import sqrt
+from ngsolve.fem.pybind11_detail_function_record_v1_system_libcpp_abi1 import tan
 from ngsolve.la import BaseMatrix
 from ngsolve.la import BaseVector
 from ngsolve.la import BlockMatrix
 from ngsolve.la import BlockVector
-from ngsolve.la import CGSolver
 from ngsolve.la import ConstEBEMatrix
-from ngsolve.la import CreateVVector
 from ngsolve.la import DiagonalMatrix
 from ngsolve.la import Embedding
-from ngsolve.la import GMRESSolver
 from ngsolve.la import IdentityMatrix
 from ngsolve.la import MultiVector
 from ngsolve.la import PARALLEL_STATUS
 from ngsolve.la import ParallelMatrix
 from ngsolve.la import PermutationMatrix
 from ngsolve.la import Projector
-from ngsolve.la import QMRSolver
+from ngsolve.la.pybind11_detail_function_record_v1_system_libcpp_abi1 import ArnoldiSolver
+from ngsolve.la.pybind11_detail_function_record_v1_system_libcpp_abi1 import CGSolver
+from ngsolve.la.pybind11_detail_function_record_v1_system_libcpp_abi1 import CreateVVector
+from ngsolve.la.pybind11_detail_function_record_v1_system_libcpp_abi1 import GMRESSolver
+from ngsolve.la.pybind11_detail_function_record_v1_system_libcpp_abi1 import QMRSolver
 from ngsolve.ngstd import IntRange
-from ngsolve.solve import Draw
-from ngsolve.solve import SetVisualization
+from ngsolve.solve.pybind11_detail_function_record_v1_system_libcpp_abi1 import Draw
+from ngsolve.solve.pybind11_detail_function_record_v1_system_libcpp_abi1 import SetVisualization
 from ngsolve.solve_implementation.Application import Solve
 from ngsolve.timing import Timing
 from ngsolve.utils import Deviator
@@ -168,10 +168,10 @@ from ngsolve.utils import printonce
 import os as os
 from pyngcore.pyngcore import BitArray
 from pyngcore.pyngcore import PajeTrace
-from pyngcore.pyngcore import SetNumThreads
 from pyngcore.pyngcore import TaskManager
 from pyngcore.pyngcore import Timer
-from pyngcore.pyngcore import Timers
+from pyngcore.pyngcore.pybind11_detail_function_record_v1_system_libcpp_abi1 import SetNumThreads
+from pyngcore.pyngcore.pybind11_detail_function_record_v1_system_libcpp_abi1 import Timers
 import sys as sys
 from . import bla
 from . import bvp
@@ -219,7 +219,7 @@ TET: fem.ET  # value = <ET.TET: 20>
 TRIG: fem.ET  # value = <ET.TRIG: 10>
 VERTEX: fem.NODE_TYPE  # value = <NODE_TYPE.VERTEX: 0>
 VOL: comp.VorB  # value = <VorB.VOL: 0>
-__version__: str = '6.2.2506-60-g1d0a21976'
+__version__: str = '6.2.2506-74-gd430159c1'
 ds: comp.DifferentialSymbol  # value = <ngsolve.comp.DifferentialSymbol object>
 dx: comp.DifferentialSymbol  # value = <ngsolve.comp.DifferentialSymbol object>
 ngsglobals: comp.GlobalVariables  # value = <ngsolve.comp.GlobalVariables object>