ngsolve 6.2.2501.post83.dev1__cp313-cp313-macosx_10_15_universal2.whl → 6.2.2502.post11.dev1__cp313-cp313-macosx_10_15_universal2.whl

netgen/include/ngblas.hpp

@@ -655,6 +655,20 @@ namespace ngbla
   template <bool ADD, ORDERING ord>
   extern NGS_DLL_HEADER  
   void NgGEMV (Complex s, BareSliceMatrix<double,ord> a, FlatVector<const Complex> x, FlatVector<Complex> y) NETGEN_NOEXCEPT;
+  
+
+
+  template <bool ADD, ORDERING ord, int S>
+  void NgGEMV (Complex s, BareSliceMatrix<double,ord> a, SliceVector<Vec<S,Complex>> x, SliceVector<Vec<S,Complex>> y) NETGEN_NOEXCEPT
+  {
+    SliceMatrix<Complex> mx(x.Size(), S, S*x.Dist(), x.Data()->Data());
+    SliceMatrix<Complex> my(y.Size(), S, S*y.Dist(), y.Data()->Data());
+    if constexpr (ADD)
+      my += s * a.AddSize(y.Size(), x.Size()) * mx;
+    else
+      my = s * a.AddSize(y.Size(), x.Size()) * mx;
+  }
+
 
   
   template <bool ADD, ORDERING ord>

netgen/include/potentialtools.hpp

@@ -0,0 +1,20 @@
+#ifndef FILE_POTENTIALS
+#define FILE_POTENTIALS
+
+/*
+  tools for computing with potentials using multipoles
+ */
+
+
+#include <mptools.hpp>
+#include <meshaccess.hpp>
+
+
+namespace ngcomp
+{
+  
+  extern void AddCurrentDensity (SingularMLMultiPole<Vec<3,Complex>> & mp, shared_ptr<CoefficientFunction> current, ngcomp::Region reg);
+ 
+}
+
+#endif

netgen/include/prolongation.hpp

@@ -8,7 +8,7 @@
 /*********************************************************************/
 
 namespace ngcomp { class NedelecFESpace; }
-
+namespace ngla { class SparseFactorization; } 
 namespace ngmg
 {
 
@@ -229,6 +229,148 @@ namespace ngmg
     ///
     virtual void RestrictInline (int finelevel, BaseVector & v) const override;
   };
