ngsolve 6.2.2501.post21.dev1__cp313-cp313-macosx_10_15_universal2.whl → 6.2.2501.post34.dev1__cp313-cp313-macosx_10_15_universal2.whl

ngsolve/comp/pml.pyi

@@ -1,89 +0,0 @@
-"""
-module for perfectly matched layers
-"""
-from __future__ import annotations
-import ngsolve.bla
-import ngsolve.fem
-import typing
-__all__ = ['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: float = 1, alpha: complex = 1j) -> PML:
-    """
-    radial pml transformation
-    
-    origin is a list/tuple with as many entries as dimenson
-    """

ngsolve/config/__init__.py

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

ngsolve/config/__init__.pyi

@@ -1,43 +0,0 @@
-from __future__ import annotations
-from ngsolve.config.config import get_cmake_dir
-from . import config
-__all__ = ['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.13/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.13/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.2501-21-g62c380ea8'
-NGSOLVE_VERSION_GIT: str = 'v6.2.2501-21-g62c380ea8'
-NGSOLVE_VERSION_HASH: str = 'g62c380ea8'
-NGSOLVE_VERSION_MAJOR: str = '6'
-NGSOLVE_VERSION_MINOR: str = '2'
-NGSOLVE_VERSION_PATCH: str = '2501'
-NGSOLVE_VERSION_PYTHON: str = '6.2.2501.post21.dev1'
-NGSOLVE_VERSION_TWEAK: str = '21'
-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.2501-21-g62c380ea8'

ngsolve/config/__main__.py

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

ngsolve/config/config.py

@@ -1,60 +0,0 @@
-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.13/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.2501-21-g62c380ea8"
-NGSOLVE_VERSION_GIT = "v6.2.2501-21-g62c380ea8"
-NGSOLVE_VERSION_PYTHON = "6.2.2501.post21.dev1"
-
-NGSOLVE_VERSION_MAJOR = "6"
-NGSOLVE_VERSION_MINOR = "2"
-NGSOLVE_VERSION_TWEAK = "21"
-NGSOLVE_VERSION_PATCH = "2501"
-NGSOLVE_VERSION_HASH = "g62c380ea8"
-
-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.13/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

@@ -1,45 +0,0 @@
-from __future__ import annotations
-__all__ = ['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.13/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.13/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.2501-21-g62c380ea8'
-NGSOLVE_VERSION_GIT: str = 'v6.2.2501-21-g62c380ea8'
-NGSOLVE_VERSION_HASH: str = 'g62c380ea8'
-NGSOLVE_VERSION_MAJOR: str = '6'
-NGSOLVE_VERSION_MINOR: str = '2'
-NGSOLVE_VERSION_PATCH: str = '2501'
-NGSOLVE_VERSION_PYTHON: str = '6.2.2501.post21.dev1'
-NGSOLVE_VERSION_TWEAK: str = '21'
-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.2501-21-g62c380ea8'

ngsolve/demos/TensorProduct/tp_dg_1d_1d.py

@@ -1,80 +0,0 @@
-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

@@ -1,73 +0,0 @@
-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

@@ -1,72 +0,0 @@
-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

@@ -1,66 +0,0 @@
-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

@@ -1,44 +0,0 @@
-#
-# 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

@@ -1,53 +0,0 @@
-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

@@ -1,30 +0,0 @@
-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")