ngsolve 6.2.2506.post60.dev0__cp312-cp312-macosx_10_15_universal2.whl → 6.2.2506.post74.dev0__cp312-cp312-macosx_10_15_universal2.whl

ngsolve/fem.pyi

@@ -2,10 +2,11 @@
 pybind fem
 """
 from __future__ import annotations
+import collections.abc
 import ngsolve.bla
 import ngsolve.ngstd
 import numpy
-import os
+import numpy.typing
 import pyngcore.pyngcore
 import typing
 __all__ = ['BFI', 'BSpline', 'BSpline2D', 'BaseMappedIntegrationPoint', 'BlockBFI', 'BlockLFI', 'CELL', 'CacheCF', 'CoefficientFunction', 'Cof', 'CompilePythonModule', 'CompoundBFI', 'CompoundLFI', 'Conj', 'CoordCF', 'CoordinateTrafo', 'Cross', 'Det', 'DifferentialOperator', 'EDGE', 'ELEMENT', 'ET', 'Einsum', 'ElementTopology', 'ElementTransformation', 'FACE', 'FACET', 'FiniteElement', 'GenerateL2ElementCode', 'H1FE', 'HCurlFE', 'HDivDivFE', 'HDivFE', 'HEX', 'Id', 'IfPos', 'IntegrationPoint', 'IntegrationRule', 'Inv', 'L2FE', 'LFI', 'LeviCivitaSymbol', 'LoggingCF', 'MeshPoint', 'MinimizationCF', 'MixedFE', 'NODE_TYPE', 'NewtonCF', 'POINT', 'PRISM', 'PYRAMID', 'Parameter', 'ParameterC', 'PlaceholderCF', 'PointEvaluationFunctional', 'QUAD', 'SEGM', 'ScalarFE', 'SetPMLParameters', 'Skew', 'SpecialCFCreator', 'Sym', 'TET', 'TRIG', 'Trace', 'VERTEX', 'VoxelCoefficient', 'Zero', 'acos', 'asin', 'atan', 'atan2', 'ceil', 'cos', 'cosh', 'erf', 'exp', 'floor', 'log', 'pow', 'sin', 'sinh', 'specialcf', 'sqrt', 'tan']
@@ -41,7 +42,10 @@ class BFI:
     @staticmethod
     def __special_treated_flags__() -> dict:
         ...
-    def ApplyElementMatrix(self, fel: FiniteElement, vec: ngsolve.bla.FlatVectorD, trafo: ElementTransformation, heapsize: int = 10000) -> typing.Any:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
+    def ApplyElementMatrix(self, fel: FiniteElement, vec: ngsolve.bla.FlatVectorD, trafo: ElementTransformation, heapsize: typing.SupportsInt = 10000) -> typing.Any:
         """
         Apply element matrix of a specific element.
         
@@ -59,7 +63,7 @@ class BFI:
         heapsize : int
           input heapsize
         """
-    def CalcElementMatrix(self, fel: FiniteElement, trafo: ElementTransformation, heapsize: int = 10000, complex: bool = False) -> typing.Any:
+    def CalcElementMatrix(self, fel: FiniteElement, trafo: ElementTransformation, heapsize: typing.SupportsInt = 10000, complex: bool = False) -> typing.Any:
         """
         Calculate element matrix of a specific element.
         
@@ -77,7 +81,7 @@ class BFI:
         complex : bool
           input complex
         """
-    def CalcLinearizedElementMatrix(self, fel: FiniteElement, vec: ngsolve.bla.FlatVectorD, trafo: ElementTransformation, heapsize: int = 10000) -> typing.Any:
+    def CalcLinearizedElementMatrix(self, fel: FiniteElement, vec: ngsolve.bla.FlatVectorD, trafo: ElementTransformation, heapsize: typing.SupportsInt = 10000) -> typing.Any:
         """
         Calculate (linearized) element matrix of a specific element.
         
@@ -129,7 +133,7 @@ class BFI:
         intrule : ngsolve.fem.Integrationrule
           input integration rule
         """
-    def __init__(self, name: str = '', coef: typing.Any, dim: int = -1, imag: bool = False, filename: str = '', **kwargs) -> None:
+    def __init__(self, name: str = '', coef: typing.Any, dim: typing.SupportsInt = -1, imag: bool = False, filename: str = '', **kwargs) -> None:
         ...
     def __initialize__(self, **kwargs) -> None:
         ...
@@ -160,6 +164,9 @@ class BSpline:
       list of float
     
     """
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def Differentiate(self) -> BSpline:
         """
         Differentiate the BSpline
@@ -169,12 +176,12 @@ class BSpline:
         Integrate the BSpline
         """
     @typing.overload
-    def __call__(self, pt: float) -> float:
+    def __call__(self, pt: typing.SupportsFloat) -> float:
         ...
     @typing.overload
     def __call__(self, cf: CoefficientFunction) -> CoefficientFunction:
         ...
-    def __init__(self, order: int, knots: list, vals: list) -> None:
+    def __init__(self, order: typing.SupportsInt, knots: list, vals: list) -> None:
         """
         B-Spline of a certain order, provide knot and value vectors
         """
