MatMatMult_MPIDense_MPIDense() works currently?

[Yujie]

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>