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

[Matthew Knepley]

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

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151015/257353e4/attachment.html>