<div dir="ltr">On Wed, Oct 30, 2013 at 9:19 PM, Rhodri Davies <span dir="ltr"><<a href="mailto:rhodri.davies@anu.edu.au" target="_blank">rhodri.davies@anu.edu.au</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">



<div style="word-wrap:break-word">
Dear all,
<div><br>
</div>
<div>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 suggestions?</div>
</div></blockquote><div><br></div><div>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</div><div>of the MPIAIJ code goes straight through, it just a matter of handling the blocking in a few places.</div>
<div><br></div><div>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,</div><div>the test assembly would be more. If you give us a nice test example for it, writing it could happen pretty quickly.</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"><div style="word-wrap:break-word">
<div>Best wishes,</div>
<div><br>
</div>
<div>Rhodri</div>
<div><br>
<div>
<div>---------------------------------------------------------------------------------------------------<br>
Research Fellow<br>
Research School of Earth Sciences<br>
The Australian National University<br>
Canberra, ACT, 0200, Australia<br>
<br>
T: <a href="tel:%2B61%20%280%292%206125%203643" value="+61261253643" target="_blank">+61 (0)2 6125 3643</a> <br>
E: <a href="mailto:rhodri.davies@anu.edu.au" target="_blank">rhodri.davies@anu.edu.au</a><br>
W: <a href="http://rses.anu.edu.au/people/dr-rhodri-davies" target="_blank">http://rses.anu.edu.au/people/dr-rhodri-davies</a> <br>
---------------------------------------------------------------------------------------------------<br>
<br>
<br>
</div>
</div>
<div><br>
<div>Begin forwarded message:</div>
<br>
<blockquote type="cite">
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px">
<span style="font-family:'Helvetica';font-size:medium;color:rgba(0,0,0,1.0)"><b>From:
</b></span><span style="font-family:'Helvetica';font-size:medium">Stephan Kramer <<a href="mailto:s.kramer@imperial.ac.uk" target="_blank">s.kramer@imperial.ac.uk</a>><br>
</span></div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px">
<span style="font-family:'Helvetica';font-size:medium;color:rgba(0,0,0,1.0)"><b>Subject:
</b></span><span style="font-family:'Helvetica';font-size:medium"><b>Fwd: [petsc-users] MatZeroRowsColumns for mpibaij</b><br>
</span></div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px">
<span style="font-family:'Helvetica';font-size:medium;color:rgba(0,0,0,1.0)"><b>Date:
</b></span><span style="font-family:'Helvetica';font-size:medium">31 October 2013 13:16:05 GMT+11:00<br>
</span></div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px">
<span style="font-family:'Helvetica';font-size:medium;color:rgba(0,0,0,1.0)"><b>To:
</b></span><span style="font-family:'Helvetica';font-size:medium">Rhodri Davies <<a href="mailto:rhodri.davies@anu.edu.au" target="_blank">rhodri.davies@anu.edu.au</a>><br>
</span></div>
<br>
<div><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>
<br>
<br>
<br>
<br>
<br>
</div>
</blockquote>
</div>
<br>
</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>