about MatMat*() functions

Dave May dave.mayhem23 at gmail.com
Sun Feb 24 16:00:19 CST 2008


Hey,
    MatMatMult() will work for MPIAIJ matrices. So will MatPtAP(). If you
are
ever in doubt, the easiest way (I find) to check whether I certain operation
is
supported is to just look at the source and see which ops. are defined.

I can usually find the answer with the online docs. In this case, starting
with the
type (MATMPIAIJ), and then searching mpiaij.c for MatMatMult. If you see
something like
your desired op. in  struct _MatOps (i.e. MatMatMult_MPIAIJ_MPIAIJ) then
it's supported.
Now you know the function defining your operation and you can search for it
(with grep as it might be in another file not online) to find out exactly
what it does.


Cheers,
    Dave





On Mon, Feb 25, 2008 at 7:07 AM, Yujie <recrusader at gmail.com> wrote:

> hi,
>
> I am wondering whether all the MatMat*() only are suitable for sequential
> matrix. I know MatMatSolve() is for sequential matrix. How about
> MatMatMult()? Thanks a lot.
>
> Regards,
> Yujie
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080225/f5440a5e/attachment.htm>


More information about the petsc-users mailing list