[petsc-users] MatMatMul with mpidense

[Matthew Knepley]

On Sat, Jul 23, 2011 at 9:59 PM, Kyle Wendt <kylewendt at gmail.com> wrote:

> Thanks, I'll look into what going wrong over the next few days.


I wrote the initial PLAPACK implementation, and it is very simple. If you
are truly going to spend
time on the code, your time is better spent replacing PLAPACK with its
successor, Elemental

  http://code.google.com/p/elemental/

than with debugging old PLAPACK code. It will take the same amount of time
(or less), perform
better, and provide every other user a benefit.

  Thanks,

     Matt


> On Sat, Jul 23, 2011 at 5:18 PM, Satish Balay <balay at mcs.anl.gov> wrote:
> > 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?
> >> >>>>
> >> >>>
> >> >>>
> >> >
> >> >
> >>
> >
>



-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110723/76e4f536/attachment-0001.htm>