ngsolve 6.2.2506.post74.dev0__cp314-cp314-macosx_10_15_universal2.whl

ngsolve/comp/pml.pyi

@@ -0,0 +1,89 @@
+"""
+module for perfectly matched layers
+"""
+from __future__ import annotations
+import ngsolve.bla
+import ngsolve.fem
+import typing
+__all__: list[str] = ['BrickRadial', 'Cartesian', 'Compound', 'Custom', 'HalfSpace', 'PML', 'Radial']
+class PML:
+    """
+    Base PML object
+    
+    can only be created by generator functions. Use PML(x, [y, z]) to evaluate the scaling.
+    """
+    @staticmethod
+    def __call__(*args) -> ngsolve.bla.VectorC:
+        """
+        map a point
+        """
+    @staticmethod
+    def call_jacobian(*args) -> ngsolve.bla.MatrixC:
+        """
+        evaluate PML jacobian at point x, [y, z]
+        """
+    def __add__(self, pml: PML) -> PML:
+        ...
+    def __str__(self) -> str:
+        ...
+    @property
+    def Det_CF(self) -> ngsolve.fem.CoefficientFunction:
+        """
+        the determinant of the jacobian as coefficient function
+        """
+    @property
+    def JacInv_CF(self) -> ngsolve.fem.CoefficientFunction:
+        """
+        the inverse of the jacobian as coefficient function
+        """
+    @property
+    def Jac_CF(self) -> ngsolve.fem.CoefficientFunction:
+        """
+        the jacobian of the PML as coefficient function
+        """
+    @property
+    def PML_CF(self) -> ngsolve.fem.CoefficientFunction:
+        """
+        the scaling as coefficient function
+        """
+    @property
+    def dim(self) -> int:
+        """
+        dimension
+        """
+def BrickRadial(mins: typing.Any, maxs: typing.Any, origin: typing.Any = (0.0, 0.0, 0.0), alpha: complex = 1j) -> PML:
+    """
+    radial pml on a brick
+    
+    mins, maxs and origin are given as tuples/lists
+    """
+def Cartesian(mins: typing.Any, maxs: typing.Any, alpha: complex = 1j) -> PML:
+    """
+    cartesian pml transformation
+    
+    mins and maxs are tuples/lists determining the dimension
+    """
+def Compound(pml1: PML, pml2: PML, dims1: typing.Any = None, dims2: typing.Any = None) -> PML:
+    """
+    tensor product of two pml transformations
+    
+            dimensions are optional, given as tuples/lists and start with 1
+    """
+def Custom(trafo: ngsolve.fem.CoefficientFunction, jac: ngsolve.fem.CoefficientFunction) -> PML:
+    """
+    custom pml transformation
+    
+    trafo and jac are coefficient functions of the scaling and the jacobian
+    """
+def HalfSpace(point: typing.Any, normal: typing.Any, alpha: complex = 1j) -> PML:
+    """
+    half space pml
+    
+    scales orthogonal to specified plane in direction of normal point and normal are given as tuples/lists determining the dimension
+    """
+def Radial(origin: typing.Any, rad: typing.SupportsFloat = 1, alpha: complex = 1j) -> PML:
+    """
+    radial pml transformation
+    
+    origin is a list/tuple with as many entries as dimenson
+    """

ngsolve/config/__init__.py

@@ -0,0 +1 @@
+from .config import *

ngsolve/config/__init__.pyi

