  Local scale invariance (LSI) has been recently proposed as a possible
extension of the dynamical scaling in systems at the critical point and during
phase ordering. LSI has been applied inter alia to provide predictions for the
scaling properties of the response function of non-equilibrium critical systems
in the aging regime following a quench from the high-temperature phase to the
critical point. These predictions have been confirmed by Monte Carlo
simulations and analytical results for some specific models, but they are in
disagreement with field-theoretical predictions. By means of Monte Carlo
simulations of the critical two- and three-dimensional Ising model with Glauber
dynamics, we study the intermediate integrated response, finding deviations
from the corresponding LSI predictions that are in qualitative agreement with
the field-theoretical computations. This result casts some doubts on the
general applicability of LSI to critical dynamics.
