<p>Dear Barry:</p><p>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(). <br>
</p><p>The errors I got is "Caught signal number 11 SEGV: Segmentation Violation". It takes place in "PLA_Obj_set_to_zero(lu->A);" in </p><p>MatMPIDenseCopyToPlapack(). </p><p>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.</p>
<p>Regards,</p><p>Yujie</p><br><div class="gmail_quote">On Thu, Apr 2, 2009 at 8:02 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><br>
On Apr 2, 2009, at 9:19 PM, Yujie wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi, PETSc Developers<br>
<br>
I am wondering whether MatMatMult_MPIDense_MPIDense() works currently based on PLAPACK? Thanks a lot.<br>
<br>
</blockquote>
<br></div>
No, if you run it you will see it print an error message.<br>
<br>
I tried to debug PLAPACK to determine the problem but it was awfully complicated and had to give up. Certainly someone else<br>
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.<br>
<br>
Barry<br>
<br>
<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Regards,<br>
<br>
Yujie<br>
<br>
</blockquote>
<br>
</blockquote></div><br>