Hi Yujie,
as a work around have you tried converting your dense matrices to aij format and using MatMatMult_MPIAIJ_MPIAIJ()?? df On Fri, 3 Apr 2009, Yujie wrote: > > Dear Barry: > > I am trying to debug the codes you have written with ex123.c. After > commenting the error output > (SETERRQ(PETSC_ERR_LIB,"Due to aparent bugs in PLAPACK,this is not currently > supported");) in > MatMatMultSymbolic_MPIDense_MPIDense().? > > The errors I got is "Caught signal number 11 SEGV: Segmentation Violation". > It takes place in > ?"PLA_Obj_set_to_zero(lu->A);" in? > > MatMPIDenseCopyToPlapack().? > > To my understanding, if you want to set "lu->A" to zero, you first assign > memory to "lu->A". However, I can't find > which function you do this in? Could you give me some advice? thanks a lot. > > Regards, > > Yujie > > > On Thu, Apr 2, 2009 at 8:02 PM, Barry Smith <bsmith at mcs.anl.gov> wrote: > > On Apr 2, 2009, at 9:19 PM, Yujie wrote: > > Hi, PETSc Developers > > I am wondering whether MatMatMult_MPIDense_MPIDense() works > currently based on PLAPACK? > Thanks a lot. > > > ?No, if you run it you will see it print an error message. > > ?I tried to debug PLAPACK to determine the problem but it was awfully > complicated and had to give up. > Certainly someone else > could try to debug PLAPACK to determine the problem. PLAPACK is not supported > so unfortunately there is no one > to complain to about it and you'd have to fix it yourself. > > ?Barry > > > > Regards, > > Yujie > > > > >
