<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=gb2312">
<META content="MSHTML 6.00.2800.1400" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial>Hi,</FONT></DIV>
<DIV><FONT face=Arial>I am used to use the Ofelib(one fem lib) to resolve the 
fem problem, but when I find that itk also has the fem lib I&nbsp;intend&nbsp;to 
use it in my app. I create a .fem file and&nbsp;use the itk::fem::solver by 
reading the .fem file. However, two problems make me very puzzled.</FONT><FONT 
face=Arial></FONT></DIV>
<DIV><FONT face=Arial></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial>1.if the size of the mesh is large, e.g. the number of the 
nodes &gt; 1000, the solver will take a very very long time(1-2 hours) to finish 
the computation.&nbsp;I once use other fem libs such as Ofeli which only take 15 
minutes.</FONT></DIV>
<DIV><FONT face=Arial></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial>2.I fix some Dof when&nbsp;applying BC&nbsp;, However, the 
solution is not same with the values I set&nbsp;on these&nbsp;DOF.</FONT></DIV>
<DIV><FONT face=Arial></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial>I look at the source code of itk::fem and I find it use 
the iterative least&nbsp;square to solve the linear system. I do not know 
if&nbsp;this is the reason&nbsp;causing the very long&nbsp;time computation. 
Becaus&nbsp;the&nbsp;Ofelib&nbsp;use the CG to&nbsp;solve the linear system and 
only need very little time.</FONT></DIV>
<DIV><FONT face=Arial></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial>I only need solve the static linear elastic problem, use 
the tetrahedron as the element of the mesh. </FONT></DIV>
<DIV><FONT face=Arial></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial>&lt;MaterialLinearElasticity&gt;<BR>&nbsp;0&nbsp;%Global 
material number<BR>&nbsp;E: 3000 % E<BR>&nbsp;A: 0&nbsp; % A<BR>&nbsp;I: 0&nbsp; 
% I<BR>&nbsp;nu: 0.45&nbsp; % nu<BR>&nbsp;h: 1<BR>&nbsp;RhoC: 
1<BR>&nbsp;</FONT></DIV>
<DIV><FONT face=Arial>&lt;Element3DC0LinearTetrahedronStrain&gt;</FONT></DIV>
<DIV><FONT face=Arial></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial>thnk you in advance!</FONT></DIV>
<DIV><FONT face=Arial></FONT>&nbsp;</DIV>
<DIV><FONT size=2><FONT face=Arial size=3>Yixun Liu</FONT></DIV>
<DIV><BR></DIV></FONT></BODY></HTML>