  In this paper, the performance characteristics of different solution
techniques and program packages to solve a large number of independent ordinary
differential equation systems is examined. The employed hardware are an Intel
Core i7-4820K CPU with 30.4 GFLOPS peak double-precision performance per cores
and an Nvidia GeForce Titan Black GPU that has a total of 1707 GFLOPS peak
double-precision performance. The tested systems (Lorenz equation,
Keller--Miksis equation and a pressure relief valve model) are non-stiff and
have low dimension. Thus, the performance of the codes are not limited by
memory bandwidth, and Runge--Kutta type solvers are efficient and suitable
choices. The tested program packages are MPGOS written in C++ and specialised
only for GPUs; ODEINT implemented in C++, which supports execution on both CPUs
and GPUs; finally, DifferentialEquations.jl written in Julia that also supports
execution on both CPUs and GPUs. Using GPUs, the program package MPGOS is
superior. For CPU computations, the ODEINT program package has the best
performance.
