ngsolve 6.2.2501.post83.dev1__cp39-cp39-win_amd64.whl → 6.2.2502__cp39-cp39-win_amd64.whl

netgen/include/facetfespace.hpp

@@ -88,6 +88,8 @@ namespace ngcomp
     using FESpace::GetDofNrs;
     virtual void GetDofNrs (ElementId ei, Array<DofId> & dnums) const override;
 
+    const auto & GetFirstFacetDof() const { return first_facet_dof; }
+    
     IntRange GetFacetDofs (int nr) const
     { 
       return IntRange (first_facet_dof[nr], first_facet_dof[nr+1]); 

netgen/include/fespace.hpp

@@ -92,6 +92,7 @@ ANY                  1 1 1 1 | 15
 
 
   class FESpace;
+  class BilinearForm;
 
   // will be size_t some day 
   typedef int DofId;
@@ -639,7 +640,7 @@ ANY                  1 1 1 1 | 15
     /// Set multigrid prolongation
     // void SetProlongation (ngmg::Prolongation * aprol)
     // { prol = aprol; }
-
+    virtual void SetHarmonicProlongation (shared_ptr<BilinearForm> bfa);
 
     /// returns function-evaluator
     shared_ptr<DifferentialOperator> GetEvaluator (VorB vb = VOL) const

netgen/include/mptools.hpp

@@ -14,23 +14,21 @@ namespace ngfem
 {
 
 
-  
+  template <typename entry_type = Complex>
   class NGS_DLL_HEADER SphericalHarmonics
   {
     int order;
-    Vector<Complex> coefs;
+    Vector<entry_type> coefs;
 
   public:
     SphericalHarmonics (int aorder)
       : order(aorder), coefs(sqr(order+1)) { coefs=0.0; }
 
     int Order() const { return order; }
-    FlatVector<Complex> Coefs() const { return coefs; }
-
-    Complex & Coef(int n, int m)
-    {
-      return coefs(n*(n+1) + m);  
-    }
+    FlatVector<entry_type> Coefs() const { return coefs; }
+    
+    entry_type & Coef(int n, int m) { return coefs(n*(n+1) + m); }
+    entry_type Coef(int n, int m) const { return coefs(n*(n+1) + m); }    
 
     auto CoefsN (int n) const
     {
@@ -53,29 +51,29 @@ namespace ngfem
       return { theta, phi };
     }
     
-    Complex Eval (Vec<3> x) const
+    entry_type Eval (Vec<3> x) const
     {
       auto [theta, phi] = Polar(x);
       return Eval(theta, phi);
     }
   
-    Complex Eval (double theta, double phi) const;
+    entry_type Eval (double theta, double phi) const;
     
-    Complex EvalOrder (int n, Vec<3> x) const
+    entry_type EvalOrder (int n, Vec<3> x) const
     {
       auto [theta, phi] = Polar (x);
       return EvalOrder(n, theta, phi);
     }
   
-    Complex EvalOrder (int n, double theta, double phi) const;
+    entry_type EvalOrder (int n, double theta, double phi) const;
 
-    void EvalOrders (Vec<3> x, FlatVector<Complex> vals) const
+    void EvalOrders (Vec<3> x, FlatVector<entry_type> vals) const
     {
       auto [theta, phi] = Polar(x);
       return EvalOrders(theta, phi, vals);
     }
   
-    void EvalOrders (double theta, double phi, FlatVector<Complex> vals) const;
+    void EvalOrders (double theta, double phi, FlatVector<entry_type> vals) const;
     
     void Calc (Vec<3> x, FlatVector<Complex> shapes);
 
@@ -106,7 +104,7 @@ namespace ngfem
     
     // Nail A. Gumerov and Ramani Duraiswami book, formula (2.2.12)
     // add directional derivative divided by kappa to res, both multipoles need same scaling
-    void DirectionalDiffAdd (Vec<3> d, SphericalHarmonics & res, double scale = 1);
+    void DirectionalDiffAdd (Vec<3> d, SphericalHarmonics<entry_type> & res, double scale = 1) const;
 
   };
 
@@ -234,24 +232,24 @@ namespace ngfem
   
 
 
-  template <typename RADIAL>
+  template <typename RADIAL, typename entry_type=Complex>
   class NGS_DLL_HEADER MultiPole
   {
-    SphericalHarmonics sh;
+    SphericalHarmonics<entry_type> sh;
     double kappa;
     double scale;
   public:
     MultiPole (int aorder, double akappa, double ascale = 1) 
       : sh(aorder), kappa(akappa), scale(ascale) { }
 
-    Complex & Coef(int n, int m) { return sh.Coef(n,m); }
+    entry_type & Coef(int n, int m) { return sh.Coef(n,m); }
     auto & SH() { return sh; }
     const auto & SH() const { return sh; }
     double Kappa() const { return kappa; }
     double Scale() const { return scale; }    
     int Order() const { return sh.Order(); }
     
-    MultiPole<RADIAL> Truncate(int neworder) const
+    MultiPole Truncate(int neworder) const
     {
       if (neworder > sh.Order()) neworder=sh.Order();
       MultiPole nmp(neworder, kappa);
@@ -266,10 +264,12 @@ namespace ngfem
       return *this;
     }
     
-    Complex Eval (Vec<3> x) const;
+    entry_type Eval (Vec<3> x) const;
+    entry_type EvalDirectionalDerivative (Vec<3> x, Vec<3> d) const;
 
-    void AddCharge (Vec<3> x, Complex c);
-    void AddDipole (Vec<3> x, Vec<3> d, Complex c);
+    void AddCharge (Vec<3> x, entry_type c);
+    void AddDipole (Vec<3> x, Vec<3> d, entry_type c);
+    void AddCurrent (Vec<3> ap, Vec<3> ep, Complex j, int num=100);
     
     
     void ChangeScaleTo (double newscale)
@@ -295,7 +295,7 @@ namespace ngfem
 
     
     template <typename TARGET>
-    void Transform (MultiPole<TARGET> & target, Vec<3> dist) const
+    void Transform (MultiPole<TARGET,entry_type> & target, Vec<3> dist) const
     {
       if (target.SH().Order() < 0) return;
       if (SH().Order() < 0)
@@ -321,7 +321,7 @@ namespace ngfem
         phi = atan2(dist(1), dist(0));
         
       
-      MultiPole<RADIAL> tmp(*this);
+      MultiPole<RADIAL,entry_type> tmp(*this);
       tmp.SH().RotateZ(phi);
       tmp.SH().RotateY(theta);
 
@@ -332,18 +332,18 @@ namespace ngfem
     }
     
     template <typename TARGET>
-    void TransformAdd (MultiPole<TARGET> & target, Vec<3> dist) const
+    void TransformAdd (MultiPole<TARGET,entry_type> & target, Vec<3> dist) const
     {
       if (SH().Order() < 0) return;
       if (target.SH().Order() < 0) return;      
       
-      MultiPole<TARGET> tmp{target};
+      MultiPole<TARGET,entry_type> tmp{target};
       Transform(tmp, dist);
       target.SH().Coefs() += tmp.SH().Coefs();
     }
 
     template <typename TARGET>
-    void ShiftZ (double z, MultiPole<TARGET> & target);
+    void ShiftZ (double z, MultiPole<TARGET,entry_type> & target);
     
   };
   
@@ -351,12 +351,14 @@ namespace ngfem
 
   // ***************** parameters ****************
 
-  static int MPOrder (double rho_kappa)
+  static constexpr int MPOrder (double rho_kappa)
   {
     return max (20, int(2*rho_kappa));
   }
   static constexpr int maxdirect = 100;
 
+
+  template <typename entry_type=Complex>
   class SingularMLMultiPole
   {
     static Array<size_t> nodes_on_level;    
@@ -367,10 +369,11 @@ namespace ngfem
       double r;
       int level;
       std::array<unique_ptr<Node>,8> childs;
-      MultiPole<MPSingular> mp;
+      MultiPole<MPSingular, entry_type> mp;
 
-      Array<tuple<Vec<3>, Complex>> charges;
-      Array<tuple<Vec<3>, Vec<3>, Complex>> dipoles;
+      Array<tuple<Vec<3>, entry_type>> charges;
+      Array<tuple<Vec<3>, Vec<3>, entry_type>> dipoles;
+      Array<tuple<Vec<3>, Vec<3>, Complex,int>> currents;
       int total_sources;
       
       Node (Vec<3> acenter, double ar, int alevel, int order, double kappa)
@@ -397,7 +400,7 @@ namespace ngfem
       }
       
 
-      void AddCharge (Vec<3> x, Complex c)
+      void AddCharge (Vec<3> x, entry_type c)
       {
         if (childs[0])
           {
@@ -422,13 +425,16 @@ namespace ngfem
           AddCharge (x,c);
         for (auto [x,d,c] : dipoles)
           AddDipole (x,d,c);
+        for (auto [sp,ep,j,num] : currents)
+          AddCurrent (sp,ep,j,num);
         
         charges.SetSize0();
-        dipoles.SetSize0();        
+        dipoles.SetSize0();
+        currents.SetSize0();
       }
 
 
-      void AddDipole (Vec<3> x, Vec<3> d, Complex c)
+      void AddDipole (Vec<3> x, Vec<3> d, entry_type c)
       {
         if (childs[0])
           {
@@ -453,15 +459,70 @@ namespace ngfem
           AddCharge (x,c);
         for (auto [x,d,c] : dipoles)
           AddDipole (x,d,c);
+        for (auto [sp,ep,j,num] : currents)
+          AddCurrent (sp,ep,j,num);
 
         charges.SetSize0();
         dipoles.SetSize0();        
+        currents.SetSize0();        
+      }
+
+      void AddCurrent (Vec<3> sp, Vec<3> ep, Complex j, int num)
+      {
+        if (childs[0])
+          {
+            // split line and send to childs
+            Array<double> split;
+            split.Append(0);
+            for (int i = 0; i < 3; i++)
+              if (sp(i) < center(i) != ep(i) < center(i))
+                split += (center(i)-sp(i)) / (ep(i)-sp(i));  // segment cuts i-th coordinate plane
+            split.Append(1);
+            BubbleSort(split);
+
+            for (int i = 0; i < split.Size()-1; i++)
+              if (split[i+1] > split[i])
+                {
+                  Vec<3> spi = sp + split[i]*(ep-sp);
+                  Vec<3> epi = sp + split[i+1]*(ep-sp);
+                  
+                  Vec<3> x = 0.5*(spi+epi);
+                  
+                  int childnum  = 0;
+                  if (x(0) > center(0)) childnum += 1;
+                  if (x(1) > center(1)) childnum += 2;
+                  if (x(2) > center(2)) childnum += 4;
+                  childs[childnum] -> AddCurrent(spi, epi, j, num);
+                }
+            return;
+          }
+        
+        currents.Append (tuple{sp,ep,j,num});
+
+        // if (currents.Size() < maxdirect || r < 1e-8)
+        if (currents.Size() < 4 || r < 1e-8)        
+          return;
+        
+        CreateChilds();
+
+        for (auto [x,c] : charges)
+          AddCharge (x,c);
+        for (auto [x,d,c] : dipoles)
+          AddDipole (x,d,c);
+        for (auto [sp,ep,j,num] : currents)
+          AddCurrent (sp,ep,j,num);
+
+        charges.SetSize0();
+        dipoles.SetSize0();
+        currents.SetSize0();
       }
 
+
       
-      Complex Evaluate(Vec<3> p) const
+      
+      entry_type Evaluate(Vec<3> p) const
       {
-        Complex sum = 0;
+        entry_type sum{0.0};
         if (childs[0])
           {
             for (auto & child : childs)
@@ -471,17 +532,63 @@ namespace ngfem
 
         for (auto [x,c] : charges)
           if (double rho = L2Norm(p-x); rho > 0)
-            sum += c*(1/(4*M_PI))*exp(Complex(0,rho*mp.Kappa())) / rho;
+            sum += (1/(4*M_PI))*exp(Complex(0,rho*mp.Kappa())) / rho * c;
         
         for (auto [x,d,c] : dipoles)
           if (double rho = L2Norm(p-x); rho > 0)
             {
               Vec<3> drhodp = 1.0/rho * (p-x);
-              Complex dGdrho = c*(1/(4*M_PI))*exp(Complex(0,rho*mp.Kappa())) *
+              Complex dGdrho = (1/(4*M_PI))*exp(Complex(0,rho*mp.Kappa())) *
                 (Complex(0, mp.Kappa())/rho - 1.0/sqr(rho));
-              sum += dGdrho * InnerProduct(drhodp, d);
+              sum += dGdrho * InnerProduct(drhodp, d) * c;
             }
 
+        for (auto [sp,ep,j,num] : currents)
+          {
+            // should use explizit formula instead ...
+            
+            Vec<3> tau = ep-sp;
+            Vec<3> tau_num = 1.0/num *  tau;
+            for (int i = 0; i < num; i++)
+              {
+                Vec<3> x = sp+(i+0.5)*tau_num;
+                
+                if (double rho = L2Norm(p-x); rho > 0)
+                  {
+                    Vec<3> drhodp = 1.0/rho * (p-x);
+                    Complex dGdrho = (1/(4*M_PI))*exp(Complex(0,rho*mp.Kappa())) *
+                      (Complex(0, mp.Kappa())/rho - 1.0/sqr(rho));
+
+                    if constexpr (std::is_same<entry_type, Vec<3,Complex>>())
+                      sum += j*dGdrho * Cross(drhodp, tau_num);
+                  }
+              }
+          }
+        
+        return sum;
+      }
+
+      entry_type EvaluateDeriv(Vec<3> p, Vec<3> d) const
+      {
+        entry_type sum = 0;
+        if (childs[0])
+        {
+            for (auto & child : childs)
+              sum += child->EvaluateDeriv(p, d);
+            return sum;
+        }
+
+        if (dipoles.Size())
+            throw Exception("EvaluateDeriv not implemented for dipoles in SingularMLMultiPole");
+
+        for (auto [x,c] : charges)
+          if (double rho = L2Norm(p-x); rho > 0)
+          {
+            Vec<3> drhodp = 1.0/rho * (p-x);
+            Complex dGdrho = (1/(4*M_PI))*exp(Complex(0,rho*mp.Kappa())) *
+            (Complex(0, mp.Kappa())/rho - 1.0/sqr(rho));
+            sum += dGdrho * InnerProduct(drhodp, d) * c;
+          }
         return sum;
       }
 
@@ -516,9 +623,9 @@ namespace ngfem
           }
         else
           {
-            if (charges.Size()+dipoles.Size() == 0)
+            if (charges.Size()+dipoles.Size()+currents.Size() == 0)
               {
-                mp = MultiPole<MPSingular> (-1, mp.Kappa());
+                mp = MultiPole<MPSingular,entry_type> (-1, mp.Kappa());
                 return;
               }
 
@@ -527,10 +634,13 @@ namespace ngfem
             
             for (auto [x,d,c] : dipoles)
               mp.AddDipole (x-center, d, c);
+
+            for (auto [sp,ep,j,num] : currents)
+              mp.AddCurrent (sp-center, ep-center, j, num);
           }
       }
       
-      Complex EvaluateMP(Vec<3> p) const
+      entry_type EvaluateMP(Vec<3> p) const
       {
         if (charges.Size() || dipoles.Size())
           return Evaluate(p);
@@ -541,22 +651,43 @@ namespace ngfem
         if (!childs[0]) //  || level==1)
           return Evaluate(p);
           
-        Complex sum = 0.0;
+        entry_type sum{0.0};
         for (auto & child : childs)
           sum += child->EvaluateMP(p);
         return sum;
       }
-      
-      
-      void Print (ostream & ost) const
+
+      entry_type EvaluateMPDeriv(Vec<3> p, Vec<3> d) const
       {
-        ost << "c = " << center << ", r = " << r << endl;
+        // cout << "EvaluateMPDeriv Singular, p = " << p << ", d = " << d << ", r = " << r << ", center = " << center <<  endl;
+        // cout << "Norm: " << L2Norm(p-center) << " > " << 3*r << endl;
+        // cout << "charges.Size() = " << charges.Size() << ", dipoles.Size() = " << dipoles.Size() << endl;
+        if (charges.Size() || dipoles.Size() || !childs[0])
+          return EvaluateDeriv(p, d);
+
+        if (L2Norm(p-center) > 3*r)
+          return mp.EvalDirectionalDerivative(p-center, d);
+
+        entry_type sum{0.0};
+        for (auto & child : childs)
+          sum += child->EvaluateMPDeriv(p, d);
+        return sum;
+      }
+
+      void Print (ostream & ost, size_t childnr = -1) const
+      {
+        if (childnr == -1)
+          ost << "c = " << center << ", r = " << r << ", level = " << level << endl;
+        else
+          ost << "c = " << center << ", r = " << r << ", level = " << level << ", childnr = " << childnr << endl;
         // for (int i = 0; i < loc_pnts.Size(); i++)
         for (auto [x,c] : charges)
           ost << "xi = " << x << ", ci = " << c << endl;
+        for (auto [x,d,c] : dipoles)
+          ost << "xi = " << x << ", di = " << d << ", ci = " << c << endl;
 
         for (int i = 0; i < 8; i++)
-          if (childs[i]) childs[i] -> Print (ost);
+          if (childs[i]) childs[i] -> Print (ost, i);
       }
 
       double Norm () const
@@ -601,6 +732,31 @@ namespace ngfem
       root.AddDipole(x, d, c);
     }
 
+    void AddCurrent (Vec<3> sp, Vec<3> ep, Complex j, int num)
+    {
+      if constexpr (!std::is_same<entry_type, Vec<3,Complex>>())
+        throw Exception("AddCurrent needs a singular vectorial MP");
+      
+      root.AddCurrent (sp, ep, j, num);
+      /*
+         // for testing
+      Vec<3> tau = ep-sp;
+      Vec<3> tau_num = 1.0/num *  tau;
+      for (int i = 0; i < num; i++)
+        {
+          for (int k = 0; k < 3; k++)
+            {
+              Vec<3> ek{0.0}; ek(k) = 1;
+              Vec<3> cp = Cross(tau, ek);
+              Vec<3,Complex> source{0.0};
+              source(k) = j/double(num);
+              if constexpr (std::is_same<entry_type, Vec<3,Complex>>())
+                root.AddDipole (sp+(i+0.5)*tau_num, cp, source);
+            }
+        }
+      */
+    }
+    
     void Print (ostream & ost) const
     {
       root.Print(ost);
@@ -634,7 +790,7 @@ namespace ngfem
       havemp = true;
     }
 
-    Complex Evaluate (Vec<3> p) const
+    entry_type Evaluate (Vec<3> p) const
     {
       if (havemp)
         return root.EvaluateMP(p);
@@ -642,17 +798,20 @@ namespace ngfem
         return root.Evaluate(p);
     }
 
+    template <typename entry_type2>
     friend class RegularMLMultiPole;
   };
 
 
-  inline ostream & operator<< (ostream & ost, const SingularMLMultiPole & mlmp)
+  template <typename entry_type>
+  inline ostream & operator<< (ostream & ost, const SingularMLMultiPole<entry_type> & mlmp)
   {
     mlmp.Print(ost);
     return ost;
   }
 
 
+  template <typename elem_type=Complex>
   class NGS_DLL_HEADER RegularMLMultiPole
   {
     static Array<size_t> nodes_on_level;
@@ -663,11 +822,11 @@ namespace ngfem
       double r;
       int level;
       std::array<unique_ptr<Node>,8> childs;
-      MultiPole<MPRegular> mp;
+      MultiPole<MPRegular,elem_type> mp;
       Array<Vec<3>> targets;
       int total_targets;
 
-      Array<const SingularMLMultiPole::Node*> singnodes;
+      Array<const typename SingularMLMultiPole<elem_type>::Node*> singnodes;
 
       Node (Vec<3> acenter, double ar, int alevel, int order, double kappa)
         : center(acenter), r(ar), level(alevel), mp(MPOrder(ar*kappa), kappa, 1.0/min(1.0, 0.25*r*kappa))
@@ -692,7 +851,7 @@ namespace ngfem
           }
       }
 
