And use the same sequential communicator when you create the KSP on each rank<br><br>On Monday, 23 June 2014, Dave May <<a href="mailto:dave.mayhem23@gmail.com">dave.mayhem23@gmail.com</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
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="javascript:_e(%7B%7D,'cvml','bogdan@lmn.pub.ro');" target="_blank">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>bogdan@lmn.pub.ro</a><br>
</blockquote></div>
</blockquote>