@@ -186,15 +193,18 @@ class BSpline2D:
     Bilinear intepolation of data given on a regular grid
     
     """
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     @typing.overload
-    def __call__(self, x: numpy.ndarray[numpy.float64], y: numpy.ndarray[numpy.float64]) -> typing.Any:
+    def __call__(self, x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64], y: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
         ...
     @typing.overload
     def __call__(self, cx: CoefficientFunction, cy: CoefficientFunction) -> CoefficientFunction:
         ...
     def __getstate__(self) -> tuple:
         ...
-    def __init__(self, x: list, y: list, vals: list, order: int = 1, extrapolate: bool = True) -> None:
+    def __init__(self, x: list, y: list, vals: list, order: typing.SupportsInt = 1, extrapolate: bool = True) -> None:
         """
         x : list, y: list
           sorted list of grid coordinates
@@ -213,6 +223,9 @@ class BSpline2D:
     def __str__(self) -> str:
         ...
 class BaseMappedIntegrationPoint:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def __init__(self, arg0: MeshPoint) -> None:
         ...
     def __str__(self) -> str:
@@ -261,8 +274,10 @@ class CoefficientFunction:
         Creates a domain-wise CF, use with generator expressions and mesh.GetMaterials()
         and mesh.GetBoundaries()
     """
-    spacedim: int
-    def Compile(self, realcompile: bool = False, maxderiv: int = 2, wait: bool = False, keep_files: bool = False) -> CoefficientFunction:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
+    def Compile(self, realcompile: bool = False, maxderiv: typing.SupportsInt = 2, wait: bool = False, keep_files: bool = False) -> CoefficientFunction:
         """
         Compile list of individual steps, experimental improvement for deep trees
         
@@ -288,7 +303,7 @@ class CoefficientFunction:
         """
         Compute directional derivative with respect to variable
         """
-    def DiffShape(self, direction: CoefficientFunction = 1.0, Eulerian: list[CoefficientFunction] = []) -> CoefficientFunction:
+    def DiffShape(self, direction: CoefficientFunction = 1.0, Eulerian: collections.abc.Sequence[CoefficientFunction] = []) -> CoefficientFunction:
         """
         Compute shape derivative in direction
         """
@@ -327,7 +342,7 @@ class CoefficientFunction:
         """
         Evaluate on other element, as needed for DG jumps
         """
-    def Replace(self, arg0: dict[CoefficientFunction, CoefficientFunction]) -> CoefficientFunction:
+    def Replace(self, arg0: collections.abc.Mapping[CoefficientFunction, CoefficientFunction]) -> CoefficientFunction:
         ...
     def Reshape(self, arg0: tuple) -> CoefficientFunction:
         """
@@ -335,13 +350,13 @@ class CoefficientFunction:
         """
     def TensorTranspose(self, arg0: tuple) -> CoefficientFunction:
         ...
-    def _BuildFieldLines(self, mesh: ..., start_points: list[tuple[float, float, float]], num_fieldlines: int = 100, length: float = 0.5, max_points: float = 500, thickness: float = 0.0015, tolerance: float = 0.0005, direction: int = 0, randomized: bool = True, critical_value: float = -1) -> dict:
+    def _BuildFieldLines(self, mesh: ..., start_points: collections.abc.Sequence[tuple[typing.SupportsFloat, typing.SupportsFloat, typing.SupportsFloat]], num_fieldlines: typing.SupportsInt = 100, length: typing.SupportsFloat = 0.5, max_points: typing.SupportsFloat = 500, thickness: typing.SupportsFloat = 0.0015, tolerance: typing.SupportsFloat = 0.0005, direction: typing.SupportsInt = 0, randomized: bool = True, critical_value: typing.SupportsFloat = -1) -> dict:
         ...
     @typing.overload
     def __add__(self, cf: CoefficientFunction) -> CoefficientFunction:
         ...
     @typing.overload
-    def __add__(self, value: float) -> CoefficientFunction:
+    def __add__(self, value: typing.SupportsFloat) -> CoefficientFunction:
         ...
     @typing.overload
     def __add__(self, value: complex) -> CoefficientFunction:
