<div dir="auto">If it's not specialized, as fast as one matrix vector multiplication (with all ones)</div><br><div class="gmail_quote"><div dir="ltr">Il Ven 27 Apr 2018, 19:58 Derek Gaston <<a href="mailto:friedmud@gmail.com" target="_blank" rel="noreferrer">friedmud@gmail.com</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">I was wondering about the comment on MatGetRowSum that says:<div><br></div><div>"This code is slow since it is not currently specialized for different formats"</div><div><br></div><div>How slow are we talking about here?  Is it actually going to loop over i,j from 0 to N on a sparse matrix?</div><div><br></div><div>I need to get a vector that represents the "lumped" mass matrix: but I would like to get it fairly quickly.  I would just compute the lumped matrix directly and then use MatGetDiag - but I actually want both (unlumped/consistent and lumped) and MatGetRowSum seemed like the best direction.  Willing to take advice though!</div><div><br></div><div>Derek</div></div>
</blockquote></div>