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