@@ -352,7 +367,7 @@ class CoefficientFunction:
         evaluate CF at a mapped integrationpoint mip. mip can be generated by calling mesh(x,y,z)
         """
     @typing.overload
-    def __call__(self, x: float, y: float | None = None, z: float | None = None) -> ...:
+    def __call__(self, x: typing.SupportsFloat, y: typing.SupportsFloat | None = None, z: typing.SupportsFloat | None = None) -> ...:
         ...
     @typing.overload
     def __call__(self, arg0: CoordinateTrafo) -> CoefficientFunction:
@@ -361,10 +376,10 @@ class CoefficientFunction:
     def __call__(self, arg0: MeshPoint) -> typing.Any:
         ...
     @typing.overload
-    def __call__(self, arg0: numpy.ndarray[MeshPoint]) -> numpy.ndarray:
+    def __call__(self, arg0: typing.Annotated[numpy.typing.ArrayLike, MeshPoint]) -> numpy.ndarray:
         ...
     @typing.overload
-    def __getitem__(self, comp: int) -> CoefficientFunction:
+    def __getitem__(self, comp: typing.SupportsInt) -> CoefficientFunction:
         """
         returns component comp of vectorial CF
         """
@@ -392,7 +407,7 @@ class CoefficientFunction:
     def __mul__(self, cf: CoefficientFunction) -> CoefficientFunction:
         ...
     @typing.overload
-    def __mul__(self, value: float) -> CoefficientFunction:
+    def __mul__(self, value: typing.SupportsFloat) -> CoefficientFunction:
         ...
     @typing.overload
     def __mul__(self, value: complex) -> CoefficientFunction:
@@ -409,32 +424,32 @@ class CoefficientFunction:
     def __or__(self, cf: ...) -> typing.Any:
         ...
     @typing.overload
-    def __pow__(self, exponent: int) -> CoefficientFunction:
+    def __pow__(self, exponent: typing.SupportsInt) -> CoefficientFunction:
         ...
     @typing.overload
-    def __pow__(self, arg0: float) -> typing.Any:
+    def __pow__(self, arg0: typing.SupportsFloat) -> typing.Any:
         ...
     @typing.overload
     def __pow__(self, arg0: CoefficientFunction) -> typing.Any:
         ...
     @typing.overload
-    def __radd__(self, value: float) -> CoefficientFunction:
+    def __radd__(self, value: typing.SupportsFloat) -> CoefficientFunction:
         ...
     @typing.overload
     def __radd__(self, value: complex) -> CoefficientFunction:
         ...
     @typing.overload
-    def __rmul__(self, value: float) -> CoefficientFunction:
+    def __rmul__(self, value: typing.SupportsFloat) -> CoefficientFunction:
         ...
     @typing.overload
     def __rmul__(self, value: complex) -> CoefficientFunction:
         ...
-    def __rpow__(self, arg0: float) -> typing.Any:
+    def __rpow__(self, arg0: typing.SupportsFloat) -> typing.Any:
         ...
-    def __rsub__(self, value: float) -> CoefficientFunction:
+    def __rsub__(self, value: typing.SupportsFloat) -> CoefficientFunction:
         ...
     @typing.overload
-    def __rtruediv__(self, value: float) -> CoefficientFunction:
+    def __rtruediv__(self, value: typing.SupportsFloat) -> CoefficientFunction:
         ...
     @typing.overload
     def __rtruediv__(self, value: complex) -> CoefficientFunction:
@@ -447,13 +462,13 @@ class CoefficientFunction:
     def __sub__(self, cf: CoefficientFunction) -> CoefficientFunction:
         ...
     @typing.overload
-    def __sub__(self, value: float) -> CoefficientFunction:
+    def __sub__(self, value: typing.SupportsFloat) -> CoefficientFunction:
         ...
     @typing.overload
     def __truediv__(self, cf: CoefficientFunction) -> CoefficientFunction:
         ...
     @typing.overload
-    def __truediv__(self, value: float) -> CoefficientFunction:
+    def __truediv__(self, value: typing.SupportsFloat) -> CoefficientFunction:
         ...
     @typing.overload
     def __truediv__(self, value: complex) -> CoefficientFunction:
@@ -495,14 +510,26 @@ class CoefficientFunction:
         shape of CF
         """
     @property
+    def spacedim(self) -> int:
+        ...
+    @spacedim.setter
+    def spacedim(self, arg1: typing.SupportsInt) -> None:
+        ...
+    @property
     def trans(self) -> CoefficientFunction:
         """
         transpose of matrix-valued CF
         """
 class CoordinateTrafo:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def __init__(self, arg0: ..., arg1: ...) -> None:
         ...
 class DifferentialOperator:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def __call__(self, arg0: FiniteElement, arg1: MeshPoint) -> ngsolve.bla.MatrixD:
         ...
     def __timing__(self, arg0: FiniteElement, arg1: ElementTransformation, arg2: IntegrationRule) -> list[tuple[str, float]]:
@@ -538,6 +565,9 @@ class ET:
     TET: typing.ClassVar[ET]  # value = <ET.TET: 20>
     TRIG: typing.ClassVar[ET]  # value = <ET.TRIG: 10>
     __members__: typing.ClassVar[dict[str, ET]]  # value = {'POINT': <ET.POINT: 0>, 'SEGM': <ET.SEGM: 1>, 'TRIG': <ET.TRIG: 10>, 'QUAD': <ET.QUAD: 11>, 'TET': <ET.TET: 20>, 'PRISM': <ET.PRISM: 22>, 'PYRAMID': <ET.PYRAMID: 21>, 'HEX': <ET.HEX: 24>}
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def __eq__(self, other: typing.Any) -> bool:
         ...
     def __getstate__(self) -> int:
@@ -546,7 +576,7 @@ class ET:
         ...
     def __index__(self) -> int:
         ...
-    def __init__(self, value: int) -> None:
+    def __init__(self, value: typing.SupportsInt) -> None:
         ...
     def __int__(self) -> int:
         ...
@@ -554,7 +584,7 @@ class ET:
         ...
     def __repr__(self) -> str:
         ...
-    def __setstate__(self, state: int) -> None:
+    def __setstate__(self, state: typing.SupportsInt) -> None:
         ...
     def __str__(self) -> str:
         ...
@@ -575,6 +605,9 @@ class ElementTopology:
       input element type
     
     """
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def __init__(self, et: ET) -> None:
         ...
     @property
@@ -588,14 +621,17 @@ class ElementTopology:
         Vertices of the element topology
         """
 class ElementTransformation:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     @typing.overload
