C++ code for solving Incompressible Navier-Stokes