For what it's worth, I've attached the results from np=1 for MatMatTransposeMult. It is significantly slower, although it converges in 2 fewer iterations with a smaller residual.
John On Fri, Aug 17, 2012 at 11:43 AM, John Fettig <john.fettig at gmail.com> wrote: > On Fri, Aug 17, 2012 at 11:27 AM, Mark F. Adams <mark.adams at columbia.edu> > wrote: >> John: could you perhaps test MatMatMult and MatMatTransposeMult so we could >> get some concrete numbers? > > For MatMatTransposeMult I get this error: > > [0]PETSC ERROR: --------------------- Error Message > ------------------------------------ > [0]PETSC ERROR: No support for this operation for this object type! > [0]PETSC ERROR: MatMatTransposeMult not supported for A of type mpiaij! > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: Petsc Development HG revision: > f9c6cac2d69c724a2258d4e0dc2f51b0825aa875 HG Date: Thu Aug 16 08:37:18 > 2012 -0700 > > For MatMatMult, my problem is not symmetric. It doesn't look like > MatTranspose() followed by MatMatMult is a viable option. > > John -------------- next part -------------- A non-text attachment was scrubbed... Name: mmtm_np1.out Type: application/octet-stream Size: 14909 bytes Desc: not available URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120817/5177f3bb/attachment-0001.obj>
