[Trilinos-Users] Re: multiply two Epetra_CrsMatrix and get back the diagonal

[Kurt Stokbro]

I need a subroutine that can
Epetra_Vector multiplyDiag(Epetra_CrsMatrix& a, Epetra_CrsMatrix& b)
 {
   for(int j=0;j<a.numRows;j++)
   {
         v(j)=0.0
        for(int i=0;i<a.numCols;i++)
                 v(j)+=a(j,i)*b(i,j)
  }
   return v
}

Does anybody have an efficient routine for this?
Thanks,
Kurt Stokbro