1 2 /* 3 Support for the parallel AIJ matrix vector multiply 4 */ 5 #include "mpiaij.h" 6 #include "vec/vecimpl.h" 7 #include "../seq/aij.h" 8 9 int MPIAIJSetUpMultiply(Mat mat) 10 { 11 Matimpiaij *aij = (Matimpiaij *) mat->data; 12 Matiaij *B = (Matiaij *) (aij->B->data); 13 int N = aij->N,i,j,*indices; 14 15 /* For the first stab we make an array as long as the number of columns */ 16 indices = (int *) MALLOC( N*sizeof(int) ); CHKPTR(indices); 17 MEMSET(indices,0,N*sizeof(int)); 18 19 return 0; 20 } 21