On Sat, Jul 23, 2011 at 9:59 PM, Kyle Wendt <span dir="ltr"><<a href="mailto:kylewendt@gmail.com">kylewendt@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Thanks, I'll look into what going wrong over the next few days.</blockquote><div><br></div><div>I wrote the initial PLAPACK implementation, and it is very simple. If you are truly going to spend</div><div>time on the code, your time is better spent replacing PLAPACK with its successor, Elemental</div>
<div><br></div><div> <a href="http://code.google.com/p/elemental/">http://code.google.com/p/elemental/</a></div><div><br></div><div>than with debugging old PLAPACK code. It will take the same amount of time (or less), perform</div>
<div>better, and provide every other user a benefit.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div class="h5">
On Sat, Jul 23, 2011 at 5:18 PM, Satish Balay <<a href="mailto:balay@mcs.anl.gov">balay@mcs.anl.gov</a>> wrote:<br>
> the error is from petsc code saying this codepath to plapack is buggy.<br>
><br>
> If you wish to debug plapack - you can comment out the following line<br>
> in src/mat/impls/dense/mpi/mpidense.c<br>
><br>
> SETERRQ(PETSC_ERR_LIB,"Due to apparent bugs in PLAPACK,this is not currently supported");<br>
><br>
> and then figureout if/why you are getting wrong numbers out of this call..<br>
><br>
> Satish<br>
><br>
> On Sat, 23 Jul 2011, Kyle Wendt wrote:<br>
><br>
>> Ok, is there a way to get more detailed information as to what went wrong?<br>
>><br>
>> On Sat, Jul 23, 2011 at 4:09 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br>
>> ><br>
>> > On Jul 23, 2011, at 2:12 PM, Kyle Wendt wrote:<br>
>> ><br>
>> >> I also tried that and when i run my code i get this error:<br>
>> >> [0]PETSC ERROR: --------------------- Error Message<br>
>> >> ------------------------------------<br>
>> >> [0]PETSC ERROR: Error in external library!<br>
>> >> [0]PETSC ERROR: Due to apparent bugs in PLAPACK,this is not currently supported!<br>
>> ><br>
>> > 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.<br>
>> ><br>
>> > Barry<br>
>> ><br>
>> >><br>
>> >> On Sat, Jul 23, 2011 at 2:44 PM, Satish Balay <<a href="mailto:balay@mcs.anl.gov">balay@mcs.anl.gov</a>> wrote:<br>
>> >>> you need --download-plapack for this.<br>
>> >>><br>
>> >>> [petsc does not use scalapack. Its a dependency of mumps - so you need<br>
>> >>> it only if using --download-mumps]<br>
>> >>><br>
>> >>> Satish<br>
>> >>><br>
>> >>> On Sat, 23 Jul 2011, Kyle Wendt wrote:<br>
>> >>><br>
>> >>>> I am trying to multiple two mpidense matrices with matmatmil and I get<br>
>> >>>> the following runtime error:<br>
>> >>>> [0]PETSC ERROR: --------------------- Error Message<br>
>> >>>> ------------------------------------<br>
>> >>>> [0]PETSC ERROR: No support for this operation for this object type!<br>
>> >>>> [0]PETSC ERROR: MatMatMult not supported for B of type mpidense!<br>
>> >>>> [0]PETSC ERROR:<br>
>> >>>> ------------------------------------------------------------------------<br>
>> >>>> [0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 8, Thu Mar 17<br>
>> >>>> 13:37:48 CDT 2011<br>
>> >>>> [0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
>> >>>> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
>> >>>><br>
>> >>>> I have built and linked petsc again scalapack and blacs but this have<br>
>> >>>> not helped, is this operation really not supported?<br>
>> >>>><br>
>> >>><br>
>> >>><br>
>> ><br>
>> ><br>
>><br>
><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>