[petsc-users] MatPtAP for involving MPIDENSE and MPIAIJ matrices

[Barry Smith]

> On Oct 15, 2015, at 8:40 AM, Matthew Knepley <knepley at gmail.com> wrote:
> 
> On Thu, Oct 15, 2015 at 8:20 AM, Bikash Kanungo <bikash at umich.edu> wrote:
> HI,
> 
> I'm trying to evaluate B = P^T * A * P, where P is an MPIDENSE matrix and A is an MPIAIJ matrix. Neither MatTransposeMatMult nor MatPtAP support such an operation for the given type of matrices. I tried to look into MATELEMENTAL class and didn't find any MatTransposeMatMult function for a pair of Elemental matrices. 
> 
> What would be a possible and efficient way to perform such the above operation?

   Hong will see about providing the missing piece.  What are you sizes of the matrix P? Presumably it has many rows but few columns?

  Barry

> 
> Jed, is TAIJ the way to do this?
> 
>   Matt
>  
> Thanks,
> Bikash 
> 
> -- 
> Bikash S. Kanungo
> PhD Student
> Computational Materials Physics Group
> Mechanical Engineering 
> University of Michigan
> 
> 
> 
> 
> -- 
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
> -- Norbert Wiener