@@ -0,0 +1,43 @@
+from __future__ import annotations
+from ngsolve.config.config import get_cmake_dir
+from . import config
+__all__: list[str] = ['BUILD_STUB_FILES', 'BUILD_UMFPACK', 'CMAKE_CUDA_COMPILER', 'CMAKE_CXX_COMPILER', 'CMAKE_CXX_COMPILER_LAUNCHER', 'CMAKE_C_COMPILER', 'CMAKE_INSTALL_PREFIX', 'CMAKE_LINKER', 'ENABLE_UNIT_TESTS', 'INSTALL_DEPENDENCIES', 'MKL_LINK', 'NETGEN_DIR', 'NGSOLVE_COMPILE_DEFINITIONS', 'NGSOLVE_COMPILE_DEFINITIONS_PRIVATE', 'NGSOLVE_COMPILE_INCLUDE_DIRS', 'NGSOLVE_COMPILE_OPTIONS', 'NGSOLVE_INSTALL_DIR_BIN', 'NGSOLVE_INSTALL_DIR_CMAKE', 'NGSOLVE_INSTALL_DIR_INCLUDE', 'NGSOLVE_INSTALL_DIR_LIB', 'NGSOLVE_INSTALL_DIR_PYTHON', 'NGSOLVE_INSTALL_DIR_RES', 'NGSOLVE_VERSION', 'NGSOLVE_VERSION_GIT', 'NGSOLVE_VERSION_HASH', 'NGSOLVE_VERSION_MAJOR', 'NGSOLVE_VERSION_MINOR', 'NGSOLVE_VERSION_PATCH', 'NGSOLVE_VERSION_PYTHON', 'NGSOLVE_VERSION_TWEAK', 'USE_CCACHE', 'USE_HYPRE', 'USE_LAPACK', 'USE_MKL', 'USE_MUMPS', 'USE_PARDISO', 'USE_UMFPACK', 'config', 'get_cmake_dir', 'is_python_package', 'version']
+BUILD_STUB_FILES: bool = True
+BUILD_UMFPACK: bool = False
+CMAKE_CUDA_COMPILER: str = ''
+CMAKE_CXX_COMPILER: str = '/Library/Developer/CommandLineTools/usr/bin/c++'
+CMAKE_CXX_COMPILER_LAUNCHER: str = '/usr/local/bin/ccache'
+CMAKE_C_COMPILER: str = '/Library/Developer/CommandLineTools/usr/bin/cc'
+CMAKE_INSTALL_PREFIX: str = '/Users/gitlab-runner/builds/builds/rL7WHzyj/0/ngsolve/ngsolve/_skbuild/macosx-10.15-universal2-3.14/cmake-install'
+CMAKE_LINKER: str = '/Library/Developer/CommandLineTools/usr/bin/ld'
+ENABLE_UNIT_TESTS: bool = False
+INSTALL_DEPENDENCIES: bool = False
+MKL_LINK: str = ''
+NETGEN_DIR: str = '/Users/gitlab-runner/Library/Python/3.14/lib/python/site-packages'
+NGSOLVE_COMPILE_DEFINITIONS: str = 'HAVE_NETGEN_SOURCES;HAVE_DLFCN_H;HAVE_CXA_DEMANGLE;USE_TIMEOFDAY;MSG_NOSIGNAL=0;TCL;LAPACK;NGS_PYTHON;USE_UMFPACK'
+NGSOLVE_COMPILE_DEFINITIONS_PRIVATE: str = ''
+NGSOLVE_COMPILE_INCLUDE_DIRS: str = ''
+NGSOLVE_COMPILE_OPTIONS: str = '$<$<COMPILE_LANGUAGE:CXX>:-std=c++17>;$<$<COMPILE_LANGUAGE:CXX>:-Wno-undefined-var-template;-Wno-vla-extension>;-DMAX_SYS_DIM=3'
+NGSOLVE_INSTALL_DIR_BIN: str = 'bin'
+NGSOLVE_INSTALL_DIR_CMAKE: str = 'ngsolve/cmake'
+NGSOLVE_INSTALL_DIR_INCLUDE: str = 'netgen/include'
+NGSOLVE_INSTALL_DIR_LIB: str = 'netgen'
+NGSOLVE_INSTALL_DIR_PYTHON: str = '.'
+NGSOLVE_INSTALL_DIR_RES: str = 'share'
+NGSOLVE_VERSION: str = '6.2.2506-74-gd430159c1'
+NGSOLVE_VERSION_GIT: str = 'v6.2.2506-74-gd430159c1'
+NGSOLVE_VERSION_HASH: str = 'gd430159c1'
+NGSOLVE_VERSION_MAJOR: str = '6'
+NGSOLVE_VERSION_MINOR: str = '2'
+NGSOLVE_VERSION_PATCH: str = '2506'
+NGSOLVE_VERSION_PYTHON: str = '6.2.2506.post74.dev0'
+NGSOLVE_VERSION_TWEAK: str = '74'
+USE_CCACHE: bool = True
+USE_HYPRE: bool = False
+USE_LAPACK: bool = True
+USE_MKL: bool = False
+USE_MUMPS: bool = False
+USE_PARDISO: bool = False
+USE_UMFPACK: bool = True
+is_python_package: bool = True
+version: str = 'v6.2.2506-74-gd430159c1'

ngsolve/config/__main__.py

@@ -0,0 +1,4 @@
+from .config import get_cmake_dir
+
+if __name__ == '__main__':
+    print(get_cmake_dir())

ngsolve/config/config.py

