ngsolve 6.2.2506.post85.dev0__cp313-cp313-manylinux_2_28_x86_64.whl → 6.2.2506.post98.dev0__cp313-cp313-manylinux_2_28_x86_64.whl

{ngsolve-6.2.2506.post85.dev0.data → ngsolve-6.2.2506.post98.dev0.data}/data/bin/ngs_nvcc

@@ -11,4 +11,4 @@ if [ -f /usr/local/cuda/bin/nvcc ]
   else NVCC=nvcc
 fi
 
-$NVCC  -DHAVE_NETGEN_SOURCES -DHAVE_DLFCN_H -DHAVE_CXA_DEMANGLE -DUSE_TIMEOFDAY -DTCL -DLAPACK -DUSE_PARDISO -DNGS_PYTHON -D_GLIBCXX_USE_CXX11_ABI=1 -DNETGEN_PYTHON -DNG_PYTHON -DPYBIND11_SIMPLE_GIL_MANAGEMENT -DPARALLEL -DNG_MPI_WRAPPER -arch=native  -Xcompiler -std=c++17 -Xcompiler -DMAX_SYS_DIM=3 -Xcompiler -march=core-avx2 -Xcompiler -fabi-version=17 -Xcompiler -fPIC  -I/opt/_internal/cpython-3.13.8/include -I/opt/_internal/cpython-3.13.8/include/python3.13 -I$NGSCXX_DIR/../include/netgen -I$NGSCXX_DIR/../include/netgen/include $PY_INCLUDE_FLAGS $*
+$NVCC  -DHAVE_NETGEN_SOURCES -DHAVE_DLFCN_H -DHAVE_CXA_DEMANGLE -DUSE_TIMEOFDAY -DTCL -DLAPACK -DUSE_PARDISO -DNGS_PYTHON -D_GLIBCXX_USE_CXX11_ABI=1 -DNETGEN_PYTHON -DNG_PYTHON -DPYBIND11_SIMPLE_GIL_MANAGEMENT -DPARALLEL -DNG_MPI_WRAPPER -arch=native  -Xcompiler -std=c++17 -Xcompiler -DMAX_SYS_DIM=3 -Xcompiler -march=core-avx2 -Xcompiler -fabi-version=17 -Xcompiler -fPIC  -I/opt/_internal/cpython-3.13.8/include -I/opt/_internal/cpython-3.13.9/include/python3.13 -I$NGSCXX_DIR/../include/netgen -I$NGSCXX_DIR/../include/netgen/include $PY_INCLUDE_FLAGS $*

{ngsolve-6.2.2506.post85.dev0.data → ngsolve-6.2.2506.post98.dev0.data}/data/bin/ngscxx

@@ -10,6 +10,6 @@ export CCACHE_NOHASHDIR=1
 
 
 if [ -f /usr/bin/ccache ]