+
+  class NGS_DLL_HEADER HarmonicProlongation : public Prolongation
+  {
+    shared_ptr<Prolongation> baseprol;
+    weak_ptr<BilinearForm> wbfa;
+    Array<shared_ptr<BaseMatrix>> innerinverses;
+    Array<size_t> facets_on_level;
+    Array<shared_ptr<BitArray>> freedofs;
+    string inverse;
+  public:
+    HarmonicProlongation (shared_ptr<Prolongation> abaseprol,
+                          shared_ptr<BilinearForm> abfa, string ainverse);
+
+    virtual void Update (const FESpace & fes) override;
+    
+    virtual size_t GetNDofLevel (int level) override { return baseprol->GetNDofLevel(level); } 
+
+    virtual shared_ptr<SparseMatrix< double >> CreateProlongationMatrix( int finelevel ) const override;
+    virtual void ProlongateInline (int finelevel, BaseVector & v) const override;
+    virtual void RestrictInline (int finelevel, BaseVector & v) const override;
+  };
+
+
+
+
+  
+  template <int DIM>
+  class FacetProlongation : public Prolongation
+  {
+  protected:
+    shared_ptr<MeshAccess> ma;
+    
+    int dofs_per_facet = -1;
+    Array<Array<int>> first_dofs;
+    Array<int> prol_class;  // which prol to use ?
+    Array<size_t> facets_on_level;
+
+    static constexpr int NUMPROL = (DIM==2) ? 2 : 32;
+    array<Matrix<double>,NUMPROL> facetprol;
+    
+    bool haveprols = false;
+    
+  public:
+    FacetProlongation (shared_ptr<MeshAccess> ama)
+      : ma(ama) { }
+
+
+    virtual shared_ptr<SparseMatrix< double >> CreateProlongationMatrix( int finelevel ) const override
+    {
+      return nullptr;  // or NULL ? 
+    }
+
+    
+    virtual void ProlongateInline (int finelevel, BaseVector & v) const override
+    {
+      auto fv = v.FV<double>();
+      Matrix<double> tmp(facets_on_level[finelevel], dofs_per_facet);
+      tmp = 0.0;
+      
+      for (size_t i = 0; i < first_dofs[finelevel-1].Size()-1; i++)
+        {
+          IntRange r(first_dofs[finelevel-1][i], first_dofs[finelevel-1][i+1]);
+          if (r.Size() > 0)
+            tmp.Row(i) = fv.Range(r);
+        }
+
+      if constexpr (DIM==2)
+        {
+          for (int i = facets_on_level[finelevel-1]; i < facets_on_level[finelevel]; i++)
+            if (auto parents = get<1>(ma->GetParentEdges(i)); parents[1] == -1)
+              {
+                int pe = parents[0];
+                tmp.Row(i) = facetprol[prol_class[i]] * tmp.Row(pe);
+              }
+        }
+      else
+        {
+          for (int i = facets_on_level[finelevel-1]; i < facets_on_level[finelevel]; i++)
+          if (auto parents = get<1>(ma->GetParentFaces(i)); parents[1] == -1)
+            {
+              int pe = parents[0];
+              tmp.Row(i) = facetprol[prol_class[i]] * tmp.Row(pe);
+            }
+        }
+      
+      for (size_t i = 0; i < first_dofs[finelevel].Size()-1; i++)
+        {
+          IntRange r(first_dofs[finelevel][i], first_dofs[finelevel][i+1]);
+          if (r.Size() > 0)
+            fv.Range(r) = tmp.Row(i);
+        }
+    }
+    
+    virtual void RestrictInline (int finelevel, BaseVector & v) const override
+    {
+      auto fv = v.FV<double>();
+      Matrix<double> tmp(facets_on_level[finelevel], dofs_per_facet);
+      tmp = 0.0;
+      
+      for (size_t i = 0; i < first_dofs[finelevel].Size()-1; i++)
+        {
+          IntRange r(first_dofs[finelevel][i], first_dofs[finelevel][i+1]);
+          if (r.Size() > 0)
+            tmp.Row(i) = fv.Range(r);
+        }
+
+
+      if constexpr (DIM==2)
+        {
+          for (int i = facets_on_level[finelevel-1]; i < facets_on_level[finelevel]; i++)
+            if (auto parents = get<1>(ma->GetParentEdges(i)); parents[1] == -1)
+              {
+                int pe = parents[0];
+                tmp.Row(pe) += Trans(facetprol[prol_class[i]]) * tmp.Row(i);
+              }
+        }
+      else
+        {
+          // for (int i = facets_on_level[finelevel-1]; i < facets_on_level[finelevel]; i++)
+          for (int i = facets_on_level[finelevel]-1; i >= facets_on_level[finelevel-1]; i--)      
+            if (auto parents = get<1>(ma->GetParentFaces(i)); parents[1] == -1)
+              {
+                int pe = parents[0];
+                tmp.Row(pe) += Trans(facetprol[prol_class[i]]) * tmp.Row(i);
+                tmp.Row(i) = 0.0;  // optional, for testing
+              }
+        }
+      
+      fv = 0.0;
+      for (size_t i = 0; i < first_dofs[finelevel-1].Size()-1; i++)
+        {
+          IntRange r(first_dofs[finelevel-1][i], first_dofs[finelevel-1][i+1]);
+          if (r.Size() > 0)
+            fv.Range(r) = tmp.Row(i);
+        }
+
+    }
+    
+  };
+  
+  
+  
 }
 
 

netgen/include/python_ngstd.hpp

@@ -12,6 +12,7 @@
 #endif
 
 #include <pybind11/pybind11.h>
+#include <pybind11/iostream.h>
 #include <pybind11/eval.h>
 #include <pybind11/operators.h>
 #include <pybind11/complex.h>

netgen/include/sparsecholesky.hpp

@@ -34,7 +34,7 @@ namespace ngla
 
     void DoArchive(Archive& ar) override;
 
-    virtual bool IsComplex() const override { return matrix.lock()->IsComplex(); }
+    // virtual bool IsComplex() const override { return matrix.lock()->IsComplex(); }
 
     virtual void Smooth (BaseVector & u, const BaseVector & f, BaseVector & y) const;
 

netgen/include/sparsematrix_impl.hpp