-    def __call__(self, x: float, y: float = 0, z: float = 0) -> BaseMappedIntegrationPoint:
+    def __call__(self, x: typing.SupportsFloat, y: typing.SupportsFloat = 0, z: typing.SupportsFloat = 0) -> BaseMappedIntegrationPoint:
         ...
     @typing.overload
     def __call__(self, ip: IntegrationPoint) -> BaseMappedIntegrationPoint:
         ...
     @typing.overload
-    def __call__(self, arg0: IntegrationRule) -> numpy.ndarray[MeshPoint]:
+    def __call__(self, arg0: IntegrationRule) -> numpy.typing.NDArray[MeshPoint]:
         ...
     def __init__(self, et: ET = ..., vertices: list) -> None:
         ...
@@ -623,6 +659,9 @@ class FiniteElement:
     """
     any finite element
     """
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def __str__(self) -> str:
         ...
     def __timing__(self) -> list[tuple[str, float]]:
@@ -656,10 +695,13 @@ class HCurlFE(FiniteElement):
     """
     an H(curl) finite element
     """
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def CalcCurlShape(self, mip: ...) -> ngsolve.bla.MatrixD:
         ...
     @typing.overload
-    def CalcShape(self, x: float, y: float = 0.0, z: float = 0.0) -> ngsolve.bla.MatrixD:
+    def CalcShape(self, x: typing.SupportsFloat, y: typing.SupportsFloat = 0.0, z: typing.SupportsFloat = 0.0) -> ngsolve.bla.MatrixD:
         ...
     @typing.overload
     def CalcShape(self, mip: ...) -> ngsolve.bla.MatrixD:
@@ -668,19 +710,28 @@ class HDivDivFE(FiniteElement):
     """
     an H(div div) finite element
     """
-    def CalcDivShape(self, x: float, y: float = 0.0, z: float = 0.0) -> ngsolve.bla.MatrixD:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
+    def CalcDivShape(self, x: typing.SupportsFloat, y: typing.SupportsFloat = 0.0, z: typing.SupportsFloat = 0.0) -> ngsolve.bla.MatrixD:
         ...
-    def CalcShape(self, x: float, y: float = 0.0, z: float = 0.0) -> ngsolve.bla.MatrixD:
+    def CalcShape(self, x: typing.SupportsFloat, y: typing.SupportsFloat = 0.0, z: typing.SupportsFloat = 0.0) -> ngsolve.bla.MatrixD:
         ...
 class HDivFE(FiniteElement):
     """
     an H(div) finite element
     """
-    def CalcDivShape(self, x: float, y: float = 0.0, z: float = 0.0) -> ngsolve.bla.VectorD:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
+    def CalcDivShape(self, x: typing.SupportsFloat, y: typing.SupportsFloat = 0.0, z: typing.SupportsFloat = 0.0) -> ngsolve.bla.VectorD:
         ...
-    def CalcShape(self, x: float, y: float = 0.0, z: float = 0.0) -> ngsolve.bla.MatrixD:
+    def CalcShape(self, x: typing.SupportsFloat, y: typing.SupportsFloat = 0.0, z: typing.SupportsFloat = 0.0) -> ngsolve.bla.MatrixD:
         ...
 class IntegrationPoint:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     @property
     def point(self) -> tuple:
         """
@@ -725,11 +776,14 @@ class IntegrationRule:
     @typing.overload
     def __init__(*args, **kwargs) -> None:
         ...
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def Integrate(self, func: typing.Any) -> typing.Any:
         """
         Integrates a given function
         """
-    def __getitem__(self, nr: int) -> IntegrationPoint:
+    def __getitem__(self, nr: typing.SupportsInt) -> IntegrationPoint:
         """
         Return integration point at given position
         """
@@ -779,11 +833,14 @@ class LFI:
       input definedonelem
     
     """
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     @typing.overload
     def CalcElementVector(self, fel: FiniteElement, trafo: ElementTransformation, vec: ngsolve.bla.FlatVectorD, lh: ngsolve.ngstd.LocalHeap) -> None:
         ...
     @typing.overload
-    def CalcElementVector(self, fel: FiniteElement, trafo: ElementTransformation, heapsize: int = 10000, complex: bool = False) -> typing.Any:
+    def CalcElementVector(self, fel: FiniteElement, trafo: ElementTransformation, heapsize: typing.SupportsInt = 10000, complex: bool = False) -> typing.Any:
         ...
     def GetDefinedOn(self) -> pyngcore.pyngcore.BitArray:
         """
@@ -810,7 +867,7 @@ class LFI:
         ir : ngsolve.fem.IntegrationRule
           input integration rule
         """
-    def __init__(self, name: str = 'lfi', dim: int = -1, coef: typing.Any, definedon: ... | list | None = None, imag: bool = False, flags: pyngcore.pyngcore.Flags = {}, definedonelements: pyngcore.pyngcore.BitArray = None) -> None:
+    def __init__(self, name: str = 'lfi', dim: typing.SupportsInt = -1, coef: typing.Any, definedon: ... | list | None = None, imag: bool = False, flags: pyngcore.pyngcore.Flags = {}, definedonelements: pyngcore.pyngcore.BitArray = None) -> None:
         ...
     def __str__(self) -> str:
         ...
