Peter, could it be that this version was the starting point of a discussion last summer and we have then created the new class FdmTimeDepDirichletBoundary, which covers the new functionality and have left the old class untouched?

I think the reason for creating a new class instead of extending the existing one fits to Luigi's arguments. If a remeber correctly we've finally reverted the change 

        void applyBeforeSolving(operator_type&, array_type&) const {}

because this one shows a better convergence for dt->0 for barrier options.