Pe trunk interpolation

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