вычисление двойных интегралов в R быстро

Я ищу решение для двойного интеграла, которое быстрее, чем

integrate(function(y) { 
   sapply(y, function(y) {
     integrate(function(x) myfun(x,y), llim, ulim)$value
   })
 }, llim, ulim)

, например,

myfun <- function(x,y) cos(x+y)
llim <- -0.5
ulim <- 0.5

Я нашел старый документ , в котором говорилось о программе FORTRAN под названием quad2d , но я не смог найти ничего, кроме нескольких страниц справки по Matlab для остальных. Поэтому я ищу библиотеку C или FORTRAN, которая может быстро выполняйте двойные интегралы (т.е. без цикла sapply), и это может быть вызвано из R. Все идеи очень ценятся , если все они совместимы с GPL.

Если решение включает вызов других функций из библиотек, которые уже поставляются с R, я бы тоже хотел получить известие от них.