[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