def METHOD_NAME(self, linearization_implementation, num_samples):
    rng = np.random.default_rng(seed=1)
    self.smoothing_posterior.sample(
        rng=rng, t=self.dense_locations, size=num_samples
    )