-      void AddSingularNode (const SingularMLMultiPole::Node & singnode, bool allow_refine)
+      void AddSingularNode (const typename SingularMLMultiPole<elem_type>::Node & singnode, bool allow_refine)
       {
         if (mp.SH().Order() < 0) return;
         if (singnode.mp.SH().Order() < 0) return;
@@ -781,15 +940,15 @@ namespace ngfem
                   mp.TransformAdd (ch->mp, ch->center-center);
                 ch->LocalizeExpansion(allow_refine);
               }
-            mp = MultiPole<MPRegular>(-1, mp.Kappa());
+            mp = MultiPole<MPRegular,elem_type>(-1, mp.Kappa());
             //mp.SH().Coefs()=0.0;
           }
       }
       
-      Complex Evaluate (Vec<3> p) const
+      elem_type Evaluate (Vec<3> p) const
       {
         // *testout << "eval p = " << p << ", level = " << level << ", center = " << center <<  ", r = " << r << endl;
-        Complex sum = 0.0;
+        elem_type sum{0.0};
         /*
         if (childs[0])
           {
@@ -817,6 +976,28 @@ namespace ngfem
         return sum;
       }
 
+      elem_type EvaluateDirectionalDerivative (Vec<3> p, Vec<3> d) const
+      {
+        elem_type sum{0.0};
+        // cout << "EvaluateDirectionalDerivative RegularMLMP, r = " << r << ", level = " << level << ", center = " << center << endl;
+        // cout << "Singnodes: " << singnodes.Size() << ", childs: " << childs[0] << endl;
+
+        int childnum = 0;
+        if (p(0) > center(0)) childnum += 1;
+        if (p(1) > center(1)) childnum += 2;
+        if (p(2) > center(2)) childnum += 4;
+        if (childs[childnum])
+          sum = childs[childnum]->EvaluateDirectionalDerivative(p, d);
+        else
+          sum = mp.EvalDirectionalDerivative(p-center, d);
+
+        static Timer t("mptool direct evaluate deriv"); RegionTimer r(t);
+        for (auto sn : singnodes)
+          sum += sn->EvaluateMPDeriv(p, d);
+
+        return sum;
+      }
+
       double Norm() const
       {
         double norm = L2Norm(mp.SH().Coefs());
@@ -886,13 +1067,27 @@ namespace ngfem
           mp = MultiPole<MPRegular>(-1, mp.Kappa());
       }
       
+
+      void Print (ostream & ost, size_t childnr = -1) const
+      {
+        if (childnr == -1)
+          ost << "c = " << center << ", r = " << r << ", level = " << level << endl;
+        else
+          ost << "c = " << center << ", r = " << r << ", level = " << level << ", childnr = " << childnr << endl;
+        for (auto x : targets)
+          ost << "xi = " << x << endl;
+
+        for (int i = 0; i < 8; i++)
+          if (childs[i]) childs[i] -> Print (ost, i);
+      }
+
     };
     
     Node root;
-    shared_ptr<SingularMLMultiPole> singmp;
+    shared_ptr<SingularMLMultiPole<elem_type>> singmp;
     
   public:
-    RegularMLMultiPole (shared_ptr<SingularMLMultiPole> asingmp, Vec<3> center, double r, int order)
+  RegularMLMultiPole (shared_ptr<SingularMLMultiPole<elem_type>> asingmp, Vec<3> center, double r, int order)
       : root(center, r, 0, order, asingmp->Kappa()), singmp(asingmp)
     {
       if (!singmp->havemp) throw Exception("first call Calc for singular MP");
@@ -934,7 +1129,7 @@ namespace ngfem
       root.AddTarget (t);
     }
 
-    void CalcMP(shared_ptr<SingularMLMultiPole> asingmp)
+    void CalcMP(shared_ptr<SingularMLMultiPole<elem_type>> asingmp)
     {
       singmp = asingmp;
 
@@ -955,6 +1150,11 @@ namespace ngfem
       root.LocalizeExpansion(false);
     }
 
+    void Print (ostream & ost) const
+    {
+      root.Print(ost);
+    }
+
     double Norm() const
     {
       return root.Norm();
@@ -965,21 +1165,26 @@ namespace ngfem
       return root.NumCoefficients();
     }
 
-    Complex Evaluate (Vec<3> p) const
+    elem_type Evaluate (Vec<3> p) const
     {
       // static Timer t("mptool Eval MLMP regular"); RegionTimer r(t);
-      if (L2Norm(p-root.center) > root.r) return 0.0;
+      if (L2Norm(p-root.center) > root.r) return elem_type{0.0};
       return root.Evaluate(p);
     }
 
-    
+    elem_type EvaluateDirectionalDerivative (Vec<3> p, Vec<3> d) const
+    {
+        if (L2Norm(p-root.center) > root.r) return elem_type{0.0};
+        return root.EvaluateDirectionalDerivative(p, d);
+    }
+
   };
-  
 
-  inline ostream & operator<< (ostream & ost, const RegularMLMultiPole & mlmp)
+  template <typename elem_type>
+  inline ostream & operator<< (ostream & ost, const RegularMLMultiPole<elem_type> & mlmp)
   {
-    // mlmp.Print(ost);
-    ost << "RegularMLMultiPole" << endl;
+    mlmp.Print(ost);
+    // ost << "RegularMLMultiPole" << endl;
     return ost;
   }
 
@@ -993,7 +1198,7 @@ namespace ngfem
   
   class SphericalHarmonicsCF : public CoefficientFunction
   {
-    SphericalHarmonics sh;
+    SphericalHarmonics<Complex> sh;
   public:
     SphericalHarmonicsCF (int order)
       : CoefficientFunction(1, true), sh(order) { }
@@ -1021,16 +1226,16 @@ namespace ngfem
 
 
 
-  template <typename RADIAL>
+  template <typename RADIAL, typename entry_type=Complex>
   class MultiPoleCF : public CoefficientFunction
   {
-    MultiPole<RADIAL> mp;
+    MultiPole<RADIAL, entry_type> mp;
     Vec<3> center;
   public:
     MultiPoleCF (int order, double kappa, Vec<3> acenter, double scale = 1)
-      : CoefficientFunction(1, true), mp(order, kappa, scale), center(acenter) { }
+      : CoefficientFunction(sizeof(entry_type)/sizeof(Complex), true), mp(order, kappa, scale), center(acenter) { }
 
-    Complex & Coef(int n, int m) { return mp.Coef(n,m); } 
+    entry_type & Coef(int n, int m) { return mp.Coef(n,m); } 
     auto & SH() { return mp.SH(); }
     auto & MP() { return mp; }
     Vec<3> Center() const { return center; }
@@ -1040,56 +1245,72 @@ namespace ngfem
 
     virtual void Evaluate (const BaseMappedIntegrationPoint & mip, FlatVector<Complex> values) const override
     {
-      values(0) = mp.Eval(mip.GetPoint()-center);
+      if constexpr (std::is_same<entry_type, Complex>())
+        values(0) = mp.Eval(mip.GetPoint()-center);
+      else
+        values = mp.Eval(mip.GetPoint()-center);
     }
 
     template <typename TARGET>
-    void ShiftZ (double z, MultiPole<TARGET> & target) { mp.ShiftZ(z, target); }
+    void ShiftZ (double z, MultiPole<TARGET, entry_type> & target) { mp.ShiftZ(z, target); }
 
     using CoefficientFunction::Transform;        
     template <typename TARGET>
-    void Transform (MultiPoleCF<TARGET> & target)
+    void Transform (MultiPoleCF<TARGET, entry_type> & target)
     {
       mp.Transform (target.MP(), target.Center()-center);
     }
   };
 
-  
+  template <typename entry_type>
   class SingularMLMultiPoleCF : public CoefficientFunction
   {
-    shared_ptr<SingularMLMultiPole> mlmp;
+    shared_ptr<SingularMLMultiPole<entry_type>> mlmp;
   public:
     SingularMLMultiPoleCF (Vec<3> center, double r, int order, double kappa)
-      : CoefficientFunction(1, true), mlmp{make_shared<SingularMLMultiPole>(center, r, order, kappa)} { }
+      : CoefficientFunction(sizeof(entry_type)/sizeof(Complex), true), mlmp{make_shared<SingularMLMultiPole<entry_type>>(center, r, order, kappa)} { }
     
     virtual double Evaluate (const BaseMappedIntegrationPoint & ip) const override
     { throw Exception("real eval not available"); }
     
     virtual void Evaluate (const BaseMappedIntegrationPoint & mip, FlatVector<Complex> values) const override
     {
-      values(0) = mlmp->Evaluate(mip.GetPoint());
+      // values(0) = mlmp->Evaluate(mip.GetPoint());
+
+      if constexpr (std::is_same<entry_type, Complex>())
+        values(0) = mlmp->Evaluate(mip.GetPoint());
+      else
+        values = mlmp->Evaluate(mip.GetPoint());
+
+      
     }
     
-    shared_ptr<SingularMLMultiPole> MLMP() { return mlmp; }
+    shared_ptr<SingularMLMultiPole<entry_type>> MLMP() { return mlmp; }
   };
   
 
+  template <typename entry_type>
   class RegularMLMultiPoleCF : public CoefficientFunction
   {
-    shared_ptr<RegularMLMultiPole> mlmp;
+    shared_ptr<RegularMLMultiPole<entry_type>> mlmp;
   public:
-    RegularMLMultiPoleCF (shared_ptr<SingularMLMultiPoleCF> asingmp, Vec<3> center, double r, int order)
-      : CoefficientFunction(1, true), mlmp{make_shared<RegularMLMultiPole>(asingmp->MLMP(), center, r, order)} { } 
+    RegularMLMultiPoleCF (shared_ptr<SingularMLMultiPoleCF<entry_type>> asingmp, Vec<3> center, double r, int order)
+      : CoefficientFunction(sizeof(entry_type)/sizeof(Complex), true), mlmp{make_shared<RegularMLMultiPole<entry_type>>(asingmp->MLMP(), center, r, order)} { } 
     
     virtual double Evaluate (const BaseMappedIntegrationPoint & ip) const override
     { throw Exception("real eval not available"); }
 
     virtual void Evaluate (const BaseMappedIntegrationPoint & mip, FlatVector<Complex> values) const override
     {
-      values(0) = mlmp->Evaluate(mip.GetPoint());
+      // values(0) = mlmp->Evaluate(mip.GetPoint());
+
+      if constexpr (std::is_same<entry_type, Complex>())
+        values(0) = mlmp->Evaluate(mip.GetPoint());
+      else
+        values = mlmp->Evaluate(mip.GetPoint());
     }
 
-    shared_ptr<RegularMLMultiPole> MLMP() { return mlmp; }
+    shared_ptr<RegularMLMultiPole<entry_type>> MLMP() { return mlmp; }
   };
 
   

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/prolongation.hpp

@@ -8,7 +8,7 @@
 /*********************************************************************/
 
 namespace ngcomp { class NedelecFESpace; }
-
+namespace ngla { class SparseFactorization; } 
 namespace ngmg
 {
 
@@ -229,6 +229,25 @@ namespace ngmg
     ///
     virtual void RestrictInline (int finelevel, BaseVector & v) const override;
   };
+
+  class NGS_DLL_HEADER HarmonicProlongation : public Prolongation
+  {
+    shared_ptr<Prolongation> baseprol;
+    shared_ptr<BilinearForm> bfa;
+    Array<shared_ptr<BaseMatrix>> innerinverses;
+    Array<size_t> edges_on_level;
+  public:
+    HarmonicProlongation (shared_ptr<Prolongation> abaseprol,
+                          shared_ptr<BilinearForm> abfa);
+
+    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;
+  };
 }
 
 

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/tensor.hpp

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

