Yes, just assemble the same sequential matrices on each rank. To do this, create the matrix using the communicator PETSC_COMM_SELF<div><br></div><div>Cheers</div><div>  Dave<br><br>On Monday, 23 June 2014, Bogdan Dita <<a href="mailto:bogdan@lmn.pub.ro">bogdan@lmn.pub.ro</a>> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
Hello,<br>
<br>
I wanted to see how well Umfpack performs in PETSc but i encountered a<br>
problem regarding the matrix type used and the way is distributed. I'm<br>
trying to solve for multiple frequencies in parallel, where system matrix<br>
A = omega*E-C, and omega is a frequency dependent scalar. Can i send<br>
matrix E and C to each process, plus a set of frequencies to compute omega<br>
and form A_k = omega_k * E - C and then solve for every omega? Is this<br>
even possible given that Umfpack only works with SeqAij?<br>
<br>
<br>
Best regards,<br>
Bogdan<br>
<br>
<br>
----<br>
Bogdan DITA, PhD Student<br>
UPB - EE Dept.- CIEAC-LMN<br>
Splaiul Independentei 313<br>
060042 Bucharest, Romania<br>
email: <a href="javascript:;" onclick="_e(event, 'cvml', 'bogdan@lmn.pub.ro')">bogdan@lmn.pub.ro</a><br>
</blockquote></div>