Work on pdd, now it mostly works for the stationary case, also for higher order element. Just some small bugs :)