Hi Christophe, I used only second-order *line* and *triangle* elements in a 2D mesh. And I have turned the "incomplete" option off, but that didn't help. Now after reading the "Useless nodes" thread,

<span style="font-size: larger;"> I'm somehow convinced the problem is caused by these extra nodes. My program pre-sets the global sparse matrix dimensions according to the number of nodes. The filling of the global matrices are done by examining each triangles. Now, the existence of some extra nodes will leave the global matrix with zero-valued diagonal elements! And that leads to NaN solutions. I wonder if there is any way to discard, or to identify these extra nodes?

<br><br>BTW, please also confirm in v2.0 if the nodes are re-ordered by default.<br><br>Thanks for the nice package.<br><br>Min Yan<br>

<br><div><span class="gmail_quote">On 2/12/07, <b class="gmail_sendername">Christophe Geuzaine</b> <<a href="mailto:cgeuzaine@ulg.ac.be" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">cgeuzaine@ulg.ac.be

</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

The default behavior in Gmsh 2.0 is to use "incomplete" second order<br>elements (8-node quads, 20-node hexas, etc.). Maybe that's causing some<br>trouble in your code?<br><br>You can revert to using "complete" second order elements by unselecting

<br>"Use imcomplete second order elements" in Options->Mesh->General.</blockquote><div> </div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

--<br>Prof. Christophe Geuzaine<br>University of Liege, Electrical Engineering and Computer Science<br><a href="http://www.montefiore.ulg.ac.be/%7Egeuzaine" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">

http://www.montefiore.ulg.ac.be/~geuzaine</a><br></blockquote></div>

<br>