-  then /usr/bin/ccache $NGSCXX -O3 -DNDEBUG   -DHAVE_NETGEN_SOURCES -DHAVE_DLFCN_H -DHAVE_CXA_DEMANGLE -DUSE_TIMEOFDAY -DTCL -DLAPACK -DUSE_PARDISO -DNGS_PYTHON -D_GLIBCXX_USE_CXX11_ABI=1 -DNETGEN_PYTHON -DNG_PYTHON -DPYBIND11_SIMPLE_GIL_MANAGEMENT -DPARALLEL -DNG_MPI_WRAPPER -std=c++17 -DMAX_SYS_DIM=3 -march=core-avx2 -fabi-version=17 -fPIC  -I/opt/_internal/cpython-3.13.8/include -I/opt/_internal/cpython-3.13.8/include/python3.13 -I$NGSCXX_DIR/../include/netgen -I$NGSCXX_DIR/../include/netgen/include $PY_INCLUDE_FLAGS $*
-  else $NGSCXX -O3 -DNDEBUG   -DHAVE_NETGEN_SOURCES -DHAVE_DLFCN_H -DHAVE_CXA_DEMANGLE -DUSE_TIMEOFDAY -DTCL -DLAPACK -DUSE_PARDISO -DNGS_PYTHON -D_GLIBCXX_USE_CXX11_ABI=1 -DNETGEN_PYTHON -DNG_PYTHON -DPYBIND11_SIMPLE_GIL_MANAGEMENT -DPARALLEL -DNG_MPI_WRAPPER -std=c++17 -DMAX_SYS_DIM=3 -march=core-avx2 -fabi-version=17 -fPIC  -I/opt/_internal/cpython-3.13.8/include -I/opt/_internal/cpython-3.13.8/include/python3.13 -I$NGSCXX_DIR/../include/netgen -I$NGSCXX_DIR/../include/netgen/include $PY_INCLUDE_FLAGS $*
+  then /usr/bin/ccache $NGSCXX -O3 -DNDEBUG   -DHAVE_NETGEN_SOURCES -DHAVE_DLFCN_H -DHAVE_CXA_DEMANGLE -DUSE_TIMEOFDAY -DTCL -DLAPACK -DUSE_PARDISO -DNGS_PYTHON -D_GLIBCXX_USE_CXX11_ABI=1 -DNETGEN_PYTHON -DNG_PYTHON -DPYBIND11_SIMPLE_GIL_MANAGEMENT -DPARALLEL -DNG_MPI_WRAPPER -std=c++17 -DMAX_SYS_DIM=3 -march=core-avx2 -fabi-version=17 -fPIC  -I/opt/_internal/cpython-3.13.8/include -I/opt/_internal/cpython-3.13.9/include/python3.13 -I$NGSCXX_DIR/../include/netgen -I$NGSCXX_DIR/../include/netgen/include $PY_INCLUDE_FLAGS $*
+  else $NGSCXX -O3 -DNDEBUG   -DHAVE_NETGEN_SOURCES -DHAVE_DLFCN_H -DHAVE_CXA_DEMANGLE -DUSE_TIMEOFDAY -DTCL -DLAPACK -DUSE_PARDISO -DNGS_PYTHON -D_GLIBCXX_USE_CXX11_ABI=1 -DNETGEN_PYTHON -DNG_PYTHON -DPYBIND11_SIMPLE_GIL_MANAGEMENT -DPARALLEL -DNG_MPI_WRAPPER -std=c++17 -DMAX_SYS_DIM=3 -march=core-avx2 -fabi-version=17 -fPIC  -I/opt/_internal/cpython-3.13.8/include -I/opt/_internal/cpython-3.13.9/include/python3.13 -I$NGSCXX_DIR/../include/netgen -I$NGSCXX_DIR/../include/netgen/include $PY_INCLUDE_FLAGS $*
 fi

ngsolve-6.2.2506.post98.dev0.data/data/bin/ngspy

@@ -0,0 +1,2 @@
+#! /usr/bin/bash
+LD_PRELOAD=$LD_PRELOAD: /opt/_internal/cpython-3.13.9/bin/python $*

{ngsolve-6.2.2506.post85.dev0.data → ngsolve-6.2.2506.post98.dev0.data}/data/include/netgen/bdbequations.hpp

@@ -2783,9 +2783,11 @@ namespace ngfem
   extern template class NGS_DLL_HEADER T_DifferentialOperator<DiffOpGradVectorH1<2> >;
   extern template class NGS_DLL_HEADER T_DifferentialOperator<DiffOpGradVectorH1<3> >;
 
+  extern template class NGS_DLL_HEADER T_DifferentialOperator<DiffOpGradientBoundary<1> >;
   extern template class NGS_DLL_HEADER T_DifferentialOperator<DiffOpGradientBoundary<2> >;
   extern template class NGS_DLL_HEADER T_DifferentialOperator<DiffOpGradientBoundary<3> >;
 
+  extern template class NGS_DLL_HEADER T_DifferentialOperator<DiffOpGradBoundaryVectorH1<1> >;
   extern template class NGS_DLL_HEADER T_DifferentialOperator<DiffOpGradBoundaryVectorH1<2> >;
   extern template class NGS_DLL_HEADER T_DifferentialOperator<DiffOpGradBoundaryVectorH1<3> >;
 

{ngsolve-6.2.2506.post85.dev0.data → ngsolve-6.2.2506.post98.dev0.data}/data/include/netgen/gridfunction.hpp

@@ -157,7 +157,7 @@ namespace ngcomp
     int GetMultiDim () const { return multidim; }
 
     /// increase multidim and copy vec to new component
-    void AddMultiDimComponent (BaseVector & vec);
+    void AddMultiDimComponent (BaseVector & vec, bool copy = true);
   
     int GetLevelUpdated() const { return level_updated; }
     ///

{ngsolve-6.2.2506.post85.dev0.data → ngsolve-6.2.2506.post98.dev0.data}/data/include/netgen/kernels.hpp

