From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 2423 invoked by alias); 13 Dec 2012 13:22:26 -0000 Received: (qmail 2407 invoked by uid 22791); 13 Dec 2012 13:22:23 -0000 X-SWARE-Spam-Status: No, hits=-3.6 required=5.0 tests=BAYES_00,KHOP_RCVD_UNTRUST,KHOP_THREADED,RCVD_IN_DNSWL_LOW,RCVD_IN_HOSTKARMA_YE X-Spam-Check-By: sourceware.org Received: from mail-vb0-f41.google.com (HELO mail-vb0-f41.google.com) (209.85.212.41) by sourceware.org (qpsmtpd/0.43rc1) with ESMTP; Thu, 13 Dec 2012 13:22:19 +0000 Received: by mail-vb0-f41.google.com with SMTP id l22so2238044vbn.0 for ; Thu, 13 Dec 2012 05:22:19 -0800 (PST) X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=20120113; h=message-id:date:from:user-agent:mime-version:to:cc:subject :references:in-reply-to:content-type:content-transfer-encoding :x-gm-message-state; bh=IDO2zgs/oBnMrI/UoYgeSsRbgYcu/B3Gv7ECsshG3nY=; b=efKvhbyALwMZcOiwjIIs4Z1GVzwvujsAFES1Rzfr2lZ13b23ssx/UIyx4x2ATQr/1u GUd2BybPn0CuOy8yTk7sgS+0A+nQYdYmKq1br4PeHEIV7lI+9PU7FplGeh3BrArovvVX II1U4QV96WIg+xRVDzxdU7A1SHVkcu29z5cWdoRET4DB+/yhjd5ByVrwBbFenH4Hm6Zc 3pGZiZXxA3xr9b4y8uaTXjy7T3pXEujGxqd1qeVd/gRQIULs2oRNluSqo9w8U2PxtVh7 VvheJeIvBMDvCb+KDvK/gBr4/HV+ww3elDk4YVT9/DxjwwwMk++HDs5PDzRp8h3r5f0G uCxA== Received: by 10.52.67.198 with SMTP id p6mr2871385vdt.3.1355404938973; Thu, 13 Dec 2012 05:22:18 -0800 (PST) Received: from MacBook-Pro-de-Maxime.local ([2607:fa48:6e41:cb80:a0ed:53ab:c9a3:c993]) by mx.google.com with ESMTPS id dp6sm1204381vec.11.2012.12.13.05.22.16 (version=SSLv3 cipher=OTHER); Thu, 13 Dec 2012 05:22:17 -0800 (PST) Message-ID: <50C9D68A.1050903@calculquebec.ca> Date: Thu, 13 Dec 2012 13:22:00 -0000 From: Maxime Boissonneault User-Agent: Mozilla/5.0 (Macintosh; Intel Mac OS X 10.8; rv:17.0) Gecko/17.0 Thunderbird/17.0 MIME-Version: 1.0 To: Rhys Ulerich CC: gsl-discuss@sourceware.org Subject: Re: Adding OpenMP support for some of the GSL functions References: <50C791BB.4060303@calculquebec.ca> <50C8F364.2010109@calculquebec.ca> <50C914C4.2060102@calculquebec.ca> In-Reply-To: Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 8bit X-Gm-Message-State: ALoCoQkMTYKZOd802uplz1YWT5OXrZPz0MdOUOnUos2jH6iKYZOvGX15viY0oPDMmL8wRHME/D3Z Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sourceware.org X-SW-Source: 2012-q4/txt/msg00010.txt.bz2 I am doing that too, but any gain we can get is an important one, and it turns out that by parallelizing rkf45_apply, my simulation runs 30% faster on 8 cores. Maxime Le 2012-12-12 21:28, Rhys Ulerich a écrit : >> 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. > True, solving one differential equation is serial. I'd thought you > were performing a parameter sweep. > > - Rhys