<div dir="ltr"><div><div><div>Thank you for the quick answer.<br><br></div>I didn't realize that I could use PC without KSP interface. I also think that it is what I wanted. Nevertheless, as long as I figured out from the source code, PcApply for PCSOR basically do some interface and preparations and then calls MatSOR. I saw that depending on the matrix ('BAIJ, SBAIJ, and AIJ matrices with Inodes') it does SOR smoothing or block SOR smoothing. <br><br>I think that in my case the seqaij matrix with bs=4 had Inodes with size 4. That is why calling MatSOR with seqaij or calling with seqbaij converted from the seqaij seem to give the same result.<br>However, with the matrix seqaij with bs = 16 I can guess that the rows inside a block dont have the same nonzero pattern, so Inodes size are different from block size. I happened to see the follow note in the MatSOR website page: "Developer Note: We should add block SOR support for AIJ matrices with block size set to great than one and no inodes
". This may be the reason why seqaij and seqbaij are leading to different results with my matrix of bs = 16. I think that answer all may previous questions. I am sorry, I've got confused and wrote MatSolve instead of MatSOR in my previous email, what changes it completely.<br><br></div>Best Regards<br><br></div>Eduardo<br><div><div><br><br><br><br><div class="gmail_extra"><br><div class="gmail_quote">2016-07-20 0:03 GMT-03:00 Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="">On Tue, Jul 19, 2016 at 8:17 PM, Eduardo Jourdan <span dir="ltr"><<a href="mailto:eduardojourdan92@gmail.com" target="_blank">eduardojourdan92@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div>Hi all,<br><br></div>I would like to perform a specific number (for instance 4 of forward and backward sweeps with a seqaij matrix with block size 4, vectors b and x. Also, I need to do this same procedure with another matrix seqaij block size 16. I would appreciate if someone knows the best way to do it. <br></div></div></div></div></div></blockquote><div><br></div></span><div>It sounds like you want PCSOR and PCApply, not MatSolve.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><span class=""><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>1 - I've been trying to use MatSolve. For the bs=4 it seems to work, but with the other matrix with bs=16 the residue diverges. When I call matConvert to convert the later matrix for a seqbaij with bs=16 the result changes and the linear residue is reduced. It is supposed to happen or it is more possible that i am doing something wrong? <br><br></div>2 - MatSolve for seqbaij and seqaij with the same block sizes gives the same results in terms of solution (not performace, memory) ?<br><br></div>3 - Can do I do a specific number of sweeps as told before with the KSP/PC interface?<br><br></div><div>4 - I saw the manual for the MatSolve and It says that it is for factored matrix. Can I use a matrix just after the MatAssembly calls?<br></div><div><br></div>Best regards,<br><br></div>Eduardo Jourdan<br></div>
</blockquote></span></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div 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>
</font></span></div></div>
</blockquote></div><br></div></div></div></div>