MatMatMult_MPIDense_MPIDense() works currently?

Yujie recrusader at
Fri Apr 3 13:02:33 CDT 2009

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

The errors I got is "Caught signal number 11 SEGV: Segmentation Violation".
It takes place in  "PLA_Obj_set_to_zero(lu->A);" in


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.



On Thu, Apr 2, 2009 at 8:02 PM, Barry Smith <bsmith at> 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: <>

More information about the petsc-users mailing list