[petsc-users] How to efficiently represent a diagonal matrix?

[Smith, Barry F.]

  What are you doing with the matrix?

  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.

  Barry




> On Feb 14, 2018, at 3:57 PM, Fande Kong <fdkong.jd at gmail.com> wrote:
> 
> Hi All,
> 
> 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?
> 
> 
> Fande,