X-Recipient: archive-cygwin AT delorie DOT com DomainKey-Signature: a=rsa-sha1; c=nofws; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:to:from:subject:date:message-id:references :mime-version:content-type:content-transfer-encoding; q=dns; s= default; b=k40d2/dYejkBIUa8VX5TwWjxs3l+thJfPYPMAIW0Rj3+G0yJ8L7ZY +f3pGJa9swCf23oy2M5DfWOGY1oxJvr6/tIRXKbetow7SUdf5+M7g5M4JZz/Gr+5 uFobs/nNXjEkonqErizLGD0oOw+VxwB6C/iCyZhrm73V8v1FmWJTDM= DKIM-Signature: v=1; a=rsa-sha1; c=relaxed; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:to:from:subject:date:message-id:references :mime-version:content-type:content-transfer-encoding; s=default; bh=bkUbDtiBFkXKDu44zH/AxpAOomQ=; b=IpLWv5B2sPArwEfDwmKpgzYpVK79 tqGbEEySNq7CDoATmSu23ktyiiIb5hIrodwYMrdMDzeiUOQYHzFBw5Sa1P7zzaE9 WOJTqOfR+jgs5Re9TWBeKfdVbWFblGEXwhhm6L9gD9dc38bZnsa+slY1qdVDPINb lpkDrPaF0KaTR+k= Mailing-List: contact cygwin-help AT cygwin DOT com; run by ezmlm List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: cygwin-owner AT cygwin DOT com Mail-Followup-To: cygwin AT cygwin DOT com Delivered-To: mailing list cygwin AT cygwin DOT com Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=1.6 required=5.0 tests=AWL,BAYES_50,RCVD_IN_DNSWL_LOW,RP_MATCHES_RCVD,SPAM_BODY,SPF_HELO_PASS,SPF_PASS,UNPARSEABLE_RELAY autolearn=no version=3.3.2 spammy=tony, Tony, Parallelization, UD:netcologne.de X-HELO: plane.gmane.org To: cygwin AT cygwin DOT com From: Tony Kelman Subject: Re: Parallelization Date: Sun, 17 Jul 2016 11:56:55 +0000 (UTC) Lines: 33 Message-ID: References: Mime-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit User-Agent: Loom/3.14 (http://gmane.org/) X-IsSubscribed: yes Thomas Koenig netcologne.de> writes: > > I am working on a non-linear, boring system of ODEs, boring being > defined as non-stiff and without other numerical surprises. > So, CV_ADAMS works well. > > The only interesting part is that there are very many of the ODEs, > around 30000 at the moment, and that each ODE depends on around half of > the others, so a Jacobian (were I to calculate it) would be quite > densely packed. > > Sundials works well for this problem, but it takes rather a lnog time > and only utilizes a single processor on a 20-core machine. > > So, the question: Is there a way to parallelize the calculation for > this? The references in the Sundials docs that I have seen only refer > to parallelization of solving linear equations, which I do not need > to do. > > Regards > > Thomas > This is more of a general numerical ODE or Sundials usage question rather than a cygwin specific one, but I would try openmp or mpi on your function evaluations first, if that's taking most of the time. Do you know whether sundials is thread-safe? -Tony -- Problem reports: http://cygwin.com/problems.html FAQ: http://cygwin.com/faq/ Documentation: http://cygwin.com/docs.html Unsubscribe info: http://cygwin.com/ml/#unsubscribe-simple