Effect of Flux Function Order and Working Precision in Spectral Element Methods