ngsolve 6.2.2501.post60.dev1__cp310-cp310-manylinux_2_28_x86_64.whl → 6.2.2501.post77.dev1__cp310-cp310-manylinux_2_28_x86_64.whl

ngsolve-6.2.2501.post77.dev1.data/data/include/netgen/ectypes.hpp

@@ -0,0 +1,121 @@
+#include "expr.hpp"
+
+// arithmetics with error control (better say error monitoring)
+
+
+namespace ngbla
+{
+  
+  
+  class EC_double
+  {
+    double val;
+    double err;
+  public:
+    EC_double () : val(0), err(0) { }
+    EC_double (double a)
+      : val(a), err(fabs(a)*1e-16) { }
+    EC_double (double a, double aerr)
+      : val(a), err(aerr)
+    {
+      Check();
+    }
+    EC_double (EC_double&) = default;    
+
+    EC_double & operator= (EC_double&) = default;
+    EC_double & operator= (double a)
+    {
+      val = a; err = fabs(a)*1e-16;
+      return *this;
+    }
+    
+    double Val() const { return val; }
+    double Err() const { return err; }    
+    operator double() const { return val; }
+
+    EC_double & operator*= (EC_double b)
+    {
+      *this = EC_double(*this)*b;
+      return *this;
+    }
+
+    EC_double & operator+= (EC_double b)
+    {
+      val += b.Val();
+      err += b.Err();
+      Check();
+      return *this;
+    }
+    
+    EC_double & operator-= (EC_double b)
+    {
+      *this = EC_double(*this) + (-b);
+      return *this;
+    }
+
+    EC_double operator- () const
+    {
+      return EC_double(-val, err);
+    }
+
+    void Check() const
+    {
+      if (err > 1e-4 * abs(val))
+        {
+          // throw Exception("roundoff is big, val = "+ToString(val)+", err = ", ToString(err));
+          cout << "roundoff is big, val = "+ToString(val)+", err = " + ToString(err) << endl;
+        }
+    }
+  };
+
+
+  inline EC_double operator* (EC_double a, EC_double b)
+  {
+    return EC_double(a.Val()*b.Val(), abs(a.Val())*b.Err()+abs(b.Val())*a.Err());
+  }
+
+  inline EC_double operator* (EC_double a, double b)
+  {
+    return a * EC_double(b);
+  }
+
+  inline EC_double operator* (double b, EC_double a)
+  {
+    return a * EC_double(b);
+  }
+  
+  inline EC_double operator+ (EC_double a, EC_double b)
+  {
+    return EC_double(a.Val()+b.Val(), a.Err()+b.Err());
+  }
+
+  typedef std::complex<EC_double> EC_Complex;
+
+  inline auto operator* (double a, EC_Complex b)
+  {
+    return EC_double(a) * b;
+  }
+
+  inline auto operator* (EC_Complex a, const Complex b)
+  {
+    return a*EC_Complex(b);
+  }
+
+  inline double L2Norm2 (EC_Complex val)
+  {
+    return sqr(val.real().Val()) + sqr(val.imag().Val());
+  }
+
+
+  // typedef std::complex<__128float> Complex128;
+  // typedef std::complex<double> Complex128;  
+
+
+
+  
+  template<> struct is_scalar_type<ngfem::EC_Complex> { static constexpr bool value = true; };
+  template<> struct is_scalar_type<ngfem::EC_double> { static constexpr bool value = true; };
+  template <> inline constexpr bool IsComplex<ngfem::EC_Complex> () { return true; }
+
+  
+}

{ngsolve-6.2.2501.post60.dev1.data → ngsolve-6.2.2501.post77.dev1.data}/data/include/netgen/h1lofe.hpp

@@ -1095,6 +1095,9 @@ namespace ngfem
   H1LOFE_EXTERN template class T_ScalarFiniteElement<FE_TSegmL2<2>,ET_SEGM>;
   H1LOFE_EXTERN template class T_ScalarFiniteElement<FE_TSegmL2<3>,ET_SEGM>;
 
