ngstrefftz 0.3.2402.dev34__tar.gz → 0.3.2403.dev5__tar.gz

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/.github/workflows/build_pip.sh

@@ -25,7 +25,7 @@ do
 
     #rm -rf /home/app/ngstrefftz/make
     rm -rf _skbuild
-    $PYDIR/pip install pytest-check numpy wheel scikit-build mkl==2023.* mkl-devel==2023.* setuptools setuptools_scm
+    $PYDIR/pip install pytest-check numpy wheel scikit-build mkl==2023.* mkl-devel==2023.* setuptools==69.5.1 setuptools_scm==8.1.0
     $PYDIR/pip install -r ./.github/workflows/ngsolve_version.txt
 
     $PYDIR/pip wheel -vvv .

ngstrefftz-0.3.2403.dev5/.github/workflows/ngsolve_version.txt

@@ -0,0 +1 @@
+ngsolve>=6.2.2402

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/Dockerfile

@@ -2,6 +2,11 @@ FROM python:3.10-slim-bullseye
 
 WORKDIR /home/app
 
+#RUN apt --yes update && apt --yes install git pandoc sphinx
+#COPY . /home/app/ngstrefftz
+#RUN cd /home/app/ngstrefftz/docs && pip install -r requirements.txt
+#RUN cd /home/app/ngstrefftz/docs && sphinx-build -M html . _build -vvv
+
 RUN pip install ngstrefftz 
 ENV PYTHONPATH=$PYTHONPATH:/usr/local/lib/python3.10/site-packages/
 ENV LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/usr/local/lib/

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/PKG-INFO

@@ -1,13 +1,13 @@
 Metadata-Version: 2.1
 Name: ngstrefftz
-Version: 0.3.2402.dev34
+Version: 0.3.2403.dev5
 Summary: NGSTrefftz is an add-on to NGSolve for Trefftz methods.
 Home-page: https://github.com/PaulSt/ngstrefftz
 Author: Paul Stocker
 Author-email: p.stocker@math.uni-goettingen.de
 License: LGPL2.1
 License-File: LICENSE
-Requires-Dist: ngsolve>=6.2.2401.post24.dev0
+Requires-Dist: ngsolve>=6.2.2402
 Requires-Dist: ngstents>=0.0.2.dev39
 Requires-Dist: mkl
 

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/docs/requirements.txt

@@ -1,11 +1,11 @@
-
-setuptools
+setuptools>=61
 setuptools_scm 
 jupyter>=1.0.0
 
 matplotlib>=3.6.1
 matplotlib-inline>=0.1.6
 
+pandoc>=2.3
 nbconvert>=7.2.1
 nbformat>=5.7.0
 nbsphinx>=0.8.9
@@ -20,6 +20,6 @@ myst_parser>=0.15.0
 numpy>=1.23.4
 scipy>=1.11.1 
 ngsolve>=6.2.2403
-xfem>=2.1.2303.post72.dev0
-ngstrefftz>=0.2.2401
+xfem==2.1.2303.post137.dev0
+ngstrefftz>=0.2.2402
 webgui-jupyter-widgets>=0.2.31

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/external_dependencies/ngstents/.github/workflows/build.yml

@@ -24,11 +24,11 @@ jobs:
              pip3 install pytest matplotlib numpy
       - name: set path
         run: |
-             echo "PYTHONPATH=/usr/lib/python3/dist-packages/:$GITHUB_WORKSPACE/install" >> $GITHUB_ENV
+             echo "PYTHONPATH=/usr/lib/python3/dist-packages/:$GITHUB_WORKSPACE/install/lib/python3.10/dist-packages" >> $GITHUB_ENV
              echo "NETGENDIR=/usr/bin/" >> $GITHUB_ENV
-      - uses: actions/checkout@v3
+      - uses: actions/checkout@v4
       - name: cmake ngstents
