<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Nov 1, 2016 at 6:21 AM, Lukasz Kaczmarczyk <span dir="ltr"><<a href="mailto:Lukasz.Kaczmarczyk@glasgow.ac.uk" target="_blank">Lukasz.Kaczmarczyk@glasgow.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
<br>
I solve problems of identifying forces on the surface, which are results in displacements (measured) on some other surface. This technique is used in cell engineering where we can observe the movement of dots through transparent gels. That is unknown to us is tractions on the surface and stresses in the body.<br>
<br>
The body is composed of two materials, which two orders difference in stiffness.<br>
<br>
We constructed problem which gives the stable solution, works great with direct solvers like mumps or superlu_dist, however, we have the problem with scalability, for a larger number of processors and DOFs (>1m).<br>
<br>
The general system which we solving looks like this,<br>
[ K K K 0 ][u]=[0]<br>
[ K K K C ][u]=[0]<br>
[ K K K 0 ][u]=[0]<br>
[ 0 0 B 0 ][L]=[0]<br>
<br>
I am trying to use PCFIELDSPLIT however I have the problem which the best mix of preconditioners and set-up. I can converge for a very small problem, but overall efficiency is far from perfect.  Note that matrices C, B are small compared to the block of K matrices. Block of  K matrices is positive define and symmetric. Note that [0 0 B]*[ 0 C 0]^T = 0.<br></blockquote><div><br></div><div>Always start from a known position. Can you run full Schur complement with an exact solve for K and very accurate solve for the Schur block? Does it</div><div>converge in 1 iterate? That is the first step. Then you need to look at the two block solves. Generally, the hard part is solving the Schur block. Do you</div><div>know what would be a good preconditioner? It looks like it could be worth directly forming your Schur complement.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Anyone have any advise what I can do with this.<br>
<br>
The full code can bee seen here<br>
<a href="http://mofem.eng.gla.ac.uk/mofem/html/cell__forces_8cpp.html" rel="noreferrer" target="_blank">http://mofem.eng.gla.ac.uk/<wbr>mofem/html/cell__forces_8cpp.<wbr>html</a><br>
<br>
Kind regards<br>
<span class="HOEnZb"><font color="#888888">Lukasz</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>