14
Tôi đang tìm kiếm một giải pháp cho một thể thiếu đôi đó là nhanh hơn so vớitính tích phân kép trong R nhanh chóng
integrate(function(y) {
sapply(y, function(y) {
integrate(function(x) myfun(x,y), llim, ulim)$value
})
}, llim, ulim)
với ví dụ
myfun <- function(x,y) cos(x+y)
llim <- -0.5
ulim <- 0.5
Tôi tìm thấy một old paper rằng gọi một chương trình FORTRAN được gọi là quad2d, nhưng tôi không thể tìm thấy bất kỳ điều gì khác ngoài một số trang trợ giúp cho MATLAB cho phần còn lại. Vì vậy, tôi đang tìm một thư viện C hoặc FORTRAN có thể làm nhanh gấp đôi tích phân (tức là không có vòng lặp thuận), và có thể được gọi từ R. Tất cả các ý tưởng đều được đánh giá rất cao, miễn là chúng đều tương thích với GPL.
Nếu giải pháp liên quan đến việc gọi các chức năng khác từ các thư viện đã được gửi cùng với R, tôi cũng muốn nghe từ chúng.
Bạn đã xem xét: 'pracma :: dblquad',' pracma: simpson2d', hoặc các chức năng trong các gói hình khối và R2Cuba? –