You can intermix this with CUDA kernels, but you cannot call them from within the kernel itself. While things like matrix multiplication already work (no need to kernels) it is simple 46 to do to write a grid-stride loop to do that in a kernel if you want.