@@ -823,6 +880,9 @@ class LFI:
     def simd_evaluate(self, arg1: bool) -> None:
         ...
 class MeshPoint:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     @property
     def mesh(self) -> ...:
         ...
@@ -841,6 +901,9 @@ class MixedFE(FiniteElement):
     """
     pair of finite elements for trial and test-functions
     """
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def __init__(self, arg0: FiniteElement, arg1: FiniteElement) -> None:
         ...
 class NODE_TYPE:
@@ -868,6 +931,9 @@ class NODE_TYPE:
     FACET: typing.ClassVar[NODE_TYPE]  # value = <NODE_TYPE.FACET: 5>
     VERTEX: typing.ClassVar[NODE_TYPE]  # value = <NODE_TYPE.VERTEX: 0>
     __members__: typing.ClassVar[dict[str, NODE_TYPE]]  # value = {'VERTEX': <NODE_TYPE.VERTEX: 0>, 'EDGE': <NODE_TYPE.EDGE: 1>, 'FACE': <NODE_TYPE.FACE: 2>, 'CELL': <NODE_TYPE.CELL: 3>, 'ELEMENT': <NODE_TYPE.ELEMENT: 4>, 'FACET': <NODE_TYPE.FACET: 5>}
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def __eq__(self, other: typing.Any) -> bool:
         ...
     def __getstate__(self) -> int:
@@ -876,7 +942,7 @@ class NODE_TYPE:
         ...
     def __index__(self) -> int:
         ...
-    def __init__(self, value: int) -> None:
+    def __init__(self, value: typing.SupportsInt) -> None:
         ...
     def __int__(self) -> int:
         ...
@@ -884,7 +950,7 @@ class NODE_TYPE:
         ...
     def __repr__(self) -> str:
         ...
-    def __setstate__(self, state: int) -> None:
+    def __setstate__(self, state: typing.SupportsInt) -> None:
         ...
     def __str__(self) -> str:
         ...
@@ -905,11 +971,14 @@ class Parameter(CoefficientFunction):
       Parameter value
     
     """
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def Get(self) -> float:
         """
         return parameter value
         """
-    def Set(self, value: float) -> None:
+    def Set(self, value: typing.SupportsFloat) -> None:
         """
         Modify parameter value.
         
@@ -918,27 +987,27 @@ class Parameter(CoefficientFunction):
         value : double
           input scalar  
         """
-    def __ge__(self, arg0: float) -> bool:
+    def __ge__(self, arg0: typing.SupportsFloat) -> bool:
         ...
     def __getstate__(self) -> tuple:
         ...
-    def __gt__(self, arg0: float) -> bool:
+    def __gt__(self, arg0: typing.SupportsFloat) -> bool:
         ...
-    def __iadd__(self, arg0: float) -> Parameter:
+    def __iadd__(self, arg0: typing.SupportsFloat) -> Parameter:
         ...
-    def __imul__(self, arg0: float) -> Parameter:
+    def __imul__(self, arg0: typing.SupportsFloat) -> Parameter:
         ...
-    def __init__(self, value: float) -> None:
+    def __init__(self, value: typing.SupportsFloat) -> None:
         """
         Construct a ParameterCF from a scalar
         """
-    def __isub__(self, arg0: float) -> Parameter:
+    def __isub__(self, arg0: typing.SupportsFloat) -> Parameter:
         ...
-    def __itruediv__(self, arg0: float) -> Parameter:
+    def __itruediv__(self, arg0: typing.SupportsFloat) -> Parameter:
         ...
-    def __le__(self, arg0: float) -> bool:
+    def __le__(self, arg0: typing.SupportsFloat) -> bool:
         ...
-    def __lt__(self, arg0: float) -> bool:
+    def __lt__(self, arg0: typing.SupportsFloat) -> bool:
         ...
     def __setstate__(self, arg0: tuple) -> None:
         ...
@@ -953,6 +1022,9 @@ class ParameterC(CoefficientFunction):
       Parameter value
     
     """
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def Get(self) -> complex:
         """
         return parameter value
@@ -975,17 +1047,26 @@ class ParameterC(CoefficientFunction):
     def __setstate__(self, arg0: tuple) -> None:
         ...
 class PlaceholderCF(CoefficientFunction):
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def Set(self, arg0: CoefficientFunction) -> None:
         ...
     def __init__(self, arg0: CoefficientFunction) -> None:
         ...
 class PointEvaluationFunctional:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     def Assemble(self) -> ngsolve.bla.SparseVector:
         ...
 class ScalarFE(FiniteElement):
     """
     a scalar-valued finite element
     """
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
     @typing.overload
     def CalcDShape(self, mip: ...) -> ngsolve.bla.MatrixD:
         """
@@ -997,7 +1078,7 @@ class ScalarFE(FiniteElement):
           input mapped integration point
         """
     @typing.overload
-    def CalcDShape(self, x: float, y: float = 0.0, z: float = 0.0) -> ngsolve.bla.MatrixD:
+    def CalcDShape(self, x: typing.SupportsFloat, y: typing.SupportsFloat = 0.0, z: typing.SupportsFloat = 0.0) -> ngsolve.bla.MatrixD:
         """
         Computes derivative of the shape in an integration point.
         
