<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">James:<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Thanks a lot for looking into this. I'm still working on debugging this on our side. It might be an issue with us. I will keep you updated.<br></blockquote><div> </div><div>Take your time. </div><div>Hong </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="im HOEnZb">
<br>
----- Original Message -----<br>
From: "Hong" <<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>><br>
</span><span class="im HOEnZb">To: "Barry Smith" <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>><br>
Cc: "James A Charles" <<a href="mailto:charlesj@purdue.edu">charlesj@purdue.edu</a>>, "Hong" <<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>>, <a href="mailto:dlemus@purdue.edu">dlemus@purdue.edu</a>, "PETSc users list" <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>>, "Daniel Mejia" <<a href="mailto:denphi@denphi.com">denphi@denphi.com</a>>, "Tillmann Kubis" <<a href="mailto:tkubis@purdue.edu">tkubis@purdue.edu</a>><br>
Sent: Sunday, April 5, 2015 6:13:42 PM<br>
Subject: Re: [petsc-users] MatMatMult with dense matrices.<br>
<br>
<br>
</span><span class="im HOEnZb">James,<br>
<br>
<br>
My understanding is that<br>
you call MatMatMult in sequential for computing C = A*B, in which<br>
A is in aij format, and B is dense, and the product C is dense, but gives different values than what you got from Matlab.<br>
<br>
<br>
Do you have a stand-alone petsc code to reproduce this so I can investigate what goes wrong?<br>
<br>
<br>
Hong<br>
<br>
<br>
<br>
On Sun, Apr 5, 2015 at 3:44 PM, Barry Smith < <a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a> > wrote:<br>
<br>
<br>
<br>
We would need to see the PETSc side of the code to see if there is anything wrong there.<br>
<br>
<br>
<br>
> On Apr 5, 2015, at 3:35 PM, James A Charles < <a href="mailto:charlesj@purdue.edu">charlesj@purdue.edu</a> > wrote:<br>
><br>
> Hi Hong,<br>
><br>
> You can open up matmatmult_comparison.mat in MATLAB ( load matmatmult_comparison.mat ). That .mat contains A, B, C, and C_MATLAB. I called matmatmult in sequential. We use AIJ when A2 was sparse but at the time of this multiplication everything was dense.<br>
><br>
> Thanks,<br>
> James<br>
><br>
><br>
> ----- Original Message -----<br>
> From: "Hong" < <a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a> ><br>
> To: "James A Charles" < <a href="mailto:charlesj@purdue.edu">charlesj@purdue.edu</a> ><br>
</span><div class="HOEnZb"><div class="h5">> Cc: "PETSc users list" < <a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a> >, <a href="mailto:dlemus@purdue.edu">dlemus@purdue.edu</a> , "Daniel Mejia" < <a href="mailto:denphi@denphi.com">denphi@denphi.com</a> >, "Tillmann Kubis" < <a href="mailto:tkubis@purdue.edu">tkubis@purdue.edu</a> ><br>
> Sent: Sunday, April 5, 2015 4:25:39 PM<br>
> Subject: Re: [petsc-users] MatMatMult with dense matrices.<br>
><br>
><br>
><br>
><br>
> James:<br>
><br>
><br>
> I'm multiplying two dense matrices. I compared the results against MATLAB and I am getting different results. Attached you will see a zip that contains a .mat file that has the A, B, and C. C_MATLAB is the multiplication of A*B in MATLAB and C is the PETSc 3.4.3 result. I also attached the .m files A, B, and C that were saved from PETSc. Is there something I'm doing wrong or<br>
><br>
><br>
> What is the file type of matmatmult_comparison.mat?<br>
> How to read it?<br>
> I need to see how it defers from C.m.<br>
><br>
><br>
><br>
> any reason due to the shape of the matrices that this would fail? If it gives any more information we convert the previous matrix that A is formed of A2 (A = A1*A2) to dense prior to the multiplication using MatConvert.<br>
><br>
><br>
><br>
> It seems both A and B are dense, complex square matrices.<br>
> Did you call MatMatMult() in sequential or parallel? What matrix format did you use?<br>
><br>
><br>
> Hong<br>
<br>
<br>
</div></div></blockquote></div><br></div></div>