ngsolve 6.2.2505__cp310-cp310-macosx_10_15_universal2.whl → 6.2.2505.post17.dev0__cp310-cp310-macosx_10_15_universal2.whl

netgen/include/fespace.hpp

@@ -448,7 +448,8 @@ ANY                  1 1 1 1 | 15
     virtual void GetDofNrs (NodeId ni, Array<DofId> & dnums) const;
     BitArray GetDofs (const Region & reg) const;
     Table<int> CreateDofTable (VorB vorb) const;
-
+    virtual void SelectDofs (const string & name, BitArray & dofs) const;
+                  
     /// get coupling types of dofs
     virtual void GetDofCouplingTypes (int elnr, Array<COUPLING_TYPE> & dnums) const;
     
@@ -1141,7 +1142,8 @@ ANY                  1 1 1 1 | 15
     void GetEdgeDofNrs (int ednr, Array<DofId> & dnums) const override;
     void GetFaceDofNrs (int fanr, Array<DofId> & dnums) const override;
     void GetInnerDofNrs (int elnr, Array<DofId> & dnums) const override;
-
+    void SelectDofs (const string & name, BitArray & dofs) const override;
+    
     void SolveM(CoefficientFunction * rho, BaseVector & vec, Region * definedon,
                         LocalHeap & lh) const override;
     void ApplyM(CoefficientFunction * rho, BaseVector & vec, Region * definedon,

netgen/include/h1amg.hpp

@@ -9,6 +9,20 @@
 namespace ngcomp
 {
   using namespace ngla;
+
+
+
+  struct H1AMG_Parameters
+  {
+    int verbose = 0;
+    int smoothing_steps = 3;
+    int coarsenings_per_level = 3; // number of binary coarsenings per level
+    bool block_smoother = true;    // block or point smoother ?
+    bool use_smoothed_prolongation = true;
+    int max_coarse = 10;
+    int max_level = 20;
+  };
+
   
   template <class SCAL>
   class NGS_DLL_HEADER H1AMG_Matrix : public ngla::BaseMatrix
@@ -26,6 +40,7 @@ namespace ngcomp
                   ngcore::FlatArray<ngcore::IVec<2>> e2v,
                   ngcore::FlatArray<double> edge_weights,
                   ngcore::FlatArray<double> vertex_weights,
+                  const H1AMG_Parameters & param,                  
                   size_t level);
 
     virtual int VHeight() const override { return size; }
@@ -50,7 +65,9 @@ namespace ngcomp
 
     ParallelHashTable<IVec<2>,double> edge_weights_ht;
     ParallelHashTable<IVec<1>,double> vertex_weights_ht;
-
+    
+    H1AMG_Parameters param;
+    
   public:
 
     static shared_ptr<Preconditioner> CreateBF (shared_ptr<BilinearForm> bfa, const Flags & flags, const string & name)
@@ -71,6 +88,12 @@ namespace ngcomp
         cout << IM(3) << "Create H1AMG" << endl;
       else
         cout << IM(3) << "Create H1AMG, complex" << endl;
+      
+      param.verbose = int(flags.GetNumFlag("verbose", 0));
+      param.smoothing_steps = int(flags.GetNumFlag("smoothingsteps", 1));    
+      param.use_smoothed_prolongation = flags.GetDefineFlagX("smoothedprolongation").IsMaybeTrue();
+      param.max_coarse = int(flags.GetNumFlag("maxcoarse", 10));
+      param.max_level = int(flags.GetNumFlag("maxlevel", 20));
     }
 
     virtual void InitLevel (shared_ptr<BitArray> _freedofs) override

netgen/include/hdivhofespace.hpp

@@ -143,6 +143,8 @@ namespace ngcomp
 	} 
     }
 
+    void SelectDofs (const string & name, BitArray & dofs) const override;
+
     int GetFirstInnerDof(int elnr) const { return(first_inner_dof[elnr]);}; 
     // virtual int LowOrderDof() const { if(discont) return(0); else return(1);} 
 

netgen/include/mptools.hpp

@@ -877,6 +877,34 @@ namespace ngsbem
   class NGS_DLL_HEADER RegularMLMultiPole
   {
     static Array<size_t> nodes_on_level;
+
+    
+    struct RecordingRS
+    {
+      const MultiPole<MPSingular,elem_type> * mpS;
+      MultiPole<MPRegular,elem_type> * mpR;
+      Vec<3> dist;
+      double len, theta, phi;
+    public:
+      RecordingRS() = default;
+      RecordingRS (const MultiPole<MPSingular,elem_type> * ampS,
+                   MultiPole<MPRegular,elem_type> * ampR,
+                   Vec<3> adist)
+        : mpS(ampS), mpR(ampR), dist(adist)
+      {
+        len = L2Norm(dist);
+        if (len < 1e-30)
+          theta = 0;
+        else
+          theta = acos (dist(2) / len);
+        
+        if (sqr(dist(0))+sqr(dist(1)) < 1e-30)
+          phi = 0;
+        else
+          phi = atan2(dist(1), dist(0));
+      }
+    };
+
     
     struct Node
     {
@@ -912,8 +940,9 @@ namespace ngsbem
             childs[i] = make_unique<Node> (cc, r/2, level+1, mp.Kappa());
           }
       }
-
-      void AddSingularNode (const typename SingularMLMultiPole<elem_type>::Node & singnode, bool allow_refine)
+      
+      void AddSingularNode (const typename SingularMLMultiPole<elem_type>::Node & singnode, bool allow_refine,
+                            Array<RecordingRS> * recording)
       {
         if (mp.SH().Order() < 0) return;
         if (singnode.mp.SH().Order() < 0) return;
@@ -936,12 +965,15 @@ namespace ngsbem
                 singnode.childs[0]->mp.Order() < singnode.mp.Order())
               {
                 for (auto & child : singnode.childs)
-                  AddSingularNode (*child, allow_refine);
+                  AddSingularNode (*child, allow_refine, recording);
                 return;
               }
 
             // static Timer t("mptool transform Helmholtz-criterion"); RegionTimer r(t);
