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 >> >> > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090403/dca3d74d/attachment.htm>
