WIP: Gaussian quadrature #4633

Job	Run time
CI-Cuda-Host (g++, Release, cuda, serial)	37m 35s
CI-Cuda-Host (g++, Release, cuda, mpi)	40m 25s
CI-Rocm-HIP (hipcc, Release, hip, serial)	31m 28s
CI-Cuda-Host (g++, Release, host, serial)	11m 8s
CI-Rocm-HIP (hipcc, Release, hip, mpi)	31m 22s
CI-Cuda-Host (g++, Release, host, mpi)	11m 10s
CI-Rocm-HIP (hipcc, Debug, hip, serial)	53m 9s
CI-Cuda-Host (g++, Debug, cuda, serial)	1h 7m 19s
CI-Rocm-HIP (hipcc, Debug, hip, mpi)	52m 10s
CI-Cuda-Host (g++, Debug, cuda, mpi)	1h 11m 6s
CI-Cuda-Host (g++, Debug, host, serial)	9m 50s
CI-Cuda-Host (g++, Debug, host, mpi)	9m 53s
CI-Cuda-Host (clang++-15, Release, cuda, serial)	24m 30s
CI-Cuda-Host (clang++-15, Release, cuda, mpi)	25m 30s
CI-Cuda-Host (clang++-15, Release, host, serial)	8m 52s
CI-Cuda-Host (clang++-15, Release, host, mpi)	9m 20s
CI-Cuda-Host (clang++-15, Debug, host, serial)	7m 9s
CI-Cuda-Host (clang++-15, Debug, host, mpi)	7m 25s
	8h 29m 21s

Provide feedback