From 17b5f1c7a47e07ce149620bd65b8c2f2debdfc0f Mon Sep 17 00:00:00 2001
From: Matthias Hochsteger <mhochsteger@cerbsim.com>
Date: Tue, 22 Jun 2021 11:35:44 +0200
Subject: [PATCH] preserve locked points in all domains

---
 libsrc/meshing/meshfunc.cpp | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/libsrc/meshing/meshfunc.cpp b/libsrc/meshing/meshfunc.cpp
index 9c0d1135..fe75af64 100644
--- a/libsrc/meshing/meshfunc.cpp
+++ b/libsrc/meshing/meshfunc.cpp
@@ -97,6 +97,11 @@ namespace netgen
         }
       }
 
+      // mark locked/fixed points for each domain TODO: domain bounding box to add only relevant points?
+      for(auto pi : mesh.LockedPoints())
+        for(auto i : Range(ret))
+          ipmap[i][pi] = 1;
+
       // add used points to domain mesh, build point mapping
       for(auto i : Range(ret))
       {