@@ -1013,7 +1094,7 @@ class ScalarFE(FiniteElement):
           input z value
         """
     @typing.overload
-    def CalcShape(self, x: float, y: float = 0.0, z: float = 0.0) -> ngsolve.bla.VectorD:
+    def CalcShape(self, x: typing.SupportsFloat, y: typing.SupportsFloat = 0.0, z: typing.SupportsFloat = 0.0) -> ngsolve.bla.VectorD:
         """
         Parameters:
         
@@ -1035,40 +1116,43 @@ class ScalarFE(FiniteElement):
           input mapped integration point
         """
 class SpecialCFCreator:
-    def EdgeCurvature(self, dim: int) -> ...:
+    @staticmethod
+    def _pybind11_conduit_v1_(*args, **kwargs):
+        ...
+    def EdgeCurvature(self, dim: typing.SupportsInt) -> ...:
         """
         EdgeCurvature 
         space-dimension must be provided
         """
-    def EdgeFaceTangentialVectors(self, dim: int) -> ...:
+    def EdgeFaceTangentialVectors(self, dim: typing.SupportsInt) -> ...:
         """
         EdgeFaceTangentialVectors 
         space-dimension must be provided
         """
     @typing.overload
-    def JacobianMatrix(self, dim: int) -> ...:
+    def JacobianMatrix(self, dim: typing.SupportsInt) -> ...:
         """
         Jacobian matrix of transformation to physical element
         space-dimension must be provided
         """
     @typing.overload
-    def JacobianMatrix(self, dimr: int, dims: int) -> ...:
+    def JacobianMatrix(self, dimr: typing.SupportsInt, dims: typing.SupportsInt) -> ...:
         """
         Jacobian matrix of transformation to physical element
         space-dimensions dimr >= dims must be provided
         """
-    def VertexTangentialVectors(self, dim: int) -> ...:
+    def VertexTangentialVectors(self, dim: typing.SupportsInt) -> ...:
         """
         VertexTangentialVectors 
         space-dimension must be provided
         """
-    def Weingarten(self, dim: int) -> ...:
+    def Weingarten(self, dim: typing.SupportsInt) -> ...:
         """
         Weingarten tensor 
         space-dimension must be provided
         """
     @typing.overload
-    def normal(self, arg0: int) -> ...:
+    def normal(self, arg0: typing.SupportsInt) -> ...:
         """
         depending on contents: normal-vector to geometry or element
         space-dimension must be provided.
@@ -1082,12 +1166,12 @@ class SpecialCFCreator:
         """
         number of elements on facet, available for element-bnd integrals, and surface integrals
         """
-    def tangential(self, dim: int, consistent: bool = False) -> ...:
+    def tangential(self, dim: typing.SupportsInt, consistent: bool = False) -> ...:
         """
         depending on contents: tangential-vector to element
         space-dimension must be provided
         """
-    def xref(self, dim: int) -> ...:
+    def xref(self, dim: typing.SupportsInt) -> ...:
         """
         element reference-coordinates
         """
@@ -1096,7 +1180,7 @@ class SpecialCFCreator:
         """
         local mesh-size (approximate element diameter) as CF
         """
-def BlockBFI(bfi: BFI = 0, dim: int = 2, comp: int = 0) -> ...:
+def BlockBFI(bfi: BFI = 0, dim: typing.SupportsInt = 2, comp: typing.SupportsInt = 0) -> ...:
     """
     Block Bilinear Form Integrator
     
@@ -1111,7 +1195,7 @@ def BlockBFI(bfi: BFI = 0, dim: int = 2, comp: int = 0) -> ...:
     comp : int
       input comp
     """
-def BlockLFI(lfi: LFI = 0, dim: int = 2, comp: int = 0) -> LFI:
+def BlockLFI(lfi: LFI = 0, dim: typing.SupportsInt = 2, comp: typing.SupportsInt = 0) -> LFI:
     """
     Block Linear Form Integrator
     
@@ -1165,7 +1249,7 @@ def CompilePythonModule(code: str, init_function_name: str = 'init', add_header:
     }
     """
 @typing.overload
-def CompilePythonModule(file: os.PathLike, init_function_name: str = 'init') -> typing.Any:
+def CompilePythonModule(file: os.PathLike | str | bytes, init_function_name: str = 'init') -> typing.Any:
     """
     Utility function to compile a c++ file with python bindings at run-time.
     
@@ -1176,7 +1260,7 @@ def CompilePythonModule(file: os.PathLike, init_function_name: str = 'init') ->
     init_function_name (default = "init"): Function, which is called after the compiled code is loaded. The prototype must match:
     extern "C" void init_function_name(pybind11::object & res);
     """
-def CompoundBFI(bfi: BFI = 0, comp: int = 0) -> ...:
+def CompoundBFI(bfi: BFI = 0, comp: typing.SupportsInt = 0) -> ...:
     """
     Compound Bilinear Form Integrator
     
@@ -1188,7 +1272,7 @@ def CompoundBFI(bfi: BFI = 0, comp: int = 0) -> ...:
     comp : int
       input component
     """