-            singnode.mp.TransformAdd(mp, dist);
+            if (recording)
+              *recording += RecordingRS(&singnode.mp, &mp, dist);
+            else
+              singnode.mp.TransformAdd(mp, dist);
             return;
           }
 
@@ -960,21 +992,21 @@ namespace ngsbem
                   CreateChilds();
                 
                 for (auto & ch : childs)
-                  ch -> AddSingularNode (singnode, allow_refine);
+                  ch -> AddSingularNode (singnode, allow_refine, recording);
               }
             else
               {
-                if (total_targets < 1000)
+                if (total_targets < 1000 || recording)
                   {
                     for (auto & ch : childs)
                       if (ch)
-                        ch -> AddSingularNode (singnode, allow_refine);
+                        ch -> AddSingularNode (singnode, allow_refine, recording);
                   }
                 else
                   ParallelFor (8, [&] (int nr)
                                {
                                  if (childs[nr])
-                                   childs[nr] -> AddSingularNode (singnode, allow_refine);
+                                   childs[nr] -> AddSingularNode (singnode, allow_refine, recording);
                                });
                 
                 if (targets.Size())
@@ -984,7 +1016,7 @@ namespace ngsbem
         else
           {
             for (auto & childsing : singnode.childs)
-              AddSingularNode (*childsing, allow_refine);
+              AddSingularNode (*childsing, allow_refine, recording);
           }
       }
 
@@ -1158,8 +1190,8 @@ namespace ngsbem
       nodes_on_level = 0;
       nodes_on_level[0] = 1;
       {
-        static Timer t("mptool compute regular MLMP"); RegionTimer rg(t);            
-        root.AddSingularNode(singmp->root, true);
+        static Timer t("mptool compute regular MLMP"); RegionTimer rg(t);
+        root.AddSingularNode(singmp->root, true, nullptr);
         // cout << "norm after S->R conversion: " << root.Norm() << endl;
       }
 
@@ -1200,9 +1232,27 @@ namespace ngsbem
 
       root.CalcTotalTargets();
       root.RemoveEmptyTrees();
-        
-      root.AddSingularNode(singmp->root, false);
 
+      
+      root.AddSingularNode(singmp->root, false, nullptr);
+      /*
+      Array<RecordingRS> recording;      
+      root.AddSingularNode(singmp->root, false, &recording);
+      // cout << "recorded: " << recording.Size() << endl;
+      QuickSort (recording, [] (auto & a, auto & b)
+      {
+        if (a.len < (1-1e-8) * b.len) return true;
+        if (a.len > (1+1e-8) * b.len) return false;
+        return a.theta < b.theta;
+      });
+      for (auto & record : recording)
+        {
+          record.mpS->TransformAdd(*record.mpR, record.dist);
+          // *testout << record.len << ", " << record.theta << ", " << record.phi << endl;
+        }
+      */
+
+      
       /*
       int maxlevel = 0;
       for (auto [i,num] : Enumerate(RegularMLMultiPole::nodes_on_level))

netgen/include/ngblas.hpp

@@ -584,12 +584,31 @@ namespace ngbla
       }
   }
 
-  template <typename TM, typename TVX, typename TVY>
-  extern void TestFunc (TM m, TVX x, TVY y);
-
+  // template <typename TM, typename TVX, typename TVY>
+  // extern void TestFunc (TM m, TVX x, TVY y);
+  
+  
+  template <typename TS, typename T> constexpr bool IsVec = false;
+  template <typename TS, int S> constexpr bool IsVec<TS, Vec<S,TS>> = true;  
+  
   template <bool ADD, bool POS, typename TM, ORDERING ORD, typename TX, typename TY>
   INLINE void NgGEMV (BareSliceMatrix<TM,ORD> a, FlatVector<const TX> x, FlatVector<TY> y)
   {
+    if constexpr (std::is_same<TM,double>() && std::is_same<TX,TY>() && IsVec<Complex,TX>)
+      {
+        FlatMatrix<double> mx(x.Size(), sizeof(TX)/sizeof(double), (double*)(void*)x.Addr(0));
+        FlatMatrix<double> my(y.Size(), sizeof(TX)/sizeof(double), (double*)(void*)y.Addr(0));
+        NgGEMM<ADD,POS> (a.AddSize(y.Size(), x.Size()),make_SliceMatrix(mx), make_SliceMatrix(my));
+        return;
+      }
+    if constexpr (std::is_same<TM,Complex>() && std::is_same<TX,TY>() && IsVec<Complex,TX>)
+      {
+        FlatMatrix<Complex> mx(x.Size(), sizeof(TX)/sizeof(Complex), &const_cast<Complex&>(*(x.Data()->Data())));
+        FlatMatrix<Complex> my(y.Size(), sizeof(TX)/sizeof(Complex), y.Data()->Data());
+        NgGEMM<ADD,POS> (a.AddSize(y.Size(), x.Size()),make_SliceMatrix(mx), make_SliceMatrix(my));
+        return;
+      }
+
     if (!ADD)
       {
         if (!POS)
@@ -635,6 +654,16 @@ namespace ngbla
   extern NGS_DLL_HEADER  
   void NgGEMV (double s, BareSliceMatrix<double,ord> a, SliceVector<double> x, SliceVector<double> y) NETGEN_NOEXCEPT;
   */
+
+
+
+
+
+
+
+  /* *********************** GEMV - SliceVector **************************** */
+
+
   
   template <bool ADD, ORDERING ord>
   extern NGS_DLL_HEADER  
@@ -714,7 +743,30 @@ namespace ngbla
   }
 
 
