[petsc-dev] MtPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy does not set MatInfo

Jed Brown jedbrown at mcs.anl.gov
Thu Feb 7 01:06:51 CST 2013


Hong, this commit made -matptap_scalable the default, but the function does
not set MatInfo,

https://bitbucket.org/petsc/petsc-dev/commits/e0ff50333da9dbf29963b66b8316bc5638f0ff5e

Consequently, snes 'runex5' has concerning diffs like below. I think you
just need to set PetscInfo fields similar to
MatMatMultSymbolic_SeqAIJ_SeqAIJ(), but the information is a little
distributed so it didn't look super fast for me to fix. Can you do this?

324c273
<           total: nonzeros=169, allocated nonzeros=169
---
>           total: nonzeros=169, allocated nonzeros=0
330c279
<           Chebyshev: eigenvalue estimates:  min = 0.501622, max = 1.10357
---
>           Chebyshev: eigenvalue estimates:  min = 0.498899, max = 1.09758
347c296
<               total: nonzeros=625, allocated nonzeros=625
---
>               total: nonzeros=625, allocated nonzeros=0
362c311
<           total: nonzeros=625, allocated nonzeros=625
---
>           total: nonzeros=625, allocated nonzeros=0
369c318
<           Chebyshev: eigenvalue estimates:  min = 0.498945, max = 1.09768
---
>           Chebyshev: eigenvalue estimates:  min = 0.49182, max = 1.082
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130207/043079b2/attachment.html>


More information about the petsc-dev mailing list