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

ngsolve/bvp.py

@@ -0,0 +1,78 @@
+
+def BVP(bf, lf, gf, \
+            pre=None, pre_flags={}, \
+            solver=None, solver_flags={},
+            maxsteps=200, tol=1e-8, print=True, inverse="umfpack",
+            needsassembling=True):
+    """
+    Solve a linear boundary value problem A(u,v) = f(v)
+
+Parameters
+----------
+
+bf : BilinearForm
+  provides the matrix. 
+
+lf : LinearForm
+  provides the right hand side.
+
+gf : GridFunction
+  provides the solution vector
+
+pre : Basematrix or class or string = None
+  used if an iterative solver is used
+  can be one of 
+    * a preconditioner object
+    * a preconditioner class
+    * a preconditioner class name
+
+pre_flags : dictionary = { }
+  flags used to create preconditioner
+
+
+    """
+    from ngsolve import Projector, Preconditioner
+    from ngsolve.krylovspace import CG
+
+
+    if isinstance(pre,type):
+        pre = pre(bf, **pre_flags)
+        if not needsassembling:
+            pre.Update()
+
+    if isinstance(pre,str):
+        pre = Preconditioner(bf, pre, **pre_flags)
+        if not needsassembling:
+            pre.Update()
+
+        
+    if needsassembling:
+        bf.Assemble()
+        lf.Assemble()
+    
+    r = lf.vec.CreateVector()
+    r.data = lf.vec
+    
+    if bf.condense:
+        r.data += bf.harmonic_extension_trans * r
+
+        # zero local dofs
+        innerbits = gf.space.FreeDofs(False) & ~gf.space.FreeDofs(True)
+        Projector(innerbits, False).Project(gf.vec)
+
+    if solver:
+        inv = solver(mat=bf.mat, pre=pre, **solver_flags)
+        r.data -= bf.mat * gf.vec
+        gf.vec.data += inv * r
+                   
+    elif pre:
+        CG(mat = bf.mat, rhs = r, pre=pre, sol=gf.vec, tol=tol, maxsteps=maxsteps, initialize=False, printrates=print)
+    else:
+        inv = bf.mat.Inverse(gf.space.FreeDofs(bf.condense), inverse=inverse)
+        r.data -= bf.mat * gf.vec
+        gf.vec.data += inv * r
+
+    if bf.condense:
+        gf.vec.data += bf.harmonic_extension * gf.vec
+        gf.vec.data += bf.inner_solve * r
+            

ngsolve/cmake/NGSolveConfig.cmake

