Add files for high order integration with quadratures

The file is a bit complex to compile and clang hangs quite a lot
during the optimization process (>1h and not compiled yet!).