[petsc-users] The product of two MPIDENSE matrices

[Joon Hee Choi]

Hello Hong,

I tried to get Elemental but it does not support 64-bit integer now. I need to use 64-bit integer because the matrix I am using is very large. I need just MatMatMult() for MPIDENSE matrices, but it also didn't support MPIDENSE. Is there any method to multiply MPIDENSE and MPIDENSE? Thank you.

Joon

----- Original Message -----
From: "Hong Zhang" <hzhang at mcs.anl.gov>
To: "Joon Hee Choi" <choi240 at purdue.edu>
Cc: "PETSc users list" <petsc-users at mcs.anl.gov>
Sent: Saturday, September 14, 2013 9:18:44 PM
Subject: Re: [petsc-users] The product of two MPIDENSE matrices


Petsc does not support C=A x B^T for MPIDENSE_ MPIDENSE matrices, but supports it via elemental interface, 
see http://code.google.com/p/elemental/ and 
~petsc//src/mat/impls/elemental/matelem.cxx. 


Hong 



On Sat, Sep 14, 2013 at 12:47 AM, Joon Hee Choi < choi240 at purdue.edu > wrote: 


Hello, 

I want to calculate M = C x C^T. C is a MPIDENSE matrix. 
I tried the calculation with MatTransposeMatMult (C, C, MAT_INITIAL_MATRIX, PETSC_DEFAULT, &M), but it did not support MPIDENSE. 
Is there any method to calculate the product of two MPIDENSE matrices? 

Thanks, 
Joon