ngsolve 6.2.2506.post216.dev0__cp314-cp314-manylinux_2_28_x86_64.whl

ngsolve-6.2.2506.post216.dev0.data/data/lib/cmake/ngsolve/ngsolve-targets-release.cmake

@@ -0,0 +1,89 @@
+#----------------------------------------------------------------
+# Generated CMake target import file for configuration "Release".
+#----------------------------------------------------------------
+
+# Commands may need to know the format version.
+set(CMAKE_IMPORT_FILE_VERSION 1)
+
+# Import target "ngstd" for configuration "Release"
+set_property(TARGET ngstd APPEND PROPERTY IMPORTED_CONFIGURATIONS RELEASE)
+set_target_properties(ngstd PROPERTIES
+  IMPORTED_LOCATION_RELEASE "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngstd.so"
+  IMPORTED_SONAME_RELEASE "libngstd.so"
+  )
+
+list(APPEND _cmake_import_check_targets ngstd )
+list(APPEND _cmake_import_check_files_for_ngstd "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngstd.so" )
+
+# Import target "ngbla" for configuration "Release"
+set_property(TARGET ngbla APPEND PROPERTY IMPORTED_CONFIGURATIONS RELEASE)
+set_target_properties(ngbla PROPERTIES
+  IMPORTED_LOCATION_RELEASE "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngbla.so"
+  IMPORTED_SONAME_RELEASE "libngbla.so"
+  )
+
+list(APPEND _cmake_import_check_targets ngbla )
+list(APPEND _cmake_import_check_files_for_ngbla "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngbla.so" )
+
+# Import target "ngla" for configuration "Release"
+set_property(TARGET ngla APPEND PROPERTY IMPORTED_CONFIGURATIONS RELEASE)
+set_target_properties(ngla PROPERTIES
+  IMPORTED_LOCATION_RELEASE "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngla.so"
+  IMPORTED_SONAME_RELEASE "libngla.so"
+  )
+
+list(APPEND _cmake_import_check_targets ngla )
+list(APPEND _cmake_import_check_files_for_ngla "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngla.so" )
+
+# Import target "ngfem" for configuration "Release"
+set_property(TARGET ngfem APPEND PROPERTY IMPORTED_CONFIGURATIONS RELEASE)
+set_target_properties(ngfem PROPERTIES
+  IMPORTED_LOCATION_RELEASE "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngfem.so"
+  IMPORTED_SONAME_RELEASE "libngfem.so"
+  )
+
+list(APPEND _cmake_import_check_targets ngfem )
+list(APPEND _cmake_import_check_files_for_ngfem "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngfem.so" )
+
+# Import target "ngsbem" for configuration "Release"
+set_property(TARGET ngsbem APPEND PROPERTY IMPORTED_CONFIGURATIONS RELEASE)
+set_target_properties(ngsbem PROPERTIES
+  IMPORTED_LOCATION_RELEASE "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngsbem.so"
+  IMPORTED_SONAME_RELEASE "libngsbem.so"
+  )
+
+list(APPEND _cmake_import_check_targets ngsbem )
+list(APPEND _cmake_import_check_files_for_ngsbem "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngsbem.so" )
+
+# Import target "ngcomp" for configuration "Release"
+set_property(TARGET ngcomp APPEND PROPERTY IMPORTED_CONFIGURATIONS RELEASE)
+set_target_properties(ngcomp PROPERTIES
+  IMPORTED_LOCATION_RELEASE "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngcomp.so"
+  IMPORTED_SONAME_RELEASE "libngcomp.so"
+  )
+
+list(APPEND _cmake_import_check_targets ngcomp )
+list(APPEND _cmake_import_check_files_for_ngcomp "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngcomp.so" )
+
+# Import target "ngsolve" for configuration "Release"
+set_property(TARGET ngsolve APPEND PROPERTY IMPORTED_CONFIGURATIONS RELEASE)
+set_target_properties(ngsolve PROPERTIES
+  IMPORTED_LOCATION_RELEASE "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngsolve.so"
+  IMPORTED_SONAME_RELEASE "libngsolve.so"
+  )
+
+list(APPEND _cmake_import_check_targets ngsolve )
+list(APPEND _cmake_import_check_files_for_ngsolve "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngsolve.so" )
+
+# Import target "ngscudalib" for configuration "Release"
+set_property(TARGET ngscudalib APPEND PROPERTY IMPORTED_CONFIGURATIONS RELEASE)
+set_target_properties(ngscudalib PROPERTIES
+  IMPORTED_LOCATION_RELEASE "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngscudalib.so"
+  IMPORTED_SONAME_RELEASE "libngscudalib.so"
+  )
+
+list(APPEND _cmake_import_check_targets ngscudalib )
+list(APPEND _cmake_import_check_files_for_ngscudalib "${_IMPORT_PREFIX}/lib/python3.14/site-packages/netgen_mesher.libs/libngscudalib.so" )
+
+# Commands beyond this point should not need to know the version.
+set(CMAKE_IMPORT_FILE_VERSION)

