diff --git a/Geo/GFaceCompound.cpp b/Geo/GFaceCompound.cpp
index 14d50896f1ac64117b83a3932724a9dd7ce5e599..5dab88d033c11eab362d792b233511fe63765ec2 100644
--- a/Geo/GFaceCompound.cpp
+++ b/Geo/GFaceCompound.cpp
@@ -513,8 +513,8 @@ bool GFaceCompound::parametrize() const
   else if (_mapping == CONFORMAL){
     Msg::Debug("Parametrizing surface %d with 'conformal map'", tag());
     fillNeumannBCS();
-    //bool withoutFolding = parametrize_conformal_spectral() ;
-    bool withoutFolding = parametrize_conformal();
+    bool withoutFolding = parametrize_conformal_spectral() ;
+    //bool withoutFolding = parametrize_conformal();
     if ( withoutFolding == false ){
       //printStuff(); exit(1);
       Msg::Warning("$$$ Parametrization switched to harmonic map");
diff --git a/Solver/function.cpp b/Solver/function.cpp
index 2769248eea235d3feeb9f76477cac23dbd2f95b8..9b0ca8199ddb0eb6b472053e8e9272080fea7c7f 100644
--- a/Solver/function.cpp
+++ b/Solver/function.cpp
@@ -400,6 +400,12 @@ void function::registerBindings(binding *b){
   mb->setDescription("return the unique instance of this class");
   cb->setParentClass<function>();
 
+  cb = b->addClass<functionSolution>("functionNormals");
+  cb->setDescription("A function to access the face normals. This is a single-instance class, use the 'get' member to access the instance.");
+  mb = cb->addMethod("get",&functionNormals::get);
+  mb->setDescription("return the unique instance of this class");
+  cb->setParentClass<function>();
+
   cb = b->addClass<functionSolutionGradient>("functionSolutionGradient");
   cb->setDescription("A function to access the gradient of the solution. This is a single-instance class, use the 'get' member to access the instance.");
   mb = cb->addMethod("get",&functionCoordinates::get);
diff --git a/Solver/multiscaleLaplace.cpp b/Solver/multiscaleLaplace.cpp
index 1da275c77effda70946be825e69902eb817a3b70..4bb7d9b9d6f7c21fc5811bca618b22a464c8fdd4 100644
--- a/Solver/multiscaleLaplace.cpp
+++ b/Solver/multiscaleLaplace.cpp
@@ -653,8 +653,7 @@ static void printLevel(const char* fn,
       if (dx)fprintf(fp, "%d %22.15E %22.15E 0\n", (*it)->getIndex(), dx[2]*(p.x()-dx[0]), dx[2]*(p.y()-dx[1]));
       else   fprintf(fp, "%d %22.15E %22.15E 0\n", (*it)->getIndex(), p.x(), p.y());
     }
-    else fprintf(fp, "%d %g %g %g\n", (*it)->getIndex(),
-                 (*it)->x(), (*it)->y(), (*it)->z());
+    else fprintf(fp, "%d %g %g %g\n", (*it)->getIndex(),(*it)->x(), (*it)->y(), (*it)->z());
   }
   fprintf(fp, "$EndNodes\n");
   
@@ -984,7 +983,7 @@ void multiscaleLaplace::parametrize (multiscaleLaplaceLevel & level){
     MElement *e = level.elements[i];
     std::vector<SPoint2> localCoord;
     double local_size = localSize(e,solution);
-    if (local_size < 1.e-5*global_size) //1.e-5
+    if (local_size < 1.e-6*global_size) //1.e-5
       tooSmall.push_back(e);
     else  goodSize.push_back(e);
   }
@@ -1033,7 +1032,7 @@ void multiscaleLaplace::parametrize (multiscaleLaplaceLevel & level){
 //       double area = fabs(triangle_area(q0, q1, q2));   
 //       totArea  += area;
       double local_size = localSize(e,solution);
-      if (local_size < 1.e-7 * global_size) //1.e-7
+      if (local_size < 1.e-8 * global_size) //1.e-7
         really_small_elements = true;
     }
     //center *= (1./regions_[i].size());