Hi Hrvoje, To be 100% sure that the bug is fixed, we should wait for the results from Flavio's run with 128 nodes. For testing locally you can run valgrind's massif with and without the patch and look at the difference of the allocated memory. In order to spot it, you'll need a decent number of ranks. I used 1024, which produces lists of 4MB. I took the standard 2D cavity case with 1000x1000 cells. The command to run: mpirun -np 1024 --oversubscribe valgrind --tool=massif icoFoam -parallel If you...
Hi Flavio, the problem comes from the Pstream class while allocating linear and tree communication lists (discovered using valgrind's massif tool). The lists are N large and their entries are of type commsStruct which is approximately also N large, where N is the number of MPI ranks. Thus, these lists are by design of size N^2. The difference to the OpenFOAM version is that, there, the lists are not allocated at start-up, but only sized to N and each entry (commsStruct) is constructed (and therefore...
Bugfix solutionControl: initialize faceU.oldTime() from U.oldTime()
Bugfix DLB: do not account for flipped faces for surfacePatchFields other than of type scalar in fvFieldDecomposer
Bugfix DLB: add write() method in passiveProcessorFvPatchField and passiveProcessorFvsPatchField to enable their encapsulation in a dictionary