<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 14, 2018 at 4:35 PM, Smith, Barry F. <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><br>
  What are you doing with the matrix?<br></blockquote><div><br></div><div>We are doing an explicit method. PDEs are discretized using a finite element method, so there is a mass matrix. The mass matrix will be lumped, and it becomes diagonal. We want to compute the inverse of the lumped matrix, and also do a few of matrix-vector multiplications using  the lumped matrix or  its inverse.</div><div><br></div><div>The specific implementation won't make this more efficient?   </div><div><br></div><div><br></div><div>Fande,</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">
<br>
  We don't have a diagonal matrix but it would be easy to add such a beast if it was performance critical, which it probably isn't.<br>
<span class="gmail-HOEnZb"><font color="#888888"><br>
  Barry<br>
</font></span><div class="gmail-HOEnZb"><div class="gmail-h5"><br>
<br>
<br>
<br>
> On Feb 14, 2018, at 3:57 PM, Fande Kong <<a href="mailto:fdkong.jd@gmail.com">fdkong.jd@gmail.com</a>> wrote:<br>
><br>
> Hi All,<br>
><br>
> If  a matrix is always diagonal, what a good way to represent the matrix?  Still MPIAIJ, MPIBAIJ? Can we have a specific implementation for this?<br>
><br>
><br>
> Fande,<br>
<br>
</div></div></blockquote></div><br></div></div>