Create AugRieFun.hpp

Functional implementation of the Augemted Riemann solver:
- based on implementation in AugRieCuda
- improved solution of linear system
- several optimizations