@@ -0,0 +1,60 @@
+def _cmake_to_bool(s):
+    return s.upper() not in ['', '0','FALSE','OFF','N','NO','IGNORE','NOTFOUND']
+
+is_python_package    = _cmake_to_bool("TRUE")
+
+BUILD_STUB_FILES     = _cmake_to_bool("ON")
+BUILD_UMFPACK        = _cmake_to_bool("")
+ENABLE_UNIT_TESTS    = _cmake_to_bool("OFF")
+INSTALL_DEPENDENCIES = _cmake_to_bool("OFF")
+USE_CCACHE           = _cmake_to_bool("ON")
+USE_HYPRE            = _cmake_to_bool("OFF")
+USE_LAPACK           = _cmake_to_bool("ON")
+USE_MKL              = _cmake_to_bool("OFF")
+USE_MUMPS            = _cmake_to_bool("OFF")
+USE_PARDISO          = _cmake_to_bool("OFF")
+USE_UMFPACK          = _cmake_to_bool("ON")
+
+NETGEN_DIR = "/Users/gitlab-runner/Library/Python/3.14/lib/python/site-packages"
+
+NGSOLVE_COMPILE_DEFINITIONS         = "HAVE_NETGEN_SOURCES;HAVE_DLFCN_H;HAVE_CXA_DEMANGLE;USE_TIMEOFDAY;MSG_NOSIGNAL=0;TCL;LAPACK;NGS_PYTHON;USE_UMFPACK"
+NGSOLVE_COMPILE_DEFINITIONS_PRIVATE = ""
+NGSOLVE_COMPILE_INCLUDE_DIRS        = ""
+NGSOLVE_COMPILE_OPTIONS             = "$<$<COMPILE_LANGUAGE:CXX>:-std=c++17>;$<$<COMPILE_LANGUAGE:CXX>:-Wno-undefined-var-template;-Wno-vla-extension>;-DMAX_SYS_DIM=3"
+
+NGSOLVE_INSTALL_DIR_PYTHON   = "."
+NGSOLVE_INSTALL_DIR_BIN      = "bin"
+NGSOLVE_INSTALL_DIR_LIB      = "netgen"
+NGSOLVE_INSTALL_DIR_INCLUDE  = "netgen/include"
+NGSOLVE_INSTALL_DIR_CMAKE    = "ngsolve/cmake"
+NGSOLVE_INSTALL_DIR_RES      = "share"
+
+NGSOLVE_VERSION = "6.2.2506-74-gd430159c1"
+NGSOLVE_VERSION_GIT = "v6.2.2506-74-gd430159c1"
+NGSOLVE_VERSION_PYTHON = "6.2.2506.post74.dev0"
+
+NGSOLVE_VERSION_MAJOR = "6"
+NGSOLVE_VERSION_MINOR = "2"
+NGSOLVE_VERSION_TWEAK = "74"
+NGSOLVE_VERSION_PATCH = "2506"
+NGSOLVE_VERSION_HASH = "gd430159c1"
+
+CMAKE_CXX_COMPILER           = "/Library/Developer/CommandLineTools/usr/bin/c++"
+CMAKE_CUDA_COMPILER          = ""
+CMAKE_C_COMPILER             = "/Library/Developer/CommandLineTools/usr/bin/cc"
+CMAKE_LINKER                 = "/Library/Developer/CommandLineTools/usr/bin/ld"
+CMAKE_INSTALL_PREFIX         = "/Users/gitlab-runner/builds/builds/rL7WHzyj/0/ngsolve/ngsolve/_skbuild/macosx-10.15-universal2-3.14/cmake-install"
+CMAKE_CXX_COMPILER_LAUNCHER  = "/usr/local/bin/ccache"
+
+version = NGSOLVE_VERSION_GIT
+
+MKL_LINK = ""
+
+def get_cmake_dir():
+    import os.path as p
+    d_python = p.dirname(p.dirname(p.dirname(__file__)))
+    py_to_cmake = p.relpath(
+            NGSOLVE_INSTALL_DIR_CMAKE,
+            NGSOLVE_INSTALL_DIR_PYTHON
+            )
+    return p.normpath(p.join(d_python,py_to_cmake))

ngsolve/config/config.pyi

