The paper analyses the fast method for solving large scale algebraic saddle-point systems arising from a smooth variant of a fictitious domain method.