added functions for inviscid flux initialization by conservative vars, and addition/subtraction of inviscid fluxes
changed hyperplanes calculation to save i, j, k indices
starting to transform LUSGS to sweep on hyperplanes