LSODA using OpenMP

[LeighKorbel]

Has anyone acheived OpenMP threading with this algorithm? I took the test example and tried using basic OpenMP compiler pragmas to thread the for loop that calls LSODA. The example with OpenMP pragmas is here (and includes omp.h):

test_openmp.zip

The algorithm fails because the initial time step seems randomized rather than what is hard coded. Multiple runs of the test program results in different "t" values to begin each time. If anyone has a minute and doesn't mind looking I would appreciate it. Might be worth noting that I am using gcc.

[mattfidler]

I don't have time to look into your code. However, I have used open MP threading with a very close variant of this package with RxODE https://github.com/nlmxirdevelopment/RxODE. Because RxODE has been on CRAN it has been tested in multiple environments including linux, mac with various different linear algebra systems, so it should work.