@@ -203,22 +203,28 @@ namespace ngsbem
 
   /** HelmholtzSLkernel is the kernel for the double layer potential of the 
       Helmholtz equation $ -\Delta u - \kappa^2 u = 0, \; \kappa>0\,. $ */
-  template <int DIM> class HelmholtzSLKernel;
+  template <int DIM, int COMPS=1> class HelmholtzSLKernel;
 
   /** HelmholtzSLkernel in 3D reads 
       $$ G(x-y) = \frac{1 }{4\,\pi} \,\frac{e^{i\,\kappa \, |x-y| }{|x-y|} \, 
           \quad x, y \in \mathbb R^3, \; x\not=y\,. $$ */
-  template<>
-  class HelmholtzSLKernel<3> : public BaseKernel
+  template<int COMPS>
+  class HelmholtzSLKernel<3, COMPS> : public BaseKernel
   {
     double kappa;
   public:
     typedef Complex value_type;
+    using mp_type = typename std::conditional<COMPS == 1,
+                                              Complex,
+                                              Vec<COMPS, Complex>>::type;
     static string Name() { return "HelmholtzSL"; }
-    static auto Shape() { return IVec<2>(1,1); }
+    static auto Shape() { return IVec<2>(COMPS,COMPS); }
     
     /** Construction of the kernel specifies the wavenumber $\kappa$. */
-    HelmholtzSLKernel (double _kappa) : kappa(_kappa) { }
+    HelmholtzSLKernel (double _kappa) : kappa(_kappa) {
+      for (size_t i = 0; i < COMPS; i++)
+        terms += {1.0, 0, i, i};
+    }
 
     template <typename T>
     auto Evaluate (Vec<3,T> x, Vec<3,T> y, Vec<3,T> nx, Vec<3,T> ny) const
@@ -229,26 +235,32 @@ namespace ngsbem
       return Vec<1,decltype(kern)> (kern);
     }
     double GetKappa() const { return kappa; }
-    Array<KernelTerm> terms = { KernelTerm{1.0, 0, 0, 0}, };    
+    Array<KernelTerm> terms;
 
     auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<SingularMLExpansion<Complex>> (c, r, kappa, fmm_params);
+      return make_shared<SingularMLExpansion<mp_type>> (c, r, kappa, fmm_params);
     }
 
     auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
     {
-      return make_shared<RegularMLExpansion<Complex>> (c, r, kappa, fmm_params);
+      return make_shared<RegularMLExpansion<mp_type>> (c, r, kappa, fmm_params);
     }
 
