[petsc-dev] Regression in Vec assembly

Stefano Zampini stefano.zampini at gmail.com
Tue Feb 19 22:47:38 CST 2019


In the link I have sent you, ex29 was modified to invert the order of values set. There’s also a command line that shows the error

mpiexec -n 3 ./ex29 -n 4 -invert -bs 2


I have added the correct output, but if you run with bts you’ll see some of the entries are scrambled up.


> On Feb 20, 2019, at 1:27 AM, Jed Brown <jed at jedbrown.org> wrote:
> 
> Is this what I'm supposed to see?  What should I be looking at?
> 
> $ mpiexec.hydra -n 3 mpich-clang/tests/vec/vec/examples/tests/ex29 -n 126 -info | grep malloc
> [0] VecAssemblyBegin_MPI_BTS(): Stash has 252 entries, uses 2 mallocs.
> [0] VecAssemblyBegin_MPI_BTS(): Block-Stash has 0 entries, uses 0 mallocs.
> [0] VecAssemblyBegin_MPI_BTS(): Stash has 0 entries, uses 0 mallocs.
> [0] VecAssemblyBegin_MPI_BTS(): Block-Stash has 252 entries, uses 2 mallocs.
> 
> $ mpiexec.hydra -n 3 mpich-clang/tests/vec/vec/examples/tests/ex29 -n 126 -vec_assembly_legacy -info | grep malloc
> [0] VecAssemblyBegin_MPI(): Stash has 252 entries, uses 2 mallocs.
> [0] VecAssemblyBegin_MPI(): Block-Stash has 0 entries, uses 0 mallocs.
> [0] VecAssemblyBegin_MPI(): Stash has 0 entries, uses 0 mallocs.
> [0] VecAssemblyBegin_MPI(): Block-Stash has 252 entries, uses 2 mallocs.
> 
> Stefano Zampini <stefano.zampini at gmail.com> writes:
> 
>> Jed,
>> 
>> I’m resurrecting this thread since I believe this bug has not yet been fixed.
>> I have spent a couple of hours today tracking down a bug in my code, that I finally confined to be related with VecSetValuesBlocked.
>> 
>> I have pushed an example that shows the problem here https://bitbucket.org/petsc/petsc/commits/1cb7fb90b416ce1d4e4767a0d5a2ef5043656d53?at=stefano_zampini/feature-pcbddc-enablelor <https://bitbucket.org/petsc/petsc/commits/1cb7fb90b416ce1d4e4767a0d5a2ef5043656d53?at=stefano_zampini/feature-pcbddc-enablelor>
>> 
>> Could you please have a look?
>> 
>> Thanks,
>> Stefano
>> 
>>> On Apr 24, 2016, at 10:52 PM, Jed Brown <jed at jedbrown.org> wrote:
>>> 
>>> Lisandro Dalcin <dalcinl at gmail.com> writes:
>>> 
>>>> Do you remember about this regression? My PetIGA tests are still
>>>> failing with petsc/master.
>>> 
>>> Yes, I have most of a fix in a branch that I'm working to finish.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20190220/5e332a48/attachment.html>


More information about the petsc-dev mailing list