MatMatMult_MPIDense_MPIDense() works currently?

[David Fuentes]

also, MatMatMult_MPIAIJ_MPIDense() and
MatMatMult_MPIDense_MPIAIJ()
seem to work if you don't want/need to convert both matrices.




df





On Sun, 5 Apr 2009, David Fuentes wrote:

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