-    void AddSource (SingularMLExpansion<Complex> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
+    void AddSource (SingularMLExpansion<mp_type> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
     {
-      mp.AddCharge(pnt, val(0));
+      if constexpr (COMPS == 1)
+        mp.AddCharge (pnt, val(0));
+      else
+        mp.AddCharge (pnt, val);
     }
 
-    void EvaluateMP (RegularMLExpansion<Complex> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
+    void EvaluateMP (RegularMLExpansion<mp_type> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
     {
-      val(0) = mp.Evaluate (pnt);
+      if constexpr (COMPS == 1)
+        val(0) = mp.Evaluate (pnt);
+      else
+        val = mp.Evaluate (pnt);
     }
   };
 
@@ -317,56 +329,6 @@ namespace ngsbem
   };
 
 
-  template <int DIM> class HelmholtzSLVecKernel;
-
-  template<>
-  class HelmholtzSLVecKernel<3> : public BaseKernel
-  {
-    double kappa;
-  public:
-    typedef Complex value_type;
-    static string Name() { return "HelmholtzSLVec"; }
-    static auto Shape() { return IVec<2>(3,3); }
-
-    HelmholtzSLVecKernel (double _kappa) : kappa(_kappa) { }
-
-    template <typename T>
-    auto Evaluate (Vec<3,T> x, Vec<3,T> y, Vec<3,T> nx, Vec<3,T> ny) const
-    {
-      T norm = L2Norm(x-y);
-      auto kern = exp(Complex(0,kappa)*norm) / (4 * M_PI * norm);
-      return Vec<1,decltype(kern)> (kern);
-    }
-    double GetKappa() const { return kappa; }
-    Array<KernelTerm> terms =
-      {
-        KernelTerm{1.0, 0, 0, 0},
-        KernelTerm{1.0, 0, 1, 1},
-        KernelTerm{1.0, 0, 2, 2},
-      };
-
-    auto CreateMultipoleExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
-    {
-      return make_shared<SingularMLExpansion<Vec<3,Complex>>> (c, r, kappa, fmm_params);
-    }
-
-    auto CreateLocalExpansion (Vec<3> c, double r, FMM_Parameters fmm_params) const
-    {
-      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
-    {
-      mp.AddCharge(pnt, val);
-    }
-
-    void EvaluateMP (RegularMLExpansion<Vec<3,Complex>> & mp, Vec<3> pnt, Vec<3> nv, BareSliceVector<Complex> val) const
-    {
-      val = mp.Evaluate (pnt);
-    }
-  };
-
-
   template <int DIM> class HelmholtzHSKernel;
   
   template<>

{ngsolve-6.2.2506.post85.dev0.data → ngsolve-6.2.2506.post98.dev0.data}/data/include/netgen/meshaccess.hpp

@@ -1201,7 +1201,7 @@ namespace ngcomp
     shared_ptr<BitArray> MaskPtr() { return mask; }
     const shared_ptr<MeshAccess> & Mesh() const { return mesh; }
     void DoArchive(Archive& ar)
-    { ar & mesh & vb & mask; }
+    { ar.Shallow(mesh) & vb & mask; }
     Region operator+ (const Region & r2) const
     {
       return Region (mesh, vb, BitArray(*mask).Or(r2.Mask()));

{ngsolve-6.2.2506.post85.dev0.data → ngsolve-6.2.2506.post98.dev0.data}/data/lib/cmake/ngsolve/NGSolveConfig.cmake

@@ -1,4 +1,4 @@
-set(PACKAGE_VERSION "6.2.2506-85-g5e4fc20dd")
+set(PACKAGE_VERSION "6.2.2506-98-gf95a1c05a")
 find_package(Netgen CONFIG REQUIRED HINTS
   ${CMAKE_CURRENT_LIST_DIR}
   ${CMAKE_CURRENT_LIST_DIR}/..
@@ -20,7 +20,7 @@ set(NGSOLVE_CMAKE_BUILD_TYPE "Release")
 set(NGSOLVE_CMAKE_THREAD_LIBS_INIT "")
 set(NGSOLVE_MKL_LIBRARIES "MKL::MKL")
 set(NGSOLVE_PYBIND_INCLUDE_DIR "")
-set(NGSOLVE_PYTHON_INCLUDE_DIRS "/opt/_internal/cpython-3.13.8/include/python3.13")
+set(NGSOLVE_PYTHON_INCLUDE_DIRS "/opt/_internal/cpython-3.13.9/include/python3.13")
 set(NGSOLVE_PYTHON_LIBRARIES "Python3_LIBRARY-NOTFOUND")
 set(NGSOLVE_PYTHON_PACKAGES_INSTALL_DIR  "")
 set(NGSOLVE_TCL_INCLUDE_PATH "/usr/include")
@@ -51,7 +51,7 @@ set(NGSOLVE_MAX_SYS_DIM       3)
 
 set(NGSOLVE_COMPILE_FLAGS " -DHAVE_NETGEN_SOURCES -DHAVE_DLFCN_H -DHAVE_CXA_DEMANGLE -DUSE_TIMEOFDAY -DTCL -DLAPACK -DUSE_PARDISO -DNGS_PYTHON -D_GLIBCXX_USE_CXX11_ABI=1 -DNETGEN_PYTHON -DNG_PYTHON -DPYBIND11_SIMPLE_GIL_MANAGEMENT -DPARALLEL -DNG_MPI_WRAPPER" CACHE STRING "Preprocessor definitions of ngscxx")
 set(NGSOLVE_LINK_FLAGS " -L$NGSCXX_DIR/../lib/python3.13/site-packages/netgen_mesher.libs -Wl,--rpath=$NGSCXX_DIR/../lib/python3.13/site-packages/netgen_mesher.libs" CACHE STRING "Link flags set in ngsld")
-set(NGSOLVE_INCLUDE_DIRS /opt/_internal/cpython-3.13.8/include;/opt/_internal/cpython-3.13.8/include/python3.13 CACHE STRING "Include dirs set in ngscxx")
+set(NGSOLVE_INCLUDE_DIRS /opt/_internal/cpython-3.13.8/include;/opt/_internal/cpython-3.13.9/include/python3.13 CACHE STRING "Include dirs set in ngscxx")
 
 set(NGSOLVE_INSTALL_DIR_PYTHON  lib/python3.13/site-packages)
 set(NGSOLVE_INSTALL_DIR_BIN     bin)

{ngsolve-6.2.2506.post85.dev0.data → ngsolve-6.2.2506.post98.dev0.data}/data/lib/python3.13/site-packages/ngsolve/config/config.py

@@ -29,15 +29,15 @@ NGSOLVE_INSTALL_DIR_INCLUDE  = "include/netgen"
 NGSOLVE_INSTALL_DIR_CMAKE    = "lib/cmake/ngsolve"
 NGSOLVE_INSTALL_DIR_RES      = "share"
 
-NGSOLVE_VERSION = "6.2.2506-85-g5e4fc20dd"
-NGSOLVE_VERSION_GIT = "v6.2.2506-85-g5e4fc20dd"
-NGSOLVE_VERSION_PYTHON = "6.2.2506.post85.dev0"
+NGSOLVE_VERSION = "6.2.2506-98-gf95a1c05a"
+NGSOLVE_VERSION_GIT = "v6.2.2506-98-gf95a1c05a"
+NGSOLVE_VERSION_PYTHON = "6.2.2506.post98.dev0"
 
 NGSOLVE_VERSION_MAJOR = "6"
 NGSOLVE_VERSION_MINOR = "2"
-NGSOLVE_VERSION_TWEAK = "85"
+NGSOLVE_VERSION_TWEAK = "98"
 NGSOLVE_VERSION_PATCH = "2506"
-NGSOLVE_VERSION_HASH = "g5e4fc20dd"
+NGSOLVE_VERSION_HASH = "gf95a1c05a"
 
 CMAKE_CXX_COMPILER           = "/opt/rh/gcc-toolset-14/root/usr/bin/c++"
 CMAKE_CUDA_COMPILER          = ""

{ngsolve-6.2.2506.post85.dev0.data → ngsolve-6.2.2506.post98.dev0.data}/data/lib/python3.13/site-packages/ngsolve/nonlinearsolvers.py

@@ -1,4 +1,4 @@
-from ngsolve.la import InnerProduct
+from ngsolve.la import InnerProduct, SparseFactorizationInterface
 from math import sqrt
 from ngsolve import Projector, Norm
 from .utils import TimeFunction
@@ -97,7 +97,7 @@ class NewtonSolver:
             else:
                 self.inv.Update()
         else:
-            if self.inverse == "sparsecholesky" and self.inv is not None:
+            if self.inv is not None and (self.inverse == "sparsecholesky" or isinstance(self.inv, SparseFactorizationInterface)):
                 self.inv.Update()
             else:
                 self.inv = self.a.mat.Inverse(self.freedofs,

{ngsolve-6.2.2506.post85.dev0.data → ngsolve-6.2.2506.post98.dev0.data}/data/lib/python3.13/site-packages/ngsolve/solvers/cudss.py

@@ -8,8 +8,11 @@ import ngsolve.la as ngla
 import scipy.sparse as sp
 import numpy as np
 
+from ngsolve import TimeFunction
+
 
 class CudssSolver(ngla.SparseFactorizationInterface):
+    @TimeFunction
     def Analyze(self):
         self._mat = self.GetInnerMatrix()
         csr = sp.csr_matrix(self._mat.CSR())
@@ -18,9 +21,11 @@ class CudssSolver(ngla.SparseFactorizationInterface):
         self.solver = nvs.DirectSolver(csr, self._tmp, options=options)
         self.solver.plan()
 
+    @TimeFunction
     def Factor(self):
         self.solver.factorize()
 
+    @TimeFunction
     def Solve(self, b, sol):
         self.solver.reset_operands(b=b.FV().NumPy())
         sol.FV().NumPy()[:] = self.solver.solve()

{ngsolve-6.2.2506.post85.dev0.dist-info → ngsolve-6.2.2506.post98.dev0.dist-info}/METADATA

@@ -1,11 +1,11 @@
 Metadata-Version: 2.4
 Name: ngsolve
-Version: 6.2.2506.post85.dev0
+Version: 6.2.2506.post98.dev0
 Summary: NGSolve
 Author: The NGSolve team
 License: LGPL2.1
 License-File: LICENSE
-Requires-Dist: netgen-mesher==6.2.2506.post14.dev0
+Requires-Dist: netgen-mesher==6.2.2506.post24.dev0
 Requires-Dist: mkl
 Dynamic: author
 Dynamic: license