@@ -51,28 +51,94 @@ namespace ngla
                  FlatArray<TM> val, size_t h, size_t w)
   {
     static Timer t("SparseMatrix::CreateFromCOO"); RegionTimer r(t);
-    Array<int> cnt(h);
+    static Timer t1("SparseMatrix::CreateFromCOO 1");
+    static Timer t2("SparseMatrix::CreateFromCOO 2");
+    static Timer t3("SparseMatrix::CreateFromCOO 3");
 
     /*
-    cnt = 0;
-    for (auto i : indi) cnt[i]++;
-    */
+    {
+    Array<int> cnt(h);
 
+    t1.Start();
     DynamicTable<int> tab(h);
     for (size_t i = 0; i < indi.Size(); i++)
       tab.AddUnique(indi[i], indj[i]);
+    t1.Stop();
     for (size_t i = 0; i < h; i++)
       cnt[i] = tab.EntrySize(i);
     
     auto matrix = make_shared<SparseMatrix<TM>> (cnt, w);
+    t2.Start();
     for (auto k : ngstd::Range(indi))
       matrix->CreatePosition(indi[k], indj[k]);
+    t2.Stop();
     matrix->SetZero();
 
+    t3.Start();
     for (auto k : ngstd::Range(indi))
       (*matrix)(indi[k], indj[k]) += val[k];
+    t3.Stop();
+    // return matrix;
+    }
+    */
+    
+    Array<int> cnt(h);
+    cnt = 0;
+    for (auto i : indi)
+      cnt[i]++;
+    
+    Table<int> tab(cnt);
+    cnt = 0;
+    
+    for (auto [i,j] : Zip(indi, indj))
+      tab[i][cnt[i]++] = j;
+
+    cnt = 0;
+    // for (int i = 0; i < tab.Size(); i++)
+    ParallelFor (tab.Size(), [&] (size_t i)
+      {
+        QuickSort (tab[i]);
+
+        int prev = -1;
+        for (auto j : tab[i])
+          {
+            if (j != prev) cnt[i]++;
+            prev = j;
+          }
+      });
+    
+    auto matrix = make_shared<SparseMatrix<TM>> (cnt, w);
+    t2.Start();
+    // for (auto k : ngstd::Range(indi))
+    // matrix->CreatePosition(indi[k], indj[k]);
 
-    return matrix;
+    cnt = 0;
+    for (int i = 0; i < tab.Size(); i++)
+      {
+        int prev = -1;
+        for (auto j : tab[i])
+          {
+            auto cols = matrix->GetRowIndices(i);
+            if (j != prev)
+              cols[cnt[i]++] = j;
+            prev = j;
+          }
+      }
+    
+    t2.Stop();
+    matrix->SetZero();
+
+    t3.Start();
+    /*
+    for (auto k : ngstd::Range(indi))
+      (*matrix)(indi[k], indj[k]) += val[k];
+    */
+    ParallelFor (indi.Size(), [&](size_t k)
+    {
+      AtomicAdd ( (*matrix)(indi[k], indj[k]),  val[k]);
+    });
+    t3.Stop();
+    return matrix;    
   }
   
 

netgen/include/tangentialfacetfespace.hpp

