> >> 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
