We now compute LxF flux like shu, use cell average value for the wave speed.
Added new time step function for cartesian grid and general grids, uses cell average value.
Added different ssprk scheme based on chosen degree
Size of unit_support_points is equal to fe.dofs_per_cell. We dont need to use base_i.