[petsc-dev] MatMatMult crashes with MATMPIAIJ using one processor

[Hong]

This is in master now.
Hong

On Wed, Mar 22, 2017 at 9:43 PM, Ce Qin <qince168 at gmail.com> wrote:

> Dear Hong,
>
> I have tested the new feature, it works now.
> Thanks.
>
> Best,
> Ce Qin
>
>
> 2017-03-17 23:58 GMT+08:00 Hong <hzhang at mcs.anl.gov>:
>
>> Ce Qin:
>> This is fixed in branch hzhang/matmatops_mpiaij4np1.
>> I consider it as a new feature - support mpiaij for np=1, not a bugfix.
>>
>> Once it passes our regression tests, I'll merge it to petsc master branch.
>> You may test it. Let us know you see any problem.
>> Thanks for reporting it.
>>
>> Hong
>>
>> Dear all,
>>>
>>> I encountered a problem that MatMatMult crashes using only one
>>> processor. The matrix is MATMPIAIJ. I am using the master branch. The same
>>> program can run in parallel without problem.
>>>
>>> Valgrind gives the following stack trace:
>>>
>>> ==5543== Invalid read of size 8
>>> ==5543==    at 0xDC094CE: MatGetBrowsOfAoCols_MPIAIJ (mpiaij.c:4509)
>>> ==5543==    by 0xDC3FB95: MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable
>>> (mpimatmatmult.c:218)
>>> ==5543==    by 0xDC45193: MatMatMult_MPIAIJ_MPIAIJ (mpimatmatmult.c:61)
>>> ==5543==    by 0xD9C657D: MatMatMult (matrix.c:9450)
>>>
>>> ==5543==
>>> ==5543== Invalid read of size 4
>>> ==5543==    at 0xDC094D7: MatGetBrowsOfAoCols_MPIAIJ (mpiaij.c:4510)
>>> ==5543==    by 0xDC3FB95: MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable
>>> (mpimatmatmult.c:218)
>>> ==5543==    by 0xDC45193: MatMatMult_MPIAIJ_MPIAIJ (mpimatmatmult.c:61)
>>> ==5543==    by 0xD9C657D: MatMatMult (matrix.c:9450)
>>>
>>> ==5543==
>>> ==5543== Invalid read of size 4
>>> ==5543==    at 0xDC0793B: MatGetBrowsOfAoCols_MPIAIJ (mpiaij.c:4513)
>>> ==5543==    by 0xDC3FB95: MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable
>>> (mpimatmatmult.c:218)
>>> ==5543==    by 0xDC45193: MatMatMult_MPIAIJ_MPIAIJ (mpimatmatmult.c:61)
>>> ==5543==    by 0xD9C657D: MatMatMult (matrix.c:9450)
>>>
>>> ==5543== Invalid read of size 4
>>> ==5543==    at 0xDC0794E: MatGetBrowsOfAoCols_MPIAIJ (mpiaij.c:4520)
>>> ==5543==    by 0xDC3FB95: MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable
>>> (mpimatmatmult.c:218)
>>> ==5543==    by 0xDC45193: MatMatMult_MPIAIJ_MPIAIJ (mpimatmatmult.c:61)
>>> ==5543==    by 0xD9C657D: MatMatMult (matrix.c:9450)
>>>
>>> This problem may be related to this thread
>>> http://lists.mcs.anl.gov/pipermail/petsc-users/2014-Septembe
>>> r/022864.html. If you need an example that reproduces this
>>> problem, please let me know.
>>>
>>> Best,
>>> Ce Qin
>>>
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20170322/d5cc786d/attachment.html>