Replace LU decomposition with Cholesky in GPU RBPF

[THargreaves]

We currently use cuBLAS's batched LU through a strided wrapper CUDA.jl provides. We need to write a strided wrapper for the batch Cholesky provided by cuSolver. This will greatly improve numerical stability and might help with performance.