+  H1LOFE_EXTERN template class T_ScalarFiniteElement<ScalarFE<ET_POINT, 1>, ET_POINT, ScalarFiniteElement<0>>;
+  H1LOFE_EXTERN template class T_ScalarFiniteElement<ScalarFE<ET_TRIG, 2>, ET_TRIG, ScalarFiniteElement<2>>;
+
   H1LOFE_EXTERN template class T_ScalarFiniteElement<ScalarFE<ET_TRIG,0>,ET_TRIG>;
   H1LOFE_EXTERN template class T_ScalarFiniteElement<ScalarFE<ET_TRIG,1>,ET_TRIG>;
   // H1LOFE_EXTERN template class T_ScalarFiniteElement<FE_Trig2,ET_TRIG>;

{ngsolve-6.2.2501.post60.dev1.data → ngsolve-6.2.2501.post77.dev1.data}/data/include/netgen/mptools.hpp

@@ -162,14 +162,15 @@ namespace ngfem
       values(i) = (2*i-1)/rho * values(i-1) - values(i-2);
     */
     
-    if (rho < 1e-8)
+    if (rho < 1e-100)
       {
         values = Complex(0);
         return;
       }
     Vector j(n+1), y(n+1), jp(n+1), yp(n+1);
-    SBESJY (rho, n, j, y, jp, yp);
+    // SBESJY (rho, n, j, y, jp, yp);
 
+    /*
     values = j + Complex(0,1) * y;
     if (scale != 1.0)
       {
@@ -180,6 +181,28 @@ namespace ngfem
             prod *= scale;
           }
       }
+    */
+
+    
+    // the bessel-evaluation with scale
+    besseljs3d (n, rho, 1/scale,  j, jp);
+
+    // Bessel y directly with the recurrence formula for (y, yp):
+    double x = rho;
+    double xinv = 1/x;
+    y(0) = -xinv * cos(x);
+    yp(0) = j(0)-xinv*y(0);
+
+    double sl = 0;
+    for (int l = 1; l <= n; l++)
+      {
+        y(l) = scale * (sl*y(l-1) - yp(l-1));
+        sl += xinv;
+        yp(l) = scale * y(l-1) - (sl+xinv)*y(l);
+      }
+    
+    for (int i = 0; i <= n; i++)
+      values(i) = Complex (j(i), y(i));
   }
 
 
@@ -258,6 +281,19 @@ namespace ngfem
       scale = newscale;
     }
 
