<div dir="ltr">On Mon, Nov 11, 2013 at 5:33 AM, Stephan Kramer <span dir="ltr"><<a href="mailto:s.kramer@imperial.ac.uk" target="_blank">s.kramer@imperial.ac.uk</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">On 01/11/13 06:01, Matthew Knepley wrote:<div class="im"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
On Wed, Oct 30, 2013 at 9:19 PM, Rhodri Davies <<a href="mailto:rhodri.davies@anu.edu.au" target="_blank">rhodri.davies@anu.edu.au</a> <mailto:<a href="mailto:rhodri.davies@anu.edu.au" target="_blank">rhodri.davies@anu.edu.<u></u>au</a>>> wrote:<br>
<br>
Dear all,<br>
<br>
Would you mind giving us an update on your thoughts re: this question? We're very keen to be able to use such functionally and it would be good to know if you're looking into this / have any<br>
suggestions?<br>
<br>
<br>
I looked at the code. It does not look hard to make a MPIBAIJ version. Since the scatter context for BAIJ is already exploded, most<br>
of the MPIAIJ code goes straight through, it just a matter of handling the blocking in a few places.<br>
<br>
Here is what would be really helpful for us. Write a nice test. The code would probably only take an hour of being careful. However,<br>
the test assembly would be more. If you give us a nice test example for it, writing it could happen pretty quickly.<br>
<br>
Thanks,<br>
<br>
Matt<br>
</blockquote>
<br></div>
That would be great. I had a look at the existing tests and found src/mat/examples/tests/ex140.c which already claims to test MatZeroRowsColumns with MPIBAIJ - except it isn't being run in parallel. Would running this test in parallel be a sufficient test to work with? It currently fails with the expected message "No support for this operation for this object type!". I'd be happy to write another test but it would basically be creating another matrix and testing it does the right thing with the matrix and the rhs. Let me know if there's anything else we could help out with<br>
</blockquote><div><br></div><div>As far as I can tell ex140 does not really work. I am using ex12 to check my rewrite of MatZeroRows(). I would like</div><div>a companion to that which uses RowsColumns(). And the work here is to provide the output to check against. Even</div>
<div>for ex140 I would have to do that, which means more time here.</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">
Cheers<span class="HOEnZb"><font color="#888888"><br>
Stephan</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
From: Stephan Kramer <<a href="mailto:s.kramer@imperial.ac.uk" target="_blank">s.kramer@imperial.ac.uk</a>><br>
Subject: Fwd: [petsc-users] MatZeroRowsColumns for mpibaij<br>
Date: 31 October 2013 13:16:05 GMT+11:00<br>
To: Rhodri Davies <<a href="mailto:rhodri.davies@anu.edu.au" target="_blank">rhodri.davies@anu.edu.au</a>><br>
<br>
<br>
<br>
<br>
-------- Original Message --------<br>
Subject: [petsc-users] MatZeroRowsColumns for mpibaij<br>
Date: Wed, 19 Jun 2013 13:22:31 +0100<br>
From: Stephan Kramer <<a href="mailto:s.kramer@imperial.ac.uk" target="_blank">s.kramer@imperial.ac.uk</a>><br>
To: PETSc users list <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<br>
Dear all,<br>
<br>
We have found the MatZeroRowsColumns() routine to be very useful for lifting boundary conditions. Unfortunately we also found that it isn't implemented for matrix type mpibaij and we're quite keen on<br>
keeping the baij structure as it seems to give us significant savings for DG problems (in particular the blocked sor). I had a quick look at the implementation of MatZeroRowsColums_MPIAIJ but couldn't<br>
really estimate how hard it would be to do something similar for mpibaij (an implementation for seqbaij exists already). Is this something that's likely to be looked at in the future?<br>
<br>
Stephan<br>
</blockquote></blockquote>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>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>