<div class="gmail_quote">On Mon, Aug 20, 2012 at 7:07 AM, w_ang_temp <span dir="ltr"><<a href="mailto:w_ang_temp@163.com" target="_blank">w_ang_temp@163.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="line-height:1.7;font-size:14px;font-family:arial"><div>Hello,</div>
<div>    I have a problem about the ill-conditioned matrix. I deal with the soil-water coupled problem in</div>
<div>geotechnical engineering. It can be demonstrated that the stiffness matrix of soil-water coupled problem is </div>
<div>ill-conditioned due to the coupled effects of large relative differences in stiffnesses and permeabilities </div>
<div>of materials. But it should not be ill-conditioned for soil-water uncoupled problem, although there is no derivation </div>
<div>in papers.</div>
<div>    In my three models, the condition numbers are 3.5595E+05, 1.8849E+06, 9.0206E+07 for soil-water </div>
<div>coupled models and 5.2885E+04, 3.3826 E+05, 3.4565E+05 for soil-water uncoupled ones, respectively.</div>
<div>Therefore, in my opinion, the uncoupled models are also ill-conditioned due to the large condition</div>
<div>numbers.</div>
<div>    However, when using GMRES, CG and CGS, the three models of uncoupled ones show rapid convergent rate  </div>
<div>while coupled models have bad convergence performance.</div></div></blockquote><div><br></div><div>This indicates that the spectrum of the simpler model has a small number of outliers. It could be as simple as boundary condition scaling or it could be more meaningful. The coupled model does not have that nice property, so convergence is slower. You should start with:<br>
<br>1. non-dimensionalize: scale your model to use some representative units so that the solution and residual are well-scaled</div><div><br></div><div>2. find an effective preconditioner. PCFIELDSPLIT might work well, but you have to understand methods to know how to drive it. Look at the examples and read through the literature in your field to find something you want to implement, then ask here if you need suggestions for how to implement. (It can often be done using only command line options.)</div>
</div>