+    Vector<double> Spectrum (bool scaled) const
+    {
+      Vector<double> spec(Order()+1);
+      double fac = 1;
+      for (int n = 0; n <= Order(); n++)
+        {
+          spec(n) = fac * L2Norm2(sh.CoefsN(n));
+          if (!scaled) fac *= sqr(scale);
+        }
+      return spec;
+    }
+
+    
     template <typename TARGET>
     void Transform (MultiPole<TARGET> & target, Vec<3> dist) const
     {
@@ -338,7 +374,7 @@ namespace ngfem
       int total_sources;
       
       Node (Vec<3> acenter, double ar, int alevel, int order, double kappa)
-        : center(acenter), r(ar), level(alevel), mp(MPOrder(ar*kappa), kappa, min(1.0, 1*r*kappa))
+        : center(acenter), r(ar), level(alevel), mp(MPOrder(ar*kappa), kappa, min(1.0, r*kappa))
           // : center(acenter), r(ar), level(alevel), mp(MPOrder(ar*kappa), kappa, 1.0)
       {
         // cout << "singml, add node, level = " << level << endl;
@@ -940,6 +976,13 @@ namespace ngfem
   };
   
 
+  inline ostream & operator<< (ostream & ost, const RegularMLMultiPole & mlmp)
+  {
+    // mlmp.Print(ost);
+    ost << "RegularMLMultiPole" << endl;
+    return ost;
+  }
+
 
 
 
@@ -1033,20 +1076,20 @@ namespace ngfem
 
   class RegularMLMultiPoleCF : public CoefficientFunction
   {
-    RegularMLMultiPole mlmp;
+    shared_ptr<RegularMLMultiPole> mlmp;
   public:
     RegularMLMultiPoleCF (shared_ptr<SingularMLMultiPoleCF> asingmp, Vec<3> center, double r, int order)
-      : CoefficientFunction(1, true), mlmp(asingmp->MLMP(), center, r, order) { } 
+      : CoefficientFunction(1, true), mlmp{make_shared<RegularMLMultiPole>(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());
     }
 
-    RegularMLMultiPole & MLMP() { return mlmp; }
+    shared_ptr<RegularMLMultiPole> MLMP() { return mlmp; }
   };
 
   

{ngsolve-6.2.2501.post60.dev1.data → ngsolve-6.2.2501.post77.dev1.data}/data/include/netgen/periodic.hpp

@@ -81,6 +81,10 @@ namespace ngcomp
 
     auto & GetVertexMap() const { return vertex_map; }
     auto & GetDofMap() const { return dofmap; }    
+    FlatArray<VorB> GetDualShapeNodes (VorB vb) const override
+    {
+      return space->GetDualShapeNodes(vb);
+    }
   protected:
     // overload in quasiperiodic space
     virtual void DofMapped(size_t from, size_t to, size_t idnr) { ; }

{ngsolve-6.2.2501.post60.dev1.data → ngsolve-6.2.2501.post77.dev1.data}/data/include/netgen/vector.hpp

@@ -920,7 +920,7 @@ namespace ngbla
   INLINE auto Cross (const TA & a, const TB & b)
   {
     typedef decltype (a(0)*b(0)) T;
-    return Vec<3,T>({ a(1)*b(2)-a(2)*b(1), a(2)*b(0)-a(0)*b(2), a(0)*b(1)-a(1)*b(0) });
+    return Vec<3,T>{ a(1)*b(2)-a(2)*b(1), a(2)*b(0)-a(0)*b(2), a(0)*b(1)-a(1)*b(0) };
   }
 
   template <typename S>

{ngsolve-6.2.2501.post60.dev1.data → ngsolve-6.2.2501.post77.dev1.data}/data/lib/cmake/ngsolve/NGSolveConfig.cmake

@@ -1,4 +1,4 @@
-set(PACKAGE_VERSION "6.2.2501-60-g0fec9135c")
+set(PACKAGE_VERSION "6.2.2501-77-g928ec68c3")
 find_package(Netgen CONFIG REQUIRED HINTS
   ${CMAKE_CURRENT_LIST_DIR}
   ${CMAKE_CURRENT_LIST_DIR}/..

ngsolve-6.2.2501.post77.dev1.data/data/lib/python3.10/site-packages/ngsolve/__init__.pyi

@@ -0,0 +1,234 @@
+"""
+
+NGSolve
+=======
+
+A high order finite element library
+
+Modules:
+ngsolve.bla .... simple vectors and matrices
+ngsolve.fem .... finite elements and integrators
+ngsolve.comp ... function spaces, forms
+"""
+from __future__ import annotations
+import atexit as atexit
+from builtins import sum as builtin_sum
+import ctypes as ctypes
+import importlib as importlib
+import netgen as netgen
+from netgen import Redraw
+from netgen import TimeFunction
+from ngsolve.bla import InnerProduct
+from ngsolve.bla import Matrix
+from ngsolve.bla import Norm
+from ngsolve.bla import Vector
+from ngsolve.comp import APhiHCurlAMG
+from ngsolve.comp import BilinearForm
+from ngsolve.comp import BoundaryFromVolumeCF
+from ngsolve.comp import COUPLING_TYPE
+from ngsolve.comp import Compress
+from ngsolve.comp import CompressCompound
+from ngsolve.comp import ContactBoundary
+from ngsolve.comp import ConvertOperator
+from ngsolve.comp import Discontinuous
+from ngsolve.comp import ElementId
+from ngsolve.comp import FESpace
+from ngsolve.comp import FacetFESpace
+from ngsolve.comp import FacetSurface
+from ngsolve.comp import GridFunction
+from ngsolve.comp import H1
+from ngsolve.comp import H1LumpingFESpace
+from ngsolve.comp import HCurl
+from ngsolve.comp import HCurlAMG
+from ngsolve.comp import HCurlCurl
+from ngsolve.comp import HCurlDiv
+from ngsolve.comp import HDiv
+from ngsolve.comp import HDivDiv
+from ngsolve.comp import HDivDivSurface
+from ngsolve.comp import HDivSurface
+from ngsolve.comp import Hidden as PrivateSpace
+from ngsolve.comp import Hidden
+from ngsolve.comp import Integrate
+from ngsolve.comp import Interpolate
+from ngsolve.comp import L2
+from ngsolve.comp import LinearForm
+from ngsolve.comp import MatrixValued
+from ngsolve.comp import Mesh
+from ngsolve.comp import MultiGridPreconditioner
+from ngsolve.comp import NodalFESpace
+from ngsolve.comp import NodeId
+from ngsolve.comp import NormalFacetFESpace
+from ngsolve.comp import NormalFacetSurface
+from ngsolve.comp import NumberSpace
+from ngsolve.comp import ORDER_POLICY
+from ngsolve.comp import PatchwiseSolve
+from ngsolve.comp import Periodic
+from ngsolve.comp import PlateauFESpace
+from ngsolve.comp import Preconditioner
+from ngsolve.comp import ProductSpace
+from ngsolve.comp import Region
+from ngsolve.comp import SetHeapSize
+from ngsolve.comp import SetTestoutFile
+from ngsolve.comp import SurfaceL2
+from ngsolve.comp import SymbolicBFI
+from ngsolve.comp import SymbolicEnergy
+from ngsolve.comp import SymbolicLFI
+from ngsolve.comp import TangentialFacetFESpace
+from ngsolve.comp import TangentialSurfaceL2
+from ngsolve.comp import VTKOutput
+from ngsolve.comp import Variation
+from ngsolve.comp import VectorFacetFESpace
+from ngsolve.comp import VectorFacetSurface
+from ngsolve.comp import VectorH1
+from ngsolve.comp import VectorL2
+from ngsolve.comp import VectorNodalFESpace
+from ngsolve.comp import VectorSurfaceL2
+from ngsolve.comp import VectorValued
+from ngsolve.comp import pml
+from ngsolve.fem import BFI
+from ngsolve.fem import BSpline
+from ngsolve.fem import BlockBFI
+from ngsolve.fem import BlockLFI
+from ngsolve.fem import CacheCF
+from ngsolve.fem import CoefficientFunction
+from ngsolve.fem import CoefficientFunction as CF
+from ngsolve.fem import Cof
+from ngsolve.fem import CompoundBFI
+from ngsolve.fem import CompoundLFI
+from ngsolve.fem import Conj
+from ngsolve.fem import Cross
+from ngsolve.fem import Det
+from ngsolve.fem import ET
+from ngsolve.fem import Id
+from ngsolve.fem import IfPos
+from ngsolve.fem import IntegrationRule
+from ngsolve.fem import Inv
+from ngsolve.fem import LFI
+from ngsolve.fem import Parameter
+from ngsolve.fem import ParameterC
+from ngsolve.fem import PlaceholderCF
+from ngsolve.fem import Skew
+from ngsolve.fem import Sym
+from ngsolve.fem import Trace
+from ngsolve.fem import VoxelCoefficient
+from ngsolve.fem import Zero as ZeroCF
+from ngsolve.fem import acos
+from ngsolve.fem import asin
+from ngsolve.fem import atan
+from ngsolve.fem import atan2
+from ngsolve.fem import ceil
+from ngsolve.fem import cos
+from ngsolve.fem import cosh
+from ngsolve.fem import erf
+from ngsolve.fem import exp
+from ngsolve.fem import floor
+from ngsolve.fem import log
+from ngsolve.fem import pow
+from ngsolve.fem import sin
+from ngsolve.fem import sinh
+from ngsolve.fem import sqrt
+from ngsolve.fem import tan
+from ngsolve.la import ArnoldiSolver
+from ngsolve.la import BaseMatrix
+from ngsolve.la import BaseVector
+from ngsolve.la import BlockMatrix
+from ngsolve.la import BlockVector
+from ngsolve.la import CGSolver
+from ngsolve.la import ConstEBEMatrix
+from ngsolve.la import CreateVVector
+from ngsolve.la import DiagonalMatrix
+from ngsolve.la import Embedding
+from ngsolve.la import GMRESSolver
+from ngsolve.la import IdentityMatrix
+from ngsolve.la import MultiVector
+from ngsolve.la import PARALLEL_STATUS
+from ngsolve.la import ParallelMatrix
+from ngsolve.la import PermutationMatrix
+from ngsolve.la import Projector
+from ngsolve.la import QMRSolver
+from ngsolve.ngstd import IntRange
+from ngsolve.solve import Draw
+from ngsolve.solve import SetVisualization
+from ngsolve.timing import Timing
+from ngsolve.utils import Deviator
+from ngsolve.utils import Grad
+from ngsolve.utils import Normalize
+from ngsolve.utils import OuterProduct
+from ngsolve.utils import PyCof
+from ngsolve.utils import PyCross
+from ngsolve.utils import PyDet
+from ngsolve.utils import PyId
+from ngsolve.utils import PyInv
+from ngsolve.utils import PySkew
+from ngsolve.utils import PySym
+from ngsolve.utils import PyTrace
+from ngsolve.utils import curl
+from ngsolve.utils import div
+from ngsolve.utils import grad
+from ngsolve.utils import printonce
+import os as os
+from pyngcore.pyngcore import BitArray
+from pyngcore.pyngcore import PajeTrace
+from pyngcore.pyngcore import SetNumThreads
+from pyngcore.pyngcore import TaskManager
+from pyngcore.pyngcore import Timer
+from pyngcore.pyngcore import Timers
+import sys as sys
+from . import bla
+from . import bvp
+from . import comp
+from . import config
+from . import directsolvers
+from . import eigenvalues
+from . import fem
+from . import krylovspace
+from . import la
+from . import ngstd
+from . import nonlinearsolvers
+from . import preconditioners
+from . import solve
+from . import solvers
+from . import timestepping
+from . import timing
+from . import utils
+__all__ = ['APhiHCurlAMG', 'ArnoldiSolver', 'BBBND', 'BBND', 'BFI', 'BND', 'BSpline', 'BaseMatrix', 'BaseVector', 'BilinearForm', 'BitArray', 'BlockBFI', 'BlockLFI', 'BlockMatrix', 'BlockVector', 'BoundaryFromVolumeCF', 'CELL', 'CF', 'CGSolver', 'COUPLING_TYPE', 'CacheCF', 'CoefficientFunction', 'Cof', 'CompoundBFI', 'CompoundLFI', 'Compress', 'CompressCompound', 'Conj', 'ConstEBEMatrix', 'ContactBoundary', 'ConvertOperator', 'CreateVVector', 'Cross', 'Det', 'Deviator', 'DiagonalMatrix', 'Discontinuous', 'Draw', 'EDGE', 'ELEMENT', 'ET', 'ElementId', 'Embedding', 'FACE', 'FACET', 'FESpace', 'FacetFESpace', 'FacetSurface', 'GMRESSolver', 'Grad', 'GridFunction', 'H1', 'H1LumpingFESpace', 'HCurl', 'HCurlAMG', 'HCurlCurl', 'HCurlDiv', 'HDiv', 'HDivDiv', 'HDivDivSurface', 'HDivSurface', 'HEX', 'Hidden', 'Id', 'IdentityMatrix', 'IfPos', 'InnerProduct', 'IntRange', 'Integrate', 'IntegrationRule', 'Interpolate', 'Inv', 'L2', 'LFI', 'LinearForm', 'Matrix', 'MatrixValued', 'Mesh', 'MultiGridPreconditioner', 'MultiVector', 'NodalFESpace', 'NodeId', 'Norm', 'NormalFacetFESpace', 'NormalFacetSurface', 'Normalize', 'NumberSpace', 'ORDER_POLICY', 'OuterProduct', 'PARALLEL_STATUS', 'POINT', 'PRISM', 'PYRAMID', 'PajeTrace', 'ParallelMatrix', 'Parameter', 'ParameterC', 'PatchwiseSolve', 'Periodic', 'PermutationMatrix', 'PlaceholderCF', 'PlateauFESpace', 'Preconditioner', 'PrivateSpace', 'ProductSpace', 'Projector', 'PyCof', 'PyCross', 'PyDet', 'PyId', 'PyInv', 'PySkew', 'PySym', 'PyTrace', 'QMRSolver', 'QUAD', 'Redraw', 'Region', 'SEGM', 'SetHeapSize', 'SetNumThreads', 'SetTestoutFile', 'SetVisualization', 'Skew', 'SurfaceL2', 'Sym', 'SymbolicBFI', 'SymbolicEnergy', 'SymbolicLFI', 'TET', 'TRIG', 'TangentialFacetFESpace', 'TangentialSurfaceL2', 'TaskManager', 'TimeFunction', 'Timer', 'Timers', 'Timing', 'Trace', 'VERTEX', 'VOL', 'VTKOutput', 'Variation', 'Vector', 'VectorFacetFESpace', 'VectorFacetSurface', 'VectorH1', 'VectorL2', 'VectorNodalFESpace', 'VectorSurfaceL2', 'VectorValued', 'VoxelCoefficient', 'ZeroCF', 'acos', 'asin', 'atan', 'atan2', 'atexit', 'bla', 'builtin_sum', 'bvp', 'ceil', 'comp', 'config', 'cos', 'cosh', 'ctypes', 'curl', 'directsolvers', 'div', 'ds', 'dx', 'eigenvalues', 'erf', 'exp', 'f', 'fem', 'floor', 'grad', 'importlib', 'krylovspace', 'la', 'log', 'netgen', 'ngsglobals', 'ngslib', 'ngstd', 'nonlinearsolvers', 'os', 'pi', 'pml', 'pow', 'preconditioners', 'printonce', 'sin', 'sinh', 'solve', 'solvers', 'specialcf', 'sqrt', 'sum', 'sys', 'tan', 'timestepping', 'timing', 'unit_cube', 'unit_square', 'utils', 'x', 'y', 'z']
+def _add_flags_doc(module):
+    ...
+def _jupyter_nbextension_paths():
+    ...
+def sum(iterable, start = None):
+    """
+    NGSolve sum function that uses the first element of an iterable as
+    start argument if no start argument is provided.
+    """
+BBBND: comp.VorB  # value = <VorB.BBBND: 3>
+BBND: comp.VorB  # value = <VorB.BBND: 2>
+BND: comp.VorB  # value = <VorB.BND: 1>
+CELL: fem.NODE_TYPE  # value = <NODE_TYPE.CELL: 3>
+EDGE: fem.NODE_TYPE  # value = <NODE_TYPE.EDGE: 1>
+ELEMENT: fem.NODE_TYPE  # value = <NODE_TYPE.ELEMENT: 4>
+FACE: fem.NODE_TYPE  # value = <NODE_TYPE.FACE: 2>
+FACET: fem.NODE_TYPE  # value = <NODE_TYPE.FACET: 5>
+HEX: fem.ET  # value = <ET.HEX: 24>
+POINT: fem.ET  # value = <ET.POINT: 0>
+PRISM: fem.ET  # value = <ET.PRISM: 22>
+PYRAMID: fem.ET  # value = <ET.PYRAMID: 21>
+QUAD: fem.ET  # value = <ET.QUAD: 11>
+SEGM: fem.ET  # value = <ET.SEGM: 1>
+TET: fem.ET  # value = <ET.TET: 20>
+TRIG: fem.ET  # value = <ET.TRIG: 10>
+VERTEX: fem.NODE_TYPE  # value = <NODE_TYPE.VERTEX: 0>
+VOL: comp.VorB  # value = <VorB.VOL: 0>
+__version__: str = '6.2.2501-77-g928ec68c3'
+ds: comp.DifferentialSymbol  # value = <ngsolve.comp.DifferentialSymbol object>
+dx: comp.DifferentialSymbol  # value = <ngsolve.comp.DifferentialSymbol object>
+f: importlib.metadata.PackagePath  # value = PackagePath('mkl-2023.2.0.dist-info/top_level.txt')
+ngsglobals: comp.GlobalVariables  # value = <ngsolve.comp.GlobalVariables object>
+pi: float = 3.141592653589793
+specialcf: fem.SpecialCFCreator  # value = <ngsolve.fem.SpecialCFCreator object>
+unit_cube: netgen.libngpy._NgOCC.OCCGeometry  # value = <netgen.libngpy._NgOCC.OCCGeometry object>
+unit_square: netgen.libngpy._NgOCC.OCCGeometry  # value = <netgen.libngpy._NgOCC.OCCGeometry object>
+x: fem.CoefficientFunction  # value = <ngsolve.fem.CoefficientFunction object>
+y: fem.CoefficientFunction  # value = <ngsolve.fem.CoefficientFunction object>
+z: fem.CoefficientFunction  # value = <ngsolve.fem.CoefficientFunction object>
+ngslib =