@@ -0,0 +1,102 @@
+set(PACKAGE_VERSION "6.2.2506-74-gd430159c1")
+find_package(Netgen CONFIG REQUIRED HINTS
+  ${CMAKE_CURRENT_LIST_DIR}
+  ${CMAKE_CURRENT_LIST_DIR}/..
+  ${CMAKE_CURRENT_LIST_DIR}/../netgen
+)
+
+get_filename_component(NGSOLVE_DIR  "${NETGEN_DIR}"  ABSOLUTE)
+
+get_filename_component(NGSOLVE_INSTALL_DIR "${CMAKE_CURRENT_LIST_DIR}/../../" ABSOLUTE)
+get_filename_component(NGSOLVE_INCLUDE_DIR  "${NETGEN_INCLUDE_DIR}"  ABSOLUTE)
+get_filename_component(NGSOLVE_BINARY_DIR   "${NETGEN_BINARY_DIR}"   ABSOLUTE)
+get_filename_component(NGSOLVE_LIBRARY_DIR  "${NETGEN_LIBRARY_DIR}"  ABSOLUTE)
+get_filename_component(NGSOLVE_PYTHON_DIR   "${NETGEN_PYTHON_DIR}"   ABSOLUTE)
+get_filename_component(NGSOLVE_RESOURCE_DIR "${NETGEN_RESOURCE_DIR}" ABSOLUTE)
+
+set(NGSOLVE_CXX_COMPILER "C:/Program Files/Microsoft Visual Studio/2022/Community/VC/Tools/MSVC/14.44.35207/bin/Hostx64/x64/cl.exe")
+set(NGSOLVE_CMAKE_BUILD_TYPE "Release")
+
+set(NGSOLVE_CMAKE_THREAD_LIBS_INIT "")
+set(NGSOLVE_MKL_LIBRARIES "C:/gitlabci/tools/builds/3zsqG5ns/0/ngsolve/venv_ngs/Library/lib/mkl_rt.lib")
+set(NGSOLVE_PYBIND_INCLUDE_DIR "")
+set(NGSOLVE_PYTHON_INCLUDE_DIRS "C:/Python314/Include")
+set(NGSOLVE_PYTHON_LIBRARIES "C:/Python314/libs/python314.lib")
+set(NGSOLVE_PYTHON_PACKAGES_INSTALL_DIR  "")
+set(NGSOLVE_TCL_INCLUDE_PATH "C:/gitlabci/tools/builds/3zsqG5ns/0/ngsolve/netgen/_skbuild/win-amd64-3.14/cmake-build/dependencies/src/project_tcl/generic")
+set(NGSOLVE_TCL_LIBRARY "C:/Python314/tcl/tclstub86.lib")
+set(NGSOLVE_TK_DND_LIBRARY "")
+set(NGSOLVE_TK_INCLUDE_PATH "C:/gitlabci/tools/builds/3zsqG5ns/0/ngsolve/netgen/_skbuild/win-amd64-3.14/cmake-build/dependencies/src/project_tk/generic")
+set(NGSOLVE_TK_LIBRARY "C:/Python314/tcl/tkstub86.lib")
+set(NGSOLVE_X11_X11_LIB "")
+set(NGSOLVE_X11_Xmu_LIB "")
+set(NGSOLVE_ZLIB_INCLUDE_DIRS "C:/gitlabci/tools/builds/3zsqG5ns/0/ngsolve/netgen/_skbuild/win-amd64-3.14/cmake-build/dependencies/zlib/include")
+set(NGSOLVE_ZLIB_LIBRARIES "C:/gitlabci/tools/builds/3zsqG5ns/0/ngsolve/netgen/_skbuild/win-amd64-3.14/cmake-build/dependencies/zlib/lib/zlibstatic.lib")
+
+set(NGSOLVE_INTEL_MIC       OFF)
+set(NGSOLVE_USE_CCACHE ON)
+set(NGSOLVE_USE_CUDA        OFF)
+set(NGSOLVE_USE_GUI ON)
+set(NGSOLVE_USE_LAPACK      ON)
+set(NGSOLVE_USE_MKL         ON)
+set(NGSOLVE_USE_MPI OFF)
+set(NGSOLVE_USE_MUMPS       OFF)
+set(NGSOLVE_USE_NUMA        )
+set(NGSOLVE_USE_PARDISO     OFF)
+set(NGSOLVE_USE_PYTHON ON)
+set(NGSOLVE_USE_UMFPACK     OFF)
+set(NGSOLVE_USE_VT          )
+set(NGSOLVE_USE_VTUNE       OFF)
+set(NGSOLVE_MAX_SYS_DIM       3)
+
+set(NGSOLVE_COMPILE_FLAGS " /DHAVE_NETGEN_SOURCES /DUSE_TIMEOFDAY /DTCL /DLAPACK /DUSE_PARDISO /DNGS_PYTHON /DNETGEN_PYTHON /DNG_PYTHON /DPYBIND11_SIMPLE_GIL_MANAGEMENT /D_WIN32_WINNT=0x1000 /DWNT /DWNT_WINDOW /DNOMINMAX /DMSVC_EXPRESS /D_CRT_SECURE_NO_WARNINGS /DHAVE_STRUCT_TIMESPEC /DWIN32" CACHE STRING "Preprocessor definitions of ngscxx")
+set(NGSOLVE_LINK_FLAGS " -LC:/gitlabci/tools/builds/3zsqG5ns/0/ngsolve/venv_ngs/Library/lib -l_rt /LIBPATH:\"%NGSCXX_DIR%/lib\" nglib.lib ngcore.lib libngsolve.lib /LIBPATH:\"%PYTHON_LIBDIR%\"" CACHE STRING "Link flags set in ngsld")
+set(NGSOLVE_INCLUDE_DIRS C:/gitlabci/tools/builds/3zsqG5ns/0/ngsolve/venv_ngs/Library/include;C:/Python314/Include CACHE STRING "Include dirs set in ngscxx")
+
+set(NGSOLVE_INSTALL_DIR_PYTHON  .)
+set(NGSOLVE_INSTALL_DIR_BIN     netgen)
+set(NGSOLVE_INSTALL_DIR_LIB     netgen/lib)
+set(NGSOLVE_INSTALL_DIR_INCLUDE netgen/include)
+set(NGSOLVE_INSTALL_DIR_CMAKE   ngsolve/cmake)
+set(NGSOLVE_INSTALL_DIR_RES     share)
+
+include(${CMAKE_CURRENT_LIST_DIR}/ngsolve-targets.cmake)
+
+if(WIN32)
+	# Libraries like ngstd, ngbla etc. are only dummy libs on Windows (there is only one library libngsolve.dll)
+	# make sure that all those dummy libs link ngsolve when used with CMake
+	target_link_libraries(ngstd INTERFACE ngsolve)
+endif(WIN32)
+
+if(NGSOLVE_USE_PYTHON)
+  function(add_ngsolve_python_module target)
+    if(NOT CMAKE_BUILD_TYPE)
+      message(STATUS "Setting build type to NGSolve build type: ${NGSOLVE_CMAKE_BUILD_TYPE}")
+      set(CMAKE_BUILD_TYPE ${NGSOLVE_CMAKE_BUILD_TYPE} PARENT_SCOPE)
+    endif()
+    add_library(${target} SHARED ${ARGN})
+
+    find_package(PythonInterp 3 REQUIRED)
+    find_package(PythonLibs 3 REQUIRED)
+    target_include_directories(${target} PRIVATE ${PYTHON_INCLUDE_DIR})
+
+    if(NETGEN_BUILD_FOR_CONDA AND NOT WIN32)
+        if(APPLE)
+            target_link_options(${target} PUBLIC -undefined dynamic_lookup)
+        endif(APPLE)
+    else(NETGEN_BUILD_FOR_CONDA AND NOT WIN32)
+        target_link_libraries(${target} PUBLIC ${PYTHON_LIBRARY})
+    endif(NETGEN_BUILD_FOR_CONDA AND NOT WIN32)
+
+    set_target_properties(${target} PROPERTIES PREFIX "" CXX_STANDARD 17)
+    target_link_libraries(${target} PUBLIC ngsolve)
+
+    if(WIN32)
+      set_target_properties( ${target} PROPERTIES SUFFIX ".pyd" )
+    else(WIN32)
+      set_target_properties(${target} PROPERTIES SUFFIX ".so")
+    endif(WIN32)
+
+    set_target_properties(${target} PROPERTIES INSTALL_RPATH "${NGSOLVE_LIBRARY_DIR}")
+endfunction()
+endif(NGSOLVE_USE_PYTHON)

