[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