Pseudospectral multi-domain method for incompressible viscous flow computation