@@ -16,34 +16,22 @@ namespace ngcomp
   class NGS_DLL_HEADER TangentialFacetFESpace : public FESpace
   {
   protected:
-    /// Level
-    // int level;
-    /// Number of Facets
-    // int nfacets;
-    /// 
-    // int ncfacets;
-    ///
-    // int nel;
-
     Array<int> first_facet_dof;
     Array<int> first_inner_dof;  // for highest_order_dc
-    // int ndof_lo;
 
     int rel_order;
 
     Array<IVec<2> > order_facet;
     Array<bool> fine_facet;
 
-    // int ndof;
-    // Array<int> ndlevel;
     bool var_order;
     bool print;
 
     bool highest_order_dc;
     bool hide_highest_order_dc;
+    bool all_dofs_together;
     
   public:
-    ///
     TangentialFacetFESpace (shared_ptr<MeshAccess> ama, const Flags & flags, 
                             bool parseflags = false );
 
@@ -61,20 +49,11 @@ namespace ngcomp
     virtual void SetOrder (NodeId ni, int order) override;
     virtual int GetOrder (NodeId ni) const override;
 
-    
-    // virtual size_t GetNDof() const throw() override { return ndof; }
-    // virtual size_t GetNDofLevel ( int i ) const override { return ndlevel[i]; }
-
     virtual FlatArray<VorB> GetDualShapeNodes (VorB vb) const override;
 
-    // virtual int GetNDofLowOrder () const
-    // { return ndof_lo; }
 
     virtual FiniteElement & GetFE(ElementId ei, Allocator & lh) const override;
     
-    // virtual const FiniteElement & GetFE ( int elnr, LocalHeap & lh ) const;
-    // virtual const FiniteElement & GetSFE ( int selnr, LocalHeap & lh ) const;
-
     virtual void GetFacetDofNrs (int felnr, Array<DofId> & dnums) const;
 
     virtual int GetNFacetDofs (int felnr) const;
@@ -90,6 +69,12 @@ namespace ngcomp
     virtual IVec<2> GetFacetOrder(int fnr) const;
 
     virtual int GetFirstFacetDof(int fanr) const;
+    const auto & GetFirstFacetDof() const { return first_facet_dof; }
+
+    IntRange GetFacetDofs (int nr) const
+    { 
+      return IntRange (first_facet_dof[nr], first_facet_dof[nr+1]); 
+    }
 
     virtual bool UsesHighestOrderDiscontinuous() const {return highest_order_dc;};
 

netgen/include/tensor.hpp

@@ -371,7 +371,7 @@ public:
   
   ~Tensor() 
   {
-    delete this->Data();
+    delete [] this->Data();
   }
 
   FlatTensor<DIM,T> operator= (double d)

ngsolve/__init__.pyi

@@ -44,8 +44,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
 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 L2
@@ -217,7 +217,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.2501-83-gaf867070b'
+__version__: str = '6.2.2502-11-g4d6a78d51'
 ds: comp.DifferentialSymbol  # value = <ngsolve.comp.DifferentialSymbol object>
 dx: comp.DifferentialSymbol  # value = <ngsolve.comp.DifferentialSymbol object>
 ngsglobals: comp.GlobalVariables  # value = <ngsolve.comp.GlobalVariables object>

ngsolve/cmake/NGSolveConfig.cmake

@@ -1,4 +1,4 @@
-set(PACKAGE_VERSION "6.2.2501-83-gaf867070b")
+set(PACKAGE_VERSION "6.2.2502-11-g4d6a78d51")
 find_package(Netgen CONFIG REQUIRED HINTS
   ${CMAKE_CURRENT_LIST_DIR}
   ${CMAKE_CURRENT_LIST_DIR}/..

ngsolve/comp/__init__.pyi

@@ -11,7 +11,7 @@ import numpy
 import pyngcore.pyngcore
 import typing
 from . import pml
-__all__ = ['APhiHCurlAMG', 'Array_N6ngcomp13COUPLING_TYPEE_S', 'BBBND', 'BBND', 'BDDCPreconditioner', 'BDDCPreconditioner_complex', 'BDDCPreconditioner_double', 'BND', 'BilinearForm', 'BndElementId', 'BoundaryFromVolumeCF', 'COUPLING_TYPE', 'ComponentGridFunction', 'Compress', 'CompressCompound', 'ContactBoundary', 'ConvertOperator', 'DifferentialSymbol', 'Discontinuous', 'DualProxyFunction', 'ElementId', 'ElementRange', 'FESpace', 'FESpaceElement', 'FESpaceElementRange', 'FacetFESpace', 'FacetSurface', 'FlatArray_N6ngcomp13COUPLING_TYPEE_S', 'FromArchiveCF', 'FromArchiveFESpace', 'FromArchiveMesh', 'GlobalInterfaceSpace', 'GlobalSpace', 'GlobalVariables', 'GridFunction', 'GridFunctionC', 'GridFunctionCoefficientFunction', 'GridFunctionD', 'H1', 'H1LumpingFESpace', 'HCurl', 'HCurlAMG', 'HCurlCurl', 'HCurlDiv', 'HDiv', 'HDivDiv', 'HDivDivSurface', 'HDivSurface', 'Hidden', 'Integral', 'Integrate', 'IntegrationRuleSpace', 'IntegrationRuleSpaceSurface', 'Interpolate', 'InterpolateProxy', 'KSpaceCoeffs', 'L2', 'LinearForm', 'LocalPreconditioner', 'MatrixFreeOperator', 'MatrixValued', 'Mesh', 'MeshNode', 'MeshNodeRange', 'MultiGridPreconditioner', 'NGS_Object', 'Ngs_Element', 'NodalFESpace', 'NodeId', 'NodeRange', 'NormalFacetFESpace', 'NormalFacetSurface', 'NumberSpace', 'ORDER_POLICY', 'PatchwiseSolve', 'Periodic', 'PlateauFESpace', 'Preconditioner', 'ProductSpace', 'Prolongate', 'ProlongateCoefficientFunction', 'Prolongation', 'ProxyFunction', 'QuasiPeriodicC', 'QuasiPeriodicD', 'Region', 'RegisterPreconditioner', 'Reorder', 'SetHeapSize', 'SetTestoutFile', 'SumOfIntegrals', 'SurfaceL2', 'SymbolTable_D', 'SymbolTable_sp_D', 'SymbolTable_sp_N5ngfem19CoefficientFunctionE', 'SymbolTable_sp_N6ngcomp10LinearFormE', 'SymbolTable_sp_N6ngcomp12BilinearFormE', 'SymbolTable_sp_N6ngcomp12GridFunctionE', 'SymbolTable_sp_N6ngcomp14PreconditionerE', 'SymbolTable_sp_N6ngcomp7FESpaceE', 'SymbolicBFI', 'SymbolicEnergy', 'SymbolicLFI', 'SymbolicTPBFI', 'TangentialFacetFESpace', 'TangentialSurfaceL2', 'TensorProductFESpace', 'TensorProductIntegrate', 'ToArchive', 'Transfer2StdMesh', 'VOL', 'VTKOutput', 'Variation', 'VectorFacetFESpace', 'VectorFacetSurface', 'VectorH1', 'VectorL2', 'VectorNodalFESpace', 'VectorSurfaceL2', 'VectorValued', 'VorB', 'ngsglobals', 'pml']
+__all__ = ['APhiHCurlAMG', 'AddCurrentDensity', 'Array_N6ngcomp13COUPLING_TYPEE_S', 'BBBND', 'BBND', 'BDDCPreconditioner', 'BDDCPreconditioner_complex', 'BDDCPreconditioner_double', 'BND', 'BilinearForm', 'BndElementId', 'BoundaryFromVolumeCF', 'COUPLING_TYPE', 'ComponentGridFunction', 'Compress', 'CompressCompound', 'ContactBoundary', 'ConvertOperator', 'DifferentialSymbol', 'Discontinuous', 'DualProxyFunction', 'ElementId', 'ElementRange', 'FESpace', 'FESpaceElement', 'FESpaceElementRange', 'FacetFESpace', 'FacetSurface', 'FlatArray_N6ngcomp13COUPLING_TYPEE_S', 'FromArchiveCF', 'FromArchiveFESpace', 'FromArchiveMesh', 'GlobalInterfaceSpace', 'GlobalSpace', 'GlobalVariables', 'GridFunction', 'GridFunctionC', 'GridFunctionCoefficientFunction', 'GridFunctionD', 'H1', 'H1LumpingFESpace', 'HCurl', 'HCurlAMG', 'HCurlCurl', 'HCurlDiv', 'HDiv', 'HDivDiv', 'HDivDivSurface', 'HDivSurface', 'Hidden', 'Integral', 'Integrate', 'IntegrationRuleSpace', 'IntegrationRuleSpaceSurface', 'Interpolate', 'InterpolateProxy', 'KSpaceCoeffs', 'L2', 'LinearForm', 'LocalPreconditioner', 'MatrixFreeOperator', 'MatrixValued', 'Mesh', 'MeshNode', 'MeshNodeRange', 'MultiGridPreconditioner', 'NGS_Object', 'Ngs_Element', 'NodalFESpace', 'NodeId', 'NodeRange', 'NormalFacetFESpace', 'NormalFacetSurface', 'NumberSpace', 'ORDER_POLICY', 'PatchwiseSolve', 'Periodic', 'PlateauFESpace', 'Preconditioner', 'ProductSpace', 'Prolongate', 'ProlongateCoefficientFunction', 'Prolongation', 'ProxyFunction', 'QuasiPeriodicC', 'QuasiPeriodicD', 'Region', 'RegisterPreconditioner', 'Reorder', 'SetHeapSize', 'SetTestoutFile', 'SumOfIntegrals', 'SurfaceL2', 'SymbolTable_D', 'SymbolTable_sp_D', 'SymbolTable_sp_N5ngfem19CoefficientFunctionE', 'SymbolTable_sp_N6ngcomp10LinearFormE', 'SymbolTable_sp_N6ngcomp12BilinearFormE', 'SymbolTable_sp_N6ngcomp12GridFunctionE', 'SymbolTable_sp_N6ngcomp14PreconditionerE', 'SymbolTable_sp_N6ngcomp7FESpaceE', 'SymbolicBFI', 'SymbolicEnergy', 'SymbolicLFI', 'SymbolicTPBFI', 'TangentialFacetFESpace', 'TangentialSurfaceL2', 'TensorProductFESpace', 'TensorProductIntegrate', 'ToArchive', 'Transfer2StdMesh', 'VOL', 'VTKOutput', 'Variation', 'VectorFacetFESpace', 'VectorFacetSurface', 'VectorH1', 'VectorL2', 'VectorNodalFESpace', 'VectorSurfaceL2', 'VectorValued', 'VorB', 'ngsglobals', 'pml']
 class APhiHCurlAMG(HCurlAMG):
     """
     
@@ -938,6 +938,10 @@ class FESpace(NGS_Object):
         region : ngsolve.comp.Region
           input region
         """
+    def SetHarmonicProlongation(self, bf: ..., inverse: str = 'sparsecholesky') -> None:
+        """
+        use harmonic prolongation w.r.t. bilinear-form
+        """
     @typing.overload
     def SetOrder(self, element_type: ngsolve.fem.ET, order: int) -> None:
         """
@@ -3103,6 +3107,10 @@ class MultiGridPreconditioner(Preconditioner):
         'block': Block smoother
     coarsetype: string = direct
       How to solve coarse problem.
+    cycle: int = 1
+      multigrid cycle (0 only smoothing, 1..V-cycle, 2..W-cycle.
+    smoothingsteps: int = 1
+      number of (pre and post-)smoothing steps
     coarsesmoothingsteps: int = 1
       If coarsetype is smoothing, then how many smoothingsteps will be done.
     updatealways: bool = False
@@ -5082,6 +5090,8 @@ class VorB:
     @property
     def value(self) -> int:
         ...
+def AddCurrentDensity(arg0: ngsolve.fem.SingularMLMP3, arg1: ngsolve.fem.CoefficientFunction, arg2: Region) -> None:
+    ...
 def BndElementId(nr: int) -> ElementId:
     """
     Creates an element-id for a boundary element

ngsolve/config/__init__.pyi

@@ -24,14 +24,14 @@ NGSOLVE_INSTALL_DIR_INCLUDE: str = 'netgen/include'
 NGSOLVE_INSTALL_DIR_LIB: str = 'netgen'
 NGSOLVE_INSTALL_DIR_PYTHON: str = '.'
 NGSOLVE_INSTALL_DIR_RES: str = 'share'
-NGSOLVE_VERSION: str = '6.2.2501-83-gaf867070b'
-NGSOLVE_VERSION_GIT: str = 'v6.2.2501-83-gaf867070b'
-NGSOLVE_VERSION_HASH: str = 'gaf867070b'
+NGSOLVE_VERSION: str = '6.2.2502-11-g4d6a78d51'
+NGSOLVE_VERSION_GIT: str = 'v6.2.2502-11-g4d6a78d51'
+NGSOLVE_VERSION_HASH: str = 'g4d6a78d51'
 NGSOLVE_VERSION_MAJOR: str = '6'
 NGSOLVE_VERSION_MINOR: str = '2'
-NGSOLVE_VERSION_PATCH: str = '2501'
-NGSOLVE_VERSION_PYTHON: str = '6.2.2501.post83.dev1'
-NGSOLVE_VERSION_TWEAK: str = '83'
+NGSOLVE_VERSION_PATCH: str = '2502'
+NGSOLVE_VERSION_PYTHON: str = '6.2.2502.post11.dev1'
+NGSOLVE_VERSION_TWEAK: str = '11'
 USE_CCACHE: bool = True
 USE_HYPRE: bool = False
 USE_LAPACK: bool = True
@@ -40,4 +40,4 @@ USE_MUMPS: bool = False
 USE_PARDISO: bool = False
 USE_UMFPACK: bool = True
 is_python_package: bool = True
-version: str = 'v6.2.2501-83-gaf867070b'
+version: str = 'v6.2.2502-11-g4d6a78d51'

ngsolve/config/config.py

@@ -29,15 +29,15 @@ NGSOLVE_INSTALL_DIR_INCLUDE  = "netgen/include"
 NGSOLVE_INSTALL_DIR_CMAKE    = "ngsolve/cmake"
 NGSOLVE_INSTALL_DIR_RES      = "share"
 
-NGSOLVE_VERSION = "6.2.2501-83-gaf867070b"
-NGSOLVE_VERSION_GIT = "v6.2.2501-83-gaf867070b"
-NGSOLVE_VERSION_PYTHON = "6.2.2501.post83.dev1"
+NGSOLVE_VERSION = "6.2.2502-11-g4d6a78d51"
+NGSOLVE_VERSION_GIT = "v6.2.2502-11-g4d6a78d51"
+NGSOLVE_VERSION_PYTHON = "6.2.2502.post11.dev1"
 
 NGSOLVE_VERSION_MAJOR = "6"
 NGSOLVE_VERSION_MINOR = "2"
-NGSOLVE_VERSION_TWEAK = "83"
-NGSOLVE_VERSION_PATCH = "2501"
-NGSOLVE_VERSION_HASH = "gaf867070b"
+NGSOLVE_VERSION_TWEAK = "11"
+NGSOLVE_VERSION_PATCH = "2502"
+NGSOLVE_VERSION_HASH = "g4d6a78d51"
 
 CMAKE_CXX_COMPILER           = "/Library/Developer/CommandLineTools/usr/bin/c++"
 CMAKE_CUDA_COMPILER          = ""

ngsolve/config/config.pyi

@@ -26,14 +26,14 @@ NGSOLVE_INSTALL_DIR_INCLUDE: str = 'netgen/include'
 NGSOLVE_INSTALL_DIR_LIB: str = 'netgen'
 NGSOLVE_INSTALL_DIR_PYTHON: str = '.'
 NGSOLVE_INSTALL_DIR_RES: str = 'share'
-NGSOLVE_VERSION: str = '6.2.2501-83-gaf867070b'
-NGSOLVE_VERSION_GIT: str = 'v6.2.2501-83-gaf867070b'
-NGSOLVE_VERSION_HASH: str = 'gaf867070b'
+NGSOLVE_VERSION: str = '6.2.2502-11-g4d6a78d51'
+NGSOLVE_VERSION_GIT: str = 'v6.2.2502-11-g4d6a78d51'
+NGSOLVE_VERSION_HASH: str = 'g4d6a78d51'
 NGSOLVE_VERSION_MAJOR: str = '6'
 NGSOLVE_VERSION_MINOR: str = '2'
-NGSOLVE_VERSION_PATCH: str = '2501'
-NGSOLVE_VERSION_PYTHON: str = '6.2.2501.post83.dev1'
-NGSOLVE_VERSION_TWEAK: str = '83'
+NGSOLVE_VERSION_PATCH: str = '2502'
+NGSOLVE_VERSION_PYTHON: str = '6.2.2502.post11.dev1'
+NGSOLVE_VERSION_TWEAK: str = '11'
 USE_CCACHE: bool = True
 USE_HYPRE: bool = False
 USE_LAPACK: bool = True
@@ -42,4 +42,4 @@ USE_MUMPS: bool = False
 USE_PARDISO: bool = False
 USE_UMFPACK: bool = True
 is_python_package: bool = True
-version: str = 'v6.2.2501-83-gaf867070b'
+version: str = 'v6.2.2502-11-g4d6a78d51'

ngsolve/directsolvers.pyi

@@ -9,7 +9,7 @@ import typing
 __all__ = ['BaseMatrix', 'BaseVector', 'BilinearForm', 'BitArray', 'SuperLU']
 class SuperLU(ngsolve.la.BaseMatrix):
     __firstlineno__: typing.ClassVar[int] = 3
-    __static_attributes__: typing.ClassVar[tuple] = ('lu', 'fd', 'a', 'freedofs')
+    __static_attributes__: typing.ClassVar[tuple] = ('freedofs', 'lu', 'fd', 'a')
     def Mult(self, x: ngsolve.la.BaseVector, y: ngsolve.la.BaseVector):
         ...
     def Update(self):

ngsolve/fem.pyi

@@ -8,7 +8,7 @@ import numpy
 import os
 import pyngcore.pyngcore
 import typing
-__all__ = ['BFI', 'BSpline', 'BSpline2D', 'BaseMappedIntegrationPoint', 'BlockBFI', 'BlockLFI', 'CELL', 'CacheCF', 'CoefficientFunction', 'Cof', 'CompilePythonModule', 'CompoundBFI', 'CompoundLFI', 'Conj', 'CoordCF', 'CoordinateTrafo', 'Cross', 'Det', 'DifferentialOperator', 'EDGE', 'ELEMENT', 'ET', 'Einsum', 'ElementTopology', 'ElementTransformation', 'FACE', 'FACET', 'FiniteElement', 'GenerateL2ElementCode', 'H1FE', 'HCurlFE', 'HDivDivFE', 'HDivFE', 'HEX', 'Id', 'IfPos', 'IntegrationPoint', 'IntegrationRule', 'Inv', 'L2FE', 'LFI', 'LeviCivitaSymbol', 'LoggingCF', 'MeshPoint', 'MinimizationCF', 'MixedFE', 'NODE_TYPE', 'NewtonCF', 'POINT', 'PRISM', 'PYRAMID', 'Parameter', 'ParameterC', 'PlaceholderCF', 'PointEvaluationFunctional', 'QUAD', 'RegularMLMP', 'RegularMLMultiPoleCF', 'RegularMultiPoleCF', 'SEGM', 'ScalarFE', 'SetPMLParameters', 'SingularMLMP', 'SingularMLMultiPoleCF', 'SingularMultiPoleCF', 'Skew', 'SpecialCFCreator', 'SphericalHarmonicsCF', 'Sphericalharmonics', 'Sym', 'TET', 'TRIG', 'Trace', 'VERTEX', 'VoxelCoefficient', 'Zero', 'acos', 'asin', 'atan', 'atan2', 'ceil', 'cos', 'cosh', 'erf', 'exp', 'floor', 'log', 'pow', 'sin', 'sinh', 'specialcf', 'sqrt', 'tan']
+__all__ = ['BFI', 'BSpline', 'BSpline2D', 'BaseMappedIntegrationPoint', 'BiotSavartCF', 'BiotSavartRegularMLCF', 'BiotSavartSingularMLCF', 'BlockBFI', 'BlockLFI', 'CELL', 'CacheCF', 'CoefficientFunction', 'Cof', 'CompilePythonModule', 'CompoundBFI', 'CompoundLFI', 'Conj', 'CoordCF', 'CoordinateTrafo', 'Cross', 'Det', 'DifferentialOperator', 'EDGE', 'ELEMENT', 'ET', 'Einsum', 'ElementTopology', 'ElementTransformation', 'FACE', 'FACET', 'FiniteElement', 'GenerateL2ElementCode', 'H1FE', 'HCurlFE', 'HDivDivFE', 'HDivFE', 'HEX', 'Id', 'IfPos', 'IntegrationPoint', 'IntegrationRule', 'Inv', 'L2FE', 'LFI', 'LeviCivitaSymbol', 'LoggingCF', 'MeshPoint', 'MinimizationCF', 'MixedFE', 'NODE_TYPE', 'NewtonCF', 'POINT', 'PRISM', 'PYRAMID', 'Parameter', 'ParameterC', 'PlaceholderCF', 'PointEvaluationFunctional', 'QUAD', 'RegularMLMP', 'RegularMLMultiPoleCF', 'RegularMultiPoleCF', 'SEGM', 'ScalarFE', 'SetPMLParameters', 'SingularMLMP', 'SingularMLMP3', 'SingularMLMultiPoleCF', 'SingularMultiPoleCF', 'Skew', 'SpecialCFCreator', 'SphericalHarmonicsCF', 'Sphericalharmonics', 'Sym', 'TET', 'TRIG', 'Trace', 'VERTEX', 'VoxelCoefficient', 'Zero', 'acos', 'asin', 'atan', 'atan2', 'ceil', 'cos', 'cosh', 'erf', 'exp', 'floor', 'log', 'pow', 'sin', 'sinh', 'specialcf', 'sqrt', 'tan']
 class BFI:
     """
     
@@ -242,6 +242,25 @@ class BaseMappedIntegrationPoint:
         """
         Transformation of the mapped integration point
         """
+class BiotSavartCF(CoefficientFunction):
+    def AddCurrent(self, sp: ngsolve.bla.Vec3D, ep: ngsolve.bla.Vec3D, j: complex, num: int = 100) -> None:
+        ...
+    def __init__(self, order: int, kappa: float, center: ngsolve.bla.Vec3D, scale: float = 1.0) -> None:
+        ...
+class BiotSavartRegularMLCF(CoefficientFunction):
+    def __init__(self, mp: BiotSavartSingularMLCF, center: ngsolve.bla.Vec3D, r: float, order: int = -1) -> None:
+        ...
+    @property
+    def mlmp(self) -> ...:
+        ...
+class BiotSavartSingularMLCF(CoefficientFunction):
+    def CreateRegularExpansion(self, center: ngsolve.bla.Vec3D, r: float) -> ...:
+        ...
+    def __init__(self, center: ngsolve.bla.Vec3D, r: float, order: int = -1, kappa: float) -> None:
+        ...
+    @property
+    def mlmp(self) -> SingularMLMP3:
+        ...
 class CoefficientFunction:
     """
     A CoefficientFunction (CF) is some function defined on a mesh.
@@ -979,6 +998,8 @@ class PointEvaluationFunctional:
 class RegularMLMP:
     def Norm(self) -> float:
         ...
+    def Print(self) -> None:
+        ...
     def __str__(self) -> str:
         ...
 class RegularMLMultiPoleCF(CoefficientFunction):
@@ -1060,11 +1081,22 @@ class SingularMLMP:
         ...
     def Norm(self) -> float:
         ...
+    def Print(self) -> None:
+        ...
     def __init__(self, arg0: ngsolve.bla.Vec3D, arg1: float, arg2: int, arg3: float) -> None:
         ...
     def __str__(self) -> str:
         ...
+class SingularMLMP3:
+    def AddCurrent(self, sp: ngsolve.bla.Vec3D, ep: ngsolve.bla.Vec3D, j: complex, num: int = 100) -> None:
+        ...
+    def Calc(self) -> None:
+        ...
+    def __init__(self, arg0: ngsolve.bla.Vec3D, arg1: float, arg2: int, arg3: float) -> None:
+        ...
 class SingularMLMultiPoleCF(CoefficientFunction):
+    def CreateRegularExpansion(self, center: ngsolve.bla.Vec3D, r: float) -> ...:
+        ...
     def __init__(self, center: ngsolve.bla.Vec3D, r: float, order: int = -1, kappa: float) -> None:
         ...
     @property