@@ -0,0 +1,45 @@
+from __future__ import annotations
+__all__: list[str] = ['BUILD_STUB_FILES', 'BUILD_UMFPACK', 'CMAKE_CUDA_COMPILER', 'CMAKE_CXX_COMPILER', 'CMAKE_CXX_COMPILER_LAUNCHER', 'CMAKE_C_COMPILER', 'CMAKE_INSTALL_PREFIX', 'CMAKE_LINKER', 'ENABLE_UNIT_TESTS', 'INSTALL_DEPENDENCIES', 'MKL_LINK', 'NETGEN_DIR', 'NGSOLVE_COMPILE_DEFINITIONS', 'NGSOLVE_COMPILE_DEFINITIONS_PRIVATE', 'NGSOLVE_COMPILE_INCLUDE_DIRS', 'NGSOLVE_COMPILE_OPTIONS', 'NGSOLVE_INSTALL_DIR_BIN', 'NGSOLVE_INSTALL_DIR_CMAKE', 'NGSOLVE_INSTALL_DIR_INCLUDE', 'NGSOLVE_INSTALL_DIR_LIB', 'NGSOLVE_INSTALL_DIR_PYTHON', 'NGSOLVE_INSTALL_DIR_RES', 'NGSOLVE_VERSION', 'NGSOLVE_VERSION_GIT', 'NGSOLVE_VERSION_HASH', 'NGSOLVE_VERSION_MAJOR', 'NGSOLVE_VERSION_MINOR', 'NGSOLVE_VERSION_PATCH', 'NGSOLVE_VERSION_PYTHON', 'NGSOLVE_VERSION_TWEAK', 'USE_CCACHE', 'USE_HYPRE', 'USE_LAPACK', 'USE_MKL', 'USE_MUMPS', 'USE_PARDISO', 'USE_UMFPACK', 'get_cmake_dir', 'is_python_package', 'version']
+def _cmake_to_bool(s):
+    ...
+def get_cmake_dir():
+    ...
+BUILD_STUB_FILES: bool = True
+BUILD_UMFPACK: bool = False
+CMAKE_CUDA_COMPILER: str = ''
+CMAKE_CXX_COMPILER: str = '/Library/Developer/CommandLineTools/usr/bin/c++'
+CMAKE_CXX_COMPILER_LAUNCHER: str = '/usr/local/bin/ccache'
+CMAKE_C_COMPILER: str = '/Library/Developer/CommandLineTools/usr/bin/cc'
+CMAKE_INSTALL_PREFIX: str = '/Users/gitlab-runner/builds/builds/rL7WHzyj/0/ngsolve/ngsolve/_skbuild/macosx-10.15-universal2-3.14/cmake-install'
+CMAKE_LINKER: str = '/Library/Developer/CommandLineTools/usr/bin/ld'
+ENABLE_UNIT_TESTS: bool = False
+INSTALL_DEPENDENCIES: bool = False
+MKL_LINK: str = ''
+NETGEN_DIR: str = '/Users/gitlab-runner/Library/Python/3.14/lib/python/site-packages'
+NGSOLVE_COMPILE_DEFINITIONS: str = 'HAVE_NETGEN_SOURCES;HAVE_DLFCN_H;HAVE_CXA_DEMANGLE;USE_TIMEOFDAY;MSG_NOSIGNAL=0;TCL;LAPACK;NGS_PYTHON;USE_UMFPACK'
+NGSOLVE_COMPILE_DEFINITIONS_PRIVATE: str = ''
+NGSOLVE_COMPILE_INCLUDE_DIRS: str = ''
+NGSOLVE_COMPILE_OPTIONS: str = '$<$<COMPILE_LANGUAGE:CXX>:-std=c++17>;$<$<COMPILE_LANGUAGE:CXX>:-Wno-undefined-var-template;-Wno-vla-extension>;-DMAX_SYS_DIM=3'
+NGSOLVE_INSTALL_DIR_BIN: str = 'bin'
+NGSOLVE_INSTALL_DIR_CMAKE: str = 'ngsolve/cmake'
+NGSOLVE_INSTALL_DIR_INCLUDE: str = 'netgen/include'
+NGSOLVE_INSTALL_DIR_LIB: str = 'netgen'
+NGSOLVE_INSTALL_DIR_PYTHON: str = '.'
+NGSOLVE_INSTALL_DIR_RES: str = 'share'
+NGSOLVE_VERSION: str = '6.2.2506-74-gd430159c1'
+NGSOLVE_VERSION_GIT: str = 'v6.2.2506-74-gd430159c1'
+NGSOLVE_VERSION_HASH: str = 'gd430159c1'
+NGSOLVE_VERSION_MAJOR: str = '6'
+NGSOLVE_VERSION_MINOR: str = '2'
+NGSOLVE_VERSION_PATCH: str = '2506'
+NGSOLVE_VERSION_PYTHON: str = '6.2.2506.post74.dev0'
+NGSOLVE_VERSION_TWEAK: str = '74'
+USE_CCACHE: bool = True
+USE_HYPRE: bool = False
+USE_LAPACK: bool = True
+USE_MKL: bool = False
+USE_MUMPS: bool = False
+USE_PARDISO: bool = False
+USE_UMFPACK: bool = True
+is_python_package: bool = True
+version: str = 'v6.2.2506-74-gd430159c1'

