MatMatMult_MPIDense_MPIDense() works currently?

[David Fuentes]

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
> 
> 
> 
> 
>