> On Oct 15, 2015, at 8:40 AM, Matthew Knepley <[email protected]> wrote:
> 
> On Thu, Oct 15, 2015 at 8:20 AM, Bikash Kanungo <[email protected]> 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

Reply via email to