ngsolve 6.2.2501.post48.dev1__cp310-cp310-manylinux_2_28_x86_64.whl → 6.2.2501.post55.dev1__cp310-cp310-manylinux_2_28_x86_64.whl

{ngsolve-6.2.2501.post48.dev1.data → ngsolve-6.2.2501.post55.dev1.data}/data/include/netgen/mptools.hpp

@@ -249,13 +249,47 @@ namespace ngfem
       return sgn/2 * sqrt((n-m)*(n+m+1));
     }
     
+    // 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)
+    {
+      double fx = d(0);
+      double fy = d(1);
+      double fz = d(2);
+      double invscale = 1./scale;
+      
+      for (int n = 0; n < order; n++)
+        for (int m = -n; m <= n; m++)
+          {
+            double amn = CalcAmn(m,n);
+            double bmn1 = CalcBmn(-m-1, n+1);
+            double bmn2 = CalcBmn(m-1,n+1);
+            
+            res.Coef(n+1, m-1) += invscale * Complex(0.5*fx,0.5*fy)*bmn2 * Coef(n,m);
+            res.Coef(n+1, m  ) -= invscale * fz*amn * Coef(n,m);                
+            res.Coef(n+1, m+1) += invscale * Complex(0.5*fx,-0.5*fy)*bmn1 * Coef(n,m);
+            
+            res.Coef(n, m) += scale * Complex(-0.5*fx,0.5*fy)*bmn2 * Coef(n+1,m-1);
+            res.Coef(n, m) += scale * fz*amn * Coef(n+1,m);
+            res.Coef(n, m) += scale * Complex(-0.5*fx,-0.5*fy)*bmn1 * Coef(n+1,m+1);
+          }
+    }
     
     void RotateY (double alpha)
     {
       LocalHeap lh(8*6*sqr(order) + 8*15*order + 500);
       
       static Timer t("mptool sh RotateY"); RegionTimer rg(t);
-
+      /*
+      static std::map<int, unique_ptr<Timer<>>> timers;
+      static std::map<int, unique_ptr<Timer<>>> timerstrafo;      
+      if (timers.find(order)==timers.end())
+        {
+          timers[order] = make_unique<Timer<>> (string("mptools sh RotateY ")+ToString(order));
+          timerstrafo[order] = make_unique<Timer<>> (string("mptools sh RotateY trafo ")+ToString(order));
+        }
+      RegionTimer rg(*timers[order]);
+      */
       double s = sin(alpha);
       double c = cos(alpha);
 
@@ -325,6 +359,8 @@ namespace ngfem
               trafo(m, n+mp-1) = invDmn * (  Dmn(order+mp,n)*trafo(m  ,n+mp+1)
                                              +Dmn(order+m-1,n)*trafo(m-1,n+mp));
             }
+
+          // RegionTimer rgtrafo(*timerstrafo[order]);                    
           // Step 6
           // symmetries in m and mp
           for (int m = 0; m <= n; m++)
@@ -351,8 +387,8 @@ namespace ngfem
 
           FlatVector<Complex> cn = CoefsN(n);
           FlatVector<Complex> old = cn | lh;
-          
-          cn = Trans(trafo) * old.Range(n, 2*n+1);
+
+          cn.Slice(0,1) = Trans(trafo) * old.Range(n, 2*n+1);
           cn.Slice(0,1).Reversed() += Trans(trafo.Rows(1,n+1)) * old.Range(0,n).Reversed();
 
           for (int m = 1; m <= n; m+=2)
@@ -787,10 +823,11 @@ c
       return;
       */
       
-      // book, formula (2.2.20)
       if constexpr (!std::is_same<RADIAL,MPSingular>())
         throw Exception("AddCharge assumes singular MP");
 
+      /*
+      // book, formula (2.2.20)
       // dipole in origin:
       MultiPole<MPSingular> tmp(1, kappa);
       tmp.SH().Coef(1,1)  += Complex(0,1)*sqr(kappa)*sh.CalcBmn(-1,1)/(2*sqrt(4*M_PI)) * d(0)*c;
@@ -801,6 +838,11 @@ c
       
       tmp.SH().Coef(1,0) += -Complex(0,1)*kappa*kappa*sh.CalcAmn(0,0)/sqrt(4*M_PI) *d(2)*c;
       tmp.TransformAdd (*this, -x);
+      */
+
+      MultiPole<MPSingular> tmp(Order(), kappa, Scale());
+      tmp.AddCharge(x, c);
+      tmp.SH().DirectionalDiffAdd (kappa*d,  this->SH(), Scale());
     }
 
 

{ngsolve-6.2.2501.post48.dev1.data → ngsolve-6.2.2501.post55.dev1.data}/data/include/netgen/ngblas.hpp

@@ -654,9 +654,19 @@ namespace ngbla
 
   template <bool ADD, ORDERING ord>
   extern NGS_DLL_HEADER  
-  void NgGEMV (Complex s, BareSliceMatrix<Complex,ord> a, SliceVector<Complex> x, SliceVector<Complex> y) NETGEN_NOEXCEPT;
+  void NgGEMV (Complex s, BareSliceMatrix<double,ord> a, FlatVector<const Complex> x, FlatVector<Complex> y) NETGEN_NOEXCEPT;
 
+  
+  template <bool ADD, ORDERING ord>
+  extern NGS_DLL_HEADER  
+  void NgGEMV (Complex s, BareSliceMatrix<Complex,ord> a, SliceVector<Complex> x, SliceVector<Complex> y) NETGEN_NOEXCEPT;
 
+  template <bool ADD, bool POS, ORDERING ord>
+  INLINE void NgGEMV (BareSliceMatrix<double,ord> a, FlatVector<const Complex> x, FlatVector<Complex> y)
+  {
+    NgGEMV<ADD,ord> (POS ? 1.0 : -1.0, a, x, y);
+  }
+  
   
   template <> INLINE void NgGEMV<false,true> (BareSliceMatrix<double,RowMajor> a, FlatVector<const double> x, FlatVector<double> y)
   {

{ngsolve-6.2.2501.post48.dev1.data → ngsolve-6.2.2501.post55.dev1.data}/data/lib/cmake/ngsolve/NGSolveConfig.cmake

@@ -1,4 +1,4 @@
-set(PACKAGE_VERSION "6.2.2501-48-g1ba25b2ba")
+set(PACKAGE_VERSION "6.2.2501-55-ge8d7eaad9")
 find_package(Netgen CONFIG REQUIRED HINTS
   ${CMAKE_CURRENT_LIST_DIR}
   ${CMAKE_CURRENT_LIST_DIR}/..

ngsolve-6.2.2501.post55.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
+from ngsolve.comp import Hidden as PrivateSpace
+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 as CF
+from ngsolve.fem import CoefficientFunction
+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-55-ge8d7eaad9'
+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 =