From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 3620 invoked by alias); 9 Jul 2012 17:39:24 -0000 Received: (qmail 3612 invoked by uid 22791); 9 Jul 2012 17:39:23 -0000 X-SWARE-Spam-Status: No, hits=-1.4 required=5.0 tests=AWL,BAYES_00,KHOP_THREADED,RCVD_IN_HOSTKARMA_NO,T_RP_MATCHES_RCVD X-Spam-Check-By: sourceware.org Received: from ipmx2.colorado.edu (HELO ipmx2.colorado.edu) (128.138.128.232) by sourceware.org (qpsmtpd/0.43rc1) with ESMTP; Mon, 09 Jul 2012 17:39:10 +0000 X-IronPort-Anti-Spam-Filtered: true X-IronPort-Anti-Spam-Result: Ap4EAGsW+0/AqBGb/2dsb2JhbABFhWazF4IgAQEFIxVRCxQBAwICBSECAg8CRhMGAgEBEId5qQiJcokEgSCKIIJwggqBEgOISZxygn4 Received: from omr-raz-2-priv.int.colorado.edu ([192.168.17.155]) by ipmx2-priv.int.colorado.edu with ESMTP; 09 Jul 2012 11:39:09 -0600 Received: from bonanza.ngdc.noaa.gov (EHLO bonanza.ngdc.noaa.gov) ([140.172.179.41]) by omr-raz-2-priv.int.colorado.edu (MOS 4.1.10-GA FastPath queued) with ESMTP id DRR47287 (AUTH alken); Mon, 09 Jul 2012 11:39:09 -0600 (MDT) Message-ID: <4FFB173C.50905@colorado.edu> Date: Mon, 09 Jul 2012 17:39:00 -0000 From: Patrick Alken User-Agent: Mozilla/5.0 (X11; Linux i686 on x86_64; rv:13.0) Gecko/20120614 Thunderbird/13.0.1 MIME-Version: 1.0 To: gsl-discuss@sourceware.org Subject: Re: Fwd: log probability density functions (PDFs) References: In-Reply-To: Content-Type: text/plain; charset=UTF-8; format=flowed Content-Transfer-Encoding: 7bit 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-q3/txt/msg00001.txt.bz2 The best option would be to make a standalone extension (see the main gsl page for examples), so people can look at it and try it out. We can add the extension to the main web page and possibly add it to the main source code later if it looks useful. On 07/09/2012 11:29 AM, Sam Mason wrote: > Hi, > > I sent the following message to the help-gsl list a while ago, but > haven't received a reply... > > Is this sort of contribution useful to the project? > > Thanks, > > Sam > > > ---------- Forwarded message ---------- > From: Sam Mason > Date: 15 June 2012 17:59 > To: help-gsl@gnu.org > > > Hi, > > I sent a message to this a couple of months ago enquiring about the > availability of log-density functions. I've been using my own > versions of these since then, and thought it would be good to send > them back to the community. > > I've gone through most of the distributions and added in code > calculating their log-densities. I've copied the naming and layout > convention of the dirichlet distribution; which already had the > gsl_ran_dirichlet_pdf() and gsl_ran_dirichlet_lnpdf() functions. I've > not updated the docs yet, as I wanted to gage interest first. > > The code compiles and the tests run successfully on my computer. I've > also used the test code as a mini-benchmark, and the changed code > takes basically the same amount of time to execute (median = 1.0038 > times GSL 1.15, 95%CI=1.0009-1.0064). > > Let me know what would make this more useful. > > Sam