Hello Christophe, now I have a doubt about the implementation of the Netgen interface. Lc is computed with : SPoint3 pt = gr->bounds().center(); double lc = BGM_MeshSize(gr, 0, 0, pt.x(), pt.y(), pt.z()); But what happens if the region has a hole in the center and the center point does not belong to the region ? Walter