My work for the month is to optimise the matrix adressing scheme of OpenFOAM to reduce cache miss. Initial idea is to replace the LDU Matrix addressing scheme with Diagnol matrix addressing scheme using multiple arrays for structures meshes. Will check it out on a simple 2D Poisson equation using Gauss-Seidel solver to check performance benefits, and will eventually proceed to introduce a new solver as a plugin if benefits are noticeable.
Wish me luck.