-
+  
+  
+  template <bool ADD, bool POS, typename TM, ORDERING ORD, typename TX, typename TY>
+  INLINE void NgGEMV (BareSliceMatrix<TM,ORD> a, SliceVector<TX> x, SliceVector<TY> y)
+  {
+    if constexpr (std::is_same<TM,double>() && std::is_same<TX,TY>() && IsVec<Complex,TX>)
+      {
+        constexpr int VS = sizeof(TX)/sizeof(double);
+        SliceMatrix<double> mx(x.Size(), VS, x.Dist()*VS, (double*)(void*)x.Addr(0));
+        SliceMatrix<double> my(y.Size(), VS, y.Dist()*VS, (double*)(void*)y.Addr(0));
+        NgGEMM<ADD,POS> (a.AddSize(y.Size(), x.Size()),make_SliceMatrix(mx), make_SliceMatrix(my));
+        return;
+      }
+    else if constexpr (std::is_same<TM,Complex>() && std::is_same<TX,TY>() && IsVec<Complex,TX>)
+      {
+        constexpr int VS = sizeof(TX)/sizeof(Complex);
+        SliceMatrix<Complex> mx(x.Size(), VS, x.Dist()*VS, &const_cast<Complex&>(*(x.Data()->Data())));
+        SliceMatrix<Complex> my(y.Size(), VS, y.Dist()*VS, y.Data()->Data());
+        NgGEMM<ADD,POS> (a.AddSize(y.Size(), x.Size()),make_SliceMatrix(mx), make_SliceMatrix(my));
+        return;
+      }
+    else
+      NgGEMV<ADD> (POS ? 1.0 : -1.0, a, x, y);
+  }
 
 
 
@@ -880,9 +932,14 @@ namespace ngbla
                          FlatVector<const TB>(prod.View().B().Range(0,w)),
                          FlatVector<T>(self.Spec().Range(0,h)));
       else
+        NgGEMV<ADD,POS> (make_BareSliceMatrix(prod.View().A()),
+                         SliceVector<TB>(prod.View().B().Range(0,w)),
+                         SliceVector<T>(self.Spec().Range(0,h)));
+        /*
         NgGEMV<ADD> (POS ? 1.0 : -1.0, make_BareSliceMatrix(prod.View().A()),
                      SliceVector<TB>(prod.View().B().Range(0,w)),
                      SliceVector<T>(self.Spec().Range(0,h)));
+        */
       return self.Spec();
     }
   };
@@ -1101,6 +1158,47 @@ namespace ngbla
   
 
 
+  template <bool ADD, bool POS, ORDERING OA, ORDERING OB>
+  extern NGS_DLL_HEADER void NgGEMMBare (size_t ah, size_t aw, size_t bw, BareSliceMatrix<Complex,OA> a, BareSliceMatrix<Complex,OB> b, BareSliceMatrix<Complex,RowMajor> c);
+
+
+  template <bool ADD, bool POS, ORDERING OA, ORDERING OB>
+  void NgGEMM (SliceMatrix<Complex,OA> a, SliceMatrix<Complex,OB> b, SliceMatrix<Complex,RowMajor> c)
+  {
+    size_t ah = a.Height();
+    size_t aw = a.Width();
+    size_t bw = b.Width();
+    NgGEMMBare<ADD,POS>(ah, aw, bw, make_BareSliceMatrix(a), make_BareSliceMatrix(b), make_BareSliceMatrix(c));
+  }
+  
+  template <bool ADD, bool POS, ORDERING OA, ORDERING OB>
+  void NgGEMM (SliceMatrix<Complex,OA> a, SliceMatrix<Complex,OB> b, SliceMatrix<Complex,ColMajor> c)
+  {
+    NgGEMM<ADD,POS> (Trans(b), Trans(a), Trans(c));
+  }
+  
+  template <typename OP, typename T, typename TA, typename TB>
+  class assign_trait<OP, T, MultExpr<TA, TB>,
+                     enable_if_t<IsConvertibleToSliceMatrix<TA,Complex>() &&
+                                 IsConvertibleToSliceMatrix<TB,Complex>() &&
+                                 IsConvertibleToSliceMatrix<T,Complex>(), int>>
+  {
+  public:
+    static inline T & Assign (MatExpr<T> & self, const Expr<MultExpr<TA, TB>> & prod) 
+    {
+      constexpr bool ADD = std::is_same<OP,typename MatExpr<T>::AsAdd>::value || std::is_same<OP,typename MatExpr<T>::AsSub>::value;
+      constexpr bool POS = std::is_same<OP,typename MatExpr<T>::As>::value || std::is_same<OP,typename MatExpr<T>::AsAdd>::value;
+
+      size_t n = CombinedSize(prod.View().A().Height(), self.Spec().Height());
+      size_t m = CombinedSize(prod.View().B().Width(), self.Spec().Width());
+      size_t k = CombinedSize(prod.View().A().Width(), prod.View().B().Height());
+      
+      NgGEMM<ADD,POS> (make_BareSliceMatrix(prod.View().A()).AddSize(n,k).RemoveConst(),
+                       make_BareSliceMatrix(prod.View().B()).AddSize(k,m).RemoveConst(),
+                       make_BareSliceMatrix(self.Spec()).AddSize(n,m));
+      return self.Spec();
+    }
+  };
 
   
 

netgen/include/sparsematrix_dyn.hpp

@@ -50,9 +50,9 @@ namespace ngla
     virtual void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
     
     AutoVector CreateRowVector() const override
-    { throw make_unique<VVector<TSCAL>> (VWidth()); }
+    { return make_unique<VVector<TSCAL>> (VWidth()); }
     AutoVector CreateColVector() const override
-    { throw make_unique<VVector<TSCAL>> (VHeight()); }
+    { return make_unique<VVector<TSCAL>> (VHeight()); }
     
     virtual tuple<int,int> EntrySizes() const override { return { bh, bw }; }
   };

ngsolve/cmake/NGSolveConfig.cmake

