Implementation of multimode Gaussian unitaries using primitive components