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

[Brian Gough]

address@hidden wrote:
> 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.

Hello,
In GSL the ode solvers must be allocated on a per-thread basis, their 
memory should not be shared.  If you think there is a bug please can you 
send a small example program which demonstrates it.  If you still have a 
problem try the mailing list address@hidden

--
best regards,

Brian Gough
(GSL Maintainer)

Network Theory Ltd,
Publishing the GSL Manual - http://www.network-theory.co.uk/gsl/manual/