From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 27237 invoked by alias); 12 Dec 2012 21:13:21 -0000 Received: (qmail 27224 invoked by uid 22791); 12 Dec 2012 21:13:20 -0000 X-SWARE-Spam-Status: No, hits=-3.4 required=5.0 tests=AWL,BAYES_00,KHOP_RCVD_UNTRUST,KHOP_THREADED,RCVD_IN_DNSWL_LOW,RCVD_IN_HOSTKARMA_YE X-Spam-Check-By: sourceware.org Received: from mail-ie0-f171.google.com (HELO mail-ie0-f171.google.com) (209.85.223.171) by sourceware.org (qpsmtpd/0.43rc1) with ESMTP; Wed, 12 Dec 2012 21:13:14 +0000 Received: by mail-ie0-f171.google.com with SMTP id 17so2882700iea.30 for ; Wed, 12 Dec 2012 13:13:13 -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=oXxq64vN2lQWaHRKuPTZsPtjKFfgQsBv2P4X7d9Ec8o=; b=j3sHI9JoQ5w4TlsFhqlty7GqBwGbhxaLGfcZkdz1MTWncz0disQezaXXDgA4r271Hh jEB3D9Qk8gN0KQSL9Gh3D0X+UyMzhU7hoDY0pVEsfKbEbrc1DGA5dv/SvkQpIOd1UgcI jUNmHwHSbFCqX/9hCi1kobxvNkSjmalIY+28dXuCP3V67Ttt+/z9ruawNeEg6wohMxcz xlt8wDB3ODALQxT1GekQR1CTJdgbUhir4nA0w83X4MQpeAyyala8sFUjJR4SkRWYUVKd A0AmrC6BGiOmTuPFJ8t/s8+tNJxAiWQa1e9DbeH66JY2ZrwggsTDrOJlWvSb6NgFQw7l 2qQg== Received: by 10.50.236.72 with SMTP id us8mr14788717igc.28.1355346793294; Wed, 12 Dec 2012 13:13:13 -0800 (PST) Received: from MacBook-Pro-de-Maxime.local ([75.98.19.134]) by mx.google.com with ESMTPS id as6sm2625204igc.8.2012.12.12.13.13.09 (version=SSLv3 cipher=OTHER); Wed, 12 Dec 2012 13:13:12 -0800 (PST) Message-ID: <50C8F364.2010109@calculquebec.ca> Date: Wed, 12 Dec 2012 21:13: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> In-Reply-To: Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 8bit X-Gm-Message-State: ALoCoQnr3rgl6iPPiZYYNCjMt8BsiJxq+Yw69qRZ6jQGpp2PW0fdIeyY2CUkBtoH3AdHAcAvBgFP 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/msg00004.txt.bz2 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