Jim,<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="font-family:arial"><br>
<div style="line-height:1.7;font-size:14px">    In fact, I will deal with the large sparse symmetric or unsymmetric (depending on the constitutive</div>
<div style="line-height:1.7;font-size:14px">model used for the soil and the formulation of equations) linear system which usually occurs in geotechnical </div>
<div style="line-height:1.7;font-size:14px">problems. As you know, In most geotechnical problems, the coefficient matrix can be severely</div>
<div style="line-height:1.7;font-size:14px">ill-conditioned, thus calling  for the development of robust and efficient preconditioners and iterative methods.</div>
<div style="line-height:1.7;font-size:14px">    </div>
<div style="line-height:1.7;font-size:14px">    The iterative method and parallel solver are my choices, actually are required by my professor. As I know,</div>
<div style="line-height:1.7;font-size:14px">MUMPS is just an external solver like Matlab. So is it a right way to analyse the preconditioners and iterative</div></div></blockquote><div><span class="Apple-style-span" style="font-size:14px;line-height:23px">MUMPS and Superlu_dist are parallel sparse direct solvers for ill-conditioned problems.</span></div>
<div><span class="Apple-style-span" style="font-size:14px;line-height:23px">They can solve problems with 10k unknowns sequentially and scale up to</span></div><div><span class="Apple-style-span" style="font-size:14px;line-height:23px">10+ or 100+ processors, i.e. solve problems with 100k </span><span class="Apple-style-span" style="font-size:14px;line-height:23px">unknowns.</span></div>
<div><span class="Apple-style-span" style="font-size:14px;line-height:23px">Give them a try. Then, explore other preconditioners.</span></div><div><span class="Apple-style-span" style="font-size:14px;line-height:23px"><br>
</span></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="font-family:arial">
<div style="line-height:1.7;font-size:14px">methods provided by PETSc like CG, GMRES, Jacobi, sor, etc. for my work? </div>
<div style="line-height:1.7;font-size:14px">    </div>
<div style="line-height:1.7;font-size:14px">    Also I do not quite understand the second part of what Barry said ('You will not be able to use the external</div>
<div style="line-height:1.7;font-size:14px">solvers with those mode'). So, what does 'those mode' mean. Using quad precision?</div></div></blockquote><div><br></div><div>Yes, high precisions. Start from small problems with robust solvers, then slowly</div>
<div>move to large, complex problems.</div><div><br></div><div>Hong </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="font-family:arial">
<div style="line-height:1.7;font-size:14px">  </div><div style="line-height:1.7;font-size:14px"><br></div>
<div style="line-height:1.7;font-size:14px"></div>
<div style="line-height:1.7;font-size:14px"></div>
<div style="line-height:1.7;font-size:14px"><br></div><pre style="line-height:1.7;font-size:14px"><br>At <a href="tel:2012-06-28%C2%A002" value="+12012062802" target="_blank">2012-06-28 02</a>:38:25,"Barry Smith" <<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>> wrote:
>
>   Also if it the condition number is, say, larger than 10^10 you likely would benefit from using quad precision for the calculations.  <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#precision" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html#precision</a>
>
>   You will not be able to use the external solvers with those mode but since it is 1d likely running sequentially will be fine.
>
>    Barry
>
>
>On Jun 27, 2012, at 9:46 AM, Matthew Knepley wrote:
>
>> On Wed, Jun 27, 2012 at 7:55 AM, w_ang_temp <<a href="mailto:w_ang_temp@163.com" target="_blank">w_ang_temp@163.com</a>> wrote:
>> Hello,
>> 
>>     I want to deal with the ill-conditioned problems of Terzaghi's one-dimensional consolidation
>> 
>> finite element equations in geotechnical engineering. So is there special support on ill-conditioned
>> 
>> problems using ksp in PETSc.
>> 
>> 
>> Use a direct solver, like built-in LU,  or in parallel, MUMPS (--download-mumps).
>> 
>>   Matt
>>  
>>     Thanks.
>> 
>>                                                      Jim
>> 
>> 
>> 
>> 
>> 
>> 
>> -- 
>> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
>> -- Norbert Wiener
>
</pre></div><br><br><span title="neteasefooter"><span></span></span></blockquote></div><br>