Changed type of cell_average to vector< Vector<double> >
Fixed computation of dx to be dimension independent