[petsc-users] MatMatMul with mpidense

Satish Balay balay at mcs.anl.gov
Sat Jul 23 16:18:15 CDT 2011


the error is from petsc code saying this codepath to plapack is buggy.

If you wish to debug plapack - you can comment out the following line
in src/mat/impls/dense/mpi/mpidense.c

  SETERRQ(PETSC_ERR_LIB,"Due to apparent bugs in PLAPACK,this is not currently supported");

and then figureout if/why you are getting wrong numbers out of this call..

Satish

On Sat, 23 Jul 2011, Kyle Wendt wrote:

> Ok, is there a way to get more detailed information as to what went wrong?
> 
> On Sat, Jul 23, 2011 at 4:09 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
> >
> > On Jul 23, 2011, at 2:12 PM, Kyle Wendt wrote:
> >
> >> I also tried that and when i run my code i get this error:
> >> [0]PETSC ERROR: --------------------- Error Message
> >> ------------------------------------
> >> [0]PETSC ERROR: Error in external library!
> >> [0]PETSC ERROR: Due to apparent bugs in PLAPACK,this is not currently supported!
> >
> >   The error message says it all. plapack is buggy and not supported. Currently there no good options for general purpose parallel dense matrix operations in PETSc. Sorry, our focus is on sparse matrix computations because we can't do everything so we do what we do best.
> >
> >   Barry
> >
> >>
> >> On Sat, Jul 23, 2011 at 2:44 PM, Satish Balay <balay at mcs.anl.gov> wrote:
> >>> you need --download-plapack for this.
> >>>
> >>> [petsc does not use scalapack. Its a dependency of mumps - so you need
> >>> it only if using --download-mumps]
> >>>
> >>> Satish
> >>>
> >>> On Sat, 23 Jul 2011, Kyle Wendt wrote:
> >>>
> >>>> I am trying to multiple two mpidense matrices with matmatmil and I get
> >>>> the following runtime error:
> >>>> [0]PETSC ERROR: --------------------- Error Message
> >>>> ------------------------------------
> >>>> [0]PETSC ERROR: No support for this operation for this object type!
> >>>> [0]PETSC ERROR: MatMatMult not supported for B of type mpidense!
> >>>> [0]PETSC ERROR:
> >>>> ------------------------------------------------------------------------
> >>>> [0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 8, Thu Mar 17
> >>>> 13:37:48 CDT 2011
> >>>> [0]PETSC ERROR: See docs/changes/index.html for recent updates.
> >>>> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> >>>>
> >>>> I have built and linked petsc again scalapack and blacs but this have
> >>>> not helped, is this operation really not supported?
> >>>>
> >>>
> >>>
> >
> >
> 


More information about the petsc-users mailing list