summaryrefslogtreecommitdiff
path: root/vere/pkg/noun/jets/b/sort.c
blob: ea5507c01576f309714e95c3e868aa76eada8b4e (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
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
/// @file

#include "jets/q.h"
#include "jets/w.h"

#include "noun.h"


static_assert( (UINT32_MAX > u3a_cells),
               "length precision" );

static_assert(
  (UINT32_MAX < (SIZE_MAX / (2 * sizeof(u3_noun)))),
  "len_w * sizeof u3_noun overflow"
);

static void
_merge_sort(u3_noun* restrict arr_u,
  u3_noun* restrict tmp_u,
  c3_w lef_w,
  c3_w rit_w,
  u3j_site* sit_u)
{
  if ( lef_w >= rit_w ) return;
  c3_w mid_w =  (lef_w + rit_w) / 2;
  if (mid_w < lef_w)
  {
    //  addition wrapped around
    //
    u3m_bail(c3__fail);
  }
  _merge_sort(arr_u, tmp_u, lef_w, mid_w, sit_u);
  _merge_sort(arr_u, tmp_u, mid_w + 1, rit_w, sit_u);

  c3_w i_w = lef_w, j_w = mid_w + 1, k_w = lef_w;
  while (i_w <= mid_w && j_w <= rit_w)
  {
    //  reversed comparison to mimick order reversal of pivot
    //  and compared element in Hoon
    //
    u3_noun sam = u3nc(u3k(arr_u[j_w]), u3k(arr_u[i_w]));
    c3_o hoz_o = u3x_loob(u3j_gate_slam(sit_u, sam));
    if ( c3n == hoz_o )
      tmp_u[k_w++] = arr_u[i_w++];
    else
      tmp_u[k_w++] = arr_u[j_w++];
  }

  while (i_w <= mid_w) tmp_u[k_w++] = arr_u[i_w++];
  while (j_w <= rit_w) tmp_u[k_w++] = arr_u[j_w++];

  for (i_w = lef_w; i_w <= rit_w; i_w++)
  {
    arr_u[i_w] = tmp_u[i_w];
  }
}

//  RETAINS list, transfers product
//
static u3_noun
_sort(u3j_site* sit_u, u3_noun list)
{
  if (u3_nul == list) return u3_nul;

  c3_w len_w = 1;
  {
    u3_noun t = u3t(list);
    while (u3_nul != t)
    {
      ++len_w; t = u3t(t);
    }
  }

  if (1 == len_w) return u3k(list);
  u3_noun* arr_u = u3a_malloc(sizeof(u3_noun) * len_w * 2);
  u3_noun* tmp_u = arr_u + len_w;
  for (c3_w i_w = 0; i_w < len_w; i_w++)
  {
    //  inlined u3r_cell without any checks
    //  since the list was already validated and measured
    //
    u3a_cell* cel_u = u3a_to_ptr(list);
    arr_u[i_w] = u3k(cel_u->hed);
    list = cel_u->tel;
  }

  _merge_sort(arr_u, tmp_u, 0, len_w - 1, sit_u);

  u3_noun pro = u3_nul;
  for (c3_w i_w = len_w; i_w--;)
  {
    pro = u3nc(arr_u[i_w], pro);
  }
  
  u3a_free(arr_u);

  return pro;
}

u3_noun
u3qb_sort(u3_noun a,
          u3_noun b)
{
  u3_noun  pro;
  u3j_site sit_u;
  u3j_gate_prep(&sit_u, u3k(b));
  pro = _sort(&sit_u, a);
  u3j_gate_lose(&sit_u);
  return pro;
}

u3_noun
u3wb_sort(u3_noun cor)
{
  u3_noun a, b;

  if ( c3n == u3r_mean(cor, u3x_sam_2, &a, u3x_sam_3, &b, 0) )
  {
    return u3m_bail(c3__exit);
  }
  else
  {
    return u3qb_sort(a, b);
  }
}