MatMatMult_MPIDense_MPIDense() works currently?

[Barry Smith]

   Yujie

     You are are on your own on this. I spent (wasted) many hours  
trying to debug the PLAPACK problem, I'm not going to deal with it  
again.
You'll have to work through the code yourself.

    Barry

On Apr 3, 2009, at 1:02 PM, 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
>
>
>