[petsc-users] GAMG setup scalability

John Fettig john.fettig at gmail.com
Fri Aug 17 10:52:35 CDT 2012


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>


More information about the petsc-users mailing list