Add discretized penalty matrix P for iterative parabolic regularization

[AlePalu]

The implementation of method P() for SpaceTimeParabolicBase<Model, IterativeSolver> in fdaPDE/models/space_time_parabolic_base.h is missing. This prevents, e.g., the definition of the GCV index for a model using such regularization.