add a new method which averages out a dg function in order to obtain a p1 function.  This is used to pass Neumann data coming from a rod to the hierarchic error estimator

[[Imported from SVN: r1866]]