Adding OpenMP support for some of the GSL functions

Maxime Boissonneault maxime.boissonneault@calculquebec.ca
Wed Dec 12 21:13:00 GMT 2012


Hi Rhys,
I have not, since it is much simpler to add 6 pragmas to the loops of 
rk45_apply to achieve the same goal.

Maxime

Le 2012-12-12 12:05, Rhys Ulerich a écrit :
>> I am using GSL from another library of my own to perform numerical
>> integration of vectorial differential equations. After optimizing and
>> parallelizing most of my library, I ended up with the conclusion that GSL is
>> a major bottle neck in my computation, simply because it is not parallelized
>> to exploit multi-core achitectures.
> Have you tried solving many such problems in parallel using
> threadprivate GSL workspaces?  This permits you to put a parallel
> section outside of many GSL invocations by ensuring GSL's working
> storage is kept local to each thread.
>
> - Rhys




More information about the Gsl-discuss mailing list