Pe trunk interpolation

Trying a pull request. 
Check new issue for my questions.

especially about inlining in interpolant_impl.hh 
method val() and val_p1()