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 by Ankang Gao