Thank you, GMSH support for linux

Steffens-Jr, Alfred P Alfred.Steffens-Jr at West.Boeing.com
Mon May 21 15:45:45 CEST 2001


Christophe,

Thank you for your email on linux.  I tried Mesa 3.4.1 with GMSH 1.18, and
the problems are gone.  One other thing:  It would be nice to have more
documentation.  There are still aspects of the language grammar that I don't
understand.  I have problems with two line loops that have common lines.
But I don't know how else to build one region around an embedded area.  For
example, do you see something wrong with this .geo file which I have
included?  Thanks again,

Al Steffens
alfred.p.steffens-jr at boeing.com


//
//  check7.geo
//


lc = 3.0;

//
//  offsets for the inner square
//
Ox1 = 2.0;
Oy1 = 2.0;
Oz1 = 0.0;
wid1 = 2.0;
ht1 = 3.0;

Ox2 = 0.5;
Oy2 = 1.0;
Oz2 = 0.0;
wid2 = 0.5;
ht2 = 0.5;


//
//  Outer square
//
P1x = 0.0;
P1y = 0.0;
P1z = 0.0;

P2x = 10.0;
P2y = 0.0;
P2z = 0.0;

P3x = 10.0;
P3y = 10.0;
P3z = 0.0;

P4x = 0.0;
P4y = 10.0;
P4z = 0.0;

Point(1) = {P1x, P1y, P1z, lc};
Point(2) = {P2x, P2y, P2z, lc};
Point(3) = {P3x, P3y, P3z, lc};
Point(4) = {P4x, P4y, P4z, lc};

//
//  Inner square
//
P5x = Ox1;
P5y = Oy1;
P5z = 0.0;

P6x = Ox1 + wid1;
P6y = Oy1;
P6z = 0.0;

P7x = Ox1 + wid1;
P7y = Oy1 + ht1;
P7z = 0.0;

P8x = Ox1;
P8y = Oy1 + ht1;
P8z = 0.0;

Point(5) = {P5x, P5y, P5z, lc};
Point(6) = {P6x, P6y, P6z, lc};
Point(7) = {P7x, P7y, P7z, lc};
Point(8) = {P8x, P8y, P8z, lc};


//
//  smallest square
//
P9x = Ox1 + Ox2;
P9y = Oy1 + Oy2;
P9z = 0.0;

P10x = Ox1 + Ox2 + wid2;
P10y = Oy1 + Oy2;
P10z = 0.0;

P11x = Ox1 + Ox2 + wid2;
P11y = Oy1 + Oy2 + ht2;
P11z = 0.0;

P12x = Ox1 + Ox2;
P12y = Oy1 + Oy2 + ht2;
P12z = 0.0;

Point(9) = {P9x, P9y, P9z, lc};
Point(10) = {P10x, P10y, P10z, lc};
Point(11) = {P11x, P11y, P11z, lc};
Point(12) = {P12x, P12y, P12z, lc};


//
//  points filled in outside perimeter
//
Point(13) = {P5x, P1y, 0.0, lc};
Point(14) = {P6x, P1y, 0.0, lc};
Point(15) = {P2x, P5y, 0.0, lc};
Point(16) = {P2x, P7y, 0.0, lc};
Point(17) = {P7x, P3y, 0.0, lc};
Point(18) = {P8x, P3y, 0.0, lc};
Point(19) = {P1x, P8y, 0.0, lc};
Point(20) = {P1x, P5y, 0.0, lc};

//
//  points filled on perimeter of inner rectangle
//
Point(21) = {P9x, P5y, 0.0, lc};
Point(22) = {P10x, P5y, 0.0, lc};
Point(23) = {P6x, P10y, 0.0, lc};
Point(24) = {P6x, P11y, 0.0, lc};
Point(25) = {P11x, P7y, 0.0, lc};
Point(26) = {P12x, P7y, 0.0, lc};
Point(27) = {P5x, P12y, 0.0, lc};
Point(28) = {P5x, P9y, 0.0, lc};



//
//  lines around outside perimeter
//
Line(1) = {1,13};
Line(2) = {13,14};
Line(3) = {14,2};
Line(4) = {2,15};
Line(5) = {15,16};
Line(6) = {16,3};
Line(7) = {3,17};
Line(8) = {17,18};
Line(9) = {18,4};
Line(10) = {4,19};
Line(11) = {19,20};
Line(12) = {20, 1};


//
//  lines bracing the outer to inner rectangle
//
Line(13) = {13,5};
Line(14) = {6,14};
Line(15) = {15,6};
Line(16) = {16,7};
Line(17) = {17,7};
Line(18) = {18,8};
Line(19) = {19,8};
Line(20) = {20,5};

//
//  lines around inner rectangle perimeter
//
Line(21) = {5,21};
Line(22) = {21,22};
Line(23) = {22,6};
Line(24) = {6,23};
Line(25) = {23,24};
Line(26) = {24,7};
Line(27) = {7,25};
Line(28) = {25,26};
Line(29) = {26,8};
Line(30) = {8,27};
Line(31) = {27,28};
Line(32) = {28,5};

//
//  lines bracing ineer to smallest rectangles
//
Line(33) = {21,9};
Line(34) = {22,10};
Line(35) = {23,10};
Line(36) = {24,11};
Line(37) = {25,11};
Line(38) = {26,12};
Line(39) = {27,12};
Line(40) = {28,9};


//
//  lines around the smallest rectangle
//
Line(41) = {9,10};
Line(42) = {10,11};
Line(43) = {11,12};
Line(44) = {12,9};

// ------------------- done lines -------------------------




//----------------------------------------------------------
//
//  set up outer boundary BC
//
Physical Line(1) = {1,2,3,4,5,6,7,8,9,10,11,12};  // region 1
//----------------------------------------------------------




//----------------------------------------------------------
//
//  set up region 2 -- area just inside the outer perimeter
//
Line Loop(45) = {1,2,3,4,5,6,7,17,-26,-25,-24,-23,-22,-21,-20,12};
Line Loop(46) = {-17,8,9,10,11,20,-32,-31,-30,-29,-28,-27};

//------------------------------------------------------------



//-----------------------------------------------------------
//
//  set up region 3 -- area inside inner rectangle
//
Line Loop(47) = {21,22,23,24,25,26,27,37,-42,-41,-40,32};
Line Loop(48) = {-37,28,29,30,31,40,-44,-43};
//------------------------------------------------------------


//--------------------------------------------------------------
//
// setup region 4 -- innermost rectangle
//
Line Loop(49) = {41,42,43,44};
//--------------------------------------------------------------




Plane Surface(50) = {45,46};
Physical Surface(2) = {50};	// region 2

Plane Surface(51) = {47,48};
Physical Surface(3) = {51};	// region 3

Plane Surface(52) = {49};
Physical Surface(4) = {52};	// region 4

// end of the .geo file