From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <gsl-discuss-return-5752-listarch-gsl-discuss=sources.redhat.com@sourceware.org>
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 <gsl-discuss@sourceware.org>; 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 <maxime.boissonneault@calculquebec.ca>
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 <rhys.ulerich@gmail.com>
CC: gsl-discuss@sourceware.org
Subject: Re: Adding OpenMP support for some of the GSL functions
References: <50C791BB.4060303@calculquebec.ca> <CAKDqugSjvphQBsVYejauNAO8T01GFfnZe8yKfOGFd9+iQAr_qw@mail.gmail.com>
In-Reply-To: <CAKDqugSjvphQBsVYejauNAO8T01GFfnZe8yKfOGFd9+iQAr_qw@mail.gmail.com>
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: <gsl-discuss.sourceware.org>
List-Subscribe: <mailto:gsl-discuss-subscribe@sourceware.org>
List-Archive: <http://sourceware.org/ml/gsl-discuss/>
List-Post: <mailto:gsl-discuss@sourceware.org>
List-Help: <mailto:gsl-discuss-help@sourceware.org>, <http://sourceware.org/ml/#faqs>
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