Added difusive laplacian term in incomp stabilisation,

Modified FunctionEvaluator And MeshMetricBasedOnHessian that now take groups instead of solution as arguments

Geo/GRegion.cpp

@@ -310,7 +310,7 @@ double GRegion::computeSolidProperties (std::vector<double> cg,
   double surface = 0;
   cg[0] = cg[1] = cg[2] = 0.0;
   for ( ; it != l_faces.end(); ++it,++itdir){    
-    printf("face %d dir %d %d elements\n",(*it)->tag(),*itdir,(int)(*it)->triangles.size());
+    //printf("face %d dir %d %d elements\n",(*it)->tag(),*itdir,(int)(*it)->triangles.size());
     for (unsigned int i = 0; i < (*it)->triangles.size(); ++i){
       MTriangle *e = (*it)->triangles[i];
       int npt;

benchmarks/stl/skull.geo

 Mesh.RemeshParametrization=1; //(0) harmonic (1) conformal 
 Mesh.RemeshAlgorithm=1; //(0) nosplit (1) automatic (2) split metis
+
 Mesh.CharacteristicLengthFactor=0.1;
 Mesh.Algorithm3D = 4; //Frontal (4) Delaunay(1)