integrating cuda matrix inverse in r programming

Can someone pls help us write a wrapper function to call a cuda c matrix inverse program in R? Our matrix is 10000x10000. The inversed matrix need to be accessed by R