Hello,
I have tried generating a mesh that is periodic in X and Y using the
Extrude command. This seems to work using the Linux version 2.0.8 (Build
date Thu Jan 3 15:55:51 UTC 2008).
When I try to use the same file on Windows XP version 2.4.2, the program
fails to build the mesh.
The mesh definition file and the Windows message console output are in
the attached files.
What is wrong here? Any help will be greatly appreciated.
Thank you.
Chris.
--
----------------------------
Dr Christopher J.P. Newton
Information Surfaces Laboratory
Hewlett-Packard Laboratories
Long Down Avenue
Stoke Gifford
BRISTOL
BS34 8QZ
United Kingdom
Tel: +44 (0)117 312 8727
Fax: +44 (0)117 312 9937
// Generate a mesh for a box which is periodic in x and y
// Parameters used when reading the mesh later
periodicBoundaryBaseNo=80;
periodicInX=0;
periodicInY=2;
periodicInZ=4;
// pitch
p = 1.0;
// lc mesh (bulk) parameter
lcm = 0.2;
// cell thickness
d = 1.5;
Point (1) = {0, 0, 0, lcm};
Point (2) = {p, 0, 0, lcm};
Point (3) = {p, 0, d, lcm};
Point (4) = {0, 0, d, lcm};
Line (1) = {1,2};
Line (2) = {2,3};
Line (3) = {3,4};
Line (4) = {4,1};
Line Loop (1) = {1, 2, 3, 4};
Plane Surface (1) = {1};
out[]=Extrude {0,p,0} {
Surface{1}; Layers{5};
};
Printf("out[0] = %g", out[0]);
Printf("out[1] = %g", out[1]);
Printf("out[2] = %g", out[2]);
Printf("out[3] = %g", out[3]);
Printf("out[4] = %g", out[4]);
Printf("out[5] = %g", out[5]);
// Delete unwanted entities - leaving the two surfaces periodic in Y and lines
// in the YoZ plane, which will now be used to form a new Plane Surface to be
// extruded.
Delete{
Volume{1};
Surface{13, 17, 21, 25};
Line{12, 16};
}
// Form the new Plane Surface and extrude it in the X direction.
Line Loop(2) = {4, -20, -9, 11};
Plane Surface(2)={2};
out[]=Extrude {p,0,0} {
Surface{2}; Layers{1};
};
Printf("out[0] = %g", out[0]);
Printf("out[1] = %g", out[1]);
Printf("out[2] = %g", out[2]);
Printf("out[3] = %g", out[3]);
Printf("out[4] = %g", out[4]);
Printf("out[5] = %g", out[5]);
// Delete the bottom and top ruled surfaces, these will be re-defined as Plane
// Surfaces so that they are meshed at the same resolution as the others.
Delete{
Volume{1};
Surface{39, 47};
}
// Form the new bottom and top surfaces, in a real application these will be
// more complicated.
Line Loop(3) = {-3, -31, 8, -20};
Plane Surface(3) = {3};
Line Loop(4) = {6, -29, -1, 11};
Plane Surface(4) = {4};
// Now define the volume
Surface Loop(100) = {2, 26, 3, 1, 4, 48};
Volume(100) = {100};
// Define the physical entities
pbNo=periodicBoundaryBaseNo+periodicInX;
Physical Surface(pbNo) = {2}; // left periodic
Physical Surface(pbNo+1) = {48}; // right periodic
pbNo=periodicBoundaryBaseNo+periodicInY;
Physical Surface(pbNo) = {1}; // front periodic
Physical Surface(pbNo+1) = {26}; // back periodic
Physical Surface(0)={4}; // Bottom surface
Physical Surface(1)={3}; // Top surface
Physical Volume(0)={100}; // Single region
// Use the frontal algo to make sure that the periodic surface mesh is
// not modified when meshing the volume
Mesh.Algorithm3D = 4;
Info : -------------------------------------------------------
Info : Gmsh version : 2.4.2
Info : Build OS : Windows
Info : Build options : Ann Blas Chaco DIntegration FlTree Fltk FltkJpeg
FltkPng FltkZlib Gmm Lapack MathEx Med Metis NativeFileChooser Netgen
NoSocklenT OpenCascade OpenGL Tetgen
Info : Build date : 20090921
Info : Build host : christop-e38ece
Info : Packager : Administrator
Info : Home directory : C:\DOCUME~1\cjn\LOCALS~1\Temp/
Info : Launch date : Fri Sep 24 14:20:16 2010
Info : Command line : C:\_cjn\bin\gmsh.exe
Info : -------------------------------------------------------
Info : Reading '//0-uxdata-br2.labs.hpl.hp.com/cjn/testPeriodicMesh3D.geo'
out[0] = 26
out[1] = 1
out[2] = 13
out[3] = 17
out[4] = 21
out[5] = 25
out[0] = 48
out[1] = 1
out[2] = 1
out[3] = 39
out[4] = 26
out[5] = 47
Info : Read '//0-uxdata-br2.labs.hpl.hp.com/cjn/testPeriodicMesh3D.geo'
Info : Meshing 1D...
Info : Meshing curve 1 (Line)
Info : Meshing curve 2 (Line)
Info : Meshing curve 3 (Line)
Info : Meshing curve 4 (Line)
Info : Mesh 1D complete (0 s)
Info : Meshing 2D...
Info : Meshing surface 1 (Plane, MeshAdapt)
Info : Meshing surface 2 (Plane, MeshAdapt)
Info : Meshing surface 3 (Plane, MeshAdapt)
Info : Meshing surface 4 (Plane, MeshAdapt)
Info : Meshing surface 26 (extruded)
Info : Meshing surface 48 (extruded)
Info : Linear search for (1, 1, 1.071428571430256)
Error : Could not find extruded vertex (1, 1, 1.071428571430256) in surface 48
Info : Mesh 2D complete (0.15625 s)
Info : Meshing 3D...
Info : Meshing volume 100 (Frontal)
Info : Region 100 Face 2, 0 intersect
Info : Region 100 Face 26, 1 intersect
Info : Region 100 Face 3, 0 intersect
Info : Region 100 Face 1, 0 intersect
Info : Region 100 Face 4, 1 intersect
Info : CalcLocalH: 182 Points 0 Elements 338 Surface Elements
Info : Check subdomain 1 / 1
Info : 338 open elements
Error : ERROR: Edge 3 - 13 multiple times in surface mesh
Error : ERROR: Edge 4 - 18 multiple times in surface mesh
Error : ERROR: Edge 5 - 22 multiple times in surface mesh
Error : ERROR: Edge 13 - 14 multiple times in surface mesh
Error : ERROR: Edge 2 - 27 multiple times in surface mesh
Error : ERROR: Edge 14 - 15 multiple times in surface mesh
Error : ERROR: Edge 15 - 16 multiple times in surface mesh
Error : ERROR: Edge 16 - 17 multiple times in surface mesh
Error : ERROR: Edge 17 - 18 multiple times in surface mesh
Error : ERROR: Edge 7 - 31 multiple times in surface mesh
Error : ERROR: Edge 8 - 36 multiple times in surface mesh
Error : ERROR: Edge 22 - 23 multiple times in surface mesh
Error : ERROR: Edge 23 - 24 multiple times in surface mesh
Error : ERROR: Edge 2 - 46 multiple times in surface mesh
Error : ERROR: Edge 9 - 40 multiple times in surface mesh
Error : ERROR: Edge 24 - 25 multiple times in surface mesh
Error : ERROR: Edge 6 - 45 multiple times in surface mesh
Error : ERROR: Edge 25 - 26 multiple times in surface mesh
Error : ERROR: Edge 26 - 27 multiple times in surface mesh
Error : ERROR: Edge 5 - 50 multiple times in surface mesh
Error : ERROR: Edge 6 - 49 multiple times in surface mesh
Error : ERROR: Edge 3 - 54 multiple times in surface mesh
Error : ERROR: Edge 9 - 53 multiple times in surface mesh
Error : ERROR: Edge 31 - 32 multiple times in surface mesh
Error : ERROR: Edge 7 - 57 multiple times in surface mesh
Error : ERROR: Edge 4 - 61 multiple times in surface mesh
Error : ERROR: Edge 32 - 33 multiple times in surface mesh
Error : ERROR: Edge 8 - 58 multiple times in surface mesh
Error : ERROR: Edge 33 - 34 multiple times in surface mesh
Error : ERROR: Edge 34 - 35 multiple times in surface mesh
Error : ERROR: Edge 35 - 36 multiple times in surface mesh
Error : ERROR: Edge 40 - 41 multiple times in surface mesh
Error : ERROR: Edge 41 - 42 multiple times in surface mesh
Error : ERROR: Edge 42 - 43 multiple times in surface mesh
Error : ERROR: Edge 43 - 44 multiple times in surface mesh
Error : ERROR: Edge 44 - 45 multiple times in surface mesh
Error : ERROR: Edge 46 - 47 multiple times in surface mesh
Error : ERROR: Edge 47 - 48 multiple times in surface mesh
Error : ERROR: Edge 48 - 49 multiple times in surface mesh
Error : ERROR: Edge 50 - 51 multiple times in surface mesh
Error : ERROR: Edge 51 - 52 multiple times in surface mesh
Error : ERROR: Edge 52 - 53 multiple times in surface mesh
Error : ERROR: Edge 54 - 55 multiple times in surface mesh
Error : ERROR: Edge 55 - 56 multiple times in surface mesh
Error : ERROR: Edge 56 - 57 multiple times in surface mesh
Error : ERROR: Edge 58 - 59 multiple times in surface mesh
Error : ERROR: Edge 59 - 60 multiple times in surface mesh
Error : ERROR: Edge 60 - 61 multiple times in surface mesh
Error : ERROR: Surface mesh not consistent
Error : ERROR: Stop meshing since surface mesh not consistent
Info : Mesh 3D complete (0.015625 s)
Info : 213 vertices 410 elements
Error : ------------------------------
Error : Mesh generation error summary
Error : 0 warnings
Error : 51 errors
Error : Check the full log for details
Error : ------------------------------
Info : Writing '//0-uxdata-br2.labs.hpl.hp.com/cjn/meshError.txt'
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh
