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