def METHOD_NAME(self, timeline):
    NL = timeline.number_of_time_levels()
    gdof = self.space.number_of_global_dofs()
    uh = np.zeros((gdof, NL), dtype=np.float)
    uh[:, 0] = self.space.interpolation(lambda x:self.pde.solution(x, 0.0))
    return uh