And use the same sequential communicator when you create the KSP on each rank
On Monday, 23 June 2014, Dave May <dave.mayhe...@gmail.com> wrote: > 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 > <javascript:_e(%7B%7D,'cvml','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 >> >