--Added option to output additional 'restart' file contatining x,y,z locations of all solution points
--Added functionality to set inviscid initial conditions to boundary conditions if not provided
Added shock capturing on CPU
Fixed local time stepping bug for quads, implemented variation of local time stepping on GPU
Fixed 'bug' preventing meshes with n_eles < 10 from working.
fixing TGV data analysis script.
beginning to implement forcing for linear-elasticity method.
Added default values for many commonly-used input parameters.
Merge with master and fix merge conflicts (been a while).