Hi, the mesh I sent you last time was generated from the attached .geo file using gmsh (version 2.4.2, which comes with Ubuntu 10.04, 64bit). I am working on two machines (AMD Phenom and Intel Core2Duo (Centrino2), respectively) all running Ubuntu 10.04 (64bits each).
CU andi Am 23.03.2012 16:27, schrieb Jonathan Guyer: > > On Mar 14, 2012, at 3:04 AM, Andreas Hasenkopf wrote: > >> attached are my mesh and the script. In the gmsh options I have >> unchecked everything but triangles and tetrahedra as fipy does not >> support something else (as written in the fipy documentation). > > Your meshes have illuminated some issues with Gmsh import/export and Vtk > export that I am working on resolving, however the meshes themselves are > pathological (see my comments at > http://matforge.org/fipy/ticket/434#comment:3). They have very poorly shaped > elements (very flat tets will lead to inaccurate solutions), but worse, some > cells overlap one-another, which is not permissible. > > How did you generate these meshes? > With Gmsh, I assume? > Can you send me the .geo files you used? > What version of Gmsh? > What platform? > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] -- Andreas Hasenkopf Phone: +49 151 11728439 Homepage: http://www.hasenkopf2000.net GPG Pub Key: http://goo.gl/4mOsM
Point(1) = {-60, -60, 0, 1e+22};
Point(2) = {-60, 60, 0, 1e+22};
Point(3) = {60, 60, 0, 1e+22};
Point(4) = {60, -60, 0, 1e+22};
Point(5) = {22, 22, 0, 1e+22};
Point(6) = {22, 48.5, 0, 1e+22};
Point(7) = {22, 50.5, 0, 1e+22};
Point(8) = {22, -4.5, 0, 1e+22};
Point(9) = {21, -6.5, 0, 1e+22};
Point(10) = {23, -6.5, 0, 1e+22};
Point(11) = {-6.5, 22, 0, 1e+22};
Point(12) = {-4.5, 22, 0, 1e+22};
Point(13) = {48.5, 22, 0, 1e+22};
Point(14) = {50.5, 22, 0, 1e+22};
Point(15) = {21, -60, 0, 1e+22};
Point(16) = {23, -60, 0, 1e+22};
Point(17) = {-22, -22, 0, 1e+22};
Point(18) = {-50.5, -22, 0, 1e+22};
Point(19) = {-48.5, -22, 0, 1e+22};
Point(20) = {4.5, -22, 0, 1e+22};
Point(21) = {6.5, -22, 0, 1e+22};
Point(22) = {-22, 6.5, 0, 1e+22};
Point(23) = {-22, 4.5, 0, 1e+22};
Point(24) = {-22, -48.5, 0, 1e+22};
Point(25) = {-21, -50.5, 0, 1e+22};
Point(26) = {-23, -50.5, 0, 1e+22};
Point(27) = {-23, -60, 0, 1e+22};
Point(28) = {-21, -60, 0, 1e+22};
Point(29) = {-60, -60, 10, 1e+22};
Point(30) = {-60, 60, 10, 1e+22};
Point(34) = {60, 60, 10, 1e+22};
Point(38) = {60, -60, 10, 1e+22};
Point(42) = {23, -60, 10, 1e+22};
Point(46) = {23, -6.5, 10, 1e+22};
Point(50) = {22, 22, 10, 1e+22};
Point(51) = {50.5, 22, 10, 1e+22};
Point(56) = {22, 50.5, 10, 1e+22};
Point(61) = {-6.5, 22, 10, 1e+22};
Point(66) = {21, -6.5, 10, 1e+22};
Point(70) = {21, -60, 10, 1e+22};
Point(74) = {-21, -60, 10, 1e+22};
Point(78) = {-21, -50.5, 10, 1e+22};
Point(82) = {-22, -22, 10, 1e+22};
Point(83) = {6.5, -22, 10, 1e+22};
Point(88) = {-22, 6.5, 10, 1e+22};
Point(93) = {-50.5, -22, 10, 1e+22};
Point(98) = {-23, -50.5, 10, 1e+22};
Point(102) = {-23, -60, 10, 1e+22};
Line(1) = {2, 3};
Line(2) = {3, 4};
Line(3) = {4, 16};
Line(4) = {16, 15};
Line(5) = {15, 28};
Line(6) = {28, 27};
Line(7) = {27, 1};
Line(8) = {1, 2};
Line(9) = {27, 26};
Line(10) = {25, 28};
Line(11) = {15, 9};
Line(12) = {10, 16};
Circle(13) = {23, 17, 20};
Circle(14) = {20, 17, 24};
Circle(15) = {24, 17, 19};
Circle(16) = {19, 17, 23};
Circle(17) = {22, 17, 21};
Circle(18) = {21, 17, 25};
Circle(19) = {26, 17, 18};
Circle(20) = {18, 17, 22};
Circle(21) = {6, 5, 13};
Circle(22) = {13, 5, 8};
Circle(23) = {8, 5, 12};
Circle(24) = {12, 5, 6};
Circle(25) = {7, 5, 14};
Circle(26) = {14, 5, 10};
Circle(27) = {9, 5, 11};
Circle(28) = {11, 5, 7};
Line(36) = {29, 30};
Line(37) = {30, 34};
Line(38) = {34, 38};
Line(39) = {38, 42};
Line(40) = {42, 46};
Circle(41) = {46, 50, 51};
Circle(42) = {51, 50, 56};
Circle(43) = {56, 50, 61};
Circle(44) = {61, 50, 66};
Line(45) = {66, 70};
Line(46) = {70, 74};
Line(47) = {74, 78};
Circle(48) = {78, 82, 83};
Circle(49) = {83, 82, 88};
Circle(50) = {88, 82, 93};
Circle(51) = {93, 82, 98};
Line(52) = {98, 102};
Line(53) = {102, 29};
Line(55) = {1, 29};
Line(56) = {2, 30};
Line(60) = {3, 34};
Line(64) = {4, 38};
Line(68) = {16, 42};
Line(72) = {10, 46};
Line(76) = {14, 51};
Line(80) = {7, 56};
Line(84) = {11, 61};
Line(88) = {9, 66};
Line(92) = {15, 70};
Line(96) = {28, 74};
Line(100) = {25, 78};
Line(104) = {21, 83};
Line(108) = {22, 88};
Line(112) = {18, 93};
Line(116) = {26, 98};
Line(120) = {27, 102};
Line(133) = {42, 70};
Line(167) = {74, 102};
Line Loop(30) = {8, 1, 2, 3, -12, -26, -25, -28, -27, -11, 5, -10, -18, -17,
-20, -19, -9, 7};
Plane Surface(30) = {30};
Line Loop(32) = {18, 10, 6, 9, 19, 20, 17};
Plane Surface(32) = {32};
Line Loop(34) = {27, 28, 25, 26, 12, 4, 11};
Plane Surface(34) = {34};
Line Loop(57) = {8, 56, -36, -55};
Ruled Surface(57) = {57};
Line Loop(61) = {1, 60, -37, -56};
Ruled Surface(61) = {61};
Line Loop(65) = {2, 64, -38, -60};
Ruled Surface(65) = {65};
Line Loop(69) = {3, 68, -39, -64};
Ruled Surface(69) = {69};
Line Loop(73) = {-12, 72, -40, -68};
Ruled Surface(73) = {73};
Line Loop(77) = {-26, 76, -41, -72};
Ruled Surface(77) = {77};
Line Loop(81) = {-25, 80, -42, -76};
Ruled Surface(81) = {81};
Line Loop(85) = {-28, 84, -43, -80};
Ruled Surface(85) = {85};
Line Loop(89) = {-27, 88, -44, -84};
Ruled Surface(89) = {89};
Line Loop(93) = {-11, 92, -45, -88};
Ruled Surface(93) = {93};
Line Loop(97) = {5, 96, -46, -92};
Ruled Surface(97) = {97};
Line Loop(101) = {-10, 100, -47, -96};
Ruled Surface(101) = {101};
Line Loop(105) = {-18, 104, -48, -100};
Ruled Surface(105) = {105};
Line Loop(109) = {-17, 108, -49, -104};
Ruled Surface(109) = {109};
Line Loop(113) = {-20, 112, -50, -108};
Ruled Surface(113) = {113};
Line Loop(117) = {-19, 116, -51, -112};
Ruled Surface(117) = {117};
Line Loop(121) = {-9, 120, -52, -116};
Ruled Surface(121) = {121};
Line Loop(125) = {7, 55, -53, -120};
Ruled Surface(125) = {125};
Line Loop(126) = {36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50,
51, 52, 53};
Plane Surface(126) = {126};
Line Loop(158) = {4, 92, -133, -68};
Ruled Surface(158) = {158};
Line Loop(163) = {-44, -43, -42, -41, -40, 133, -45};
Plane Surface(163) = {163};
Line Loop(183) = {6, 120, -167, -96};
Ruled Surface(183) = {183};
Line Loop(200) = {-48, -47, 167, -52, -51, -50, -49};
Plane Surface(200) = {200};
Surface Loop(1) = {30, 126, 57, 61, 65, 69, 73, 77, 81, 85, 89, 93, 97, 101,
105, 109, 113, 117, 121, 125};
Volume(1) = {1};
Surface Loop(2) = {34, 163, 89, 85, 81, 77, 73, 158, 93};
Volume(2) = {2};
Surface Loop(3) = {32, 200, 105, 101, 183, 121, 117, 113, 109};
Volume(3) = {3};
signature.asc
Description: OpenPGP digital signature
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
