Date: Fri, 16 May 2003 12:18:03 -0700
From: roon8505 AT uidaho DOT edu
Subject: Long Doubles and Non-uniform Random Numbers
To: djgpp AT delorie DOT com
Cc: 
Message-id: <4f85ed4fcbff.4fcbff4f85ed@uidaho.edu>
MIME-version: 1.0
X-Mailer: iPlanet Messenger Express 5.2 HotFix 1.12 (built Feb 13 2003)
Content-type: multipart/mixed; boundary="Boundary_(ID_diQP5OVHhvr4F6QbLJIL1A)"
Content-language: en
X-Accept-Language: en
Reply-To: djgpp AT delorie DOT com

This is a multi-part message in MIME format.

--Boundary_(ID_diQP5OVHhvr4F6QbLJIL1A)
Content-type: text/plain; charset=us-ascii
Content-transfer-encoding: 7BIT
Content-disposition: inline

Cheers.

I did manage to write a program which simulates Gamma and Beta rv's. 

There's one hitch with it at the moment. The process it iterative, requiring the calculation of some pretty miniscule values. 

When I run, (using long doubles) the values- apparently- get re-set to zero as soon as they fall below ~ 1e-324. I thought that long doubles could handle much smaller values than this- am I missing something?

Thanks-

David R 

--Boundary_(ID_diQP5OVHhvr4F6QbLJIL1A)
Content-type: message/rfc822

Return-path: <lnobody AT delorie DOT com>
Received: from conversion-daemon.mailB.its.uidaho.edu by mailB.its.uidaho.edu
 (Go Vandals!) id <0HEY00M01HY8BQ AT mailB DOT its DOT uidaho DOT edu> for
 roon8505 AT ims-ms-daemon; Thu, 15 May 2003 18:41:15 -0700 (PDT)
Received: from rain ([129.101.155.252]) by mailB.its.uidaho.edu (Go Vandals!)
 with ESMTPP id <0HEY003I1I0QTU AT mailB DOT its DOT uidaho DOT edu> for
 roon8505 AT ims-ms-daemon (ORCPT roon8505 AT uidaho DOT edu); Thu,
 15 May 2003 18:41:14 -0700 (PDT)
Received: from conversion-daemon.mailA.its.uidaho.edu by mailA.its.uidaho.edu
 (Go Vandals!) id <0HEY00101I0GOX AT mailA DOT its DOT uidaho DOT edu> for roon8505 AT uidaho DOT edu
 (ORCPT roon8505 AT uidaho DOT edu); Thu, 15 May 2003 18:41:14 -0700 (PDT)
Received: from delorie.com ([207.22.48.162])
 by mailA.its.uidaho.edu (Go Vandals!)
 with ESMTPP id <0HEY00E9LI0PNN AT mailA DOT its DOT uidaho DOT edu> for roon8505 AT uidaho DOT edu
 (ORCPT roon8505 AT uidaho DOT edu); Thu, 15 May 2003 18:41:14 -0700 (PDT)
Received: (from mailnull AT localhost)	by delorie.com (8.11.6/8.9.1)
 id h4G1R3S32215	for djgpp-list; Thu, 15 May 2003 21:27:03 -0400
Received: from sccmmhc01.mchsi.com (root AT localhost)
	by delorie.com (8.11.6/8.9.1) with SMTP id h4G1R2232203	for
 <djgpp AT delorie DOT com>; Thu, 15 May 2003 21:27:02 -0400
Received: from sccmmhc01.mchsi.com (sccmmhc01.mchsi.com [204.127.203.183])
	by delorie.com (antispam)
Received: from cmc3092889b (12-218-74-57.client.mchsi.com[12.218.74.57])
 by sccmmhc01.mchsi.com (sccmmhc01) with SMTP id <20030516012659mm1004b7kae>;
 Fri, 16 May 2003 01:26:59 +0000
Date: Thu, 15 May 2003 20:27:02 -0500
From: Wayne <dockeen AT mchsi DOT com>
Subject: Non-uniform Random Numbers
To: roon8505 AT uidaho DOT edu
Cc: djgpp AT delorie DOT com
Errors-to: nobody AT delorie DOT com
Message-id: <IHEAIGJEFNKKKDPFCPECMEGOCBAA DOT dockeen AT mchsi DOT com>
MIME-version: 1.0
X-MIMEOLE: Produced By Microsoft MimeOLE V6.00.2800.1165
X-Mailer: Microsoft Outlook IMO, Build 9.0.2416 (9.0.2910.0)
Content-type: text/plain; charset=iso-8859-1
Content-transfer-encoding: 7BIT
Importance: Normal
X-Priority: 3 (Normal)
X-MSMail-priority: Normal
Precedence: bulk
X-Unsubscribes-To: listserv AT delorie DOT com
X-Mailing-List: djgpp AT delorie DOT com
Original-recipient: rfc822;roon8505 AT uidaho DOT edu

Check out GSL, the GNU Scientific Library.
I had a link for GSL for djgpp on my other
computer, I'll see if I can hunt it up
again..

Wayne


--Boundary_(ID_diQP5OVHhvr4F6QbLJIL1A)--