-        run: cmake -DCMAKE_INSTALL_PREFIX=$GITHUB_WORKSPACE/install -B$GITHUB_WORKSPACE/make -S$GITHUB_WORKSPACE/src
+        run: cmake -DCMAKE_INSTALL_PREFIX=$GITHUB_WORKSPACE/install -DPY_INSTALL_DIR="" -B$GITHUB_WORKSPACE/make -S$GITHUB_WORKSPACE/src
       - name: compile ngstents
         run: |
              make -C$GITHUB_WORKSPACE/make
@@ -39,27 +39,32 @@ jobs:
              pytest .
 
   macos-nightly:
-    runs-on: macos-latest
+    runs-on: macos-13
     steps:     
-      - name: install ngsolve
-        run: |
-             wget https://www.asc.tuwien.ac.at/~mhochsteger/ngsolve/macos/NGSolve_nightly.dmg
-             hdiutil attach NGSolve_nightly.dmg
-             cp -R /Volumes/NGSolve/* /Applications
+      - name: Get Python 
+        uses: actions/setup-python@v4
+        with:
+          python-version: '3.12'
       - name: set path
         run: |
-             echo "PYTHONPATH=$PYTHONPATH:/Applications/Netgen.app/Contents/Resources/lib/python3.8/site-packages:." >> $GITHUB_ENV
-             echo "NETGENDIR=/Applications/Netgen.app/Contents/MacOS" >> $GITHUB_ENV
-             echo "DYLD_LIBRARY_PATH=$DYLD_LIBRARY_PATH:$NETGENDIR" >> $GITHUB_ENV
-             echo "DYLD_FRAMEWORK_PATH=$DYLD_FRAMEWORK_PATH:$NETGENDIR/../Frameworks" >> $GITHUB_ENV
-             echo "$NETGENDIR" >> $GITHUB_PATH
-      - uses: actions/checkout@v3
+             echo "PYDIR=${{ env.Python3_ROOT_DIR }}/bin" >> $GITHUB_ENV
+             echo "NETGEN_Dir=${{ env.Python3_ROOT_DIR }}/bin/../lib/python3.12/site-packages/netgen/cmake" >> $GITHUB_ENV
+             echo "NGSolve_Dir=${{ env.Python3_ROOT_DIR }}/bin/../lib/python3.12/site-packages/ngsolve/cmake" >> $GITHUB_ENV
+             echo "/Applications/CMake.app/Contents/bin" >> $GITHUB_PATH
+      - run: echo "CMAKE_PREFIX_PATH=${{ env.CMAKE_PREFIX_PATH }}:${{ env.NGSolve_Dir }}:${{ env.NETGEN_Dir }}" >> $GITHUB_ENV
+      - name: install ngsolve
+        run: pip3 install ngsolve pytest matplotlib
+      - uses: actions/checkout@v4
       - name: cmake ngstents
         run: cmake -B$GITHUB_WORKSPACE/make -S$GITHUB_WORKSPACE/src
       - name: compile ngstents
         run: |
              make -C$GITHUB_WORKSPACE/make
              sudo make -C$GITHUB_WORKSPACE/make install  
+      - name: test ngstents
+        run: |
+             cd $GITHUB_WORKSPACE/tests
+             pytest .
 
   pypi-linux:
     name: Build wheels for linux
@@ -68,7 +73,7 @@ jobs:
     steps:
       - name: clear cache
         run: rm -rf /opt/hostedtoolcache
-      - uses: actions/checkout@v3
+      - uses: actions/checkout@v4
         with:
           submodules: 'True'
           path: 'ngstents'
@@ -78,9 +83,9 @@ jobs:
       - name: build pip
         run: docker run -v $GITHUB_WORKSPACE:/workspace sameli/manylinux2014_x86_64_cuda_11.7 /bin/sh /workspace/ngstents/.github/workflows/build_pip.sh 
       - name: Store the binary wheel
-        uses: actions/upload-artifact@v3
+        uses: actions/upload-artifact@v4
         with:
-          name: python-package-distributions
+          name: python-package-distributions-linux
           path: ngstents/wheelhouse
 
   pypi-mac:
@@ -88,13 +93,13 @@ jobs:
     needs: [ubuntu-nightly, macos-nightly]
     strategy:
       matrix:
-        py: ['3.8', '3.9', '3.10', '3.11']
+        py: ['3.8', '3.9', '3.10', '3.11', '3.12']
     runs-on: macos-11
     env:
       MACOSX_DEPLOYMENT_TARGET: '10.15'
     steps:
       - name: Checkout
-        uses: actions/checkout@v3
+        uses: actions/checkout@v4
         with:
           fetch-depth: 0
           submodules: recursive
@@ -106,9 +111,9 @@ jobs:
       - name: Run build script
         run: cd ngstents/.github/workflows && bash build_pip_mac.sh ${{ matrix.py }}
       - name: Store the binary wheel
-        uses: actions/upload-artifact@v3
+        uses: actions/upload-artifact@v4
         with:
-          name: python-package-distributions
+          name: python-package-distributions-mac-${{ matrix.py }}
           path: ngstents/wheelhouse
 
   pypi-win:
@@ -116,11 +121,11 @@ jobs:
     needs: [ubuntu-nightly, macos-nightly]
     strategy:
       matrix:
-        py: ['3.8', '3.9', '3.10', '3.11']
+        py: ['3.8', '3.9', '3.10', '3.11', '3.12']
     runs-on: windows-2019
     steps:
       - name: Checkout
-        uses: actions/checkout@v3
+        uses: actions/checkout@v4
         with:
           fetch-depth: 0
           submodules: recursive
@@ -133,9 +138,9 @@ jobs:
       - name: Run build script
         run: Set-Location ngstents/.github/workflows; ./build_pip.ps1 ${{ env.pythonLocation }}
       - name: Store the binary wheel
-        uses: actions/upload-artifact@v3
+        uses: actions/upload-artifact@v4
         with:
-          name: python-package-distributions
+          name: python-package-distributions-win-${{ matrix.py }}
           path: ngstents/wheelhouse
 
   publish-pypi:
@@ -144,9 +149,10 @@ jobs:
     runs-on: ubuntu-latest
     steps:
     - name: Download all the dists
-      uses: actions/download-artifact@v2
+      uses: actions/download-artifact@v4
       with:
-        name: python-package-distributions
+        pattern: python-package-distributions-*
+        merge-multiple: true
         path: ngstents/wheelhouse/
     - name: Publish binary distributions to PyPI
       if: github.ref == 'refs/heads/master'
@@ -161,7 +167,7 @@ jobs:
     needs: [publish-pypi]
     runs-on: ubuntu-latest
     steps:
-      - uses: actions/checkout@v3
+      - uses: actions/checkout@v4
         with:
           fetch-depth: 0 # otherwise, you will failed to push refs to dest repo
           submodules: 'True'

ngstrefftz-0.3.2403.dev5/external_dependencies/ngstents/.github/workflows/build_pip.ps1

@@ -0,0 +1,24 @@
+$ErrorActionPreference = "Stop"
+
+if (test-path _skbuild) {
+    cmd.exe /c rd /s /q _skbuild
+}
+if (test-path dist) {
+    cmd.exe /c rd /s /q dist
+}
+if (test-path venv_ngs) {
+    cmd.exe /c rd /s /q venv_ngs
+}
+
+Set-Location ../..
+python.exe -m venv .\venv_ngs
+.\venv_ngs\scripts\Activate.ps1
+$env:PATH += Join-Path ";" (Get-Item .).FullName "venv_ngs\bin"
+
+pip3 install scikit-build wheel numpy twine mkl-devel==2022.* mkl==2022.* setuptools==69.5.1 setuptools_scm==8.1.0
+pip3 install -r .github/workflows/ngsolve_version.txt
+
+#$env:NGSolve_DIR = "$env:Python3_ROOT_DIR\lib\site-packages\ngsolve\cmake"
+#$env:Netgen_DIR = "$env:Python3_ROOT_DIR\lib\site-packages\netgen\cmake"
+
+python setup.py bdist_wheel -G"Visual Studio 16 2019" -d wheelhouse

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/external_dependencies/ngstents/.github/workflows/build_pip.sh

@@ -15,18 +15,18 @@ git config --global --add safe.directory '*'
 
 export ORIGINAL_PATH="$PATH"
 
-for pyversion in 38 39 310 311
+for pyversion in 38 39 310 311 312
 do
     export PYDIR="/opt/python/cp${pyversion}-cp${pyversion}/bin"
     export PATH="$ORIGINAL_PATH:$PYDIR"
 
     rm -rf _skbuild
-    $PYDIR/pip install pytest-check numpy wheel scikit-build mkl==2023.* mkl-devel==2023.* setuptools
-    $PYDIR/pip install ngsolve 
+    $PYDIR/pip install pytest-check numpy wheel scikit-build mkl==2023.* mkl-devel==2023.* setuptools==69.5.1 setuptools_scm==8.1.0
+    $PYDIR/pip install -r ./.github/workflows/ngsolve_version.txt
 
     $PYDIR/pip wheel -vvv .
     rename linux_ manylinux_2_17_x86_64.manylinux2014_ ngstents*.whl
     mv ngstents*.whl wheelhouse/
     rm -rf *.whl
-    $PYDIR/pip uninstall -y ngsolve netgen-mesher setuptools pytest-check numpy wheel scikit-build mkl mkl-devel 
+    $PYDIR/pip uninstall -y ngsolve netgen-mesher setuptools setuptools_scm pytest-check numpy wheel scikit-build mkl mkl-devel 
 done

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/external_dependencies/ngstents/.github/workflows/build_pip_mac.sh

@@ -11,9 +11,9 @@ export PATH=/Applications/CMake.app/Contents/bin:$PATH
 export NETGEN_Dir=$PYDIR/../lib/python$1/site-packages/netgen/cmake
 export NGSolve_Dir=$PYDIR/../lib/python$1/site-packages/ngsolve/cmake
 export CMAKE_PREFIX_PATH=$CMAKE_PREFIX_PATH:$NGSolve_Dir:$NETGEN_Dir
-$PYDIR/pip3 install scikit-build wheel
+$PYDIR/pip3 install scikit-build wheel setuptools==69.5.1 setuptools_scm==8.1.0
 
 export CMAKE_OSX_ARCHITECTURES='arm64;x86_64'
-$PYDIR/pip3 install ngsolve
+$PYDIR/pip3 install -r ./.github/workflows/ngsolve_version.txt
 
 $PYDIR/python3 setup.py bdist_wheel --plat-name macosx-10.15-universal2 -d wheelhouse

ngstrefftz-0.3.2403.dev5/external_dependencies/ngstents/.github/workflows/ngsolve_version.txt

@@ -0,0 +1 @@
+ngsolve>=6.2.2402

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/external_dependencies/ngstents/README.md

@@ -16,6 +16,8 @@ python >= 3.9).
 
 * `pip install ngstents`
 
+But note that the binaries are not generally kept up to date. The
+recommended option is to build from source (see below).
 
 ## Build
 

ngstrefftz-0.3.2403.dev5/external_dependencies/ngstents/setup.py

@@ -0,0 +1,27 @@
+#!/usr/bin/env python3
+# -*- coding: utf-8 -*-
+from skbuild import setup
+from setuptools_scm import get_version
+import sys
+
+_cmake_args = []
+version = get_version().split('+')[0]
+
+ngsolve_version = '6.2.2402'
+install_requires = [ 'ngsolve >= '+ngsolve_version, 'matplotlib' ] 
+
+setup(
+    name='ngstents',
+    version=version,
+    author='Jay Gopalakrishnan',
+    author_email='gjay@pdx.edu',
+    description='Spacetime tent facilities for solving hyperbolic equations',
+    long_description='This ngstents package is a c++ extension of the NGSolve finite element library, designed to ease experimentation with solvers based on spacetime tents for hyperbolic systems. A python front-end allows new equations (linear or nonlinear conservation laws) to be solved by easily defining required fluxes and numerical fluxes in a few lines of code.',
+    url="https://github.com/jayggg/ngstents",
+    license="LGPL2.1",
+    install_requires=install_requires,
+    packages=["ngstents"],
+    package_dir={"ngstents": "src"},
+    cmake_args=_cmake_args,
+    cmake_source_dir='src',
+)

ngstrefftz-0.3.2403.dev5/external_dependencies/ngstents/src/CMakeLists.txt

@@ -0,0 +1,97 @@
+cmake_minimum_required(VERSION 3.1)
+project(ngstents)
+set(CMAKE_CXX_STANDARD 17)
+set(CMAKE_CXX_STANDARD_REQUIRED ON)
+
+# Find NGSolve and Netgen using python
+if(CMAKE_VERSION VERSION_LESS "3.18")
+  find_package(Python3 REQUIRED COMPONENTS Interpreter Development)
+else()
+  find_package(Python3 REQUIRED COMPONENTS Interpreter Development.Module)
+endif()
+set(Netgen_DIR "" CACHE PATH "Path to directory containing NetgenConfig.cmake")
+set(NGSolve_DIR "" CACHE PATH "Path to directory containing NGSolveConfig.cmake")
+execute_process(COMMAND ${Python3_EXECUTABLE} -m netgen.config OUTPUT_STRIP_TRAILING_WHITESPACE OUTPUT_VARIABLE Netgen_DIR)
+execute_process(COMMAND ${Python3_EXECUTABLE} -m ngsolve.config OUTPUT_STRIP_TRAILING_WHITESPACE OUTPUT_VARIABLE NGSolve_DIR)
+find_package(NGSolve CONFIG REQUIRED)
+
+### check if CMAKE_INSTALL_PREFIX is set by user, if not install in NGSolve python dir
+if(CMAKE_INSTALL_PREFIX_INITIALIZED_TO_DEFAULT)
+    execute_process(COMMAND ${Python3_EXECUTABLE} -m site --user-site OUTPUT_STRIP_TRAILING_WHITESPACE OUTPUT_VARIABLE install_dir RESULT_VARIABLE ret)
+    if (NOT ret EQUAL 0)
+        # user site directory is disabled (are we in a virtual environment?)
+        set(install_dir ${Python3_SITEARCH})
+    endif()
+    message("The python module ${module_name} will be installed to: ${install_dir}")
+    set(CMAKE_INSTALL_PREFIX ${install_dir} CACHE PATH "Install dir" FORCE)
+else(CMAKE_INSTALL_PREFIX_INITIALIZED_TO_DEFAULT)
+    #set(PY_INSTALL_DIR ${NGSOLVE_INSTALL_DIR_PYTHON} CACHE STRING "subdir of install dir for python modules")
+    execute_process(COMMAND ${Python3_EXECUTABLE} -c "import os.path, sysconfig;print(os.path.relpath(sysconfig.get_path('platlib'), sysconfig.get_path('data')))" OUTPUT_VARIABLE PY_INSTALL_DIR OUTPUT_STRIP_TRAILING_WHITESPACE)
+endif(CMAKE_INSTALL_PREFIX_INITIALIZED_TO_DEFAULT)
+set(CMAKE_LIBRARY_OUTPUT_DIRECTORY ${CMAKE_BINARY_DIR} CACHE PATH "build dir")
+
+add_library(_pytents SHARED
+  python_tents.cpp
+  tents.cpp
+  )
+add_library(_pyconslaw SHARED
+  python_conslaw.cpp
+  burgers.cpp
+  euler.cpp
+  wave.cpp
+  advection.cpp
+  maxwell.cpp
+  symbolic.cpp
+  vis3d.cpp
+  )
+target_link_libraries(_pytents PUBLIC ngsolve Python3::Module)
+set_target_properties(_pytents PROPERTIES PREFIX "" CXX_STANDARD 17)
+target_compile_definitions(_pytents PRIVATE NGSTENT_EXPORTS)
+target_link_libraries(_pyconslaw PRIVATE _pytents)
+set_target_properties(_pyconslaw PROPERTIES PREFIX "" CXX_STANDARD 17)
+if(WIN32)
+  set_target_properties(_pytents PROPERTIES SUFFIX ".pyd" )
+  set_target_properties(_pyconslaw PROPERTIES SUFFIX ".pyd" )
+else(WIN32)
+  set_target_properties(_pytents PROPERTIES SUFFIX ".so")
+  set_target_properties(_pyconslaw PROPERTIES SUFFIX ".so")
+endif(WIN32)
+set_target_properties(_pytents PROPERTIES LIBRARY_OUTPUT_DIRECTORY "${CMAKE_LIBRARY_OUTPUT_DIRECTORY}/ngstents")
+set_target_properties(_pyconslaw PROPERTIES LIBRARY_OUTPUT_DIRECTORY "${CMAKE_LIBRARY_OUTPUT_DIRECTORY}/ngstents")
+
+option(BUILD_STUB_FILES "Generate stub files for better autocompletion support" OFF)
+if(BUILD_STUB_FILES)
+  execute_process(COMMAND ${PYTHON_EXECUTABLE} -c "import pybind11_stubgen; print(pybind11_stubgen.__file__)" OUTPUT_VARIABLE stubgen_path RESULT_VARIABLE pybind11_stubgen)
+  if(pybind11_stubgen AND NOT ${pybind11_stubgen} EQUAL 0)    
+    message(WARNING "pybind11-stubgeN not found, if you want to create stub files
+for better autocompletion support install it with pip.")
+  else()
+    message("-- Found pybind11-stubgen: ${stubgen_path}")
+    set(FOUND_STUBGEN ON)
+  endif()
+   message(STATUS
+      "This project is only compatible with\n\t\tpybind11-stubgen==0.5")
+endif(BUILD_STUB_FILES)
+
+message("With 'make install' the python package will be installed to: ${CMAKE_INSTALL_PREFIX}/${PY_INSTALL_DIR}")
+install(TARGETS _pytents DESTINATION ${CMAKE_INSTALL_PREFIX}/${PY_INSTALL_DIR}/ngstents)
+install(TARGETS _pyconslaw DESTINATION ${CMAKE_INSTALL_PREFIX}/${PY_INSTALL_DIR}/ngstents/conslaw)
+
+install(FILES
+  ../py/__init__.py
+  DESTINATION ${CMAKE_INSTALL_PREFIX}/${PY_INSTALL_DIR}/ngstents)
+
+install(FILES
+  ../py/utils/__init__.py
+  ../py/utils/_drawtents.py
+  ../py/utils/_drawtents2d.py
+  DESTINATION ${CMAKE_INSTALL_PREFIX}/${PY_INSTALL_DIR}/ngstents/utils)
+
+install(FILES
+  ../py/conslaw/__init__.py
+  DESTINATION ${CMAKE_INSTALL_PREFIX}/${PY_INSTALL_DIR}/ngstents/conslaw)
+
+if(FOUND_STUBGEN)
+  install(CODE "execute_process(COMMAND ${PYTHON_EXECUTABLE} -m pybind11_stubgen --no-setup-py  ngstents)")
+  install(DIRECTORY ${CMAKE_CURRENT_BINARY_DIR}/stubs/ngstents-stubs/ DESTINATION ${CMAKE_INSTALL_PREFIX}/${PY_INSTALL_DIR}/ngstents)
+endif(FOUND_STUBGEN)

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/external_dependencies/ngstents/src/tents.cpp

@@ -565,7 +565,7 @@ std::tuple<Table<int>,Table<int>> TentSlabPitcher::InitializeMeshData(LocalHeap
 // Get the periodic vertex associated with a primary vertex in periodic case
 void TentSlabPitcher::GetVertexElements(int vnr_main, Array<int> & elems) const
 {
-  ma->GetVertexElements (vnr_main, elems);
+  elems = ma->GetVertexElements(vnr_main); 
   if(per_verts[vnr_main].Size()==0)
     return;
   else

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/external_dependencies/ngstents/src/tents.hpp

@@ -34,6 +34,7 @@ using namespace std;
 /// vertices.
 ///
 
+
 class NGSTENT_API Tent {
 
 public:

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/src/embtrefftz.cpp

@@ -147,6 +147,7 @@ namespace ngcomp
 
     vector<shared_ptr<Matrix<SCAL>>> ETmats (ne);
     auto lfvec = make_shared<VVector<SCAL>> (ndof);
+    lfvec->operator= (0.0);
 
     size_t active_elements = 0;
     size_t test_local_ndof = test_fes->GetNDof () / ne;
@@ -249,14 +250,15 @@ namespace ngcomp
         nz = tndof;
       else
         {
-          nz = dofs.Size () - test_dofs.Size ();
-          nz = max (nz, 0);
+          nz = max<size_t> (dofs.Size () - test_dofs.Size (), 0);
           for (int i = 0; i < min (elmat.Width (), elmat.Height ()); i++)
             if (abs (elmat (i, i)) < eps)
               nz++;
         }
-      if (dofs.Size () - test_dofs.Size () > nz)
-        throw Exception ("test fes not large enough for given tndof");
+      if (dofs.Size () < nz)
+        throw Exception ("tndof too large, nz: " + to_string (nz)
+                         + ", dofs:" + to_string (dofs.Size ())
+                         + ", test_dofs:" + to_string (test_dofs.Size ()));
 
       if (getrange)
         ETmats[ei.Nr ()]

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/src/monomialfespace.cpp

@@ -98,9 +98,15 @@ namespace ngcomp
           case ET_QUAD:
           case ET_TRIG:
             {
+              Vec<2> scale = 1.0;
+              if (usescale == 2)
+                scale = { 1.0 / ElSize<2> (ei, { 1.0, 0 }),
+                          1.0 / ElSize<2> (ei, { 0, 1.0 }) };
+              else if (usescale != 0)
+                scale = 1.0 / ElSize<2> (ei);
               return *(new (alloc) ScalarMappedElement<2> (
                   local_ndof, order, basismat, eltype, ElCenter<2> (ei),
-                  1.0 / Adiam<2> (ei)));
+                  scale));
               break;
             }
           case ET_HEX:
@@ -108,9 +114,18 @@ namespace ngcomp
           case ET_PYRAMID:
           case ET_TET:
             {
+              Vec<3> scale = 1.0;
+              if (usescale == 2)
+                {
+                  double hx = ElSize<3> (ei, { 1.0, 1.0, 0 });
+                  double ht = ElSize<3> (ei, { 0, 0, 1.0 });
+                  scale = { 1.0 / hx, 1.0 / hx, 1.0 / ht };
+                }
+              else if (usescale != 0)
+                scale = 1.0 / ElSize<3> (ei);
               return *(new (alloc) ScalarMappedElement<3> (
                   local_ndof, order, basismat, eltype, ElCenter<3> (ei),
-                  1.0 / Adiam<3> (ei)));
+                  scale));
               break;
             }
           }

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/src/monomialfespace.hpp

@@ -51,47 +51,29 @@ namespace ngcomp
       MatToCSR (basis, tb);
       return tb;
     }
-
-    template <int D> double Adiam (ElementId ei) const
+    template <int D>
+    double ElSize (ElementId ei, Vec<D> coeff_const = 1.0) const
     {
-      if (usescale == 0)
-        return 1.0;
-      LocalHeap lh (1000 * 1000);
       double anisotropicdiam = 0.0;
       auto vertices_index = ma->GetElVertices (ei);
-
       for (auto vertex1 : vertices_index)
         {
           for (auto vertex2 : vertices_index)
             {
-              Vec<D> v1 = ma->GetPoint<D> (vertex1);
-              Vec<D> v2 = ma->GetPoint<D> (vertex2);
-              IntegrationRule ir (ma->GetElType (ei), 0);
-              ElementTransformation &trafo = ma->GetTrafo (ei, lh);
-              MappedIntegrationPoint<D, D> mip (ir[0], trafo);
-              mip.Point () = v1;
-              double c1 = coeff_cf ? coeff_cf->Evaluate (mip) : 1.0;
-              mip.Point () = v2;
-              double c2 = coeff_cf ? coeff_cf->Evaluate (mip) : 1.0;
-
-              anisotropicdiam = max (
-                  anisotropicdiam,
-                  sqrt (L2Norm2 (v1.Range (0, D - 1) - v2.Range (0, D - 1))
-                        + pow (c1 * v1 (D - 1) - c2 * v2 (D - 1), 2)));
+              Vec<D> v = ma->GetPoint<D> (vertex2) - ma->GetPoint<D> (vertex1);
+              vtimes (v, coeff_const);
+              anisotropicdiam = max (anisotropicdiam, sqrt (L2Norm2 (v)));
             }
         }
       return anisotropicdiam;
     }
-
     template <int D> Vec<D> ElCenter (ElementId ei) const
     {
-      if (useshift == 0)
-        return 1;
       Vec<D> center = 0;
       auto vertices_index = ma->GetElVertices (ei);
       for (auto vertex : vertices_index)
         center += ma->GetPoint<D> (vertex);
-      center *= (1.0 / vertices_index.Size ());
+      center *= (1.0 / vertices_index.Size ()) * useshift;
       return center;
     }
   };

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/src/ngstrefftz.egg-info/PKG-INFO

@@ -1,13 +1,13 @@
 Metadata-Version: 2.1
 Name: ngstrefftz
-Version: 0.3.2402.dev34
+Version: 0.3.2403.dev5
 Summary: NGSTrefftz is an add-on to NGSolve for Trefftz methods.
 Home-page: https://github.com/PaulSt/ngstrefftz
 Author: Paul Stocker
 Author-email: p.stocker@math.uni-goettingen.de
 License: LGPL2.1
 License-File: LICENSE
-Requires-Dist: ngsolve>=6.2.2401.post24.dev0
+Requires-Dist: ngsolve>=6.2.2402
 Requires-Dist: ngstents>=0.0.2.dev39
 Requires-Dist: mkl
 

{ngstrefftz-0.3.2402.dev34 → ngstrefftz-0.3.2403.dev5}/src/ngstrefftz.egg-info/SOURCES.txt

@@ -48,6 +48,7 @@ external_dependencies/ngstents/.github/workflows/build_pip.ps1
 external_dependencies/ngstents/.github/workflows/build_pip.sh
 external_dependencies/ngstents/.github/workflows/build_pip_mac.sh
 external_dependencies/ngstents/.github/workflows/fix_auditwheel_policy.py
+external_dependencies/ngstents/.github/workflows/ngsolve_version.txt
 external_dependencies/ngstents/demo/advection/Advection_Periodic_Clipping.ipynb
 external_dependencies/ngstents/demo/advection/advection2d.py
 external_dependencies/ngstents/demo/advection/advection2d_periodic.py

ngstrefftz-0.3.2403.dev5/src/ngstrefftz.egg-info/requires.txt

@@ -0,0 +1,3 @@
+ngsolve>=6.2.2402
+ngstents>=0.0.2.dev39
+mkl