Implemented canGetClosestPoint and isOnModel

Author: mortezah

apf/apfMesh.cc

@@ -175,6 +175,11 @@ bool Mesh::canSnap()
   return gmi_can_eval(getModel());
 }
 
+bool Mesh::canGetClosestPoint()
+{
+  return gmi_can_get_closest_point(getModel());
+}
+
 bool Mesh::canGetModelNormal()
 {
   return gmi_has_normal(getModel());
@@ -252,6 +257,16 @@ bool Mesh::isInClosureOf(ModelEntity* g, ModelEntity* target){
   return (res == 1) ? true : false;
 }
 
+bool Mesh::isOnModel(ModelEntity* g, Vector3 p, double scale)
+{
+  Vector3 to;
+  double param[2];
+  gmi_ent* c = (gmi_ent*)g;
+  gmi_closest_point(getModel(), c, &p[0], &to[0], param);
+  double ratio = (to - p).getLength() / scale;
+  return ratio < 0.001;
+}
+
 void Mesh::getPoint(MeshEntity* e, int node, Vector3& p)
 {
   getVector(coordinateField,e,node,p);

apf/apfMesh.h

@@ -285,6 +285,8 @@ class Mesh
     ModelEntity* findModelEntity(int type, int tag);
     /** \brief return true if the geometric model supports snapping */
     bool canSnap();
+    /** \brief return true if the geometric model supports get closest point */
+    bool canGetClosestPoint();
     /** \brief return true if the geometric model supports normal computation */
     bool canGetModelNormal();
     /** \brief evaluate parametric coordinate (p) as a spatial point (x) */
@@ -310,6 +312,8 @@ class Mesh
       Vector3 const& param, Vector3& x);
     /** \brief checks if g is in the closure of the target */
     bool isInClosureOf(ModelEntity* g, ModelEntity* target);
+    /** \brief checks if p is on model g */
+    bool isOnModel(ModelEntity* g, Vector3 p, double scale);
     /** \brief get the distribution of the mesh's coordinate field */
     FieldShape* getShape() const;
     /** \brief get the mesh's coordinate field */

gmi/gmi.c

@@ -85,6 +85,11 @@ int gmi_can_eval(struct gmi_model* m)
   return m->ops->eval != NULL;
 }
 
+int gmi_can_get_closest_point(struct gmi_model* m)
+{
+  return m->ops->closest_point != NULL;
+}
+
 int gmi_has_normal(struct gmi_model* m)
 {
   return m->ops->normal != NULL;

gmi/gmi.h

@@ -145,6 +145,8 @@ struct gmi_ent* gmi_find(struct gmi_model* m, int dim, int tag);
 struct gmi_set* gmi_adjacent(struct gmi_model* m, struct gmi_ent* e, int dim);
 /** \brief check whether the model implements gmi_eval */
 int gmi_can_eval(struct gmi_model* m);
+/** \brief check whether the model gmi_closest_point */
+int gmi_can_get_closest_point(struct gmi_model* m);
 /** \brief check whether the model implements gmi_normal */
 int gmi_has_normal(struct gmi_model* m);
 /** \brief evaluate the parametric definition of a model boundary entity