diff --git a/tutorial/t5.geo b/tutorial/t5.geo index e23dff7086a2ac0238095276caf32660f80faef1..ba2de1f0d616dc4ab0ccea9c43d60ebc896901c5 100644 --- a/tutorial/t5.geo +++ b/tutorial/t5.geo @@ -27,9 +27,10 @@ lcar3 = .055; // > gmsh t5 -clscale 0.2 // // (i.e. with all characteristic lengths divided by 5), the mesh -// counts approximately 170,000 nodes and one million tetrahedra. +// counts approximately 170,000 nodes and one million tetrahedra (and +// the computation takes 16 minutes on the same machine). -// Let us proceed by defining some elementary entities, describing a +// Let us proceed by defining some elementary entities describing a // truncated cube: Point(1) = {0.5,0.5,0.5,lcar2}; Point(2) = {0.5,0.5,0,lcar1}; @@ -66,10 +67,10 @@ Function CheeseHole // In the following commands we use the reserved variable name // `newp', which automatically selects a new point number. This // number is chosen as the highest current point number, plus - // one. (Note that, analogously to `newp', there also exists - // variables `newc', `news', `newv' and `newreg' which select the - // highest number of amongst curves, surfaces, volumes or - // any entities other than points, respectively.) + // one. (Note that, analogously to `newp', the variables `newc', + // `news', `newv' and `newreg' select the highest number amongst + // currently defined curves, surfaces, volumes and `any entities + // other than points', respectively.) p1 = newp; Point(p1) = {x, y, z, lcar3} ; p2 = newp; Point(p2) = {x+r,y, z, lcar3} ; @@ -94,7 +95,7 @@ Function CheeseHole // We need non-plane surfaces to define the spherical cheese // holes. Here we use ruled surfaces, which can have 3 or 4 - // borders: + // sides: l1 = newreg; Line Loop(l1) = {c5,c10,c4}; Ruled Surface(newreg) = {l1}; l2 = newreg; Line Loop(l2) = {c9,-c5,c1}; Ruled Surface(newreg) = {l2}; @@ -107,10 +108,10 @@ Function CheeseHole // Please note that all surface meshes are generated by projecting a // 2D planar mesh onto the surface, and that this method gives nice - // results only if the surface's curvature is relatively enough. - // If not, you will have to cut the surface in pieces. + // results only if the surface's curvature is small enough. If not, + // you will have to cut the surface in pieces. - // We then use an array of variables to store the surface loop's + // We then use an array of variables to store the surface loops // identification numbers for later reference (we will need these to // define the final volume): @@ -135,7 +136,7 @@ For t In {1:5} Call CheeseHole ; - // We define a physical volume for each cheese hole: + // We define a physical volume for each hole: Physical Volume (t) = thehole ;