diff --git a/Mesh/meshGFace.cpp b/Mesh/meshGFace.cpp index d3db3beffb1f261a34f9bad50558744d66d3f5b1..1a10eac9d2d3fa689b333caf63f131af979c2856 100644 --- a/Mesh/meshGFace.cpp +++ b/Mesh/meshGFace.cpp @@ -1,4 +1,4 @@ -// $Id: meshGFace.cpp,v 1.73 2007-04-16 09:08:28 remacle Exp $ +// $Id: meshGFace.cpp,v 1.74 2007-04-21 22:08:30 geuzaine Exp $ // // Copyright (C) 1997-2007 C. Geuzaine, J.-F. Remacle // @@ -725,8 +725,9 @@ bool gmsh2DMeshGenerator ( GFace *gf , bool debug = true) } - /// Increase the size of the bounding box by 20 % - bbox *= 1.5; + /// Increase the size of the bounding box + bbox *= 2.5; + /// add 4 points than encloses the domain /// Use negative number to distinguish thos fake vertices MVertex *bb[4]; diff --git a/benchmarks/2d/bounding_box_pb.geo b/benchmarks/2d/bounding_box_pb.geo new file mode 100644 index 0000000000000000000000000000000000000000..63a70ad370ac275e81f3f2f5504aca7180da6b3a --- /dev/null +++ b/benchmarks/2d/bounding_box_pb.geo @@ -0,0 +1,12 @@ + +// this is a simple example for which computing a bounding box with +// bbox*=1.5 fails... + +Point (1) = {0, 0, 0, 1}; +Point (2) = {0.1, 0.1, 0, 1}; +Point (3) = {0.1, 0.2, 0, 1}; +Line (1) = {1, 2}; +Line (2) = {2, 3}; +Line (3) = {3, 1}; +Line Loop(1) = {1, 2, 3}; +Plane Surface(1) = {1};