[Bug-gsl] ode solver function in multiple mpi environment

[caoj]

Hi:

My mpi program works fine using IMSL C library. But after I use gsl instead of 
IMSL, my mpi program has different results using different number of CPUs. 
These two programs are same except using different ode solver library function 
call(one is IMSL, the other is gsl). 

Do you think it is possible ode solver function in gsl may not be correct in 
the multiple mpi environment? Thanks.

Jun