ngstrefftz 0.3.2507.dev7__tar.gz → 0.3.2507.dev8__tar.gz

{ngstrefftz-0.3.2507.dev7 → ngstrefftz-0.3.2507.dev8}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: ngstrefftz
-Version: 0.3.2507.dev7
+Version: 0.3.2507.dev8
 Summary: NGSTrefftz is an add-on to NGSolve for Trefftz methods.
 Author-Email: Paul Stocker <paul.stocker@univie.ac.at>
 Classifier: License :: OSI Approved :: GNU Lesser General Public License v3 (LGPLv3)

{ngstrefftz-0.3.2507.dev7 → ngstrefftz-0.3.2507.dev8}/_version.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '0.3.2507.dev7'
-__version_tuple__ = version_tuple = (0, 3, 2507, 'dev7')
+__version__ = version = '0.3.2507.dev8'
+__version_tuple__ = version_tuple = (0, 3, 2507, 'dev8')
 
-__commit_id__ = commit_id = 'gc1672b4a2'
+__commit_id__ = commit_id = 'g25f082f1b'

{ngstrefftz-0.3.2507.dev7 → ngstrefftz-0.3.2507.dev8}/src/embtrefftz.cpp

@@ -761,6 +761,14 @@ namespace ngcomp
           size_t ndof_test = dofs_test.Size ();
           size_t ndof_conforming = dofs_conforming.Size ();
 
+          if (ndof_test == 0 && ndof_conforming == 0)
+            {
+              Identity Id (ndof);
+              etmats[element_id.Nr ()] = make_optional<Matrix<SCAL>> (Id);
+              local_ndofs_trefftz[element_id.Nr ()] = ndof;
+              return;
+            }
+
           FlatMatrix<SCAL> elmat_A (ndof_test + ndof_conforming, ndof, lh);
           FlatMatrix<SCAL> elmat_B (ndof_test + ndof_conforming,
                                     ndof_conforming, lh);

{ngstrefftz-0.3.2507.dev7 → ngstrefftz-0.3.2507.dev8}/src/tp0fespace.cpp

@@ -200,13 +200,12 @@ namespace ngcomp
 
   FiniteElement &TP0FESpace ::GetFE (ElementId ei, Allocator &alloc) const
   {
-
     Ngs_Element ngel = ma->GetElement (ei);
     ELEMENT_TYPE eltype = ngel.GetType ();
     int D = ma->GetDimension ();
     int order = order_inner.Size () > 0 ? order_inner[ei.Nr ()] : this->order;
 
-    if (ei.IsVolume ())
+    if (ei.IsVolume () && order > 1)
       {
         if (D == 2)
           switch (eltype)

{ngstrefftz-0.3.2507.dev7 → ngstrefftz-0.3.2507.dev8}/src/tp0fespace.hpp

@@ -38,6 +38,8 @@ namespace ngcomp
     {
       HeapReset hr (lh);
       // Cast (fel).CalcShape (mip.IP (), mat.Row (0));
+      if (fel.Order () < 2)
+        return;
       switch (fel.ElementType ())
         {
         case ET_QUAD: