1 /* $Id: closure.c,v 1.10 2014/02/19 00:45:42 Tom.Shields Exp $ */
9 static unsigned *first_derives;
13 static void print_closure(int);
14 static void print_EFF(void);
15 static void print_first_derives(void);
28 rowsize = WORDSIZE(nvars);
29 EFF = NEW2(nvars * rowsize, unsigned);
32 for (i = start_symbol; i < nsyms; i++)
35 for (rule = *sp; rule > 0; rule = *++sp)
37 symbol = ritem[rrhs[rule]];
40 symbol -= start_symbol;
47 reflexive_transitive_closure(EFF, nvars);
55 set_first_derives(void)
69 rulesetsize = WORDSIZE(nrules);
70 varsetsize = WORDSIZE(nvars);
71 first_derives = NEW2(nvars * rulesetsize, unsigned) - ntokens * rulesetsize;
75 rrow = first_derives + ntokens * rulesetsize;
76 for (i = start_symbol; i < nsyms; i++)
78 vrow = EFF + ((i - ntokens) * varsetsize);
80 for (j = start_symbol; j < nsyms; k++, j++)
82 if (k >= BITS_PER_WORD)
88 if (cword & (unsigned)(1 << k))
91 while ((rule = *rp++) >= 0)
102 print_first_derives();
109 closure(Value_t *nucleus, int n)
124 rulesetsize = WORDSIZE(nrules);
125 rsend = ruleset + rulesetsize;
126 for (rsp = ruleset; rsp < rsend; rsp++)
130 for (csp = nucleus; csp < csend; ++csp)
132 symbol = ritem[*csp];
135 dsp = first_derives + symbol * rulesetsize;
143 itemsetend = itemset;
145 for (rsp = ruleset; rsp < rsend; ++rsp)
150 for (i = 0; i < BITS_PER_WORD; ++i)
152 if (word & (unsigned)(1 << i))
154 itemno = rrhs[ruleno + i];
155 while (csp < csend && *csp < itemno)
156 *itemsetend++ = *csp++;
157 *itemsetend++ = itemno;
158 while (csp < csend && *csp == itemno)
163 ruleno += BITS_PER_WORD;
167 *itemsetend++ = *csp++;
175 finalize_closure(void)
179 FREE(first_derives + ntokens * WORDSIZE(nrules));
189 printf("\n\nn = %d\n\n", n);
190 for (isp = itemset; isp < itemsetend; isp++)
191 printf(" %d\n", *isp);
202 printf("\n\nEpsilon Free Firsts\n");
204 for (i = start_symbol; i < nsyms; i++)
206 printf("\n%s", symbol_name[i]);
207 rowp = EFF + ((i - start_symbol) * WORDSIZE(nvars));
211 for (j = 0; j < nvars; k++, j++)
213 if (k >= BITS_PER_WORD)
220 printf(" %s", symbol_name[start_symbol + j]);
226 print_first_derives(void)
234 printf("\n\n\nFirst Derives\n");
236 for (i = start_symbol; i < nsyms; i++)
238 printf("\n%s derives\n", symbol_name[i]);
239 rp = first_derives + i * WORDSIZE(nrules);
241 for (j = 0; j <= nrules; k++, j++)
243 if (k >= BITS_PER_WORD)
249 if (cword & (1 << k))