ngsolve/cmake/NGSolveConfig.cmake

@@ -1,4 +1,4 @@
-set(PACKAGE_VERSION "6.2.2501-83-gaf867070b")
+set(PACKAGE_VERSION "6.2.2502")
 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.2501-83-gaf867070b"
-NGSOLVE_VERSION_GIT = "v6.2.2501-83-gaf867070b"
-NGSOLVE_VERSION_PYTHON = "6.2.2501.post83.dev1"
+NGSOLVE_VERSION = "6.2.2502"
+NGSOLVE_VERSION_GIT = "v6.2.2502-0-g984bf8ae2"
+NGSOLVE_VERSION_PYTHON = "6.2.2502"
 
 NGSOLVE_VERSION_MAJOR = "6"
 NGSOLVE_VERSION_MINOR = "2"
-NGSOLVE_VERSION_TWEAK = "83"
-NGSOLVE_VERSION_PATCH = "2501"
-NGSOLVE_VERSION_HASH = "gaf867070b"
+NGSOLVE_VERSION_TWEAK = "0"
+NGSOLVE_VERSION_PATCH = "2502"
+NGSOLVE_VERSION_HASH = "g984bf8ae2"
 
 CMAKE_CXX_COMPILER           = "C:/Program Files/Microsoft Visual Studio/2022/Community/VC/Tools/MSVC/14.44.35207/bin/Hostx64/x64/cl.exe"
 CMAKE_CUDA_COMPILER          = ""

