Computational efficiency of the hybrid penalty-pseudocompressibility method for incompressible flow (Q1327277)

We implement the title method to within a specified error, and also a method of subcycling the rapidly-varying terms. This makes the hybrid method directly comparable to standard implicit methods of solving the incompressible flow equations. We find that the hybrid method is still inferior with regard to computational work to methods based on conjugate gradients, and this conclusion appears to hold independently of mesh resolution.