ngsolve-6.2.2506.post216.dev0.data/data/lib/cmake/ngsolve/ngsolve-targets.cmake

@@ -0,0 +1,180 @@
+# Generated by CMake
+
+if("${CMAKE_MAJOR_VERSION}.${CMAKE_MINOR_VERSION}" LESS 2.8)
+   message(FATAL_ERROR "CMake >= 3.0.0 required")
+endif()
+if(CMAKE_VERSION VERSION_LESS "3.0.0")
+   message(FATAL_ERROR "CMake >= 3.0.0 required")
+endif()
+cmake_policy(PUSH)
+cmake_policy(VERSION 3.0.0...3.31)
+#----------------------------------------------------------------
+# 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_lib netgen_core_lib ngs_lapack ngstd ngbla ngla ngfem ngsbem ngcomp ngsolve ngscudalib)
+  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)
+get_filename_component(_IMPORT_PREFIX "${_IMPORT_PREFIX}" PATH)
+if(_IMPORT_PREFIX STREQUAL "/")
+  set(_IMPORT_PREFIX "")
+endif()
+
+# Create imported target netgen_lib
+add_library(netgen_lib INTERFACE IMPORTED)
+
+set_target_properties(netgen_lib PROPERTIES
+  INTERFACE_LINK_LIBRARIES "nglib"
+)
+
+# Create imported target netgen_core_lib
+add_library(netgen_core_lib INTERFACE IMPORTED)
+
+set_target_properties(netgen_core_lib PROPERTIES
+  INTERFACE_LINK_LIBRARIES "ngcore"
+)
+
+# Create imported target ngs_lapack
+add_library(ngs_lapack INTERFACE IMPORTED)
+
+set_target_properties(ngs_lapack PROPERTIES
+  INTERFACE_INCLUDE_DIRECTORIES "/opt/_internal/cpython-3.14.0/include"
+  INTERFACE_LINK_LIBRARIES "MKL::MKL"
+)
+
+# Create imported target ngstd
+add_library(ngstd SHARED IMPORTED)
+
+set_target_properties(ngstd PROPERTIES
+  INTERFACE_COMPILE_DEFINITIONS "HAVE_NETGEN_SOURCES;HAVE_DLFCN_H;HAVE_CXA_DEMANGLE;USE_TIMEOFDAY;TCL;LAPACK;USE_PARDISO;NGS_PYTHON"
+  INTERFACE_COMPILE_OPTIONS "\$<\$<COMPILE_LANGUAGE:CXX>:-std=c++17>;-DMAX_SYS_DIM=3"
+  INTERFACE_LINK_LIBRARIES "netgen_core_lib"
+)
+
+# Create imported target ngbla
+add_library(ngbla SHARED IMPORTED)
+
+set_target_properties(ngbla PROPERTIES
+  INTERFACE_LINK_LIBRARIES "ngstd"
+)
+
+# Create imported target ngla
+add_library(ngla SHARED IMPORTED)
+
+set_target_properties(ngla PROPERTIES
+  INTERFACE_LINK_LIBRARIES "ngbla;ngstd"
+)
+
+# Create imported target ngfem
+add_library(ngfem SHARED IMPORTED)
+
+set_target_properties(ngfem PROPERTIES
+  INTERFACE_LINK_LIBRARIES "ngbla;ngstd"
+)
+
+# Create imported target ngsbem
+add_library(ngsbem SHARED IMPORTED)
+
+set_target_properties(ngsbem PROPERTIES
+  INTERFACE_LINK_LIBRARIES "ngcomp;ngfem;ngla;ngbla;ngstd"
+)
+
+# Create imported target ngcomp
+add_library(ngcomp SHARED IMPORTED)
+
+set_target_properties(ngcomp PROPERTIES
+  INTERFACE_LINK_LIBRARIES "ngfem;ngla;ngbla;ngstd;netgen_lib"
+)
+
+# Create imported target ngsolve
+add_library(ngsolve SHARED IMPORTED)
+
+set_target_properties(ngsolve PROPERTIES
+  INTERFACE_LINK_LIBRARIES "netgen_lib;netgen_core_lib;ngsolve;ngcomp;ngfem;ngsbem;ngbla;ngla;ngstd"
+)
+
+# Create imported target ngscudalib
+add_library(ngscudalib SHARED IMPORTED)
+
+set_target_properties(ngscudalib PROPERTIES
+  INTERFACE_COMPILE_OPTIONS "\$<\$<COMPILE_LANGUAGE:CUDA>:--expt-relaxed-constexpr;--extended-lambda;--diag-suppress=611;--diag-suppress=20011;--diag-suppress=20012;--diag-suppress=20013;--diag-suppress=20014;--diag-suppress=20015;-rdc=true>"
+  INTERFACE_LINK_LIBRARIES "CUDA::cusparse;CUDA::cublas;CUDA::cudart;ngsolve"
+)
+
+# 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)
+  if(CMAKE_VERSION VERSION_LESS "3.28"
+      OR NOT DEFINED _cmake_import_check_xcframework_for_${_cmake_target}
+      OR NOT IS_DIRECTORY "${_cmake_import_check_xcframework_for_${_cmake_target}}")
+    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()
+  endif()
+  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-6.2.2506.post216.dev0.data/data/lib/python3.14/site-packages/ngsolve/TensorProductTools.py

