Voyagers' Pseudorandom Number Generator on HP41CX

Programs for HP-41

Voyagers' Pseudorandom Number Generator on HP41CX

Postby mike-stgt » Thu Jun 25, 2020 11:00 pm

I ported the Voyagers' Pseudorandom Number Generator to the HP41CX. It could be useful when exactly the same sequence of random numbers is necessary for tests.

Warnings first
Indispensable prerequisite is an HP-41CX either as real machine or running under a so-called "emulator" (one of the CPU-simulating kind). Any former OS of 41C or 41CV will probably fail, highly likely to produce a MEMORY LOST or even worse, alter memory silently. There are no safety measures to prevent you to carry out this misuse.

Intended use
VORANOGE.zip contains two "4k ROM files" in different formats. The .MOD file may be used with V41 and the .BINJFG may be used with Emu41. See the emulators' description for details how to "plug-in application packs".

What else?
This module comes with only two commands, SEED and RN. SEED stores X in a buffer as sequence start for RN. In case there is no seed when executing RN time serves as seed (some fraction of seconds to be precise). RN computes with the seed a pseudorandom number and sets it as new seed.
Note: the seed buffer is volatile, PACKing or a power cycle (OFF-ON) will flush it.

Few more details
XROM: 18
FUNC: 3
VERS: FX20
IRVC: no interrupt vector
CKSM: 091 (145 dec) ok
FAT:
XROM 18,00 -VORANOGE B013 * NP
XROM 18,01 RN B11D
XROM 18,02 SEED B10E
Buffer-ID: 9

The fine print
VORANOGE is published under QPL, the Q Public License V 1.0 -- find it also in VORANOGE.zip


Final remark: this is not a request for comments or suggestions. But feel free to report errors.

/M.
Attachments
VORANOGE.zip
The a. m. ZIP
(4.17 KiB) Downloaded 42 times
mike-stgt
.......
.......
 
Posts: 38
Joined: Tue Dec 24, 2019 12:12 pm

Re: Voyagers' Pseudorandom Number Generator on HP41CX

Postby mike-stgt » Wed Jul 08, 2020 11:50 pm

Addendum
A very simple application using random numbers is the Monte Carlo method. Find below how to estimate the volume of a circle with r=1 (result should be Pi),
i) on an HP-11C
Code: Select all
 1-42.21.11  LBL A
 2-       4  4
 3-      15  1/x
 4-      12  e^x
 5-   44 36  SEED
 6-      26  EEX
 7-       4  4
 8-      36  ENTER^
 9-   44 25  STO I
10-42.21. 1  LBL 1
11-   42 36  RAN#
12-   43 11  x^2
13-   42 36  RAN#
14-   43 11  x^2
15-      40  +
16-   43 44  INT
17-      30  -
18-   42  5  DSE
19-   22  1  GTO 1
20-      34  x<>y
21-      10  /
22-       4  4
23-      20  *
24-   43 32  RTN
result = 3,140000

ii) same on an HP-41CX
PiMCds.png
Pi by Monte Carlo with distinct seed
PiMCds.png (870 Bytes) Viewed 558 times

result = 3,140000 (as expected -- when using the VoRaNoGe ROM)

iii) same once more on a 41CX running under Emu41
Pi by Monte Carlo.png

Note: I did all examples on virtual calculators like Emu41, V41 and -- for the 11C -- with HP's PC version of the HP-15C (ROM just replaced). The later is fast enough to experiment with larger iterations and different seeds.

BTW -- RN returns random numbers in the range 0 <= rn < 1. So in the a. m. example r is not 1, only close to 1. I leave this for you or the experts, to investigate if a correction is necessary.
mike-stgt
.......
.......
 
Posts: 38
Joined: Tue Dec 24, 2019 12:12 pm

Re: Voyagers' Pseudorandom Number Generator on HP41CX

Postby mike-stgt » Tue Aug 25, 2020 3:10 pm

One more use case for RN -- replacement of GN (Gaussian rN generator) found in PPC ROM.

Besides using my seed-in-buffer-RN instead of the register-based RN of PPC ROM I also like to get rid of the 'nasty surprise' (soft-crash) related to incorrect trig setting. The technical details in PPC ROM USR MNL show 'angular mode not used' what is only true in so far it is neither changed nor queried. But you have to set DEG before using GN, otherwise the results are (silently!) wrong. For this I first planned to set the "circular factor" according the set angular mode like this:
Code: Select all
 01 PI
 02 FC? 43
 03 R-D
 04 ,9
 05 FC? 42
 06 SIGN
 07 /

where lines 4..7 may be omitted by those who never-ever set their HP41 to GRAD.

Here a one byte shorter solution respecting angular setting including GRAD mode:
gn.png
gn.png (588 Bytes) Viewed 227 times

Used formula is the same as PPC ROM. The hindmost two lines, VIEW X and VIEW Y, are to ease transfer of the result by edit/copy - edit/paste when running on virtual equipment. It's up to you to delete them or set flag 21 according your preference. On an HP-42S I would replace R06 and R07 by named ones, 'standard deviation' and 'mean' or alike.
mike-stgt
.......
.......
 
Posts: 38
Joined: Tue Dec 24, 2019 12:12 pm


Return to HP-41 Software

Who is online

Users browsing this forum: No registered users and 1 guest

cron