added boundary points to UG data

program/source/grid/examples_ug_optics.cc

@@ -8,6 +8,7 @@
 #include "elements.h"
 #include "parti.h"
 #include "ug.h"
+#include <algorithm>
 #include "examples_ug_optics.h"
 
 ////////////////////////////////////////////////////////////////////
@@ -1580,6 +1581,8 @@ Lens_Geometry_Quad::Lens_Geometry_Quad(double Radius, double thickness, double c
     constructionParameters[10] = 0;
     constructionParameters[11] = 0;
 
+    boundaryPoints.push_back(D3vector{-Radius-1,-Radius-1,constructionParameters[8]-1});
+    boundaryPoints.push_back(D3vector{+Radius+1,+Radius+1,constructionParameters[9]+1});
 
 
 
@@ -1920,7 +1923,8 @@ Lens_Geometry_cutted_edges::Lens_Geometry_cutted_edges(double RadiusLeft, double
     pointer_global_data[11] = z_left_edge;
     pointer_global_data[12] = MechanicalRadiusRight;
     pointer_global_data[13] = MechanicalRadiusLeft;
-
+    boundaryPoints.push_back(D3vector{-MechanicalRadiusRight-1,-MechanicalRadiusRight-1,constructionParameters[8]-1});
+    boundaryPoints.push_back(D3vector{+MechanicalRadiusRight+1,+MechanicalRadiusRight+1,constructionParameters[9]+1});
 
 //    cout << "fixed lens param., channge again " << endl;
 //    radius = 1;

program/source/grid/ug.h

@@ -25,6 +25,7 @@
 #include <list>
 #include <vector>
 #include <algorithm>
+#include <math_lib.h>
 
 
 ////////////////////////////////////////////////////////////////
@@ -288,10 +289,12 @@ class Unstructured_grid : public Partitioning {
     Zordering* Give_Zordering() { return zordering; }
 
     double* Get_pointer_global_data() { return pointer_global_data; }
+    std::vector<D3vector> getBoundaryPoints(){return boundaryPoints;}
     int sizeGlobalPointerData() { return size_pointer_global_data; } 
   protected:
 	 // pointer for global objects used in coordinate transformations
      double* pointer_global_data;
+     std::vector<D3vector> boundaryPoints;
      int size_pointer_global_data;
      
      TypeOfUG typFuerSlice;