@@ -0,0 +1,210 @@
+from netgen.meshing import *
+from netgen.csg import *
+#from ngsolve import *
+from ngsolve.comp import TensorProductFESpace, Transfer2StdMesh, SymbolicTPBFI, ProlongateCoefficientFunction, TensorProductIntegrate
+import netgen.meshing as ngmeshing
+from netgen.geom2d import SplineGeometry
+import math
+
+
+
+def AddEdgeEls (p1, dx, facenr, n, pids,mesh):
+    for i in range(n):
+        base = p1 + i*dx
+        pnum = [base, base+dx]
+        elpids = [pids[p] for p in pnum]
+        mesh.Add (Element1D(elpids,index=1))
+
+def SegMesh(n,x0,x1,periodic=False):
+    mesh = Mesh(dim=1)
+    pids = []
+    for i in range(n+1):
+        pids.append (mesh.Add (MeshPoint(Pnt(x0+(x1-x0)*i/n, 0, 0))))
+    AddEdgeEls(0,1,1,n,pids,mesh)
+    mesh.Add (Element0D( pids[0], index=1))
+    mesh.Add (Element0D( pids[n], index=2))
+    mesh.SetBCName(0,"left")
+    mesh.SetBCName(1,"right")
+    if periodic == True:
+        mesh.AddPointIdentification(pids[0],pids[n],1,2)
+    return mesh
+
+# create a hexagonal mesh with corner points on the unit circle
+def MakeHexagonalMesh2D(maxh=0.1):
+    geo = SplineGeometry()
+    pnums = [ geo.AddPoint(math.cos(phi),math.sin(phi)) for phi in [xx*math.pi/3 for xx in range(6)] ]
+    l1 = geo.Append(["line", 0, 1], leftdomain=1, rightdomain=0, bc="upperRight")
+    geo.Append(["line", 4, 3], leftdomain=0, rightdomain=1, bc="lowerLeft", copy = l1)
+    l2 = geo.Append(["line", 1, 2], leftdomain=1, rightdomain=0, bc="upperCenter")
+    geo.Append(["line", 5, 4], leftdomain=0, rightdomain=1, bc="lowerCenter", copy = l2)
+    l3 = geo.Append(["line", 2, 3], leftdomain=1, rightdomain=0, bc="upperLeft")
+    geo.Append(["line", 0, 5], leftdomain=0, rightdomain=1, bc="lowerRight", copy = l3)
+    return geo.GenerateMesh(maxh=maxh)
+
+def MakeTensorProductMesh(mesh1,mesh2):
+    if mesh1.dim + mesh2.dim >3:
+        print('MakeMesh only possible if dim<=3!!')
+        return
+    if mesh1.dim + mesh2.dim == 3:
+        tpmesh = MakeMesh3D(mesh1,mesh2)
+    if mesh1.dim + mesh2.dim == 2:
+        tpmesh = MakeMesh2D(mesh1,mesh2)
+    AddSurfElements(tpmesh,mesh1,mesh2)
+    return tpmesh
+
+def AddSurfElements(tpmesh,mesh1,mesh2):
+    if mesh1.dim + mesh2.dim == 3:
+        tpmesh = AddSurfElements2D(tpmesh,mesh1,mesh2)
+    if mesh1.dim + mesh2.dim == 2:
+        tpmesh = AddSurfElements1D(tpmesh,mesh1,mesh2)
+    return tpmesh
+
+    
+def MakeMesh2D(mesh1,mesh2):
+    tpmesh = ngmeshing.Mesh(dim = mesh1.dim + mesh2.dim)
+    ngm1 = mesh1.ngmesh
+    ngm2 = mesh2.ngmesh
+
+    els1 = ngm1.Elements1D()
+    els2 = ngm2.Elements1D()
+    vert1 = ngm1.Points()
+    vert2 = ngm2.Points()
+    pids = []
+    for i in range(len(vert1)):
+        for j in range(len(vert2)):
+            pids.append(tpmesh.Add(MeshPoint(Pnt(vert1[ngmeshing.PointId(i+1)].p[0],vert2[ngmeshing.PointId(j+1)].p[0],0     ))))
+    tpmesh.Add (FaceDescriptor(surfnr=1,domin=1,bc=1))
+    for elx in els1:
+        for ely in els2:
+            pnum = [(elx.vertices[1].nr-1) * len(vert2) + ely.vertices[1].nr-1,
+                    (elx.vertices[0].nr-1) * len(vert2) + ely.vertices[1].nr-1,
+                    (elx.vertices[0].nr-1) * len(vert2) + ely.vertices[0].nr-1,
+                    (elx.vertices[1].nr-1) * len(vert2) + ely.vertices[0].nr-1]
+                    
+            elpids = [pids[p] for p in pnum]
+            tpmesh.Add(Element2D(1,elpids))
+    return tpmesh
+
+def MakeMesh3D(mesh1,mesh2):
+    tpmesh = ngmeshing.Mesh(dim = mesh1.dim + mesh2.dim)
+    ngm1 = mesh1.ngmesh
+    ngm2 = mesh2.ngmesh
+    if mesh1.dim==2:
+        els1 = ngm1.Elements2D()
+        els2 = ngm2.Elements1D()
+        vert1 = ngm1.Points()
+        vert2 = ngm2.Points()
+        pids = []
+        for i in range(len(vert1)):
+            for j in range(len(vert2)):
+                pids.append(tpmesh.Add(MeshPoint(Pnt(vert1[ngmeshing.PointId(i+1)].p[0],vert1[ngmeshing.PointId(i+1)].p[1],vert2[ngmeshing.PointId(j+1)].p[0]))))
+        for elx in els1:
+            for ely in els2:
+                pnum = []
+                for j in reversed(ely.vertices):
+                    for i in elx.vertices:
+                        pnum.append((i.nr-1) * len(vert2) + j.nr-1)
+                elpids = [pids[p] for p in pnum]
+                tpmesh.Add( Element3D(1,elpids) )
+        return tpmesh
+    else:
+        els1 = ngm1.Elements1D()
+        els2 = ngm2.Elements2D()
+        vert1 = ngm1.Points()
+        vert2 = ngm2.Points()
+        pids = []
+        for i in range(len(vert1)):
+            for j in range(len(vert2)):
+                pids.append(tpmesh.Add(MeshPoint(Pnt(vert1[ngmeshing.PointId(i+1)].p[0],vert2[ngmeshing.PointId(j+1)].p[0],vert2[ngmeshing.PointId(j+1)].p[1]))))
+        for elx in els1:
+            for ely in els2:
+                pnum = []
+                for i in reversed(elx.vertices):
+                    for j in (ely.vertices):
+                        pnum.append((i.nr-1) * len(vert2) + j.nr-1)
+                elpids = [pids[p] for p in pnum]
+                tpmesh.Add( Element3D(1,elpids) )
+        return tpmesh
+
+def AddSurfElements1D(tpmesh,mesh1,mesh2):
+    ngm1 = mesh1.ngmesh;
+    ngm2 = mesh2.ngmesh;
+    els1 = ngm1.Elements1D()
+    els2 = ngm2.Elements1D()
+    for ely in els2:
+        elpids = ely.vertices
+        elpids1=[]
+        for i in elpids:
+            elpids1.append(PointId((i.nr-1)+(len(ngm1.Points()) -1 )*(len(ngm2.Points())) + 1 ))
+        tpmesh.Add(Element1D(elpids))
+        tpmesh.Add(Element1D(elpids1))
+    for elx in els1:
+        elpids = elx.vertices
+        elpids1=[]
+        elpids2=[]
+        for i in elpids:
+            elpids1.append(PointId( (i.nr-1)*(len(ngm2.Points()))+1) )
+            elpids2.append(PointId( (i.nr-1)*(len(ngm2.Points()))+(len(ngm2.Points()))) )
+        tpmesh.Add(Element1D(elpids1))
+        tpmesh.Add(Element1D(elpids2))       
+
+def AddSurfElements2D(tpmesh,mesh1,mesh2):
+    ngm1 = mesh1.ngmesh;
+    ngm2 = mesh2.ngmesh;
+    if mesh1.dim==2:
+        els1 = ngm1.Elements2D()
+        els2 = ngm2.Elements1D()
+        tpmesh.Add (FaceDescriptor(surfnr=1,domin=1,bc=1))
+        for elx in els1:
+            vert_loc = elx.vertices
+            vert_glob = []
+            for vx in vert_loc:
+                vert_glob.append(PointId((vx.nr-1)*len(ngm2.Points())+len(ngm2.Points())))
+            tpmesh.Add(Element2D(1,vert_glob))
+        for elx in els1:
+            vert_loc = elx.vertices
+            vert_glob = []
+            for vx in vert_loc:
+                vert_glob = [PointId((vx.nr-1)*len(ngm2.Points())+1)] + vert_glob
+            tpmesh.Add(Element2D(1,vert_glob))
+        els1 = ngm1.Elements1D()
+        for elx in els1:
+            for ely in els2:
+                vert_glob=[]
+#            for vy in ely.vertices:
+#                for vx in elx.vertices:
+                vx = elx.vertices
+                vy = ely.vertices
+                vert_glob = [PointId((vx[1].nr-1)*len(ngm2.Points())+vy[0].nr),
+                            PointId((vx[1].nr-1)*len(ngm2.Points())+vy[1].nr),
+                            PointId((vx[0].nr-1)*len(ngm2.Points())+vy[1].nr),
+                            PointId((vx[0].nr-1)*len(ngm2.Points())+vy[0].nr)]
+                tpmesh.Add(Element2D(1,vert_glob))
+    else:
+        els1 = ngm1.Elements1D()
+        els2 = ngm2.Elements2D()
+        tpmesh.Add (FaceDescriptor(surfnr=1,domin=1,bc=1))
+        for ely in els2:
+            vert_loc = ely.vertices
+            vert_glob = []
+            for vy in vert_loc:
+                vert_glob.append(PointId((vy.nr)+(len(ngm1.Points())-1)*(len(ngm2.Points()))))
+            tpmesh.Add(Element2D(1,vert_glob))
+        for ely in els2:
+            vert_loc = ely.vertices
+            vert_glob = []
+            for vy in reversed(vert_loc):
+                vert_glob.append(PointId( vy.nr))
+            tpmesh.Add(Element2D(1,vert_glob))
+        els2 = ngm2.Elements1D()
+        for elx in els1:
+            for ely in els2:
+                vert_glob=[]
+                vx = elx.vertices
+                vy = ely.vertices
+                vert_glob = [PointId((vx[0].nr-1)*len(ngm2.Points())+vy[0].nr),
+                             PointId((vx[0].nr-1)*len(ngm2.Points())+vy[1].nr),
+                             PointId((vx[1].nr-1)*len(ngm2.Points())+vy[1].nr),
+                             PointId((vx[1].nr-1)*len(ngm2.Points())+vy[0].nr)]
+                tpmesh.Add(Element2D(1,vert_glob))
+    return tpmesh

