krige in gstat() package

HI,
  I wonder does anyone have experience with doing sequential gaussian
simulation with krige() function in gstat?

I find it VERY slow compared to use krige() to achieve kriging function
itself..  I wonder why, is that because it has to model the variogram, and
do the kriging separately for each point to be simulated?

so it would be N times slower to achieve the simulation than the kriging
if the number of points to be estimated is N??

Any hints is welcome and appreciated!
yan