[petsc-users] Multiple Sparse Matrix Vector products

[Nun ion]

Hello i have a conceptual idea of a sparse matvec implementation where i
have multiple matrices, how would i go about implementing something such as

for i = ...
   for k = ...
      w_{ik} = K_i * u_k
   end
end

Where each of the K_i are sparse matrices... the K_i are various stiffness
matrices whose size can range (although they are all the same size).  The
u_k are reused

Thanks!

Mark
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120219/bd9db64c/attachment.htm>