-def CompoundLFI(lfi: LFI = 0, comp: int = 0) -> LFI:
+def CompoundLFI(lfi: LFI = 0, comp: typing.SupportsInt = 0) -> LFI:
     """
     Compound Linear Form Integrator
     
@@ -1204,7 +1288,7 @@ def Conj(arg0: CoefficientFunction) -> CoefficientFunction:
     """
     complex-conjugate
     """
-def CoordCF(direction: int) -> ...:
+def CoordCF(direction: typing.SupportsInt) -> ...:
     """
     CoefficientFunction for x, y, z.
     
@@ -1236,9 +1320,9 @@ def Einsum(einsum_signature: str, *args, **kwargs) -> CoefficientFunction:
       "use_legacy_ops" (false) -- fall back to existing CFs implementing certain blas operations where possible
       "sparse_evaluation" (true) -- exploit sparsity of tensors
     """
-def GenerateL2ElementCode(arg0: int) -> str:
+def GenerateL2ElementCode(arg0: typing.SupportsInt) -> str:
     ...
-def H1FE(et: ET, order: int) -> None:
+def H1FE(et: ET, order: typing.SupportsInt) -> None:
     """
     Creates an H1 finite element of given geometric shape and polynomial order.
     
@@ -1251,7 +1335,7 @@ def H1FE(et: ET, order: int) -> None:
       input polynomial order
     """
 @typing.overload
-def Id(dim: int) -> CoefficientFunction:
+def Id(dim: typing.SupportsInt) -> CoefficientFunction:
     """
     Identity matrix of given dimension
     """
@@ -1279,11 +1363,11 @@ def IfPos(c1: ..., then_obj: typing.Any, else_obj: typing.Any) -> ...:
       ngsolve.CoefficientFunction. See help(CoefficientFunction ) for information.
     """
 @typing.overload