ngsolve-6.2.2506.post216.dev0.data/data/lib/python3.14/site-packages/ngsolve/__console.py

@@ -0,0 +1,94 @@
+import code
+# import ngfemcf
+# import ngmpi
+
+def mydir(x=None):
+    if x==None:
+        return []
+    else:
+        return [i for i in dir(x) if not '__' in i]
+
+def execfile(fname):
+    exec(open(fname).read())
+
+class InteractiveMPIConsole(code.InteractiveConsole):
+# Function copied from /usr/lib/python3.4/code.py line 38 
+    def runsource(self, source, filename="<input>", symbol="single"):
+        try:
+            compiled_code = self.compile(source, filename, symbol)
+        except (OverflowError, SyntaxError, ValueError):
+            # Case 1
+            self.showsyntaxerror(filename)
+            return False
+
+        if compiled_code is None:
+            # Case 2
+            return True
+
+        # Case 3 -- first send code to other mpi processes
+        ngmpi.SendCommand('ngs_py '+source)
+        # then run it on master
+        code.InteractiveConsole.runcode(self, compiled_code)
+
+        # Avoid the prompt to show up before other processes' output
+        self.Barrier()
+        return False
+
+    def interact(self):
+        import sys
+        self.write("MPI Shell\n")
+        self.write("================\n")
+        self.write("Use pprint(str) to print with MPI ranks\n\n")
+        self.runsource("from ngmpi import *\n")
+        self.runsource("pprint=lambda p='':print('Process ' + str(ngmpi.Rank()) + '\\n'+str(p)+'\\n')\n")
+        self.locals.update(locals())
+        self.locals.update(globals())
+        sys.ps1 = "MPI >>> "
+        sys.ps2 = "MPI ... "
+        code.InteractiveConsole.interact(self,'')
+        sys.ps1 = ">>> "
+        sys.ps2 = "... "
+    def Barrier(self):
+        source = 'ngmpi.Barrier()'
+        ngmpi.SendCommand('ngs_py '+source)
+        code.InteractiveConsole.runsource(self, source)
+
+
+
+def MpiShell():
+    import code
+    try:
+        import readline
+        import rlcompleter
+        readline.parse_and_bind("tab:complete") # autocomplete
+    except:
+        try:
+            import pyreadline as readline
+            import rlcompleter
+            readline.parse_and_bind("tab:complete") # autocomplete
+        except:
+            print('readline not found')
+    vars2 = globals()
+    vars2.update(locals())
+    mshell = InteractiveMPIConsole(vars2)
+    mshell.interact()
+
+def startConsole(vars2):
+    try:
+        import readline
+        import rlcompleter
+        if 'libedit' in readline.__doc__:
+            readline.parse_and_bind("bind ^I rl_complete")
+        else:
+            readline.parse_and_bind("tab: complete")
+    except:
+        try:
+            import pyreadline as readline
+            import rlcompleter
+            readline.parse_and_bind("tab:complete") # autocomplete
+        except:
+            print('readline not found')
+    shell = code.InteractiveConsole(vars2)
+    shell.push('from netgen import *')
+    shell.push('from ngsolve import *')
+    shell.interact(banner="NGS-Python console is up ...")