ngsolve/demos/TensorProduct/tp_dg_1d_1d.py

@@ -0,0 +1,80 @@
+from ngsolve.TensorProductTools import *
+from ngsolve.comp import *
+from ngsolve import *
+import netgen.gui
+
+mesh1 = Mesh(SegMesh(20,0,1,periodic=True) )
+mesh2 = Mesh(SegMesh(20,0,1,periodic=True) )
+
+tpmesh = Mesh(MakeTensorProductMesh(mesh1,mesh2))
+Draw(tpmesh)
+
+n=5
+m=5
+
+fesx = L2(mesh1,order=n)
+fesy = L2(mesh2,order=m)
+
+tpfes = TensorProductFESpace([fesx,fesy])
+
+fes = L2(tpmesh,order=n)
+u = tpfes.TrialFunction()
+v = tpfes.TestFunction()
+vx = v.Operator("gradx")
+vy = v.Operator("grady")
+b = (-1)*CoefficientFunction( (0.25,0.5) )
+
+
+uin = ProlongateCoefficientFunction(IfPos((1.0-x)*(x-0.5),sin((x-0.5)*2*3.14159),0) + IfPos((x-0.0)*(0.5-x),sin((x)*2*3.14159),0) , 0 , tpfes)
+
+gradv = CoefficientFunction((vx,vy))
+
+a = BilinearForm(tpfes)
+
+n = CoefficientFunction((ProlongateCoefficientFunction(specialcf.normal(1),1,tpfes),ProlongateCoefficientFunction(specialcf.normal(1),0,tpfes)))
+bn = b*n
+
+a += SymbolicTPBFI ( (-u * b*gradv).Compile() )
+a += SymbolicTPBFI ( (bn) *IfPos(bn, u, u.Other( )) * (v-v.Other()).Compile(), VOL, skeleton=True)
+a += SymbolicTPBFI ( (bn) *IfPos(bn, u, u.Other(bnd = CoefficientFunction(uin) )) * (v).Compile(), BND, skeleton=True)
+
+
+
+u = GridFunction(tpfes)
+v = GridFunction(tpfes)
+
+uu = GridFunction(fes)
+print("Setting")
+u.Set(exp( ProlongateCoefficientFunction(-200*(x-0.4)*(x-0.4), 1, tpfes)+ ProlongateCoefficientFunction(-200*(x-0.4)*(x-0.4),0,tpfes) ) )
+print("Done")
+Transfer2StdMesh(u,uu)
+
+Draw(uu,sd=2,autoscale=False)
+
+h = u.vec.CreateVector()
+#u.vec[:] = 1.0
+print('To start the simulation type Run(n_steps)!')
+
+def Step():
+    a.Apply(u.vec,v.vec)
+    h.data = 0.00125*v.vec.data
+    tpfes.SolveM(rho=CoefficientFunction(1), vec = h)
+    u.vec.data-=h.data
+    #Redraw()
+
+def Run(nsteps):
+    count = 0
+    with TaskManager():
+        for i in range(nsteps):
+            print("Step ",i+1, "/",nsteps)
+            Step()
+            count += 1
+            if count % 1 == 0:
+                Transfer2StdMesh(u,uu)
+                count = 0
+                Redraw()
+
+Transfer2StdMesh(u,uu)
+Run(100)
+for t in Timers():
+    print(t["counts"], t["time"], t["name"])

ngsolve/demos/TensorProduct/tp_dg_1d_2d.py

