Optimize compute_cell_average()

[cpraveen]

Since we use QGauss nodal basis, we can optimize this function. No need to call get_function_values.

[cpraveen]

In src we cannot do this since it also has Pk basis. It has been implemented in src_gll.