<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 11 (filtered medium)">
<!--[if !mso]>
<style>
v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style>
<![endif]-->
<style>
<!--
/* Font Definitions */
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman";}
h1
{margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:14.2pt;
margin-bottom:.0001pt;
text-indent:-14.2pt;
page-break-after:avoid;
mso-list:l0 level1 lfo3;
font-size:14.0pt;
font-family:"Times New Roman";}
h2
{margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:1.0cm;
margin-bottom:.0001pt;
text-indent:-1.0cm;
page-break-after:avoid;
mso-list:l0 level2 lfo3;
font-size:12.0pt;
font-family:"Times New Roman";}
h3
{margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:35.45pt;
margin-bottom:.0001pt;
text-indent:-35.45pt;
page-break-after:avoid;
mso-list:l0 level3 lfo3;
font-size:12.0pt;
font-family:"Times New Roman";
font-style:italic;}
h4
{margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:42.55pt;
margin-bottom:.0001pt;
text-indent:-42.55pt;
page-break-after:avoid;
mso-list:l0 level4 lfo3;
font-size:12.0pt;
font-family:"Times New Roman";
font-weight:normal;
font-style:italic;}
h5
{margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:49.65pt;
margin-bottom:.0001pt;
text-indent:-49.65pt;
mso-list:l0 level5 lfo3;
font-size:12.0pt;
font-family:"Times New Roman";
font-weight:normal;}
h6
{margin-top:12.0pt;
margin-right:0cm;
margin-bottom:3.0pt;
margin-left:0cm;
text-indent:0cm;
mso-list:l0 level6 lfo3;
font-size:12.0pt;
font-family:"Times New Roman";
font-weight:normal;
font-style:italic;}
p.MsoHeading7, li.MsoHeading7, div.MsoHeading7
{margin-top:12.0pt;
margin-right:0cm;
margin-bottom:3.0pt;
margin-left:0cm;
text-indent:0cm;
mso-list:l0 level7 lfo3;
font-size:10.0pt;
font-family:Arial;}
p.MsoHeading8, li.MsoHeading8, div.MsoHeading8
{margin-top:12.0pt;
margin-right:0cm;
margin-bottom:3.0pt;
margin-left:0cm;
text-indent:0cm;
mso-list:l0 level8 lfo3;
font-size:10.0pt;
font-family:Arial;
font-style:italic;}
p.MsoHeading9, li.MsoHeading9, div.MsoHeading9
{margin-top:12.0pt;
margin-right:0cm;
margin-bottom:3.0pt;
margin-left:0cm;
text-indent:0cm;
mso-list:l0 level9 lfo3;
font-size:9.0pt;
font-family:Arial;
font-weight:bold;
font-style:italic;}
p.MsoCaption, li.MsoCaption, div.MsoCaption
{margin-top:6.0pt;
margin-right:0cm;
margin-bottom:6.0pt;
margin-left:42.55pt;
text-indent:-42.55pt;
font-size:9.0pt;
font-family:"Times New Roman";
font-weight:bold;}
p.MsoTitle, li.MsoTitle, div.MsoTitle
{margin:0cm;
margin-bottom:.0001pt;
text-align:center;
font-size:16.0pt;
font-family:"Times New Roman";
text-transform:uppercase;
font-weight:bold;}
p.MsoSubtitle, li.MsoSubtitle, div.MsoSubtitle
{margin:0cm;
margin-bottom:.0001pt;
text-align:center;
font-size:12.0pt;
font-family:"Times New Roman";
font-weight:bold;}
a:link, span.MsoHyperlink
{color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{color:blue;
text-decoration:underline;}
p.Autor, li.Autor, div.Autor
{margin:0cm;
margin-bottom:.0001pt;
text-align:center;
font-size:11.0pt;
font-family:"Times New Roman";
text-transform:uppercase;}
p.Ecuacin, li.Ecuacin, div.Ecuacin
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman";
font-weight:bold;}
p.Organizacin, li.Organizacin, div.Organizacin
{margin:0cm;
margin-bottom:.0001pt;
text-align:center;
font-size:12.0pt;
font-family:"Times New Roman";}
p.Sangra1, li.Sangra1, div.Sangra1
{margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:35.45pt;
margin-bottom:.0001pt;
text-indent:-21.25pt;
font-size:12.0pt;
font-family:"Times New Roman";}
p.Sangra2, li.Sangra2, div.Sangra2
{margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:49.65pt;
margin-bottom:.0001pt;
text-indent:-14.2pt;
font-size:12.0pt;
font-family:"Times New Roman";}
p.Sangra3, li.Sangra3, div.Sangra3
{margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:63.8pt;
margin-bottom:.0001pt;
text-indent:-14.15pt;
font-size:12.0pt;
font-family:"Times New Roman";}
p.Tablas, li.Tablas, div.Tablas
{margin-top:3.0pt;
margin-right:0cm;
margin-bottom:3.0pt;
margin-left:0cm;
text-align:center;
font-size:12.0pt;
font-family:"Times New Roman";}
p.autor0, li.autor0, div.autor0
{margin:0cm;
margin-bottom:.0001pt;
text-align:center;
font-size:11.0pt;
font-family:"Times New Roman";
text-transform:uppercase;}
p.ecuacin0, li.ecuacin0, div.ecuacin0
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman";
font-weight:bold;}
p.organizacin0, li.organizacin0, div.organizacin0
{margin:0cm;
margin-bottom:.0001pt;
text-align:center;
font-size:12.0pt;
font-family:"Times New Roman";}
p.sangra10, li.sangra10, div.sangra10
{margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:35.45pt;
margin-bottom:.0001pt;
text-indent:-21.25pt;
font-size:12.0pt;
font-family:"Times New Roman";}
p.sangra20, li.sangra20, div.sangra20
{margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:49.65pt;
margin-bottom:.0001pt;
text-indent:-14.2pt;
font-size:12.0pt;
font-family:"Times New Roman";}
p.sangra30, li.sangra30, div.sangra30
{margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:63.8pt;
margin-bottom:.0001pt;
text-indent:-14.15pt;
font-size:12.0pt;
font-family:"Times New Roman";}
p.tablas0, li.tablas0, div.tablas0
{margin-top:3.0pt;
margin-right:0cm;
margin-bottom:3.0pt;
margin-left:0cm;
text-align:center;
font-size:12.0pt;
font-family:"Times New Roman";}
span.EstiloCorreo34
{mso-style-type:personal;
font-family:Arial;
color:navy;}
span.EstiloCorreo35
{mso-style-type:personal-reply;
font-family:Arial;
color:navy;}
@page Section1
{size:595.3pt 841.9pt;
margin:70.9pt 3.0cm 70.9pt 3.0cm;}
div.Section1
{page:Section1;}
/* List Definitions */
@list l0
{mso-list-id:-5;
mso-list-template-ids:-1;}
@list l0:level1
{mso-level-style-link:"T\00EDtulo 1";
mso-level-text:%1;
mso-level-tab-stop:none;
mso-level-number-position:left;
mso-level-legacy:yes;
mso-level-legacy-indent:0cm;
mso-level-legacy-space:7.2pt;
margin-left:0cm;
text-indent:0cm;}
@list l0:level2
{mso-level-style-link:"T\00EDtulo 2";
mso-level-text:"%1\.%2";
mso-level-tab-stop:none;
mso-level-number-position:left;
mso-level-legacy:yes;
mso-level-legacy-indent:0cm;
mso-level-legacy-space:7.2pt;
margin-left:0cm;
text-indent:0cm;}
@list l0:level3
{mso-level-style-link:"T\00EDtulo 3";
mso-level-text:"%1\.%2\.%3";
mso-level-tab-stop:none;
mso-level-number-position:left;
mso-level-legacy:yes;
mso-level-legacy-indent:0cm;
mso-level-legacy-space:7.2pt;
margin-left:0cm;
text-indent:0cm;}
@list l0:level4
{mso-level-style-link:"T\00EDtulo 4";
mso-level-text:"%1\.%2\.%3\.%4";
mso-level-tab-stop:none;
mso-level-number-position:left;
mso-level-legacy:yes;
mso-level-legacy-indent:0cm;
mso-level-legacy-space:7.2pt;
margin-left:0cm;
text-indent:0cm;}
@list l0:level5
{mso-level-style-link:"T\00EDtulo 5";
mso-level-text:"%1\.%2\.%3\.%4\.%5";
mso-level-tab-stop:none;
mso-level-number-position:left;
mso-level-legacy:yes;
mso-level-legacy-indent:0cm;
mso-level-legacy-space:7.2pt;
margin-left:0cm;
text-indent:0cm;}
@list l0:level6
{mso-level-style-link:"T\00EDtulo 6";
mso-level-text:"%1\.%2\.%3\.%4\.%5\.%6";
mso-level-tab-stop:none;
mso-level-number-position:left;
mso-level-legacy:yes;
mso-level-legacy-indent:0cm;
mso-level-legacy-space:7.2pt;
margin-left:0cm;
text-indent:0cm;}
@list l0:level7
{mso-level-style-link:"T\00EDtulo 7";
mso-level-text:"%1\.%2\.%3\.%4\.%5\.%6\.%7";
mso-level-tab-stop:none;
mso-level-number-position:left;
mso-level-legacy:yes;
mso-level-legacy-indent:0cm;
mso-level-legacy-space:7.2pt;
margin-left:0cm;
text-indent:0cm;}
@list l0:level8
{mso-level-style-link:"T\00EDtulo 8";
mso-level-text:"%1\.%2\.%3\.%4\.%5\.%6\.%7\.%8";
mso-level-tab-stop:none;
mso-level-number-position:left;
mso-level-legacy:yes;
mso-level-legacy-indent:0cm;
mso-level-legacy-space:7.2pt;
margin-left:0cm;
text-indent:0cm;}
@list l0:level9
{mso-level-style-link:"T\00EDtulo 9";
mso-level-text:"%1\.%2\.%3\.%4\.%5\.%6\.%7\.%8\.%9";
mso-level-tab-stop:none;
mso-level-number-position:left;
mso-level-legacy:yes;
mso-level-legacy-indent:0cm;
mso-level-legacy-space:7.2pt;
margin-left:0cm;
text-indent:0cm;}
ol
{margin-bottom:0cm;}
ul
{margin-bottom:0cm;}
-->
</style>
</head>
<body lang=ES-TRAD link=blue vlink=blue>
<div class=Section1>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
style='font-size:10.0pt;font-family:"Courier New"'>Hello, <o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>You can insert
points or lines in any surface. Please see the example below.<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>I hope this help
you.<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>Regards,<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>Francisco<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>//****************************************************************<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
style='font-size:10.0pt;font-family:"Courier New"'>len=0.1;<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
style='font-size:10.0pt;font-family:"Courier New"'>Point(1) = {0, 0, 0, len};<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
style='font-size:10.0pt;font-family:"Courier New"'>ext1[]=Extrude {1, 0, 0} {<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
style='font-size:10.0pt;font-family:"Courier New"'> Point{1};<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
style='font-size:10.0pt;font-family:"Courier New"'>};<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>ext2[]=Extrude
{0, 1, 0} {<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>
Line{ext1[1]};<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>};<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>my_point = newp;<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>Point(my_point) =
{0.34, 0.45, 0, len};<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>// Inserting a
point in a surface:<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>Point{my_point}
In Surface{ext2[1]};<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>my_p1=newp;
Point(my_p1) = {0.2, 0.78, 0, len};<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>my_p2=newp;
Point(my_p2) = {0.6, 0.67, 0, len};<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>my_line=newl;<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>Line(my_line)={my_p1,my_p2};<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>// Inserting a
line in surface<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'>Line{my_line} In
Surface{ext2[1]};<o:p></o:p></span></font></p>
<p class=MsoNormal style='text-autospace:none'><font size=2 face="Courier New"><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=2 face="Courier New"><span lang=EN-GB
style='font-size:10.0pt;font-family:"Courier New"'>//****************************************************************</span></font><font
size=2 color=navy face=Arial><span lang=EN-GB style='font-size:10.0pt;
font-family:Arial;color:navy'><o:p></o:p></span></font></p>
<p class=MsoNormal><font size=2 color=navy face=Arial><span lang=EN-GB
style='font-size:10.0pt;font-family:Arial;color:navy'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=2 color=navy face=Arial><span lang=EN-GB
style='font-size:10.0pt;font-family:Arial;color:navy'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=2 color=navy face=Arial><span lang=EN-GB
style='font-size:10.0pt;font-family:Arial;color:navy'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=2 color=navy face=Arial><span lang=EN-GB
style='font-size:10.0pt;font-family:Arial;color:navy'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=2 color=navy face=Arial><span lang=EN-GB
style='font-size:10.0pt;font-family:Arial;color:navy'><o:p> </o:p></span></font></p>
<div>
<div class=MsoNormal align=center style='text-align:center'><font size=3
face="Times New Roman"><span lang=ES style='font-size:12.0pt'>
<hr size=2 width="100%" align=center tabindex=-1>
</span></font></div>
<p class=MsoNormal><b><font size=2 face=Tahoma><span lang=ES style='font-size:
10.0pt;font-family:Tahoma;font-weight:bold'>De:</span></font></b><font size=2
face=Tahoma><span lang=ES style='font-size:10.0pt;font-family:Tahoma'> gmsh-bounces@ace20.montefiore.ulg.ac.be
[mailto:gmsh-bounces@ace20.montefiore.ulg.ac.be] <b><span style='font-weight:
bold'>En nombre de </span></b>jordi poblet<br>
<b><span style='font-weight:bold'>Enviado el:</span></b> lunes, 11 de octubre
de 2010 14:10<br>
<b><span style='font-weight:bold'>Para:</span></b> Bart Vandewoestyne<br>
<b><span style='font-weight:bold'>CC:</span></b> gmsh@geuz.org<br>
<b><span style='font-weight:bold'>Asunto:</span></b> Re: [Gmsh] pointsources and
Nodes and Elements lists</span></font><span lang=ES><o:p></o:p></span></p>
</div>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='margin-bottom:12.0pt'><font size=3
face="Times New Roman"><span style='font-size:12.0pt'>Hello Bart,<br>
<br>
I also need so often points in an specific position of the mesh. What I so is
to construct the mesh by parts in order to be sure that there will be a node in
the desired position (I copy an example for a single rectangle with a node in
X=0.49; Y = 0.48).<br>
I suppose that this is not the answer that you were expecting... sometimes it
is hard to build the meshes. In any case, I would also be interested in a
better procedure to specify a priori the positions of some nodes of the mesh.
But I guess that this is not natural for most of the meshing algorithms.<br>
<br>
Best regards,<br>
<br>
JORDI<br>
<br>
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------<br>
<br>
<br>
lc = 0.13 ; <br>
<br>
//Coords <br>
x0 = 0.0 ; <br>
y0 = 0.0 ; <br>
xs = 0.49 ; //X Point source <br>
ys = 0.48 ; //Y Point source <br>
x1 = 2.0 ; <br>
y1 = 3.0 ; <br>
<br>
//Points: <br>
Point(1) = {x0, y0, 0, lc}; <br>
Point(2) = {xs, y0, 0, lc} ; <br>
Point(3) = {x1, y0, 0, lc} ; <br>
Point(4) = {x0, ys, 0, lc} ; <br>
Point(5) = {xs, ys, 0, lc}; <br>
Point(6) = {x1, ys, 0, lc} ; <br>
Point(7) = {x0, y1, 0, lc} ; <br>
Point(8) = {xs, y1, 0, lc} ; <br>
Point(9) = {x1, y1, 0, lc}; <br>
<br>
//Lines: <br>
Line(1) = {1,2} ; <br>
Line(2) = {2,3} ; <br>
<br>
Line(3) = {4,5} ; <br>
Line(4) = {5,6} ; <br>
<br>
Line(5) = {7,8} ; <br>
Line(6) = {8,9} ; <br>
<br>
Line(7) = {1,4} ; <br>
Line(8) = {4,7} ; <br>
<br>
Line(9) = {2,5} ; <br>
Line(10) = {5,8} ; <br>
<br>
Line(11) = {3,6} ; <br>
Line(12) = {6,9} ; <br>
<br>
//Loops: <br>
Line Loop(10) = {1,9,-3,-7} ; <br>
Line Loop(11) = {2,11,-4,-9} ; <br>
Line Loop(12) = {3,10,-5,-8} ; <br>
Line Loop(13) = {4,12,-6,-10} ; <br>
<br>
// Surfaces <br>
Plane Surface(10) = {10} ; <br>
Plane Surface(11) = {11} ; <br>
Plane Surface(12) = {12} ; <br>
Plane Surface(13) = {13} ; <br>
<br>
//Physical properties (final physical surface): <br>
Physical Surface(100) = {10,11,12,13} ; <br>
<br>
<br>
//-------------------------------------------------------------------------------------------------------------------------------------
<o:p></o:p></span></font></p>
<div>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'>On Mon, Oct 11, 2010 at 10:40 AM, Bart Vandewoestyne <<a
href="mailto:Bart.Vandewoestyne@telenet.be">Bart.Vandewoestyne@telenet.be</a>>
wrote:<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'>Hello list,<br>
<br>
For an application I'm working on, we use gmsh as our mesher, we parse the .msh<br>
file and then run our solver in Matlab. From the .msh file, we extract
the X<br>
and Y coordinates of the nodes (from $Nodes) in the mesh. We also have a<br>
matrix representing the triangular elements in the mesh (extracted from<br>
$Elements).<br>
<br>
We would like our users to be able to specify pointsources in our computational<br>
domain. The location of the pointsource is an arbitrary point that a user<br>
specifies and has nothing to do with the mesh. In gmsh, we currently
specify<br>
it as follows:<br>
<br>
<br>
----------------------- small example ------------------<br>
lc = 0.6;<br>
<br>
// Corner points.<br>
Point(1) = {0, 0, 0, lc};<br>
Point(2) = {0, 1, 0, lc};<br>
Point(3) = {1, 1, 0, lc};<br>
Point(4) = {1, 0, 0, lc};<br>
<br>
// Point sources.<br>
Point(5) = {0.6, 0.6, 0, lc};<br>
Point(6) = {0.7, 0.7, 0, lc};<br>
<br>
Line(1) = {1, 2};<br>
Line(2) = {2, 3};<br>
Line(3) = {3, 4};<br>
Line(4) = {4, 1};<br>
<br>
Line Loop(5) = {1, 2, 3, 4};<br>
Plane Surface(6) = {5};<br>
<br>
Physical Point("Pointsource") = {5, 6};<br>
Physical Surface("My Surface") = {6};<br>
----------------------- end small example ------------------<br>
<br>
<br>
And now here is the problem: if I mesh this .geo file with<br>
<br>
gmsh -2 pointsource_test.geo<br>
<br>
then this is the result:<br>
<br>
<br>
------------------ small example mesh ---------------------<br>
<br>
$MeshFormat<br>
2.1 0 8<br>
$EndMeshFormat<br>
$PhysicalNames<br>
2<br>
0 1 "Pointsource"<br>
2 2 "My Surface"<br>
$EndPhysicalNames<br>
$Nodes<br>
7<br>
1 0 0 0<br>
2 0 1 0<br>
3 1 1 0<br>
4 1 0 0<br>
5 0.6 0.6 0<br>
6 0.7 0.7 0<br>
7 0.5 0.5 0<br>
$EndNodes<br>
$Elements<br>
6<br>
1 15 3 1 5 0 5<br>
2 15 3 1 6 0 6<br>
3 2 3 2 6 0 1 7 4<br>
4 2 3 2 6 0 3 7 2<br>
5 2 3 2 6 0 7 3 4<br>
6 2 3 2 6 0 7 1 2<br>
$EndElements<br>
<br>
----------------------- end small example mesh ---------<br>
<br>
<br>
The 'problem' is that the coordinates of the pointsource appear in the list<br>
of nodal coordinates. This has two consequences:<br>
<br>
1) Our solver cannot have the pointsource coordinates in its list of x and y<br>
coordinates. It should only have the x and y coordinates of the
*real*<br>
nodes in our mesh.<br>
<br>
2) Using the tags in the $Elements section, i can of course detect which<br>
nodes are pointsources, and remove them from the $Nodes section in my<br>
Matlab parser, but removing nodes from the $Nodes section also means
that<br>
a lot of indices in the $Elements section need to change! For
example,<br>
if I remove nodes 5 and 6, then every 7 in the $Elements section has to
be<br>
replaced by a 5 (because node 5 and 6 were removed, and node 7 has
become<br>
node 5).<br>
<br>
<br>
Is there a way to circumvent this troublesome bookkeeping? Maybe somehow<br>
prevent gmsh from putting our pointsources in the list of nodes, and returning<br>
the coordinates of our pointsources in another section somehow? (I'm just<br>
guessing here...)<br>
<br>
Kind regards,<br>
Bart<br>
<font color="#888888"><span style='color:#888888'><br>
--<br>
"Share what you know. Learn what you
don't."<br>
<br>
_______________________________________________<br>
gmsh mailing list<br>
<a href="mailto:gmsh@geuz.org">gmsh@geuz.org</a><br>
<a href="http://www.geuz.org/mailman/listinfo/gmsh" target="_blank">http://www.geuz.org/mailman/listinfo/gmsh</a></span></font><o:p></o:p></span></font></p>
</div>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'><o:p> </o:p></span></font></p>
</div>
</body>
</html>