@@ -1,4 +1,4 @@
-set(PACKAGE_VERSION "6.2.2505")
+set(PACKAGE_VERSION "6.2.2505-17-g2593d4c47")
 find_package(Netgen CONFIG REQUIRED HINTS
   ${CMAKE_CURRENT_LIST_DIR}
   ${CMAKE_CURRENT_LIST_DIR}/..

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.2505"
-NGSOLVE_VERSION_GIT = "v6.2.2505-0-g8fbbbc389"
-NGSOLVE_VERSION_PYTHON = "6.2.2505"
+NGSOLVE_VERSION = "6.2.2505-17-g2593d4c47"
+NGSOLVE_VERSION_GIT = "v6.2.2505-17-g2593d4c47"
+NGSOLVE_VERSION_PYTHON = "6.2.2505.post17.dev0"
 
 NGSOLVE_VERSION_MAJOR = "6"
 NGSOLVE_VERSION_MINOR = "2"
-NGSOLVE_VERSION_TWEAK = "0"
+NGSOLVE_VERSION_TWEAK = "17"
 NGSOLVE_VERSION_PATCH = "2505"
-NGSOLVE_VERSION_HASH = "g8fbbbc389"
+NGSOLVE_VERSION_HASH = "g2593d4c47"
 
 CMAKE_CXX_COMPILER           = "/Library/Developer/CommandLineTools/usr/bin/c++"
 CMAKE_CUDA_COMPILER          = ""

{ngsolve-6.2.2505.dist-info → ngsolve-6.2.2505.post17.dev0.dist-info}/METADATA

@@ -1,13 +1,13 @@
 Metadata-Version: 2.1
 Name: ngsolve
-Version: 6.2.2505
+Version: 6.2.2505.post17.dev0
 Summary: NGSolve
 Home-page: UNKNOWN
 Author: The NGSolve team
 License: LGPL2.1
 Platform: UNKNOWN
 License-File: LICENSE
-Requires-Dist: netgen-mesher (==6.2.2505)
+Requires-Dist: netgen-mesher (==6.2.2505.post3.dev0)
 
 UNKNOWN
 

{ngsolve-6.2.2505.dist-info → ngsolve-6.2.2505.post17.dev0.dist-info}/RECORD

@@ -1,10 +1,10 @@
-netgen/libngbla.dylib,sha256=LSO6sC7ms_Opr0pOPN54dLBj1RrvO6SsAWu1IkVGYCw,5183872
-netgen/libngcomp.dylib,sha256=za2Vreqy_oLc1fAzqWi2cmYOGkc54pRRpPAKBfsL2hk,42847792
-netgen/libngfem.dylib,sha256=vH2E_9Uvi6VtTCQqDIwzcOL7yPZIvaJICdYz_-jojhg,51449008
-netgen/libngla.dylib,sha256=LTyxtZzcw6aBZiUFhdQOywK_tVKvquaIzbLf5oerna0,25323904
-netgen/libngsbem.dylib,sha256=P3gt1MCnN9nNYUl0cMJ11rju9NTI2qEkhLY3D4KMyeQ,5000736
+netgen/libngbla.dylib,sha256=iSmy-_zc_CUppaRSGYB4slyD9RCeeUay-iF7piOZkkQ,5343120
+netgen/libngcomp.dylib,sha256=bxMalulee1YjIgYIjuWSGddT76c1TZISj0U-RApX-1c,42901584
+netgen/libngfem.dylib,sha256=sJLGft0YHgBOpoTQyaLileCldH26_D9n2Topc0OdeJI,51417376
+netgen/libngla.dylib,sha256=hzPi82Ei6ReXdK-DkDekfuncTolEz7j2QBVkOOHZyeY,25308000
+netgen/libngsbem.dylib,sha256=Je2j8_IyyAYEsBdZMsqKjd7SFKehuTewpMMqyh63QJo,5500000
 netgen/libngsolve.dylib,sha256=1Csyic1b0bLPAl45sYsHSKVB20L48Auh4mwQLbLoot8,939136
-netgen/libngstd.dylib,sha256=WEB_mrt0iefNn2fiSfzC71L4FhWNoEGVhLgCwD_lnG0,922576
+netgen/libngstd.dylib,sha256=nRbe88OAQLD47nGY9XLTTGusXXB8wiKpXA5dHC2uVqg,922576
 netgen/include/arnoldi.hpp,sha256=b8mRu9wGyaLa7z6mXCFIhxfW_xcQvsOMD6TqbfYZ05g,1437
 netgen/include/bandmatrix.hpp,sha256=zwKz89QY53MPsWgNsl0cXEAFKka_W4sXPIO54M3eT80,7778
 netgen/include/basematrix.hpp,sha256=6YGW07HMH7F2JWXUKOwlGZ2ZvEoP7uCSGihdrzzKPag,28634
@@ -54,12 +54,12 @@ netgen/include/facetsurffespace.hpp,sha256=3FTy--cwUqVFZqF9MAaK78aRSIHWQUgOvSjPr
 netgen/include/fe_interfaces.hpp,sha256=gWxTH8E-g3RgDfLG7kkggKUIBC_zWefOaJPpYFBiiDI,822
 netgen/include/fem.hpp,sha256=xXxfiwpi3EW3jJZtWE9WxRJm9htSaivfKhCbD2_ToBo,2094
 netgen/include/fesconvert.hpp,sha256=WL-icGOUDl2edBG6QoI1MCohKibLqPVFuaTzyAenIb0,560
-netgen/include/fespace.hpp,sha256=tRyAbXR5q3-X-gYt-Sgy46ZiSA3n3bBJeQBCoR-x-1I,47056
+netgen/include/fespace.hpp,sha256=xWyWuCXziIiAnV2Gh7DhwLJJ_yQfr1_yVssiPV6h74A,47227
 netgen/include/finiteelement.hpp,sha256=yGHVHO5GJqHCnEcvyNwtXZHFRhebVTdaBfM96jlIppg,9145
 netgen/include/globalinterfacespace.hpp,sha256=JMRw_B8kfWRAlhZglT0IH9fFvcQbHZgF1T0tnIuszsM,2719
 netgen/include/globalspace.hpp,sha256=Hjzmz9zaUuKVSkduBYVr6zbRHV503wK2vdFL6Vxr_fY,3964
 netgen/include/gridfunction.hpp,sha256=em4M_5oUhofbD5jyoaAinKns4v4G9TCzs49Q2qANy0A,18041
-netgen/include/h1amg.hpp,sha256=O_fYedZwGOU3QxaYeyWuOnGbRzMXYVT3yDdJI-txRgM,3002
+netgen/include/h1amg.hpp,sha256=ygadM9_MFv0Qc47jP8G__rhnyTtuj_QNETwpm-8WYw8,3801
 netgen/include/h1hofe.hpp,sha256=CSVONOnzWpt6bf-D2dacEheU-3OmID3QUeRD64Uoqz8,5844
 netgen/include/h1hofe_impl.hpp,sha256=ruuxMbKB_jNOCviiPhj5PF0vR_idTSQMvb9jfU2_ylo,35456
 netgen/include/h1hofefo.hpp,sha256=9wQDD_5zXMhlptPHZZzOUImFLkF0dZwYuYbLX1Ai2dA,4465
@@ -87,7 +87,7 @@ netgen/include/hdivfes.hpp,sha256=MrPtDM7GXnozM-PbOqyZX1X7CsDNuJxtTGE5tgtAXOE,21
 netgen/include/hdivhofe.hpp,sha256=fTcC2CaT3fghZOPwsmJUFI-XMFfZTfBgyJfCPxwu3pU,13703
 netgen/include/hdivhofe_impl.hpp,sha256=G59J6oeF5AS51G9TQBMFf-_1smPdjYFcPB5Vcux6dTI,37278
 netgen/include/hdivhofefo.hpp,sha256=_RbPV_e--FvyQheAkvGYIKCmKys-ludl26XbrD07qeo,6390
-netgen/include/hdivhofespace.hpp,sha256=31ozqB31tlWIQN3FtnQxNJrrof5OijewEeA1ITgx5ts,4955
+netgen/include/hdivhofespace.hpp,sha256=F6lNi2Yd4m_34scgyvrcvT5Eshs92NkO5nrzDC6YBEM,5031
 netgen/include/hdivhosurfacefespace.hpp,sha256=Dp8EmF_Bj-jNvPQS_KEN87R4iqQgaYZnjCysRaruGw0,2843
 netgen/include/hdivlofe.hpp,sha256=QESj6iJ8UGLboWlCE3EHmPcD4Z5sMlaf6Uo7f7KFe7M,14603
 netgen/include/hidden.hpp,sha256=nYQrLY1eO8MiVZWfYLi5Sp8bOwiK4C7ZGPjQToLx_Os,2646
@@ -112,13 +112,13 @@ netgen/include/memusage.hpp,sha256=H-8rWQ8UQDCnOiYQnbwvILPuPoJD8bqwiNGB2QAUa0o,1
 netgen/include/meshaccess.hpp,sha256=j5dinD-8bVm18MkVqG9j-IhhZUHHXw3erdZy-LHxl-8,47133
 netgen/include/mgpre.hpp,sha256=zdNozWRSyluOF8E4GLeOR290Yw2wqGYGc1sTTSB5n3U,4825
 netgen/include/mp_coefficient.hpp,sha256=lzFP-XHYnFKDdAbb_4cWqEpCIgk4gwUWxOmwvJvSN3w,4736
-netgen/include/mptools.hpp,sha256=Wh2wxGYmO7iJYCH-B4iiZ5nrwR2bepwjLQ1k-4swCpU,34474
+netgen/include/mptools.hpp,sha256=I0w37lAM8u2qZGAngoQj60HR9jKy4INUVK0Wp4LT80k,36032
 netgen/include/multigrid.hpp,sha256=ovl1GXcTfXuYLVVMcFfCKP5sg6diSGlqdjCMMPCgFUo,856
 netgen/include/multivector.hpp,sha256=nuLxMiq2y6HoH5w9j6ka2BhJjRyREPeEnNFx_x0SMU0,12635
 netgen/include/mumpsinverse.hpp,sha256=Jaa5NPQzwPo2kSbHijI_LX3tr0qUmVhqVlRM14KeD18,4062
 netgen/include/mycomplex.hpp,sha256=8xKKByk_qMMCvDAl-Wg5wN_Y8A1eCGxe8sWowWdG0co,7662
 netgen/include/ng_lapack.hpp,sha256=1sH75a5q-ja2lbGU450tiX-z1KPkerMV5ftjzUg4Inc,48408
-netgen/include/ngblas.hpp,sha256=ArfhQ60ahMUgo6voSK8Z3EHxPJhfT0BZT-mBgNGGNus,42769
+netgen/include/ngblas.hpp,sha256=4ZsAHD5j0ulPIn2atd5mxrh6zf78xniy7DTGB_HWkHQ,47303
 netgen/include/ngs_defines.hpp,sha256=bGB_MsAq5tvNml55sklzLjGkzLuH3ZHCPCriS8zM7yM,686
 netgen/include/ngs_stdcpp_include.hpp,sha256=RUjXe5YfnAIBDg8jInLBw7DXu5AYvrkja3EFhMBa5fY,2551
 netgen/include/ngs_utils.hpp,sha256=qoXJBajujUW35GXix_olD7C7S3xWUeOcQ23kIF_hjXw,3144
@@ -161,7 +161,7 @@ netgen/include/smoother.hpp,sha256=hZAe15R6BU88hMfhcTM3-eclRgzmHJnPKgYYWol9Sx8,6
 netgen/include/solve.hpp,sha256=OZQfkg6HN20eudQ4kkRSecr-miwV35NPEjAtoMcunG4,2034
 netgen/include/sparsecholesky.hpp,sha256=2HKM6aKO02YX32jxGna6CAjBBfeXvCsmPk4hoWZWIQk,8888
 netgen/include/sparsematrix.hpp,sha256=7YN50okcY5Lxipc43njHPKzzKZ_gLtMhwQB_Jsj99rE,32706
-netgen/include/sparsematrix_dyn.hpp,sha256=XGTvSgtG3sOj6ho_Bu3VcZzUxO49VWztlxyuW--3Z9M,2654
+netgen/include/sparsematrix_dyn.hpp,sha256=u9aWEYl6PhrevL9SPSXVE22_EIwIduDDvzY8xgIsw4Q,2656
 netgen/include/sparsematrix_impl.hpp,sha256=Z45m3piI8zPvBkRgdQiQLHmSRuL82gO6WPK1b2frW6M,29580
 netgen/include/special_matrix.hpp,sha256=1moUjxQ_uaS-GWrSbrbN4c8N3pDma9UpSaFqrNUXV80,16526
 netgen/include/specialelement.hpp,sha256=iiQaIE-cHbwGX014jET8k-rfb8k1GT6Gs_vsMIBWaOA,3226
@@ -213,12 +213,12 @@ ngsolve/timestepping.py,sha256=SLbsgE5pY5tXYTLGlVypk4kS2Ui8GxOwWmKVxqneZRU,8350
 ngsolve/timing.py,sha256=Xqkq2byJjNNK9vJjHG2O3cMf6xuKIyySjk8yv6Le2T8,3985
 ngsolve/utils.py,sha256=V2yVF4K3uv1007-DEs_GrzJ3BOOh09c4lDGGR1Xw-PU,4557
 ngsolve/webgui.py,sha256=EUhaOtWHZy80aLXhXbHwbCQVNjRSw-BAAE50cz7Ydaw,24313
-ngsolve/cmake/NGSolveConfig.cmake,sha256=nQ1upVJl8qYicatuXmaNsJ_OsJWTtkG92Gyn1d4omAs,5052
+ngsolve/cmake/NGSolveConfig.cmake,sha256=AioHU74ZxtvojoLmtn0zgTyfC2g_VoDaXvjKCYVRPmg,5066
 ngsolve/cmake/ngsolve-targets-release.cmake,sha256=2dUHzaaMQX08eT7Tjti0crzp_Ah3oqW86gp5UgSufmw,3426
 ngsolve/cmake/ngsolve-targets.cmake,sha256=1cQTFW-4aJfY356XRcuslSbFClmwPOSS03bj3WL15iM,5637
 ngsolve/config/__init__.py,sha256=A00lpnUKg-q3y_b-cN3jKF9urkGTGjb2Yz_gGc9WCLs,22
 ngsolve/config/__main__.py,sha256=m8k2RrXnXR5oYaUna1dwkKNhOZBiM7EyFSE9kIcZKvo,89
-ngsolve/config/config.py,sha256=y5lsmOBrzLQpM2wAaXBIwpVX3Kq4otVzz6X92m5qFJo,2439
+ngsolve/config/config.py,sha256=1aI8x1EvjzC-n4t_MV6QxJqW8jTwhb1KcNjA1C8JRtU,2467
 ngsolve/demos/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 ngsolve/demos/TensorProduct/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 ngsolve/demos/TensorProduct/tp_dg_1d_1d.py,sha256=CrJCkbJNinVYUL3DO-4DnCAMGSAhegNc0QVWiJTk3Js,2109
@@ -250,41 +250,41 @@ ngsolve/demos/mpi/mpi_cmagnet.py,sha256=HuN4lttA2q2kNWss6LZNUKraHqskRQNVEBRjkrU2
 ngsolve/demos/mpi/mpi_navierstokes.py,sha256=Sct0G6EziYVfeIeSa7_b6r3WiYXcxZekoOGsZQALx-c,2911
 ngsolve/demos/mpi/mpi_poisson.py,sha256=UvXODjpWtqIgqq7lvMJf3ZEoJZLY1Dck-82R7auiW6I,2468
 ngsolve/demos/mpi/mpi_timeDG.py,sha256=TPJZ-DwmzOpuAol-fFKeU-R2rJ74MbjOm3DDkQI9BcY,1881
-ngsolve-6.2.2505.data/data/Netgen.icns,sha256=Ebqct4LKxLzxdJLwGg0tTUQb1VjLwkKhtmdYhpRhKjg,291475
-ngsolve-6.2.2505.data/data/bin/ngscxx,sha256=HFcVHkblKlDxFT-obJQYZU1pGqopDwzoO9kU6MMGS5w,1755
-ngsolve-6.2.2505.data/data/bin/ngsld,sha256=50UErkFWLTocn8GtZDTkUWmPTIqsoK5lcfBGwyP0Vnc,557
-ngsolve-6.2.2505.data/data/bin/ngsolve.tcl,sha256=SSBYLyp4zgl6eM17LXIiLoj1Doy-h2s4EZg6oKi6ucA,20995
-ngsolve-6.2.2505.data/data/bin/ngspy,sha256=KglThT4AqnzC-4q3Q_UdoOrCEvgISINlJFEChpY8GF8,106
-ngsolve-6.2.2505.data/data/share/ngsolve/beam.geo,sha256=EtxklGAIo6SxZlHHC_FZHoQGZqZUDGJQSheEMFO-bJ8,387
-ngsolve-6.2.2505.data/data/share/ngsolve/beam.vol,sha256=aotG4VA7sIgAsPQsYkzbMP1Jta_TddG2gMCblkSPkqc,13594
-ngsolve-6.2.2505.data/data/share/ngsolve/chip.in2d,sha256=g6C-mcoD8-q9JtZpuh4mf4iFEZUIGvvivp4C5YYMlAw,551
-ngsolve-6.2.2505.data/data/share/ngsolve/chip.vol,sha256=if0yLn8E3P4mZ_X5sBcrEIp2C4ZhMysGROtB1YCkJHI,27897
-ngsolve-6.2.2505.data/data/share/ngsolve/coil.geo,sha256=PVD4I9ztZowKwNOUrjnoOP7Tj2KG9J0UCJFRS9EhlxU,364
-ngsolve-6.2.2505.data/data/share/ngsolve/coil.vol,sha256=IOZe4oPxMUa9vBGJWN4LG9M79XAruY23S2rtAh_2_98,144724
-ngsolve-6.2.2505.data/data/share/ngsolve/coilshield.geo,sha256=lCUDCAhfowBttSGEY_giyY0xEFInAbgVfIPwSWBb504,605
-ngsolve-6.2.2505.data/data/share/ngsolve/coilshield.vol,sha256=2Dq5VCX7dqGxo_ExtUhI716I-SnxhLGPawW29IR9nYY,162397
-ngsolve-6.2.2505.data/data/share/ngsolve/cube.geo,sha256=kbBNVUYVSf1ogV7-Xf31KXju3LCUIPaDLfvwZHMgddI,386
-ngsolve-6.2.2505.data/data/share/ngsolve/cube.vol,sha256=54108a3RvrTUOJSFltxUY8iiNzyWoBTtNrWM0kufj7c,90295
-ngsolve-6.2.2505.data/data/share/ngsolve/d10_DGdoubleglazing.pde,sha256=bPQAvbYlyvCfLzk9cOzjXDn94UFC_vBCONl7fEIaMVs,1304
-ngsolve-6.2.2505.data/data/share/ngsolve/d11_chip_nitsche.pde,sha256=Mq3Ap63mWqjQyXsASFAlkRQ0AWxBlY3BYYejTYd-lwg,1101
-ngsolve-6.2.2505.data/data/share/ngsolve/d1_square.pde,sha256=ysOroxhoAOePmB_Aswps9iGjA2qV5435iu3hOnSvN8M,1029
-ngsolve-6.2.2505.data/data/share/ngsolve/d2_chip.pde,sha256=AowiGbuAWEdyFUB3d1LlZeQ7rarbkqIYsj1QcWc89-A,950
-ngsolve-6.2.2505.data/data/share/ngsolve/d3_helmholtz.pde,sha256=nG0HxrxboI2qGoR5_LNZt7RYB7ji2o5jtkyIqTFR4y8,568
-ngsolve-6.2.2505.data/data/share/ngsolve/d4_cube.pde,sha256=-C4RgoVtBXZu7XHCwZy7F7pr_PGNGcRwWCFzBppq1mQ,1033
-ngsolve-6.2.2505.data/data/share/ngsolve/d5_beam.pde,sha256=QOYC7YN7LdcjSFeLdJu1Tlfi-lmsVsa7b9G6qoCsGo8,1985
-ngsolve-6.2.2505.data/data/share/ngsolve/d6_shaft.pde,sha256=Iua3INfuIUeaBk21qNg4u5KZt__H1B8BZUsy5frjhdY,1901
-ngsolve-6.2.2505.data/data/share/ngsolve/d7_coil.pde,sha256=-ELH2obynMJcvSOqoYT1OThXcllIcdss0ZrklLSimZw,1200
-ngsolve-6.2.2505.data/data/share/ngsolve/d8_coilshield.pde,sha256=WUUulvd3dsedB25sFVQgU7FgZ3B7BkbxgkXTEBglk1w,1141
-ngsolve-6.2.2505.data/data/share/ngsolve/d9_hybridDG.pde,sha256=drLxnz1GvCEeme-FV772cGbS-YTgeUB4G5DrbkLdapw,1479
-ngsolve-6.2.2505.data/data/share/ngsolve/doubleglazing.in2d,sha256=bKn7m7eCcoLimGgoO3nxhiu7Mm0lIRnbx0TJh5dC1I0,506
-ngsolve-6.2.2505.data/data/share/ngsolve/doubleglazing.vol,sha256=Eg-H6nl2Hd5Nr593SFCCLVTDl25kIXuddrWpd7xlrP8,34987
-ngsolve-6.2.2505.data/data/share/ngsolve/piezo2d40round4.vol.gz,sha256=RuxM0vtYljoAImOH6U6avt0LmWaUbGNao37n8p_Yd1w,496477
-ngsolve-6.2.2505.data/data/share/ngsolve/shaft.geo,sha256=aPbk4EItJSLN_PjpOejnijk74y53QBW-IPCXHvLNlVs,2434
-ngsolve-6.2.2505.data/data/share/ngsolve/shaft.vol,sha256=fT8qu1L-p0MkI-l3WVisKG52W5FCzO_NzXaCvMJT1V4,270450
-ngsolve-6.2.2505.data/data/share/ngsolve/square.in2d,sha256=evpVgNNzReABv6bIs4J_r1OFeedMppoS-TEzE-tfKgY,192
-ngsolve-6.2.2505.data/data/share/ngsolve/square.vol,sha256=IXIHYmq85hkFsMVZ_vEVthUnd0YTQWE4jRegjW8_ImI,9162
-ngsolve-6.2.2505.dist-info/LICENSE,sha256=oZDcnIBDdV2Q-LCnX6ZrnkLUr0yYC_XdxjPwEk2zzuc,26430
-ngsolve-6.2.2505.dist-info/METADATA,sha256=D0Nkep9naTmDimDfGndBhGKGK2uyREojZbVpf1KlJxQ,224
-ngsolve-6.2.2505.dist-info/WHEEL,sha256=vYGgz4rYivjGPkIS_2Kj6g97rQxH536HeP_jo4e1Oq8,110
-ngsolve-6.2.2505.dist-info/top_level.txt,sha256=1llU44DxNIP85uPT4RChsIuV3h2d0_5j0Ui_Df5vV8E,15
-ngsolve-6.2.2505.dist-info/RECORD,,
+ngsolve-6.2.2505.post17.dev0.data/data/Netgen.icns,sha256=Ebqct4LKxLzxdJLwGg0tTUQb1VjLwkKhtmdYhpRhKjg,291475
+ngsolve-6.2.2505.post17.dev0.data/data/bin/ngscxx,sha256=HFcVHkblKlDxFT-obJQYZU1pGqopDwzoO9kU6MMGS5w,1755
+ngsolve-6.2.2505.post17.dev0.data/data/bin/ngsld,sha256=50UErkFWLTocn8GtZDTkUWmPTIqsoK5lcfBGwyP0Vnc,557
+ngsolve-6.2.2505.post17.dev0.data/data/bin/ngsolve.tcl,sha256=SSBYLyp4zgl6eM17LXIiLoj1Doy-h2s4EZg6oKi6ucA,20995
+ngsolve-6.2.2505.post17.dev0.data/data/bin/ngspy,sha256=KglThT4AqnzC-4q3Q_UdoOrCEvgISINlJFEChpY8GF8,106
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/beam.geo,sha256=EtxklGAIo6SxZlHHC_FZHoQGZqZUDGJQSheEMFO-bJ8,387
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/beam.vol,sha256=aotG4VA7sIgAsPQsYkzbMP1Jta_TddG2gMCblkSPkqc,13594
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/chip.in2d,sha256=g6C-mcoD8-q9JtZpuh4mf4iFEZUIGvvivp4C5YYMlAw,551
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/chip.vol,sha256=if0yLn8E3P4mZ_X5sBcrEIp2C4ZhMysGROtB1YCkJHI,27897
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/coil.geo,sha256=PVD4I9ztZowKwNOUrjnoOP7Tj2KG9J0UCJFRS9EhlxU,364
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/coil.vol,sha256=IOZe4oPxMUa9vBGJWN4LG9M79XAruY23S2rtAh_2_98,144724
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/coilshield.geo,sha256=lCUDCAhfowBttSGEY_giyY0xEFInAbgVfIPwSWBb504,605
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/coilshield.vol,sha256=2Dq5VCX7dqGxo_ExtUhI716I-SnxhLGPawW29IR9nYY,162397
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/cube.geo,sha256=kbBNVUYVSf1ogV7-Xf31KXju3LCUIPaDLfvwZHMgddI,386
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/cube.vol,sha256=54108a3RvrTUOJSFltxUY8iiNzyWoBTtNrWM0kufj7c,90295
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/d10_DGdoubleglazing.pde,sha256=bPQAvbYlyvCfLzk9cOzjXDn94UFC_vBCONl7fEIaMVs,1304
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/d11_chip_nitsche.pde,sha256=Mq3Ap63mWqjQyXsASFAlkRQ0AWxBlY3BYYejTYd-lwg,1101
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/d1_square.pde,sha256=ysOroxhoAOePmB_Aswps9iGjA2qV5435iu3hOnSvN8M,1029
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/d2_chip.pde,sha256=AowiGbuAWEdyFUB3d1LlZeQ7rarbkqIYsj1QcWc89-A,950
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/d3_helmholtz.pde,sha256=nG0HxrxboI2qGoR5_LNZt7RYB7ji2o5jtkyIqTFR4y8,568
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/d4_cube.pde,sha256=-C4RgoVtBXZu7XHCwZy7F7pr_PGNGcRwWCFzBppq1mQ,1033
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/d5_beam.pde,sha256=QOYC7YN7LdcjSFeLdJu1Tlfi-lmsVsa7b9G6qoCsGo8,1985
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/d6_shaft.pde,sha256=Iua3INfuIUeaBk21qNg4u5KZt__H1B8BZUsy5frjhdY,1901
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/d7_coil.pde,sha256=-ELH2obynMJcvSOqoYT1OThXcllIcdss0ZrklLSimZw,1200
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/d8_coilshield.pde,sha256=WUUulvd3dsedB25sFVQgU7FgZ3B7BkbxgkXTEBglk1w,1141
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/d9_hybridDG.pde,sha256=drLxnz1GvCEeme-FV772cGbS-YTgeUB4G5DrbkLdapw,1479
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/doubleglazing.in2d,sha256=bKn7m7eCcoLimGgoO3nxhiu7Mm0lIRnbx0TJh5dC1I0,506
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/doubleglazing.vol,sha256=Eg-H6nl2Hd5Nr593SFCCLVTDl25kIXuddrWpd7xlrP8,34987
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/piezo2d40round4.vol.gz,sha256=RuxM0vtYljoAImOH6U6avt0LmWaUbGNao37n8p_Yd1w,496477
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/shaft.geo,sha256=aPbk4EItJSLN_PjpOejnijk74y53QBW-IPCXHvLNlVs,2434
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/shaft.vol,sha256=fT8qu1L-p0MkI-l3WVisKG52W5FCzO_NzXaCvMJT1V4,270450
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/square.in2d,sha256=evpVgNNzReABv6bIs4J_r1OFeedMppoS-TEzE-tfKgY,192
+ngsolve-6.2.2505.post17.dev0.data/data/share/ngsolve/square.vol,sha256=IXIHYmq85hkFsMVZ_vEVthUnd0YTQWE4jRegjW8_ImI,9162
+ngsolve-6.2.2505.post17.dev0.dist-info/LICENSE,sha256=oZDcnIBDdV2Q-LCnX6ZrnkLUr0yYC_XdxjPwEk2zzuc,26430
+ngsolve-6.2.2505.post17.dev0.dist-info/METADATA,sha256=m_QtjpNkzWhWCiabUf-08HoW2sXYVcMhVM6ct2l8hHY,247
+ngsolve-6.2.2505.post17.dev0.dist-info/WHEEL,sha256=vYGgz4rYivjGPkIS_2Kj6g97rQxH536HeP_jo4e1Oq8,110
+ngsolve-6.2.2505.post17.dev0.dist-info/top_level.txt,sha256=1llU44DxNIP85uPT4RChsIuV3h2d0_5j0Ui_Df5vV8E,15
+ngsolve-6.2.2505.post17.dev0.dist-info/RECORD,,