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

[Ce Qin]

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/piper
>> mail/petsc-users/2014-September/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/20170323/f9efdf95/attachment.html>