[petsc-users] Reshaping a vector into a matrix

S V N Vishwanathan vishy at stat.purdue.edu
Thu Jan 6 17:24:51 CST 2011


> >> 2. Since each vec1,vec2,... is the same size, this is likely to
> >> produce poor memory performance.  If you care about speed, I suggest
> >> interlacing the values in vec1,...,veck.  In that case, you can create
> >> an MAIJ matrix that acts on this "multi-vector".

>    You most definitely want to use the MAIJ.  MAIJ does not "repeat" X, it uses the 
> original matrix passed in but does efficient multiple matrix-vector products at the same time.

Is there any place to read up about MAIJ matrices. I read the manual
pages but they were rather cryptic. Perhaps some example code?

vishy





More information about the petsc-users mailing list