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