@@ -0,0 +1,73 @@
+from ngsolve.TensorProductTools import *
+from ngsolve.comp import *
+from ngsolve import *
+from netgen.geom2d import unit_square
+import netgen.gui
+
+mesh1 = Mesh(SegMesh(20,0,1))
+mesh2 = Mesh(unit_square.GenerateMesh(maxh=0.15))
+
+tpmesh = Mesh(MakeTensorProductMesh(mesh1,mesh2))
+Draw(tpmesh)
+
+n=3
+m=3
+
+fesx = L2(mesh1,order=n)
+fesy = L2(mesh2,order=m)
+tpfes = TensorProductFESpace([fesx,fesy])
+
+fes = L2(tpmesh,order=n)
+
+u = tpfes.TrialFunction()
+v = tpfes.TestFunction()
+
+vx = v.Operator("gradx")
+vy = v.Operator("grady")
+
+b = CoefficientFunction( (ProlongateCoefficientFunction(x-0.5,0,tpfes),ProlongateCoefficientFunction(0.5-x,1,tpfes),0) )
+
+uin = CoefficientFunction(0.0)
+
+
+gradv = CoefficientFunction((vx,vy))
+
+a = BilinearForm(tpfes)
+
+n = CoefficientFunction(( ProlongateCoefficientFunction(specialcf.normal(1)[0],1,tpfes),ProlongateCoefficientFunction(specialcf.normal(2)[0],0,tpfes),ProlongateCoefficientFunction(specialcf.normal(2)[1],0,tpfes)))
+bn = b*n
+
+a += SymbolicTPBFI ( -u * b*gradv )
+a += SymbolicTPBFI ( (bn) *IfPos(bn, u, u.Other(bnd = uin )) * (v-v.Other(bnd = 0.0)), VOL, skeleton=True)
+a += SymbolicTPBFI ( (bn) *IfPos(bn, u, u.Other(bnd = uin )) * (v), BND, skeleton=True)
+
+
+u = GridFunction(tpfes)
+v = GridFunction(tpfes)
+
+
+uu = GridFunction(fes)
+
+u.Set(exp(ProlongateCoefficientFunction( -70*(x-0.25)*(x-0.25),1,tpfes) + ProlongateCoefficientFunction(-70*(x-0.25)*(x-0.25)-70*(y-0.75)*(y-0.75),0,tpfes) ))
+Transfer2StdMesh(u,uu)
+Draw(uu,sd=3,autoscale=False)
+
+h = u.vec.CreateVector()
+print('To start the simulation type Run(n_steps)!')
+def Step():
+    a.Apply(u.vec,v.vec)
+    h.data = 0.001*v.vec.data
+    tpfes.SolveM(rho=CoefficientFunction(1), vec = h)
+    u.vec.data-=h.data
+    Transfer2StdMesh(u,uu)
+    Redraw()
+
+def Run(nsteps):
+    with TaskManager():
+        for i in range(nsteps):
+            print("Step ",i+1, "/",nsteps)
+            Step()
+
+#Run(100)
+for t in Timers():
+    print(t["counts"], t["time"], t["name"])

ngsolve/demos/TensorProduct/tp_dg_2d_1d.py

@@ -0,0 +1,72 @@
+from ngsolve.TensorProductTools import *
+from ngsolve.comp import *
+from ngsolve import *
+from netgen.geom2d import unit_square
+import netgen.gui
+
+mesh1 = Mesh(unit_square.GenerateMesh(maxh=0.15))
+mesh2 = Mesh(SegMesh(20,0,1))
+
+tpmesh = Mesh(MakeTensorProductMesh(mesh1,mesh2))
+Draw(tpmesh)
+
+n=3
+m=3
+
+fesx = L2(mesh1,order=n)
+fesy = L2(mesh2,order=m)
+tpfes = TensorProductFESpace([fesx,fesy])
+
+fes = L2(tpmesh,order=n)
+
+u = tpfes.TrialFunction()
+v = tpfes.TestFunction()
+
+vx = v.Operator("gradx")
+vy = v.Operator("grady")
+
+b = CoefficientFunction( (0,0,1) )
+
+uin = CoefficientFunction(0.0)
+
+gradv = CoefficientFunction((vx,vy))
+
+a = BilinearForm(tpfes)
+
+n = CoefficientFunction((ProlongateCoefficientFunction(specialcf.normal(2)[0],1,tpfes),ProlongateCoefficientFunction(specialcf.normal(2)[1],1,tpfes),ProlongateCoefficientFunction(specialcf.normal(1)[0],0,tpfes)))
+bn = b*n
+
+a += SymbolicTPBFI ( -u * b*gradv )
+a += SymbolicTPBFI ( (bn) *IfPos(bn, u, u.Other(bnd = uin )) * (v-v.Other(bnd = 0.0)), VOL, skeleton=True)
+a += SymbolicTPBFI ( (bn) *IfPos(bn, u, u.Other(bnd = uin )) * (v), BND, skeleton=True)
+
+
+u = GridFunction(tpfes)
+v = GridFunction(tpfes)
+
+
+uu = GridFunction(fes)
+
+u.Set(exp(ProlongateCoefficientFunction(-70*(x-0.125)*(x-0.125)-70*(y-0.125)*(y-0.125),1,tpfes)+ProlongateCoefficientFunction(-70*(x-0.75)*(x-0.75),0,tpfes) ))
+Transfer2StdMesh(u,uu)
+Draw(uu,sd=3,autoscale=False)
+
+h = u.vec.CreateVector()
+print('To start the simulation type Run(n_steps)!')
+def Step():
+    a.Apply(u.vec,v.vec)
+    h.data = 0.001*v.vec.data
+    tpfes.SolveM(rho=CoefficientFunction(1), vec = h)
+    u.vec.data-=h.data
+    Transfer2StdMesh(u,uu)
+    Redraw()
+
+def Run(nsteps):
+    with TaskManager():
+        for i in range(nsteps):
+            print("Step ",i+1, "/",nsteps)
+            Step()
+
+Run(100)
+for t in Timers():
+    print(t["counts"], t["time"], t["name"])

