<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Feb 14, 2018 at 7:11 PM, Fande Kong <span dir="ltr"><<a href="mailto:fdkong.jd@gmail.com" target="_blank">fdkong.jd@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 class="gmail_extra"><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></div></div></blockquote><div><br></div><div>I am doing this for Pylith. I think you should just do Vec operations, and pull the inverse mass matrix to the rhs.</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 dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><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="m_-4001369547178759607gmail-HOEnZb"><font color="#888888"><br>
  Barry<br>
</font></span><div class="m_-4001369547178759607gmail-HOEnZb"><div class="m_-4001369547178759607gmail-h5"><br>
<br>
<br>
<br>
> On Feb 14, 2018, at 3:57 PM, Fande Kong <<a href="mailto:fdkong.jd@gmail.com" target="_blank">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>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>