Gather gets stuck in ContField::GlobalSolve

In my parallel test wing_15.xml using 2 cores, after the first time step, one thread satisfies (contNcoeffs - NumDirBcs == 0) in function ContField::GlobalSolve. So, there is only one thread runs GlobalLinSysSharedPtr LinSys = GetGlobalLinSys(key);. This makes the code stuck in Gs::Gather (in PreconditionerDiagonal::StaticCondDiagonalPreconditionerSum: asmMap->UniversalAssembleBnd(vOutput);).

Removing if(contNcoeffs - NumDirBcs > 0) solves this problem.

Edited Jul 22, 2020 by Ankang Gao

To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information