ngsolve/demos/TensorProduct/tp_dg_2d_2d.py

@@ -0,0 +1,66 @@
+from ngsolve.TensorProductTools import *
+from ngsolve.comp import *
+from ngsolve import *
+from netgen.geom2d import unit_square
+import netgen.gui
+
+mesh1 = Mesh(unit_square.GenerateMesh(maxh=0.15))
+mesh2 = Mesh(MakeHexagonalMesh2D(maxh=0.2))
+n=4
+m=4
+SetHeapSize(1000000000)
+
+fesx = L2(mesh1,order=n)
+fesy = L2(mesh2,order=m)
+tpfes = TensorProductFESpace([fesx,fesy])
+
+u = tpfes.TrialFunction()
+v = tpfes.TestFunction()
+
+vx = v.Operator("gradx")
+vy = v.Operator("grady")
+
+b = CoefficientFunction( ( ProlongateCoefficientFunction(y-0.5,1,tpfes),ProlongateCoefficientFunction(0.5-x,1,tpfes),1,1 ) ) 
+
+uin = CoefficientFunction(0.0)
+
+gradv = CoefficientFunction((vx,vy))
+
+a = BilinearForm(tpfes)
+
+n = CoefficientFunction((ProlongateCoefficientFunction(specialcf.normal(2)[0],1,tpfes),ProlongateCoefficientFunction(specialcf.normal(2)[1],1,tpfes),ProlongateCoefficientFunction(specialcf.normal(2)[0],0,tpfes),ProlongateCoefficientFunction(specialcf.normal(2)[1],0,tpfes)))
+bn = b*n
+
+a += SymbolicTPBFI ( -u * b*gradv )
+a += SymbolicTPBFI ( (bn) *IfPos(bn, u, u.Other(bnd = uin )) * (v-v.Other(bnd = 0.0)), VOL, skeleton=True)
+a += SymbolicTPBFI ( (bn) *IfPos(bn, u, u.Other(bnd = uin )) * (v), BND, skeleton=True)
+
+
+u = GridFunction(tpfes)
+v = GridFunction(tpfes)
+print(1)
+with TaskManager():
+    u.Set(exp(ProlongateCoefficientFunction(-90*(x-0.25)*(x-0.25)-90*(y-0.5)*(y-0.5),1,tpfes)+ProlongateCoefficientFunction(-70*(x-0.75)*(x-0.75)-70*(y-0.75)*(y-0.75),0,tpfes) ))
+print(2)
+rho = GridFunction(fesx,name="rho")
+Draw(rho)
+h = u.vec.CreateVector()
+print('To start the simulation type Run(n_steps)!')
+
+def Step():
+    a.Apply(u.vec,v.vec)
+    h.data = 0.005*v.vec.data
+    tpfes.SolveM(rho=CoefficientFunction(1), vec = h)
+    u.vec.data-=h.data
+    TensorProductIntegrate(u,rho)
+    Redraw()
+
+def Run(nsteps):
+    with TaskManager():
+        for i in range(nsteps):
+            print("Step ",i+1, "/",nsteps)
+            Step()
+
+Run(1000)
+for t in Timers():
+    print(t["counts"], t["time"], t["name"])

ngsolve/demos/howto/hhj.py

