efficient stencil storage, multiple stencil option, rearranged initial and boundary condition number, modification in exact solution to cover whole domain and 1D output for density profile