ngsolve/cmake/ngsolve-targets-release.cmake

@@ -0,0 +1,19 @@
+#----------------------------------------------------------------
+# Generated CMake target import file for configuration "Release".
+#----------------------------------------------------------------
+
+# Commands may need to know the format version.
+set(CMAKE_IMPORT_FILE_VERSION 1)
+
+# Import target "ngsolve" for configuration "Release"
+set_property(TARGET ngsolve APPEND PROPERTY IMPORTED_CONFIGURATIONS RELEASE)
+set_target_properties(ngsolve PROPERTIES
+  IMPORTED_IMPLIB_RELEASE "${_IMPORT_PREFIX}/netgen/lib/libngsolve.lib"
+  IMPORTED_LOCATION_RELEASE "${_IMPORT_PREFIX}/netgen/libngsolve.dll"
+  )
+
+list(APPEND _cmake_import_check_targets ngsolve )
+list(APPEND _cmake_import_check_files_for_ngsolve "${_IMPORT_PREFIX}/netgen/lib/libngsolve.lib" "${_IMPORT_PREFIX}/netgen/libngsolve.dll" )
+
+# Commands beyond this point should not need to know the version.
+set(CMAKE_IMPORT_FILE_VERSION)

ngsolve/cmake/ngsolve-targets.cmake

