From: Maxime Boissonneault <maxime.boissonneault@calculquebec.ca>
To: Rhys Ulerich <rhys.ulerich@gmail.com>
Cc: gsl-discuss@sourceware.org
Subject: Re: Adding OpenMP support for some of the GSL functions
Date: Wed, 12 Dec 2012 23:35:00 -0000 [thread overview]
Message-ID: <50C914C4.2060102@calculquebec.ca> (raw)
In-Reply-To: <CAKDqugQKdtuQJrgXGSFXZ2_mqQ+0e_8O24JaqmkGdjR4WYKy8g@mail.gmail.com>
Hi Rhys,
I am not sure if you mean parallelizing the for loop in this code, but
this would not work, as each iteration is not independent from the next
one. Differential equation solving is intrinsically serial.
I can not either just split the problem into smaller chunks (lets say
100k istead of 1M), because they are coupled, i.e. func can only be
defined for the whole system, not for some of its parts.
Thanks,
Maxime
Le 2012-12-12 16:36, Rhys Ulerich a écrit :
>> I have not, since it is much simpler to add 6 pragmas to the loops of
>> rk45_apply to achieve the same goal.
> I suspect, from glancing at the first example at
> http://www.gnu.org/software/gsl/manual/html_node/ODE-Example-programs.html,
> you could achieve parallelization atop stock GSL builds with fewer
> than 6 pragmas.
>
> I also suspect that by using OpenMP tasks over the coarser granularity
> concept (1 ODE) instead of the parallel-fors over the finer
> granularity method (1 step) you may see better overall performance.
>
> - Rhys
next prev parent reply other threads:[~2012-12-12 23:35 UTC|newest]
Thread overview: 16+ messages / expand[flat|nested] mbox.gz Atom feed top
2012-12-11 20:04 Maxime Boissonneault
2012-12-12 16:35 ` Frank Reininghaus
2012-12-12 21:11 ` Maxime Boissonneault
2012-12-12 21:41 ` Rhys Ulerich
2012-12-12 23:32 ` Maxime Boissonneault
2012-12-12 17:05 ` Rhys Ulerich
2012-12-12 21:13 ` Maxime Boissonneault
2012-12-12 21:36 ` Rhys Ulerich
2012-12-12 23:35 ` Maxime Boissonneault [this message]
2012-12-13 2:29 ` Rhys Ulerich
2012-12-13 13:22 ` Maxime Boissonneault
2012-12-13 15:53 ` Rhys Ulerich
2012-12-13 16:44 ` Rhys Ulerich
2012-12-13 21:07 ` Maxime Boissonneault
2012-12-13 21:05 ` Maxime Boissonneault
2012-12-13 21:14 ` Maxime Boissonneault
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=50C914C4.2060102@calculquebec.ca \
--to=maxime.boissonneault@calculquebec.ca \
--cc=gsl-discuss@sourceware.org \
--cc=rhys.ulerich@gmail.com \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).