Thank you Matt,<br><br>Any suggestion to solve the problem I am trying to tackle. I want to solve a linear system: <br><br>Sum(A_i) u= Sum(f_i) , i=1.... to No. of CPUs. <br><br>Where&nbsp; A_i is a sparse sequential matrix and f_i is a sequential vector. Each CPU has one matrix and one vector of the same size. Now I want to sum up and solve the system in parallel.<br><br>Thanks again<br><br>Waad<br><br><b><i>Matthew Knepley &lt;knepley@gmail.com&gt;</i></b> wrote:<blockquote class="replbq" style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"> On Tue, May 20, 2008 at 2:12 PM, Waad Subber <w_subber@yahoo.com> wrote:<br>&gt; Hi,<br>&gt;<br>&gt; I am trying to construct a sparse parallel matrix (MPIAIJ) by adding up<br>&gt; sparse sequential matrices (SeqAIJ) from each CPU. I am using<br>&gt;<br>&gt; MatMerge_SeqsToMPI(MPI_Comm comm,Mat seqmat,PetscInt m,PetscInt n,MatReuse<br>&gt; scall,Mat *mpimat)<br>&gt;<br>&gt; to do that. However, when I
 compile the code I get the following<br>&gt;<br>&gt;  undefined reference to `matmerge_seqstompi_'<br>&gt; collect2: ld returned 1 exit status<br>&gt; make: *** [all] Error 1<br>&gt;<br>&gt; Am I using this function correctly ?<br><br>These have no Fortran bindings right now.<br><br>  Matt<br><br>&gt; Thanks<br>&gt;<br>&gt; Waad<br>&gt;<br><br><br><br>-- <br>What most experimenters take for granted before they begin their<br>experiments is infinitely more interesting than any results to which<br>their experiments lead.<br>-- Norbert Wiener<br><br></w_subber@yahoo.com></blockquote><br><p>&#32;