[petsc-users] The product of two MPIDENSE matrices

[S V N Vishwanathan]

Hi

> > I am wondering whether I can multiply MATELEMENTAL and MATAIJ
> > matrices. If I cannot multiply those, could you build the
> > multiplication?
> 
> What are the dimensions of the matrices?
> 
> Jack, does Clique support multiplication of sparse matrices with
> [MC,MR] dense matrices?

The dense matrix is very flat and long (e.g. 10 x few 1000s) and the
sparse matrix is also flat and very very long (e.g. few 1000s x few
millions). I saw the code which is used to multiply a MPIAIJ matrix with
a dense matrix. We need something very similar, but which uses elemental
matrices. 

There seems to be some inconsistency right now in PETSc. I can multiple
a MPIAIJ matrix with MPIDENSE, but I cannot multiply MPIDENSE with
MPIDENSE since there is no PLAPACK support anymore. On other hand, I can
multiply a MATELEMENTAL with MATELEMENTAL but I cannot multiply a
MPIDENSE with MATELEMENTAL. Is the eventual plan to remove MPIDENSE
completely and replace it with MATELEMENTAL?

thanks

vishy