[petsc-users] MatMatMult with dense matrices.

[James A Charles]

Hello,

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 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. 

Thanks for your help,
James Charles
-------------- next part --------------
A non-text attachment was scrubbed...
Name: matmatmult_comparison.zip
Type: application/zip
Size: 1804789 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150404/49c8dcce/attachment-0001.zip>