{ngsolve-6.2.2501.post83.dev1.dist-info → ngsolve-6.2.2502.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: ngsolve
-Version: 6.2.2501.post83.dev1
+Version: 6.2.2502
 Summary: NGSolve
 Home-page: UNKNOWN
 Author: The NGSolve team

{ngsolve-6.2.2501.post83.dev1.dist-info → ngsolve-6.2.2502.dist-info}/RECORD

@@ -1,4 +1,4 @@
-netgen/libngsolve.dll,sha256=y0FOIn2p2t4DDhUuHltHK1xqB2s1pd33_CLCSbmYRs0,35603456
+netgen/libngsolve.dll,sha256=s-5HInrZFbSomtZ6iehm1wfCNV-Ao9-4x8QbP_-znt8,35779584
 netgen/ngscxx.bat,sha256=DjwBtuQlVFSBYmdu08f4H88Li7YN0uCu2LRLkd_MsLI,697
 netgen/ngsld.bat,sha256=YD1I1rK70K_YvLev7aYQMD9mCcwF55k1B_8pRWPpob0,428
 netgen/ngspy,sha256=RJeJY2bsimN2kO97US752sBKFTV-yv7sniRlz7Fy1j8,69
@@ -45,13 +45,13 @@ netgen/include/elementtransformation.hpp,sha256=H8eUIkDAxgeyeFL_cMJgDbQp00hRcwQP
 netgen/include/evalfunc.hpp,sha256=-I7SJFWJvLp8sFWYeOnNzkdDLaPLnWqAZSjRZHkLXfg,10424
 netgen/include/expr.hpp,sha256=8LlNBC9vk0duLRoa93ub3I60c2n_4Zbijr0VHNl0Yqo,46784
 netgen/include/facetfe.hpp,sha256=5weP0JK-l_UhPgHnRAx8G1y7arKZ_dKqo1qLUysPVjw,5095
-netgen/include/facetfespace.hpp,sha256=G3DQMNkkdACRNQqkTdj6X9bpoFUETDUtxkeYh7aM7ps,4842
+netgen/include/facetfespace.hpp,sha256=h4oBcysnGy8qKrWvg0K69VH-tmbcNxJ_pGH2gyCuzsA,4919
 netgen/include/facethofe.hpp,sha256=n2gCSrlazvCG0hs3MlNVfzIS83V0nYMAEbqholRwNxQ,3812
 netgen/include/facetsurffespace.hpp,sha256=qO4tLszTL3YHmWCsWZn9QKTW45GiIJ4UFnf5ySPWPuM,2521
 netgen/include/fe_interfaces.hpp,sha256=gGLCSJtPBuOCQZuo2t0VGvXrRwNt0OOPLRag--qXfS8,854
 netgen/include/fem.hpp,sha256=DPZpvRbOezchnVcUGlOE5R6xKFvmAqUeRza4GC285rw,2181
 netgen/include/fesconvert.hpp,sha256=hbGst9rksfUJX3BguHcgiF62myv3D3a2zkbjv2Ccl4g,574
-netgen/include/fespace.hpp,sha256=Y3VtWPInfrbvXzSad58Wkv4ybun3LOkFk42uhwxbtTQ,48313
+netgen/include/fespace.hpp,sha256=T0ZIfr2k7yj4wtXLJrqIgFSfx4C4R-1eS0KlBAvQeyI,48408
 netgen/include/finiteelement.hpp,sha256=ULds8D80WXlpVqgCppjHKrXhFuo-WewzIr3IalsIv8g,9431
 netgen/include/globalinterfacespace.hpp,sha256=atR009fWSlJpYjRI3OgMcvInOajz1W5egPLdfPAPFu8,2796
 netgen/include/globalspace.hpp,sha256=ZYni-ie6O2DD0N3u-fdmVfwygpV5lwJ3SOAnSSgWMNo,4079
@@ -107,13 +107,13 @@ netgen/include/matrix.hpp,sha256=0UT4qQMFBMf3eq8FS7XTkMbQaVSgz--EiNdVMzGq-Sc,579
 netgen/include/memusage.hpp,sha256=m_-fuupqepSve37WkRkoFHHyTtE9PdGnOVPIdQHlDiw,1209
 netgen/include/meshaccess.hpp,sha256=Yk5VSmQx4cp5WEHnfsZXUOijVlooEoozTAaUu75zla8,48419
 netgen/include/mgpre.hpp,sha256=5Zqt3pzrvm3EOfr84iU6rQx23LTkntUTd_s4_i6wEzI,5029
-netgen/include/mptools.hpp,sha256=PU8gcIVRawb1oDORD6xUzopSNel0hNIikqUlPPRpr5M,29607
+netgen/include/mptools.hpp,sha256=LUuPHuRv06i7UxENbUJcp3HQRFJuwoBfGT5nzw2C3EQ,38375
 netgen/include/multigrid.hpp,sha256=7YgJfeGTBnb-aFWHrbcVk3dos6XYU0iVCgiY8UvRopY,898
 netgen/include/multivector.hpp,sha256=5rkNSqyPhIXAb1lGKLcy908D9Uw9oqY4wpUXZEXbEPE,13082
 netgen/include/mumpsinverse.hpp,sha256=8YiUlQZMFVsJ0kuwaSphGZvxyEZJAZZApH8XahiTeEo,4249
 netgen/include/mycomplex.hpp,sha256=dj031oIJSDeCEB22UoP2FVv3FcNRsi6ssf_nTpHsVaw,8023
 netgen/include/ng_lapack.hpp,sha256=VFbQIzdfon96GEtyviYaqM8ul1tTiH3QyQM5U0Tviow,50069
-netgen/include/ngblas.hpp,sha256=JIkg45wZLeceZcpDCjumBBMjPDQjXF27CqCZ_nnhvQw,43407
+netgen/include/ngblas.hpp,sha256=ZwrjI7PqRMAr-Vg5CaswCMgVkxP6JjholA_E54zZd-g,43892
 netgen/include/ngs_defines.hpp,sha256=t7TFqQPQtm83a2NDH69TFOOtSQvFyJPbjaHXIP_XsxY,716
 netgen/include/ngs_stdcpp_include.hpp,sha256=GV74VWQLG5Z2unDAnjjGAWAheFg5ZtMH584pTGSjQAk,2657
 netgen/include/ngs_utils.hpp,sha256=yO0yJKC1q5BZP39tE2LabiG8pFV9mtNbrzE2b6g3-uU,3265
@@ -138,11 +138,11 @@ netgen/include/pmltrafo.hpp,sha256=j6Sk5tHbzKPT1hj_jI6n6-cRUgV5b3_o3kRNev95wDE,2
 netgen/include/postproc.hpp,sha256=TiEqDkANZJtxl-aOYMac-FKhtSkokXS3AuxAScqUDBQ,5334
 netgen/include/precomp.hpp,sha256=mx_zGsnBxqRVArwgkDpvUcX8EzP8v_8QwVdOp5m1nWQ,1337
 netgen/include/preconditioner.hpp,sha256=3UCrQqZ49YUHneg24-JoYVO40B4Di1DR-A3TQLZme2A,14630
-netgen/include/prolongation.hpp,sha256=eguG1dvLQTuo89qv2CDHP7bWx1TlYOQnv6QgVIaxFew,6615
+netgen/include/prolongation.hpp,sha256=JGW4mvgtPsqANsNZpTUeEMJuSp9eh07M_Whmoa2qxAQ,7468
 netgen/include/python_comp.hpp,sha256=FgujlcVk9FKWwMBEo9CfJ-q_Kb4TNcCziotVTFIuXvY,3571
 netgen/include/python_fem.hpp,sha256=pdSB1RLhv2Vi52TaRoGpjQ21EHOA3RFAF0jeiWDp6EQ,3161
 netgen/include/python_linalg.hpp,sha256=HdWTpdddyVmv3ZO73sG0BoHleKiREg9emIcyjFl3au8,1761
-netgen/include/python_ngstd.hpp,sha256=Z-6gEclm77s7Yzn7cFLAWv0cd_BsCam5aFq1zx6iLbo,10771
+netgen/include/python_ngstd.hpp,sha256=8NptP1VGDLDZ6gRC2kbayi-fHVg1IRNCZhgalJ4LUQ4,10803
 netgen/include/recursive_pol.hpp,sha256=u7hH83vY9w-F2zXh5jpqYuMkTNwMxOwp7xUy0f2j_KY,123587
 netgen/include/recursive_pol_tet.hpp,sha256=TKnMfmn1K7XvmNRD8kmLxPlooLRK6irBMKOwvRowsn4,9863
 netgen/include/recursive_pol_trig.hpp,sha256=UylKioGKIXXgSlSEbD4XPnzhJvP24JwD3bBYvFppHpg,11333
@@ -153,7 +153,7 @@ netgen/include/shapefunction_utils.hpp,sha256=um-JOAFzOJVP-zW-lrvYR6ojDblocO48H0
 netgen/include/simd_complex.hpp,sha256=HwMHrVpF8rqzTheKwrPUkZStgk9hTKG2-o2KUeISBkM,10010
 netgen/include/smoother.hpp,sha256=SoVZvhFynyPSIAe-x8EloVluIHUE0k78Dp3sQFu-Yfo,6896
 netgen/include/solve.hpp,sha256=lzCGHU1QOj9VXG43Q8MzvDi48-BfVH_fG8R0ttU8Xfg,2123
-netgen/include/sparsecholesky.hpp,sha256=5RYlQ6gbPTVqDy8npnDBydAF6pvkl8I4mHi_JjYJ6lM,9198
+netgen/include/sparsecholesky.hpp,sha256=hOzFoXMEelAdJROF1q4zkYv0yfZgSDhjwaYvRRmiMcE,9201
 netgen/include/sparsematrix.hpp,sha256=vyk1ViL0oXSiPZX4FfrFgy8whefffUrTmSlTS60SbrU,33744
 netgen/include/sparsematrix_dyn.hpp,sha256=RSxxEAg6Y82DDt7S46Hw_KCxvQBPvEoBvH8j-H5plz4,2792
 netgen/include/sparsematrix_impl.hpp,sha256=C2yLlIx3QPXCreg5GLBy2mV8E-1DUD1S8m3jzbWzlWQ,28732
@@ -166,7 +166,7 @@ netgen/include/symbolicintegrator.hpp,sha256=tk8SwRT2swPiafO1JevZskjFVwvAxbGd0HT
 netgen/include/symmetricmatrix.hpp,sha256=M79mq0z70ohLFoi8cxPAxxlqD3RjFZAkbCrHwQrhvY8,2954
 netgen/include/tangentialfacetfe.hpp,sha256=iEP_cU3nBtVUnH67dvzeMqIS-JM6SeILhB738YX9Spg,7491
 netgen/include/tangentialfacetfespace.hpp,sha256=DoijcYOUEG8Tud_RqYOsmOZlpsQ2swx0XTLgI6eGGG4,3198
-netgen/include/tensor.hpp,sha256=MEkTOpdqmb53A430k8ngtrtZvfd_sRhvxx9fwpDkOjo,12532
+netgen/include/tensor.hpp,sha256=nzrF84BknvEixYfMyN-RdaWVlBIXJkYw82kpjmfyqCE,12535
 netgen/include/tensorcoefficient.hpp,sha256=eA-SVbv_N89PnBLAEf-O--8SwicgiCwRPbx_NHQp-fw,17522
 netgen/include/tensorproductintegrator.hpp,sha256=Sfo_YsCL_TZnbqR2gBWcFLDr7DOlmSHAUJu-mCeYrI0,4643
 netgen/include/thcurlfe.hpp,sha256=48OMvHXQgn5QWu9AgA8ByPrIhdIrQoKD3z9uahlpx60,5099
@@ -185,7 +185,7 @@ netgen/include/voxelcoefficientfunction.hpp,sha256=zPlatX_nf2ODK2X979wvpm1HXHEzn
 netgen/include/vtkoutput.hpp,sha256=GwPaFtOVY27_1hKwSqpNHV4ppOk3CVsLdMpocm27odQ,6628
 netgen/include/vvector.hpp,sha256=vpnwWlNWsEJB26ve8cKADbUMGVBTTVCFvrBC4r_LDDQ,4785
 netgen/include/webgui.hpp,sha256=xm1S5KHockrzf_R4zHzQ9fzwfQ4H5aZaZm9ItCaR41s,2199
-netgen/lib/libngsolve.lib,sha256=auYbUlqkze8keSvy8EQrT_tLzSg5hblA945ZR6xyKM4,13366128
+netgen/lib/libngsolve.lib,sha256=NyQso40noaNdGKAII1SxjaB0ksz6KltH6Ul2ze2-C1Y,13475284
 ngsolve/TensorProductTools.py,sha256=jD4LtJdEiLt2dDsU6STMXb10yrau_i20gatwGlqbiFQ,8562
 ngsolve/__console.py,sha256=QLnfw1ONU9eVPQZ_LjzY-7b_gfZpRRRdI7imHR1lJR8,2983
 ngsolve/__expr.py,sha256=ZAtv94PVtAdvWWKFIvQGJC55aIEY6ajo-l5Hp-1V-6s,4396
@@ -200,7 +200,7 @@ ngsolve/meshes.py,sha256=nz95LP1EEIVAEJ0jyQGqq0E-62uuy2tvIRgqIrEqDHQ,27534
 ngsolve/ngs2petsc.py,sha256=PLAkvu4H-ZEXd0-bM2v0OJX7gga-OJ9d-NB7g322-wQ,12809
 ngsolve/ngscxx.py,sha256=r9TV49w0MKIxE6fgJYO4FNv-5wcM9eAFYIG9GspUy_A,1335
 ngsolve/ngslib.lib,sha256=_HzV4t-JUpCMSs1rNu4Jjtmr69v0FQUaOatAn17WF4A,1716
-ngsolve/ngslib.pyd,sha256=-_AkM3HOuVkdhB14e0dljDY5DeO56xpNeKXOMmFrQWk,121344
+ngsolve/ngslib.pyd,sha256=yJPNiC_ZYJtukc1zZwIwsaN87K4gY3bYjkqoGbBWmUw,121344
 ngsolve/nonlinearsolvers.py,sha256=nKd-SVJseTWztHZPEPKsbv0mLwQjFILZtOh7Ggft0xs,8419
 ngsolve/preconditioners.py,sha256=2To9cGdhZscFRKrIb6NfjKLBXoIfkemnhryaXrdgNC8,218
 ngsolve/solvers.py,sha256=Vbm5RfIJ83nMtQ1pYJL0uBBuys4-mcY6b63MNkyiqxU,311
@@ -208,12 +208,12 @@ ngsolve/timestepping.py,sha256=SuKsFS7cDXRixKWdMMIWW4flet12HLeqWTBC0DFhPwg,8535
 ngsolve/timing.py,sha256=j_Plbbh51pS9QwMq93hyvFDSrbywmdKZ6Z6OOZmZoe4,4093
 ngsolve/utils.py,sha256=jg6smLuKF6G9Rvg-EXjyYwjFJGyyiqSRknevdjp-9J4,4724
 ngsolve/webgui.py,sha256=MG9bzYvfZYodiC5-DR-9auveCzKX3d13AGOgKANne-A,24902
-ngsolve/cmake/NGSolveConfig.cmake,sha256=wH7-MiurL4tcjSKVmd-g8NpCQxwo7k7Jr6kHRnRLpsg,5292
+ngsolve/cmake/NGSolveConfig.cmake,sha256=CDDptKNu85-IxE0xWxQ_K1mzMAzxKrjMYFr-jzr1rI8,5278
 ngsolve/cmake/ngsolve-targets-release.cmake,sha256=knQjCmz6J-x_Egg5DosY5hU_CZdZlSQxQ4Za7pS07fA,907
 ngsolve/cmake/ngsolve-targets.cmake,sha256=By4x6oAQKohvxBjQ03BSG6FyD9dC3T1da43a5rXiM9U,5966
 ngsolve/config/__init__.py,sha256=IjGKWBHvRx572Z6m-5z0q82xuQaWlBABy1whwy2Uxak,23
 ngsolve/config/__main__.py,sha256=anI_XpgrpEnVnMyrFiH3Yon_tpQZkehnqRTNC02NiBI,93
-ngsolve/config/config.py,sha256=afyK3u85RJsxmDc_GpiJZU_Ck7lKRSE04U2n39UWfpc,2548
+ngsolve/config/config.py,sha256=ZEeqHqs5fnI3uVn-2MpnxcwgFzlJx5bGsVd3CzQ-nQw,2520
 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=3Q9ntmVJGIKQ-JLdGWgECR_kDjTG8FPVCSiTJkxyi40,2188
@@ -245,37 +245,37 @@ ngsolve/demos/mpi/mpi_cmagnet.py,sha256=59TuQvIIGoDStyJneH47NG131nppBoQpZQ9rfvG_
 ngsolve/demos/mpi/mpi_navierstokes.py,sha256=ieAWMDd9LoH99L14iBgF_YmX_xk7gltZ9MruTFN6xAY,3028
 ngsolve/demos/mpi/mpi_poisson.py,sha256=_-H0HFqlwlpasiVD9OQjUvJbWJcQ34-4y-va1WOSOzE,2557
 ngsolve/demos/mpi/mpi_timeDG.py,sha256=FUOh_1dzOGUU-uAsD4iyz0RVT382jMiMiuYkZjarhAE,1963
-ngsolve-6.2.2501.post83.dev1.data/data/Scripts/ngsolve.tcl,sha256=8KBOXZKYU_qx8785PIjdmVjcUbCxR5aYrod0BZ4UgGc,21637
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/beam.geo,sha256=HTv0A6IUoVuavifolhGKY6U44z0hni07Htonqe6rw_A,404
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/beam.vol,sha256=QZbM9JE4rBs5nVXaw7ZAbDf60yXbX19bm1OnOScF57M,13834
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/chip.in2d,sha256=561AjEIPnlzMZHIzep5lRjocQDm8w7o7IDna5u8s9JQ,592
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/chip.vol,sha256=gkEmBWlQh8BkSFNTf-XISfmuGzHXDKANV4Uoww_BwJI,28511
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/coil.geo,sha256=ylB2OiWEZztFXPx99qTNRJKmZFIDPEezEgQt4EPtZRA,376
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/coil.vol,sha256=t8cWP_YX38FMw0S-7vEVUOcSH0VK6DlaqVWG_zmBo-4,147284
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/coilshield.geo,sha256=Sq8w7i22WQ0viaiM3bU_63lZGFWp8JgyDqs3RdjeTe0,629
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/coilshield.vol,sha256=JiZOM9A529-hikvv1kyLwffTRpeatLtaVTNcu42o1fo,165576
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/cube.geo,sha256=Wd9AnPdoQgjQ-JK6zmoXjsNxPR9B8ED4pdHVwmRoe1M,405
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/cube.vol,sha256=RS8bay_KlF6IZM0QDDr7_dGaFLqoxszQU_w8TYutMKY,92127
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/d10_DGdoubleglazing.pde,sha256=lK3sXmK74eFEMk8F6hwFHe_tBYnOdRq3F-RbzKpDK0c,1354
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/d11_chip_nitsche.pde,sha256=VuEQtkZaIJA9yuvjG_jcx1UB9UtAaZSoeJYFHcHeYxk,1141
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/d1_square.pde,sha256=24ZZ0nlbr5LKl_mPyeZZMSkiCBiBk1fzNElo3yITX2A,1072
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/d2_chip.pde,sha256=yGIbzuTdCEjdoC9Ig_hbcy5qz69_wxgPZeYmhY1rN7k,985
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/d3_helmholtz.pde,sha256=gWqQTh1M9NiKGohlu5JxDJuvRz-T28WlNSVbwmMvM-0,590
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/d4_cube.pde,sha256=6kL43ehAh1PdwnFlLjadg325bFjpz02p2PpAOxlDFuc,1079
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/d5_beam.pde,sha256=VQt9WbqAVWObFH9Jy1YieE9WHnVZ_hnsKVTz6ysEBxI,2059
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/d6_shaft.pde,sha256=y10nLFA2-GwEwnAs5R-8wppRq-09vYJ_C9LJ3DSvlS0,1974
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/d7_coil.pde,sha256=KFkqtZiPD4LQsGnghwwdq1qmpWf0T5xQR4HwbRL9ZZA,1250
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/d8_coilshield.pde,sha256=9XG0TcgZ7mComh_pVntnuzDldEx5PD-C7gqK05PTWus,1190
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/d9_hybridDG.pde,sha256=vYschVRu5vC91YpyVZFCBFuno-ZlOcRznPD-OvJTLtQ,1551
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/doubleglazing.in2d,sha256=lDCukUh9iTpkiH-qJAVdeAlMpiYF0AibAA6g3jOeWDY,533
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/doubleglazing.vol,sha256=IUOhUrZi8WMnvrIJe--A8cYPN6mscslA8_4yxvKShb0,35724
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/piezo2d40round4.vol.gz,sha256=RuxM0vtYljoAImOH6U6avt0LmWaUbGNao37n8p_Yd1w,496477
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/shaft.geo,sha256=Y3OVQlWjgICe2Dwz6-Qs8r0tLanjwtWP0xFx5AwNRJk,2507
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/shaft.vol,sha256=IqViJ4hJKk3Uf66PkVJ-kHpBQ_ZpkigRuh5g1Po4euY,274741
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/square.in2d,sha256=5dewiEjTsRFiNkBlmaTjBAZejq55TTlusp8ALfTMc2M,209
-ngsolve-6.2.2501.post83.dev1.data/data/share/ngsolve/square.vol,sha256=tWdPeQ52deUNgWOsWKhqgNVsJKhkJmEAVp5-xjhzEwg,9311
-ngsolve-6.2.2501.post83.dev1.dist-info/LICENSE,sha256=iFoD9UsVeWEjb0aEPnmXKr_NaJC2y7NovH7KMo_5WhI,26934
-ngsolve-6.2.2501.post83.dev1.dist-info/METADATA,sha256=RsgS39HKdh3f0_jafN2XficDPMaZ1OPw_SaJceiKaAw,278
-ngsolve-6.2.2501.post83.dev1.dist-info/WHEEL,sha256=qDFAnWyk8KPVgFtn4Ts7GL4WTDWO6sseESMI7KyWq64,94
-ngsolve-6.2.2501.post83.dev1.dist-info/top_level.txt,sha256=1llU44DxNIP85uPT4RChsIuV3h2d0_5j0Ui_Df5vV8E,15
-ngsolve-6.2.2501.post83.dev1.dist-info/RECORD,,
+ngsolve-6.2.2502.data/data/Scripts/ngsolve.tcl,sha256=8KBOXZKYU_qx8785PIjdmVjcUbCxR5aYrod0BZ4UgGc,21637
+ngsolve-6.2.2502.data/data/share/ngsolve/beam.geo,sha256=HTv0A6IUoVuavifolhGKY6U44z0hni07Htonqe6rw_A,404
+ngsolve-6.2.2502.data/data/share/ngsolve/beam.vol,sha256=QZbM9JE4rBs5nVXaw7ZAbDf60yXbX19bm1OnOScF57M,13834
+ngsolve-6.2.2502.data/data/share/ngsolve/chip.in2d,sha256=561AjEIPnlzMZHIzep5lRjocQDm8w7o7IDna5u8s9JQ,592
+ngsolve-6.2.2502.data/data/share/ngsolve/chip.vol,sha256=gkEmBWlQh8BkSFNTf-XISfmuGzHXDKANV4Uoww_BwJI,28511
+ngsolve-6.2.2502.data/data/share/ngsolve/coil.geo,sha256=ylB2OiWEZztFXPx99qTNRJKmZFIDPEezEgQt4EPtZRA,376
+ngsolve-6.2.2502.data/data/share/ngsolve/coil.vol,sha256=t8cWP_YX38FMw0S-7vEVUOcSH0VK6DlaqVWG_zmBo-4,147284
+ngsolve-6.2.2502.data/data/share/ngsolve/coilshield.geo,sha256=Sq8w7i22WQ0viaiM3bU_63lZGFWp8JgyDqs3RdjeTe0,629
+ngsolve-6.2.2502.data/data/share/ngsolve/coilshield.vol,sha256=JiZOM9A529-hikvv1kyLwffTRpeatLtaVTNcu42o1fo,165576
+ngsolve-6.2.2502.data/data/share/ngsolve/cube.geo,sha256=Wd9AnPdoQgjQ-JK6zmoXjsNxPR9B8ED4pdHVwmRoe1M,405
+ngsolve-6.2.2502.data/data/share/ngsolve/cube.vol,sha256=RS8bay_KlF6IZM0QDDr7_dGaFLqoxszQU_w8TYutMKY,92127
+ngsolve-6.2.2502.data/data/share/ngsolve/d10_DGdoubleglazing.pde,sha256=lK3sXmK74eFEMk8F6hwFHe_tBYnOdRq3F-RbzKpDK0c,1354
+ngsolve-6.2.2502.data/data/share/ngsolve/d11_chip_nitsche.pde,sha256=VuEQtkZaIJA9yuvjG_jcx1UB9UtAaZSoeJYFHcHeYxk,1141
+ngsolve-6.2.2502.data/data/share/ngsolve/d1_square.pde,sha256=24ZZ0nlbr5LKl_mPyeZZMSkiCBiBk1fzNElo3yITX2A,1072
+ngsolve-6.2.2502.data/data/share/ngsolve/d2_chip.pde,sha256=yGIbzuTdCEjdoC9Ig_hbcy5qz69_wxgPZeYmhY1rN7k,985
+ngsolve-6.2.2502.data/data/share/ngsolve/d3_helmholtz.pde,sha256=gWqQTh1M9NiKGohlu5JxDJuvRz-T28WlNSVbwmMvM-0,590
+ngsolve-6.2.2502.data/data/share/ngsolve/d4_cube.pde,sha256=6kL43ehAh1PdwnFlLjadg325bFjpz02p2PpAOxlDFuc,1079
+ngsolve-6.2.2502.data/data/share/ngsolve/d5_beam.pde,sha256=VQt9WbqAVWObFH9Jy1YieE9WHnVZ_hnsKVTz6ysEBxI,2059
+ngsolve-6.2.2502.data/data/share/ngsolve/d6_shaft.pde,sha256=y10nLFA2-GwEwnAs5R-8wppRq-09vYJ_C9LJ3DSvlS0,1974
+ngsolve-6.2.2502.data/data/share/ngsolve/d7_coil.pde,sha256=KFkqtZiPD4LQsGnghwwdq1qmpWf0T5xQR4HwbRL9ZZA,1250
+ngsolve-6.2.2502.data/data/share/ngsolve/d8_coilshield.pde,sha256=9XG0TcgZ7mComh_pVntnuzDldEx5PD-C7gqK05PTWus,1190
+ngsolve-6.2.2502.data/data/share/ngsolve/d9_hybridDG.pde,sha256=vYschVRu5vC91YpyVZFCBFuno-ZlOcRznPD-OvJTLtQ,1551
+ngsolve-6.2.2502.data/data/share/ngsolve/doubleglazing.in2d,sha256=lDCukUh9iTpkiH-qJAVdeAlMpiYF0AibAA6g3jOeWDY,533
+ngsolve-6.2.2502.data/data/share/ngsolve/doubleglazing.vol,sha256=IUOhUrZi8WMnvrIJe--A8cYPN6mscslA8_4yxvKShb0,35724
+ngsolve-6.2.2502.data/data/share/ngsolve/piezo2d40round4.vol.gz,sha256=RuxM0vtYljoAImOH6U6avt0LmWaUbGNao37n8p_Yd1w,496477
+ngsolve-6.2.2502.data/data/share/ngsolve/shaft.geo,sha256=Y3OVQlWjgICe2Dwz6-Qs8r0tLanjwtWP0xFx5AwNRJk,2507
+ngsolve-6.2.2502.data/data/share/ngsolve/shaft.vol,sha256=IqViJ4hJKk3Uf66PkVJ-kHpBQ_ZpkigRuh5g1Po4euY,274741
+ngsolve-6.2.2502.data/data/share/ngsolve/square.in2d,sha256=5dewiEjTsRFiNkBlmaTjBAZejq55TTlusp8ALfTMc2M,209
+ngsolve-6.2.2502.data/data/share/ngsolve/square.vol,sha256=tWdPeQ52deUNgWOsWKhqgNVsJKhkJmEAVp5-xjhzEwg,9311
+ngsolve-6.2.2502.dist-info/LICENSE,sha256=iFoD9UsVeWEjb0aEPnmXKr_NaJC2y7NovH7KMo_5WhI,26934
+ngsolve-6.2.2502.dist-info/METADATA,sha256=kq2ILsbR6Nmgrpb6SqT_t_0_Soag7EziVGb9Eecarm4,266
+ngsolve-6.2.2502.dist-info/WHEEL,sha256=qDFAnWyk8KPVgFtn4Ts7GL4WTDWO6sseESMI7KyWq64,94
+ngsolve-6.2.2502.dist-info/top_level.txt,sha256=1llU44DxNIP85uPT4RChsIuV3h2d0_5j0Ui_Df5vV8E,15
+ngsolve-6.2.2502.dist-info/RECORD,,