The R package irregulAR1 provides efficient density evaluation and random number generation for an irregularly sampled stationary AR(1) process with Gaussian errors.
You can install irregulAR1 from github with:
# install.packages("devtools")
devtools::install_github("BenjaK/irregulAR1")