From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 30887 invoked by alias); 9 Jul 2012 21:05:40 -0000 Received: (qmail 30875 invoked by uid 22791); 9 Jul 2012 21:05:39 -0000 X-SWARE-Spam-Status: No, hits=-4.6 required=5.0 tests=AWL,BAYES_00,KHOP_RCVD_UNTRUST,KHOP_THREADED,RCVD_IN_DNSWL_LOW,RCVD_IN_HOSTKARMA_NO,RCVD_IN_HOSTKARMA_W,RCVD_IN_HOSTKARMA_WL,RCVD_IN_HOSTKARMA_YE,SPF_NEUTRAL X-Spam-Check-By: sourceware.org Received: from va3ehsobe005.messaging.microsoft.com (HELO va3outboundpool.messaging.microsoft.com) (216.32.180.31) by sourceware.org (qpsmtpd/0.43rc1) with ESMTP; Mon, 09 Jul 2012 21:05:25 +0000 Received: from mail100-va3-R.bigfish.com (10.7.14.242) by VA3EHSOBE011.bigfish.com (10.7.40.61) with Microsoft SMTP Server id 14.1.225.23; Mon, 9 Jul 2012 21:03:07 +0000 Received: from mail100-va3 (localhost [127.0.0.1]) by mail100-va3-R.bigfish.com (Postfix) with ESMTP id 84C26A00A6 for ; Mon, 9 Jul 2012 21:03:06 +0000 (UTC) X-Forefront-Antispam-Report: CIP:157.55.9.135;KIP:(null);UIP:(null);IPV:NLI;H:DB3PRD0104HT014.eurprd01.prod.exchangelabs.com;RD:none;EFVD:NLI X-SpamScore: -5 X-BigFish: PS-5(zf7Izbb2dI98dI9371I1432I4015Izz1202hzz8275dhz2ei2a8h668h839h8e3h93fhf0ah) Received-SPF: neutral (mail100-va3: 157.55.9.135 is neither permitted nor denied by domain of warwick.ac.uk) client-ip=157.55.9.135; envelope-from=S.A.Mason@warwick.ac.uk; helo=DB3PRD0104HT014.eurprd01.prod.exchangelabs.com ;angelabs.com ; Received: from mail100-va3 (localhost.localdomain [127.0.0.1]) by mail100-va3 (MessageSwitch) id 1341867784399106_28965; Mon, 9 Jul 2012 21:03:04 +0000 (UTC) Received: from VA3EHSMHS032.bigfish.com (unknown [10.7.14.250]) by mail100-va3.bigfish.com (Postfix) with ESMTP id 541003600AE for ; Mon, 9 Jul 2012 21:03:04 +0000 (UTC) Received: from DB3PRD0104HT014.eurprd01.prod.exchangelabs.com (157.55.9.135) by VA3EHSMHS032.bigfish.com (10.7.99.42) with Microsoft SMTP Server (TLS) id 14.1.225.23; Mon, 9 Jul 2012 21:03:04 +0000 Received: from mail-yx0-f169.google.com (209.85.213.169) by pod51002.outlook.com (10.4.156.44) with Microsoft SMTP Server (TLS) id 14.15.86.1; Mon, 9 Jul 2012 21:05:19 +0000 Received: by yenr5 with SMTP id r5so12226101yen.0 for ; Mon, 09 Jul 2012 14:05:15 -0700 (PDT) Received: by 10.60.30.35 with SMTP id p3mr41670578oeh.16.1341867915092; Mon, 09 Jul 2012 14:05:15 -0700 (PDT) MIME-Version: 1.0 Received: by 10.76.76.196 with HTTP; Mon, 9 Jul 2012 14:04:54 -0700 (PDT) In-Reply-To: References: From: Sam Mason Date: Mon, 09 Jul 2012 21:05:00 -0000 Message-ID: Subject: Re: Fwd: log probability density functions (PDFs) To: Patrick Alken CC: "gsl-discuss@sourceware.org" Content-Type: text/plain; charset="UTF-8" X-OriginatorOrg: warwick.ac.uk 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/msg00002.txt.bz2 Hi Patrick, I've put the following together: http://www.samason.me.uk/hacking-the-gsl Is that the sort of thing you'd consider linking to? Maybe with an accompanying sentence of "log-probability density functions"? Thanks, Sam On 9 July 2012 18:39, Patrick Alken wrote: > 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 > >