<div dir="ltr">On Mon, Mar 18, 2013 at 2:34 PM, Mengda Wu <span dir="ltr"><<a href="mailto:wumengda@gmail.com" target="_blank">wumengda@gmail.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi all,<br><br>   I would like to solve a system Ax=b multiple times, with A and b stay the same from time to time. The difference among each time is that a set of DOFs in x are fixed and others are free to be solved. And the set of DOFs to be fixed/ solved varies from time to time. I know I can assemble A each time and move the set of fixed DOFs to the rhs. But I don't want to do this because the size of A is changing. <br>

<br>   Is there a simpler and more automated way to deal with this issue? <br></blockquote><div><br></div><div style>You should assemble a separate system for each "active" set. This is the simplest and fastest method. You can create the system using MatGetSubMatrix() if you don't want to call your assembly routines each time.</div>
</div></div></div>