*** empty log message ***

Mesh/meshGFace.cpp

-// $Id: meshGFace.cpp,v 1.95 2007-10-11 14:44:01 remacle Exp $
+// $Id: meshGFace.cpp,v 1.96 2007-10-11 16:08:23 remacle Exp $
 //
 // Copyright (C) 1997-2007 C. Geuzaine, J.-F. Remacle
 //
@@ -493,6 +493,27 @@ void RefineMesh ( GFace *gf, BDS_Mesh &m , const int NIT)
 
   int MAXNP = m.MAXPOINTNUMBER;
 
+  if (NIT > 0)
+    {
+      std::set<BDS_Point*,PointLessThan>::iterator itp = m.points.begin();
+      while (itp != m.points.end())
+	{
+	  std::list<BDS_Edge*>::iterator it  = (*itp)->edges.begin();
+	  std::list<BDS_Edge*>::iterator ite = (*itp)->edges.end();
+	  double L = 1.e22;
+	  while(it!=ite){
+	    double l = (*it)->length();
+	    if (l<L && (*it)->g && (*it)->g->classif_degree == 1)L=l;
+	    ++it;
+	  }
+	  if(!CTX.mesh.constrained_bgmesh)
+	    (*itp)->lc() = L;
+	  (*itp)->lcBGM() = L;
+	  ++itp;
+	}
+    }
+
+
   double OLDminL=1.E22,OLDmaxL=0;
 
   double t_spl=0, t_sw=0,t_col=0,t_sm=0;