@@ -0,0 +1,44 @@
+#
+# The Hellan-Herrmann-Johnson method for a Kirchhoff plate
+# 
+# M. I. Comodi: The Hellan–Herrmann-Johnson Method
+# Some error estimates and postprocessing. Math. Comp. 52, 17–39, 1989
+# 
+
+from ngsolve import *
+from netgen.geom2d import unit_square
+
+mesh = Mesh (unit_square.GenerateMesh(maxh=0.05))
+order = 3
+
+V = HDivDiv(mesh, order=order-1)
+Q = H1(mesh, order=order, dirichlet="left|right|top|bottom")
+X = V*Q
+
+print ("ndof-V:", V.ndof, ", ndof-Q:", Q.ndof)
+
+sigma, u = X.TrialFunction()
+tau, v = X.TestFunction()
+
+n = specialcf.normal(2)
+
+def tang(u): return u-(u*n)*n
+
+a = BilinearForm(X, symmetric=True)
+a += (InnerProduct (sigma, tau) + div(sigma)*grad(v) + div(tau)*grad(u) - 1e-10*u*v)*dx
+a += (-(sigma*n) * tang(grad(v)) - (tau*n)*tang(grad(u)))*dx(element_boundary=True)
+a.Assemble()
+
+f = LinearForm(X)
+f += 1 * v * dx
+# f += Trace(tau.Trace()) * ds("bottom")
+f.Assemble()
+
+u = GridFunction(X)
+u.vec.data = a.mat.Inverse(X.FreeDofs()) * f.vec
+
+Draw (u.components[0], mesh, name="sigma")
+Draw (u.components[1], mesh, name="disp")
+
+
+

ngsolve/demos/howto/hybrid_dg.py

@@ -0,0 +1,53 @@
+from netgen.geom2d import unit_square
+from ngsolve import *
+
+
+ngsglobals.msg_level = 1
+mesh = Mesh(unit_square.GenerateMesh(maxh=0.4))
+for k in range(5):
+    mesh.Refine()
+
+
+order = 3
+fes1 = L2(mesh, order=order)
+fes2 = FacetFESpace(mesh, order=order, dirichlet="bottom|right|top")
+
+print ("element dofs: ", fes1.ndof)
+print ("facet dofs: ", fes2.ndof)
+
+fes = fes1*fes2
+
+u,uhat = fes.TrialFunction()
+v,vhat = fes.TestFunction()
+
+n = specialcf.normal(mesh.dim)
+h = specialcf.mesh_size
+
+a = BilinearForm(fes, symmetric=True, condense = True)
+a += grad(u) * grad(v) * dx
+a += (grad(u)*n*(vhat-v)+grad(v)*n*(uhat-u)+10*order*order/h*(u-uhat)*(v-vhat))*dx(element_boundary=True)
+
+c = Preconditioner(type="direct", bf=a, inverse = "sparsecholesky")
+# c = Preconditioner(type="bddc", bf=a)
+
+with TaskManager():
+    a.Assemble()
+    ainv = CGSolver(a.mat, c.mat)
+
+f = LinearForm(fes)
+f += 1*v*dx
+f.Assemble()
+
+u = GridFunction(fes)
+
+
+f.vec.data += a.harmonic_extension_trans * f.vec
+
+u.vec.data = ainv * f.vec
+
+u.vec.data += a.harmonic_extension * u.vec
+u.vec.data += a.inner_solve * f.vec
+
+Draw (u.components[0], mesh, "sol")
+
+

ngsolve/demos/howto/mixed.py

@@ -0,0 +1,30 @@
+from netgen.geom2d import unit_square
+from ngsolve import *
+
+ngsglobals.msg_level = 1
+
+mesh = Mesh(unit_square.GenerateMesh(maxh=0.1))
+
+order = 2
+fes1 = HDiv(mesh, order=order)
+fes2 = L2(mesh, order=order-1)
+
+fes = fes1*fes2
+
+sigma,u = fes.TrialFunction()
+tau,v = fes.TestFunction()
+
+a = BilinearForm(fes)
+a += (sigma*tau + div(sigma)*v + div(tau)*u - 1e-10*u*v)*dx
+# (regularization needed for direct solver)
+a.Assemble()
+
+f = LinearForm(fes)
+f += -v*dx
+f.Assemble()
+
+u = GridFunction(fes)
+u.vec.data = a.mat.Inverse(fes.FreeDofs()) * f.vec
+
+Draw (u.components[0], mesh, "flux")
+Draw (u.components[1], mesh, "sol")