<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi,<br>
Your original script works properly with my laptop. No problems for
meshing.<br>
I don't understand why you need to make a surface loop yourself and
create a volume after.<br>
<br>
I've printed the content of z1[], using: <br>
<br>
For k In {0:#z1[]-1}<br>
Printf("z1[%g] %g", k, z1[k]);<br>
EndFor<br>
As expected the list contains the top of the cuboid at index 0 and the
volume at index 1:<br>
z1[0] = 26 ==> top surface<br>
z1[1] = 1 ==> volume<br>
<br>
The rest are the other sides:<br>
z1[2] 13<br>
z1[3] 17<br>
z1[4] 21<br>
z1[5] 25<br>
<br>
Regards,<br>
Ruth<br>
<br>
<br>
<br>
Mike B. wrote:
<blockquote cite="mid:952681.62596.qm@web110004.mail.gq1.yahoo.com"
type="cite">
<table border="0" cellpadding="0" cellspacing="0">
<tbody>
<tr>
<td
style="font-family: inherit; font-style: inherit; font-variant: inherit; font-weight: inherit; font-size: inherit; line-height: inherit; font-size-adjust: inherit; font-stretch: inherit;"
valign="top">Dear Gmsh users,<br>
<br>
Can someone please confirm this:<br>
I define a rectangle (planar surface) as Surface( 2 ) and then extrude
it to a cuboid (please see full script below). However, meshing the
cuboid fails. It seems the problem is with Extrude. From the Gmsh
manual:<br>
`As explained in 4.2.1 Floating point expressions, extrude can be used
in an expression, in which case it returns a list of identification
numbers. By default, the list contains the "top" of the extruded entity
at index 0 and the extruded entity at index 1, followed by the "sides"
of the extruded entity at indices 2, 3, etc.'<br>
<br>
In the script I use<br>
Plane Surface( 2 ) = {1};<br>
z1[] = Extrude{ 0, 0, 3 }{ Surface{ 2 }; };<br>
<br>
so z1 is the variable holding the indices of the surfaces created by
extrusion. This means z1[0] should hold the index of the original
planar surface, that is, 2. However, it shows 1, which is a
non-existing surface.<br>
<br>
If the index of the original planar surface is 1 then meshing works,
that is:<br>
Plane Surface( 1 ) = {1};<br>
z1[] = Extrude{ 0, 0, 3 }{ Surface{ 1 }; };<br>
<br>
<br>
Here is the full script:<br>
<br>
Point(1) = {-1, -1, -1};<br>
Point(2) = {-1, 1, -1};<br>
Point(3) = {1, 1, -1};<br>
Point(4) = {1 ,-1, -1};<br>
Line(1) = {1,2};<br>
Line(2) = {2,3};<br>
Line(3) = {3,4};<br>
Line(4) = {4,1};<br>
Line Loop(1) = {1,2,3,4};<br>
<br>
Plane Surface( 2 ) = {1};<br>
z1[] = Extrude{ 0, 0, 3 }{ Surface{ 2 }; };<br>
<br>
// meshing works when commenting the above 2 lines and uncommenting the
following 2 lines <br>
//Plane Surface( 1 ) = {1};<br>
//z1[] = Extrude{ 0, 0, 3 }{ Surface{ 1 }; };<br>
<br>
Surface Loop( 2 ) = { z1[ 0 ], z1[ 1 ], z1[ 2 ], z1[ 3 ], z1[ 4 ], z1[
5 ] };<br>
<br>
Printf( "%g %g %g %g %g %g", z1[ 0 ], z1[ 1 ], z1[ 2 ], z1[ 3 ], z1[ 4
], z1[ 5 ] );<br>
<br>
Volume( 2 ) = { 2 };<br>
Mesh 3;<br>
<br>
Thanks,<br>
Mike.</td>
</tr>
</tbody>
</table>
<br>
<pre wrap="">
<hr size="4" width="90%">
_______________________________________________
gmsh mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gmsh@geuz.org">gmsh@geuz.org</a>
<a class="moz-txt-link-freetext" href="http://www.geuz.org/mailman/listinfo/gmsh">http://www.geuz.org/mailman/listinfo/gmsh</a>
</pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">--
Dr. Ir. Ruth V. Sabariego
University of Liege, Institut Montefiore,
Dept. of Electrical Engineering & Computer Science,
Applied & Computational Electromagnetics (ACE),
Sart Tilman Campus, Grande Traverse, 10 (B28), B-4000 LIEGE, Belgium
phone: +32-4-3663737 -- fax: +32-4-3662910 -- <a class="moz-txt-link-freetext" href="http://elap.montefiore.ulg.ac.be/">http://elap.montefiore.ulg.ac.be/</a></pre>
</body>
</html>