summaryrefslogtreecommitdiffstats
path: root/util.c
blob: a2573a4a2d4f3f5d6af19868a1ec3630e26e6130 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
/*********************************************************************
 * kd-forest                                                         *
 * Copyright (C) 2014 Tavian Barnes <tavianator@tavianator.com>      *
 *                                                                   *
 * This program is free software. It comes without any warranty, to  *
 * the extent permitted by applicable law. You can redistribute it   *
 * and/or modify it under the terms of the Do What The Fuck You Want *
 * To Public License, Version 2, as published by Sam Hocevar. See    *
 * the COPYING file or http://www.wtfpl.net/ for more details.       *
 *********************************************************************/

#include "util.h"
#include <stdlib.h>

void *
xmalloc(size_t size)
{
  void *ret = malloc(size);
  if (!ret) {
    abort();
  }
  return ret;
}

void *
xrealloc(void *ptr, size_t size)
{
  void *ret = realloc(ptr, size);
  if (!ret) {
    abort();
  }
  return ret;
}

// Based on sample rand() implementation from POSIX.1-2001

static unsigned long xrand_next = 0;

void xsrand(unsigned int seed) {
  xrand_next = seed;
}

static unsigned int xrand_simple(void) {
  xrand_next = xrand_next*1103515245 + 12345;
  return (unsigned int)(xrand_next/65536)%32768;
}

static unsigned int xrand_full(void) {
  unsigned int low = xrand_simple();
  unsigned int high = xrand_simple();
  return low | (high << 15);
}

#define XRAND_RANGE 1073741824U

unsigned int
xrand(unsigned int range)
{
  // Compensate for bias if XRAND_RANGE isn't a multiple of range
  unsigned int limit = XRAND_RANGE - XRAND_RANGE%range;
  unsigned int res;
  do {
    res = xrand_full();
  } while (res >= limit);
  return res%range;
}