Replace a single loop, and a few changes are needed for the small kernel PR #432

Job	Run time
gcc-cpu_gh / define_matrix	4s
clang-format_check / C++ formatting check	1m 13s
gcc-cuda / define_matrix	3s
manual-cpu_gh	-1s
manual-gpu_cuda	-1s
gcc-cuda / gcc-cuda / Debug - single	1s
gcc-cuda / gcc-cuda / Release - double	14h 47m 35s
gcc-cuda / gcc-cuda / Release - single	1s
gcc-cuda / gcc-cuda / Debug - double	1s
gcc-cpu_gh / gcc-cpu / Debug - double	16m 29s
gcc-cpu_gh / gcc-cpu / Debug - single	3m 1s
gcc-cpu_gh / gcc-cpu / Release - double	15m 50s
gcc-cpu_gh / gcc-cpu / Release - single	3m 20s
	15h 27m 36s

Provide feedback