@@ -0,0 +1,164 @@
+# Generated by CMake
+
+if("${CMAKE_MAJOR_VERSION}.${CMAKE_MINOR_VERSION}" LESS 2.8)
+   message(FATAL_ERROR "CMake >= 2.8.0 required")
+endif()
+if(CMAKE_VERSION VERSION_LESS "2.8.3")
+   message(FATAL_ERROR "CMake >= 2.8.3 required")
+endif()
+cmake_policy(PUSH)
+cmake_policy(VERSION 2.8.3...3.25)
+#----------------------------------------------------------------
+# Generated CMake target import file.
+#----------------------------------------------------------------
+
+# Commands may need to know the format version.
+set(CMAKE_IMPORT_FILE_VERSION 1)
+
+# Protect against multiple inclusion, which would fail when already imported targets are added once more.
+set(_cmake_targets_defined "")
+set(_cmake_targets_not_defined "")
+set(_cmake_expected_targets "")
+foreach(_cmake_expected_target IN ITEMS netgen_libs ngs_lapack ngstd ngbla ngla ngfem ngsbem ngcomp ngsolve)
+  list(APPEND _cmake_expected_targets "${_cmake_expected_target}")
+  if(TARGET "${_cmake_expected_target}")
+    list(APPEND _cmake_targets_defined "${_cmake_expected_target}")
+  else()
+    list(APPEND _cmake_targets_not_defined "${_cmake_expected_target}")
+  endif()
+endforeach()
+unset(_cmake_expected_target)
+if(_cmake_targets_defined STREQUAL _cmake_expected_targets)
+  unset(_cmake_targets_defined)
+  unset(_cmake_targets_not_defined)
+  unset(_cmake_expected_targets)
+  unset(CMAKE_IMPORT_FILE_VERSION)
+  cmake_policy(POP)
+  return()
+endif()
+if(NOT _cmake_targets_defined STREQUAL "")
+  string(REPLACE ";" ", " _cmake_targets_defined_text "${_cmake_targets_defined}")
+  string(REPLACE ";" ", " _cmake_targets_not_defined_text "${_cmake_targets_not_defined}")
+  message(FATAL_ERROR "Some (but not all) targets in this export set were already defined.\nTargets Defined: ${_cmake_targets_defined_text}\nTargets not yet defined: ${_cmake_targets_not_defined_text}\n")
+endif()
+unset(_cmake_targets_defined)
+unset(_cmake_targets_not_defined)
+unset(_cmake_expected_targets)
+
+
+# Compute the installation prefix relative to this file.
+get_filename_component(_IMPORT_PREFIX "${CMAKE_CURRENT_LIST_FILE}" PATH)
+get_filename_component(_IMPORT_PREFIX "${_IMPORT_PREFIX}" PATH)
+get_filename_component(_IMPORT_PREFIX "${_IMPORT_PREFIX}" PATH)
+if(_IMPORT_PREFIX STREQUAL "/")
+  set(_IMPORT_PREFIX "")
+endif()
+
+# Create imported target netgen_libs
+add_library(netgen_libs INTERFACE IMPORTED)
+
+set_target_properties(netgen_libs PROPERTIES
+  INTERFACE_LINK_LIBRARIES "nglib;ngcore"
+)
+
+# Create imported target ngs_lapack
+add_library(ngs_lapack INTERFACE IMPORTED)
+
+set_target_properties(ngs_lapack PROPERTIES
+  INTERFACE_INCLUDE_DIRECTORIES "C:/gitlabci/tools/builds/3zsqG5ns/0/ngsolve/venv_ngs/Library/include"
+  INTERFACE_LINK_LIBRARIES "C:/gitlabci/tools/builds/3zsqG5ns/0/ngsolve/venv_ngs/Library/lib/mkl_rt.lib"
+)
+
+# Create imported target ngstd
+add_library(ngstd INTERFACE IMPORTED)
+
+set_target_properties(ngstd PROPERTIES
+  INTERFACE_COMPILE_DEFINITIONS "HAVE_NETGEN_SOURCES;USE_TIMEOFDAY;TCL;LAPACK;USE_PARDISO;NGS_PYTHON"
+  INTERFACE_COMPILE_OPTIONS "/std:c++17;/bigobj;/wd4068;-DMAX_SYS_DIM=3"
+  INTERFACE_LINK_LIBRARIES "netgen_libs;\$<LINK_ONLY:>"
+)
+
+# Create imported target ngbla
+add_library(ngbla INTERFACE IMPORTED)
+
+set_target_properties(ngbla PROPERTIES
+  INTERFACE_LINK_LIBRARIES "ngstd;\$<LINK_ONLY:>"
+)
+
+# Create imported target ngla
+add_library(ngla INTERFACE IMPORTED)
+
+set_target_properties(ngla PROPERTIES
+  INTERFACE_LINK_LIBRARIES "ngbla;ngstd;\$<LINK_ONLY:>"
+)
+
+# Create imported target ngfem
+add_library(ngfem INTERFACE IMPORTED)
+
+set_target_properties(ngfem PROPERTIES
+  INTERFACE_LINK_LIBRARIES "ngbla;ngstd;\$<LINK_ONLY:>"
+)
+
+# Create imported target ngsbem
+add_library(ngsbem INTERFACE IMPORTED)
+
+set_target_properties(ngsbem PROPERTIES
+  INTERFACE_LINK_LIBRARIES "ngcomp;ngfem;ngla;ngbla;ngstd;\$<LINK_ONLY:>"
+)
+
+# Create imported target ngcomp
+add_library(ngcomp INTERFACE IMPORTED)
+
+set_target_properties(ngcomp PROPERTIES
+  INTERFACE_LINK_LIBRARIES "ngfem;ngla;ngbla;ngstd;\$<LINK_ONLY:>"
+)
+
+# Create imported target ngsolve
+add_library(ngsolve SHARED IMPORTED)
+
+set_target_properties(ngsolve PROPERTIES
+  INTERFACE_LINK_LIBRARIES "netgen_libs;ngsolve;ngcomp;ngfem;ngsbem;ngbla;ngla;ngstd"
+)
+
+if(CMAKE_VERSION VERSION_LESS 3.0.0)
+  message(FATAL_ERROR "This file relies on consumers using CMake 3.0.0 or greater.")
+endif()
+
+# Load information for each installed configuration.
+file(GLOB _cmake_config_files "${CMAKE_CURRENT_LIST_DIR}/ngsolve-targets-*.cmake")
+foreach(_cmake_config_file IN LISTS _cmake_config_files)
+  include("${_cmake_config_file}")
+endforeach()
+unset(_cmake_config_file)
+unset(_cmake_config_files)
+
+# Cleanup temporary variables.
+set(_IMPORT_PREFIX)
+
+# Loop over all imported files and verify that they actually exist
+foreach(_cmake_target IN LISTS _cmake_import_check_targets)
+  foreach(_cmake_file IN LISTS "_cmake_import_check_files_for_${_cmake_target}")
+    if(NOT EXISTS "${_cmake_file}")
+      message(FATAL_ERROR "The imported target \"${_cmake_target}\" references the file
+   \"${_cmake_file}\"
+but this file does not exist.  Possible reasons include:
+* The file was deleted, renamed, or moved to another location.
+* An install or uninstall procedure did not complete successfully.
+* The installation package was faulty and contained
+   \"${CMAKE_CURRENT_LIST_FILE}\"
+but not all the files it references.
+")
+    endif()
+  endforeach()
+  unset(_cmake_file)
+  unset("_cmake_import_check_files_for_${_cmake_target}")
+endforeach()
+unset(_cmake_target)
+unset(_cmake_import_check_targets)
+
+# This file does not depend on other imported targets which have
+# been exported from the same project but in a separate export set.
+
+# Commands beyond this point should not need to know the version.
+set(CMAKE_IMPORT_FILE_VERSION)
+cmake_policy(POP)

ngsolve/config/__init__.py

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

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("OFF")
+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("ON")
+USE_MUMPS            = _cmake_to_bool("OFF")
+USE_PARDISO          = _cmake_to_bool("OFF")
+USE_UMFPACK          = _cmake_to_bool("OFF")
+
+NETGEN_DIR = "C:/gitlabci/tools/builds/3zsqG5ns/0/ngsolve/venv_ngs/Lib/site-packages"
+
+NGSOLVE_COMPILE_DEFINITIONS         = "HAVE_NETGEN_SOURCES;USE_TIMEOFDAY;TCL;LAPACK;USE_PARDISO;NGS_PYTHON"
+NGSOLVE_COMPILE_DEFINITIONS_PRIVATE = "USE_MKL"
+NGSOLVE_COMPILE_INCLUDE_DIRS        = ""
+NGSOLVE_COMPILE_OPTIONS             = "/std:c++17;/bigobj;/wd4068;-DMAX_SYS_DIM=3"
+
+NGSOLVE_INSTALL_DIR_PYTHON   = "."
+NGSOLVE_INSTALL_DIR_BIN      = "netgen"
+NGSOLVE_INSTALL_DIR_LIB      = "netgen/lib"
+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           = "C:/Program Files/Microsoft Visual Studio/2022/Community/VC/Tools/MSVC/14.44.35207/bin/Hostx64/x64/cl.exe"
+CMAKE_CUDA_COMPILER          = ""
+CMAKE_C_COMPILER             = "C:/Program Files/Microsoft Visual Studio/2022/Community/VC/Tools/MSVC/14.44.35207/bin/Hostx64/x64/cl.exe"
+CMAKE_LINKER                 = "C:/Program Files/Microsoft Visual Studio/2022/Community/VC/Tools/MSVC/14.44.35207/bin/Hostx64/x64/link.exe"
+CMAKE_INSTALL_PREFIX         = "C:/gitlabci/tools/builds/3zsqG5ns/0/ngsolve/ngsolve/_skbuild/win-amd64-3.14/cmake-install"
+CMAKE_CXX_COMPILER_LAUNCHER  = ""
+
+version = NGSOLVE_VERSION_GIT
+
+MKL_LINK = "sdl"
+
+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/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"])