Multiply SparseMatrix and a Vector

[Nitish Anand]

Hello,

I wish to multiply a PETSeWrapper::MPI::SparseMatrix with a Vector<double>. I would be glad if someone can point to the right material to do this efficiently.

Thanks in advance.

--Nitish

[Nitish Anand]

Errata: PETSc 

[Jean-Paul Pelteret]

Dear Nitish,

Well, in the documentation to vmult it states that when

the current object represents a parallel distributed matrix (of type PETScWrappers::MPI::SparseMatrix), then both vectors have to be distributed vectors as well

. So I suspect that what you need to do is to create a pair of distributed vectors and then simply copy the data from the deal.II vector to the PETSc one using operator=.

Does this help?

Regards,

Jean-Paul

[Nitish Anand]

In fact, that was the issue.. Works well now. :) Thanks !