🚀 go-pugleaf

I coded this up from an article I think in EDN, by an engineer at 
Tektronix.  The code in the article was in some DSP assembler, and I
converted it to C (isn't that backwards??).  Anyway, it works.  I tested
it on a DSP32c output to a D/A and to a spectrum analyzer.

-----BEGIN RAN32.C-----

#include <math.h>

/* generate 32-bit random noise using DSP32c */

void ran32 (ioval)
unsigned long *ioval;
{

#define  C    0x234567
#define  Alo  0x007465
#define  Ahi  0x000010

  unsigned long Rlo, Rhi;
  unsigned long ar1, ar2, ar3;
  unsigned long t1, t2, t3;

  Rlo = *ioval & 0x00ffff;
  Rhi = *ioval >> 16;
  ar1 = Alo * Rlo;
  ar2 = Alo * Rhi;
  ar3 = Ahi * Rlo;
  t1 = (ar2 + ar3);
  t2 = t1 << 16;
  t3 = ar1 + t2;
  *ioval = C + t3;
}