fixed problem with points besides mesh

cardmech/src/electrophysiology/problem/ElphyProblem.cpp

@@ -123,7 +123,7 @@ bool  ElphyProblem::Includes(const cell &c,const Point P){
     }
     return mx < P[0] && P[0] < Mx && my < P[1] && P[1] < My && mz < P[2] && P[2] < Mz;
   }else if(c.Type()==TETRAHEDRON){
-    Point LocalP =c.GlobalToLocal(P-c.Corner(0));
+    Point LocalP =c.GlobalToLocal(P); //-c.Corner(0)
     return LocalP[0]>=0.0 && LocalP[1]>=0.0 &&LocalP[2]>=0.0 && LocalP[0]+LocalP[1]+LocalP[2]< 1.0;
   }
   THROW("Includes only implemented for hexahedral and tetrahedral cells.")