Subsequently the left hand side for the energy system, the momentum system, the pressure system and the turbulence system are calculated in subroutine mafillvlhs.f (all equation systems except pressure for incompressible fluids; notice that this is one and the same left hand side matrix for all these equations systems) and mafillplhs.f (pressure system for incompressible fluids). For compressible fluids all systems are lumped in subroutine lump.f. The lumped matrix is diagonal. However, the lumped matrix is not stored as such. Indeed, out of efficiency considerations the diagonal of the original matrix, stored in field adb, is replaced by itself minus the lumping matrix. The inverse of the lumping matrix is stored in adl. For incompressible fluids all equation systems except the pressure equations are lumped. For these pressure equations a LU decomposition is performed for later use in the solution phase of all systems of equations.
For compressible fluids the equation system is reformatted into a row-by-row format in subroutine convert2rowbyrow for use in the smoothing procedure.
At this point the preparation phase is finished an the major loop starts calculating the solution at the subsequent time points