From 1337786e7370f48907db4db134cd4865006187ab Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Joachim=20Sch=C3=B6berl?= <joachim.schoeberl@tuwien.ac.at>
Date: Mon, 9 May 2016 05:45:09 +0200
Subject: [PATCH] fixes for surface element point search

---
 libsrc/interface/nginterface_v2.cpp | 15 ++++++++++++---
 1 file changed, 12 insertions(+), 3 deletions(-)

diff --git a/libsrc/interface/nginterface_v2.cpp b/libsrc/interface/nginterface_v2.cpp
index 6b918e54..0e082a3a 100644
--- a/libsrc/interface/nginterface_v2.cpp
+++ b/libsrc/interface/nginterface_v2.cpp
@@ -739,9 +739,18 @@ namespace netgen
     for (int i = 0; i < numind; i++) dummy[i] = indices[i]+1;
     
     double lam3[3];
-    Point<3> p2d(p[0], p[1], 0);
-    int ind = 
-      mesh->GetElementOfPoint(p2d, lam3, &dummy, build_searchtree);
+    int ind;
+
+    if (mesh->GetDimension() == 2)
+      {
+        Point<3> p2d(p[0], p[1], 0);
+        ind = mesh->GetElementOfPoint(p2d, lam3, &dummy, build_searchtree);
+      }
+    else
+      {
+        Point3d p3d(p[0], p[1], p[2]);
+        ind = mesh->GetSurfaceElementOfPoint(p3d, lam3, &dummy, build_searchtree);
+      }
     
     if (ind > 0)
       {