Why does my Rcpp function sometimes return NaN and other times the correct answer for the same inputs?

Question

I am writing an Rcpp function to multiply a banded matrix by a vector. I have an R wrapper function that makes it easier to call the C++ function. I get the same issue regar