[petsc-dev] MtPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy does not set MatInfo

Hong Zhang hzhang at mcs.anl.gov
Thu Feb 7 10:28:43 CST 2013


Thanks for reporting it. I'm checking it.
Hong

On Thu, Feb 7, 2013 at 1:06 AM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
> 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
>



More information about the petsc-dev mailing list