From d33d38f1130acd4abd314a11ee5315ebca482455 Mon Sep 17 00:00:00 2001
From: Matthias Hochsteger <mhochsteger@cerbsim.com>
Date: Wed, 7 Sep 2022 13:50:05 +0200
Subject: [PATCH] faster 2d smoothing for certain mixed meshes

don't consider quads with fixed points as "mixed mesh" (will be skipped during
optimization anyway)
---
 libsrc/meshing/smoothing2.cpp | 15 +++++++++++----
 1 file changed, 11 insertions(+), 4 deletions(-)

diff --git a/libsrc/meshing/smoothing2.cpp b/libsrc/meshing/smoothing2.cpp
index 11f21cb8..0e6db78a 100644
--- a/libsrc/meshing/smoothing2.cpp
+++ b/libsrc/meshing/smoothing2.cpp
@@ -756,10 +756,17 @@ namespace netgen
       {
         for (auto & se : mesh.SurfaceElements())
           if (se.GetNP() != 3)
-            {
-              mixed = true;
-              break;
-            }
+          {
+              for(auto pi : se.PNums())
+                  if(mesh[pi].Type() == SURFACEPOINT)
+                  {
+                      mixed = true;
+                      break;
+                  }
+              if(mixed)
+                  break;
+          }
+
         const auto & getDofs = [&] (int i)
           {
             return elementsonpoint[i+PointIndex::BASE];