-def IfPos(arg0: numpy.ndarray[numpy.float64], arg1: numpy.ndarray[numpy.float64], arg2: numpy.ndarray[numpy.float64]) -> typing.Any:
+def IfPos(arg0: typing.Annotated[numpy.typing.ArrayLike, numpy.float64], arg1: typing.Annotated[numpy.typing.ArrayLike, numpy.float64], arg2: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     ...
 def Inv(arg0: CoefficientFunction) -> CoefficientFunction:
     ...
-def L2FE(et: ET, order: int) -> None:
+def L2FE(et: ET, order: typing.SupportsInt) -> None:
     """
     Creates an L2 finite element of given geometric shape and polynomial order.
     
@@ -1295,11 +1379,11 @@ def L2FE(et: ET, order: int) -> None:
     order : int
       input polynomial order
     """
-def LeviCivitaSymbol(arg0: int) -> CoefficientFunction:
+def LeviCivitaSymbol(arg0: typing.SupportsInt) -> CoefficientFunction:
     ...
 def LoggingCF(cf: CoefficientFunction, logfile: str = 'stdout') -> CoefficientFunction:
     ...
-def MinimizationCF(expression: CoefficientFunction, startingpoint: typing.Any, tol: float | None = 1e-06, rtol: float | None = 0.0, maxiter: int | None = 20, allow_fail: bool | None = False) -> CoefficientFunction:
+def MinimizationCF(expression: CoefficientFunction, startingpoint: typing.Any, tol: typing.SupportsFloat | None = 1e-06, rtol: typing.SupportsFloat | None = 0.0, maxiter: typing.SupportsInt | None = 20, allow_fail: bool | None = False) -> CoefficientFunction:
     """
     Creates a CoefficientFunction that returns the solution to a minimization problem.
     Convergence failure is indicated by returning NaNs. Set ngsgloals.message_level
@@ -1328,7 +1412,7 @@ def MinimizationCF(expression: CoefficientFunction, startingpoint: typing.Any, t
       Returns the result of the final Newton step, even if the tolerance is not reached.
       Otherwise NaNs are returned.
     """
-def NewtonCF(expression: CoefficientFunction, startingpoint: typing.Any, tol: float | None = 1e-06, rtol: float | None = 0.0, maxiter: int | None = 10, allow_fail: bool | None = False) -> CoefficientFunction:
+def NewtonCF(expression: CoefficientFunction, startingpoint: typing.Any, tol: typing.SupportsFloat | None = 1e-06, rtol: typing.SupportsFloat | None = 0.0, maxiter: typing.SupportsInt | None = 10, allow_fail: bool | None = False) -> CoefficientFunction:
     """
     Creates a CoefficientFunction that returns the solution to a nonlinear problem.
     By default, convergence failure is indicated by returning NaNs. Set ngsgloals.message_level
@@ -1357,7 +1441,7 @@ def NewtonCF(expression: CoefficientFunction, startingpoint: typing.Any, tol: fl
         Returns the result of the final Newton step, even if the tolerance is not reached.
         Otherwise NaNs are returned.
     """
-def SetPMLParameters(rad: float = 1, alpha: float = 1) -> None:
+def SetPMLParameters(rad: typing.SupportsFloat = 1, alpha: typing.SupportsFloat = 1) -> None:
     """
     Parameters:
     
@@ -1385,12 +1469,12 @@ def VoxelCoefficient(start: tuple, end: tuple, values: numpy.ndarray, linear: bo
 def Zero(arg0: pyngcore.pyngcore.Array_I_S) -> CoefficientFunction:
     ...
 @typing.overload
-def acos(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def acos(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Inverse cosine in radians
     """
 @typing.overload
-def acos(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def acos(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Inverse cosine in radians
     """
@@ -1400,12 +1484,12 @@ def acos(x: ...) -> ...:
     Inverse cosine in radians
     """
 @typing.overload
-def asin(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def asin(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Inverse sine in radians
     """
 @typing.overload
-def asin(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def asin(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Inverse sine in radians
     """
@@ -1415,12 +1499,12 @@ def asin(x: ...) -> ...:
     Inverse sine in radians
     """
 @typing.overload
-def atan(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def atan(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Inverse tangent in radians
     """
 @typing.overload
-def atan(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def atan(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Inverse tangent in radians
     """
@@ -1430,12 +1514,12 @@ def atan(x: ...) -> ...:
     Inverse tangent in radians
     """
 @typing.overload
-def atan2(y: numpy.ndarray[numpy.float64], x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def atan2(y: typing.Annotated[numpy.typing.ArrayLike, numpy.float64], x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Four quadrant inverse tangent in radians
     """
 @typing.overload
-def atan2(y: numpy.ndarray[numpy.complex128], x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def atan2(y: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128], x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Four quadrant inverse tangent in radians
     """
@@ -1445,12 +1529,12 @@ def atan2(y: ..., x: ...) -> ...:
     Four quadrant inverse tangent in radians
     """
 @typing.overload
-def ceil(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def ceil(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Round to next greater integer
     """
 @typing.overload
-def ceil(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def ceil(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Round to next greater integer
     """
@@ -1460,12 +1544,12 @@ def ceil(x: ...) -> ...:
     Round to next greater integer
     """
 @typing.overload
-def cos(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def cos(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Cosine of argument in radians
     """
 @typing.overload
-def cos(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def cos(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Cosine of argument in radians
     """
@@ -1475,12 +1559,12 @@ def cos(x: ...) -> ...:
     Cosine of argument in radians
     """
 @typing.overload
-def cosh(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def cosh(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Hyperbolic cosine of argument in radians
     """
 @typing.overload
-def cosh(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def cosh(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Hyperbolic cosine of argument in radians
     """
@@ -1490,12 +1574,12 @@ def cosh(x: ...) -> ...:
     Hyperbolic cosine of argument in radians
     """
 @typing.overload
-def erf(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def erf(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Error function
     """
 @typing.overload
-def erf(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def erf(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Error function
     """
@@ -1505,12 +1589,12 @@ def erf(x: ...) -> ...:
     Error function
     """
 @typing.overload
-def exp(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def exp(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Exponential function
     """
 @typing.overload
-def exp(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def exp(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Exponential function
     """
@@ -1520,12 +1604,12 @@ def exp(x: ...) -> ...:
     Exponential function
     """
 @typing.overload
-def floor(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def floor(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Round to next lower integer
     """
 @typing.overload
-def floor(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def floor(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Round to next lower integer
     """
@@ -1535,12 +1619,12 @@ def floor(x: ...) -> ...:
     Round to next lower integer
     """
 @typing.overload
-def log(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def log(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Logarithm function
     """
 @typing.overload
-def log(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def log(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Logarithm function
     """
@@ -1550,12 +1634,12 @@ def log(x: ...) -> ...:
     Logarithm function
     """
 @typing.overload
-def pow(x: numpy.ndarray[numpy.float64], y: numpy.ndarray[numpy.float64]) -> typing.Any:
+def pow(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64], y: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Power function
     """
 @typing.overload
-def pow(x: numpy.ndarray[numpy.complex128], y: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def pow(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128], y: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Power function
     """
@@ -1565,12 +1649,12 @@ def pow(x: ..., y: ...) -> ...:
     Power function
     """
 @typing.overload
-def sin(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def sin(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Sine of argument in radians
     """
 @typing.overload
-def sin(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def sin(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Sine of argument in radians
     """
@@ -1580,12 +1664,12 @@ def sin(x: ...) -> ...:
     Sine of argument in radians
     """
 @typing.overload
-def sinh(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def sinh(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Hyperbolic sine of argument in radians
     """
 @typing.overload
-def sinh(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def sinh(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Hyperbolic sine of argument in radians
     """
@@ -1595,12 +1679,12 @@ def sinh(x: ...) -> ...:
     Hyperbolic sine of argument in radians
     """
 @typing.overload
-def sqrt(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def sqrt(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Square root function
     """
 @typing.overload
-def sqrt(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def sqrt(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Square root function
     """
@@ -1610,12 +1694,12 @@ def sqrt(x: ...) -> ...:
     Square root function
     """
 @typing.overload
-def tan(x: numpy.ndarray[numpy.float64]) -> typing.Any:
+def tan(x: typing.Annotated[numpy.typing.ArrayLike, numpy.float64]) -> typing.Any:
     """
     Tangent of argument in radians
     """
 @typing.overload
-def tan(x: numpy.ndarray[numpy.complex128]) -> typing.Any:
+def tan(x: typing.Annotated[numpy.typing.ArrayLike, numpy.complex128]) -> typing.Any:
     """
     Tangent of argument in radians
     """