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
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
|
//
//
// The Nimrod Compiler
// (c) Copyright 2009 Andreas Rumpf
//
// See the file "copying.txt", included in this
// distribution, for details about the copyright.
//
unit ropes;
{ Ropes for the C code generator
Ropes are a data structure that represents a very long string
efficiently; especially concatenation is done in O(1) instead of O(N).
Ropes make use a lazy evaluation: They are essentially concatenation
trees that are only flattened when converting to a native Nimrod
string or when written to disk. The empty string is represented by a
nil pointer.
A little picture makes everything clear:
"this string" & " is internally " & "represented as"
con -- inner nodes do not contain raw data
/ \
/ \
/ \
con "represented as"
/ \
/ \
/ \
/ \
/ \
"this string" " is internally "
Note that this is the same as:
"this string" & (" is internally " & "represented as")
con
/ \
/ \
/ \
"this string" con
/ \
/ \
/ \
/ \
/ \
" is internally " "represented as"
The 'con' operator is associative! This does not matter however for
the algorithms we use for ropes.
Note that the left and right pointers are not needed for leafs.
Leafs have relatively high memory overhead (~30 bytes on a 32
bit machines) and we produce many of them. This is why we cache and
share leafs accross different rope trees.
To cache them they are inserted in another tree, a splay tree for best
performance. But for the caching tree we use the leafs' left and right
pointers.
}
interface
{$include 'config.inc'}
uses
nsystem, msgs, strutils, platform, nhashes, crc;
const
CacheLeafs = true;
countCacheMisses = False; // see what our little optimization gives
type
TFormatStr = string;
// later we may change it to CString for better
// performance of the code generator (assignments copy the format strings
// though it is not necessary)
PRope = ^TRope;
TRope = object(NObject)
left, right: PRope;
len: int;
data: string; // != nil if a leaf
end {@acyclic};
// the empty rope is represented by nil to safe space
TRopeSeq = array of PRope;
function con(a, b: PRope): PRope; overload;
function con(a: PRope; const b: string): PRope; overload;
function con(const a: string; b: PRope): PRope; overload;
function con(a: array of PRope): PRope; overload;
procedure app(var a: PRope; b: PRope); overload;
procedure app(var a: PRope; const b: string); overload;
procedure prepend(var a: PRope; b: PRope);
function toRope(const s: string): PRope; overload;
function toRopeF(const r: BiggestFloat): PRope;
function toRope(i: BiggestInt): PRope; overload;
function ropeLen(a: PRope): int;
procedure WriteRope(head: PRope; const filename: string);
function writeRopeIfNotEqual(r: PRope; const filename: string): boolean;
function ropeToStr(p: PRope): string;
function ropef(const frmt: TFormatStr; const args: array of PRope): PRope;
procedure appf(var c: PRope; const frmt: TFormatStr;
const args: array of PRope);
function getCacheStats: string;
function RopeEqualsFile(r: PRope; const f: string): Boolean;
// returns true if the rope r is the same as the contents of file f
function RopeInvariant(r: PRope): Boolean;
// exported for debugging
implementation
function ropeLen(a: PRope): int;
begin
if a = nil then result := 0
else result := a.len
end;
function newRope(const data: string = snil): PRope;
begin
new(result);
{@ignore}
fillChar(result^, sizeof(TRope), 0);
{@emit}
if data <> snil then begin
result.len := length(data);
result.data := data;
end
end;
// -------------- leaf cache: ---------------------------------------
var
cache: PRope; // the root of the cache tree
misses, hits: int;
N: PRope; // dummy rope needed for splay algorithm
function getCacheStats: string;
begin
if hits+misses <> 0 then
result := 'Misses: ' +{&} ToString(misses) +{&}
' total: ' +{&} toString(hits+misses) +{&}
' quot: ' +{&} toStringF(toFloat(misses) / toFloat(hits+misses))
else
result := ''
end;
function splay(const s: string; tree: PRope; out cmpres: int): PRope;
var
le, r, y, t: PRope;
c: int;
begin
t := tree;
N.left := nil; N.right := nil; // reset to nil
le := N;
r := N;
repeat
c := cmp(s, t.data);
if c < 0 then begin
if (t.left <> nil) and (s < t.left.data) then begin
y := t.left; t.left := y.right; y.right := t; t := y
end;
if t.left = nil then break;
r.left := t; r := t; t := t.left
end
else if c > 0 then begin
if (t.right <> nil) and (s > t.right.data) then begin
y := t.right; t.right := y.left; y.left := t; t := y
end;
if t.right = nil then break;
le.right := t; le := t; t := t.right
end
else break
until false;
cmpres := c;
le.right := t.left; r.left := t.right; t.left := N.right; t.right := N.left;
result := t
end;
function insertInCache(const s: string; tree: PRope): PRope;
// Insert i into the tree t, unless it's already there.
// Return a pointer to the resulting tree.
var
t: PRope;
cmp: int;
begin
t := tree;
if t = nil then begin
result := newRope(s);
if countCacheMisses then inc(misses);
exit
end;
t := splay(s, t, cmp);
if cmp = 0 then begin
// We get here if it's already in the Tree
// Don't add it again
result := t;
if countCacheMisses then inc(hits);
end
else begin
if countCacheMisses then inc(misses);
result := newRope(s);
if cmp < 0 then begin
result.left := t.left; result.right := t; t.left := nil
end
else begin // i > t.item:
result.right := t.right; result.left := t; t.right := nil
end
end
end;
function RopeInvariant(r: PRope): Boolean;
begin
if r = nil then
result := true
else begin
result := true
(*
if r.data <> snil then
result := true
else begin
result := (r.left <> nil) and (r.right <> nil);
if result then result := ropeInvariant(r.left);
if result then result := ropeInvariant(r.right);
end *)
end
end;
function toRope(const s: string): PRope;
begin
if s = '' then
result := nil
else if cacheLeafs then begin
result := insertInCache(s, cache);
cache := result;
end
else
result := newRope(s);
assert(RopeInvariant(result));
end;
// ------------------------------------------------------------------
procedure RopeSeqInsert(var rs: TRopeSeq; r: PRope; at: Natural);
var
len, i: int;
begin
len := length(rs);
if at > len then
SetLength(rs, at+1)
else
SetLength(rs, len+1);
// move old rope elements:
for i := len downto at+1 do
rs[i] := rs[i-1]; // this is correct, I used pen and paper to validate it
rs[at] := r
end;
function con(a, b: PRope): PRope; overload;
begin
assert(RopeInvariant(a));
assert(RopeInvariant(b));
if a = nil then // len is valid for every cord not only for leafs
result := b
else if b = nil then
result := a
else begin
result := newRope();
result.len := a.len + b.len;
result.left := a;
result.right := b
end;
assert(RopeInvariant(result));
end;
function con(a: PRope; const b: string): PRope; overload;
var
r: PRope;
begin
assert(RopeInvariant(a));
if b = '' then
result := a
else begin
r := toRope(b);
if a = nil then begin
result := r
end
else begin
result := newRope();
result.len := a.len + r.len;
result.left := a;
result.right := r;
end
end;
assert(RopeInvariant(result));
end;
function con(const a: string; b: PRope): PRope; overload;
var
r: PRope;
begin
assert(RopeInvariant(b));
if a = '' then
result := b
else begin
r := toRope(a);
if b = nil then
result := r
else begin
result := newRope();
result.len := b.len + r.len;
result.left := r;
result.right := b;
end
end;
assert(RopeInvariant(result));
end;
function con(a: array of PRope): PRope; overload;
var
i: int;
begin
result := nil;
for i := 0 to high(a) do result := con(result, a[i]);
assert(RopeInvariant(result));
end;
function toRope(i: BiggestInt): PRope;
begin
result := toRope(ToString(i))
end;
function toRopeF(const r: BiggestFloat): PRope;
begin
result := toRope(toStringF(r))
end;
procedure app(var a: PRope; b: PRope); overload;
begin
a := con(a, b);
assert(RopeInvariant(a));
end;
procedure app(var a: PRope; const b: string); overload;
begin
a := con(a, b);
assert(RopeInvariant(a));
end;
procedure prepend(var a: PRope; b: PRope);
begin
a := con(b, a);
assert(RopeInvariant(a));
end;
procedure InitStack(var stack: TRopeSeq);
begin
{@ignore}
setLength(stack, 0);
{@emit stack := @[];}
end;
procedure push(var stack: TRopeSeq; r: PRope);
var
len: int;
begin
len := length(stack);
setLength(stack, len+1);
stack[len] := r;
end;
function pop(var stack: TRopeSeq): PRope;
var
len: int;
begin
len := length(stack);
result := stack[len-1];
setLength(stack, len-1);
end;
procedure WriteRopeRec(var f: TTextFile; c: PRope);
begin
assert(RopeInvariant(c));
if c = nil then exit;
if (c.data <> snil) then begin
nimWrite(f, c.data)
end
else begin
writeRopeRec(f, c.left);
writeRopeRec(f, c.right)
end
end;
procedure newWriteRopeRec(var f: TTextFile; c: PRope);
var
stack: TRopeSeq;
it: PRope;
begin
assert(RopeInvariant(c));
initStack(stack);
push(stack, c);
while length(stack) > 0 do begin
it := pop(stack);
while it.data = snil do begin
push(stack, it.right);
it := it.left;
assert(it <> nil);
end;
assert(it.data <> snil);
nimWrite(f, it.data);
end
end;
procedure WriteRope(head: PRope; const filename: string);
var
f: TTextFile; // we use a textfile for automatic buffer handling
begin
if OpenFile(f, filename, fmWrite) then begin
if head <> nil then newWriteRopeRec(f, head);
nimCloseFile(f);
end
else
rawMessage(errCannotOpenFile, filename);
end;
procedure recRopeToStr(var result: string; var resultLen: int; p: PRope);
begin
if p = nil then exit; // do not add to result
if (p.data = snil) then begin
recRopeToStr(result, resultLen, p.left);
recRopeToStr(result, resultLen, p.right);
end
else begin
CopyMem(@result[resultLen+StrStart], @p.data[strStart], p.len);
Inc(resultLen, p.len);
assert(resultLen <= length(result));
end
end;
procedure newRecRopeToStr(var result: string; var resultLen: int;
r: PRope);
var
stack: TRopeSeq;
it: PRope;
begin
initStack(stack);
push(stack, r);
while length(stack) > 0 do begin
it := pop(stack);
while it.data = snil do begin
push(stack, it.right);
it := it.left;
end;
assert(it.data <> snil);
CopyMem(@result[resultLen+StrStart], @it.data[strStart], it.len);
Inc(resultLen, it.len);
assert(resultLen <= length(result));
end
end;
function ropeToStr(p: PRope): string;
var
resultLen: int;
begin
assert(RopeInvariant(p));
if p = nil then
result := ''
else begin
result := newString(p.len);
resultLen := 0;
newRecRopeToStr(result, resultLen, p);
end
end;
function ropef(const frmt: TFormatStr; const args: array of PRope): PRope;
var
i, j, len, start, num: int;
begin
i := strStart;
len := length(frmt);
result := nil;
num := 0;
while i <= len + StrStart - 1 do begin
if frmt[i] = '$' then begin
inc(i); // skip '$'
case frmt[i] of
'$': begin app(result, '$'+''); inc(i); end;
'#': begin inc(i); app(result, args[num]); inc(num); end;
'0'..'9': begin
j := 0;
repeat
j := (j*10) + Ord(frmt[i]) - ord('0');
inc(i);
until (i > len + StrStart - 1) or not (frmt[i] in ['0'..'9']);
num := j;
if j > high(args)+1 then
internalError('ropes: invalid format string $' + toString(j));
app(result, args[j-1]);
end;
'N', 'n': begin app(result, tnl); inc(i); end;
else InternalError('ropes: invalid format string $' + frmt[i]);
end
end;
start := i;
while (i <= len + StrStart - 1) do
if (frmt[i] <> '$') then inc(i) else break;
if i-1 >= start then begin
app(result, ncopy(frmt, start, i-1));
end
end;
assert(RopeInvariant(result));
end;
procedure appf(var c: PRope; const frmt: TFormatStr; const args: array of PRope);
begin
app(c, ropef(frmt, args))
end;
const
bufSize = 1024; // 1 KB is reasonable
function auxRopeEqualsFile(r: PRope; var bin: TBinaryFile;
buf: Pointer): Boolean;
var
readBytes: int;
begin
if (r.data <> snil) then begin
if r.len > bufSize then
// A token bigger than 1 KB? - This cannot happen in reality.
internalError('ropes: token too long');
readBytes := readBuffer(bin, buf, r.len);
result := (readBytes = r.len) // BUGFIX
and equalMem(buf, addr(r.data[strStart]), r.len);
end
else begin
result := auxRopeEqualsFile(r.left, bin, buf);
if result then
result := auxRopeEqualsFile(r.right, bin, buf);
end
end;
function RopeEqualsFile(r: PRope; const f: string): Boolean;
var
bin: TBinaryFile;
buf: Pointer;
begin
result := openFile(bin, f);
if not result then exit; // not equal if file does not exist
buf := alloc(BufSize);
result := auxRopeEqualsFile(r, bin, buf);
if result then
result := readBuffer(bin, buf, bufSize) = 0; // really at the end of file?
dealloc(buf);
CloseFile(bin);
end;
function crcFromRopeAux(r: PRope; startVal: TCrc32): TCrc32;
var
i: int;
begin
if r.data <> snil then begin
result := startVal;
for i := strStart to length(r.data)+strStart-1 do
result := updateCrc32(r.data[i], result);
end
else begin
result := crcFromRopeAux(r.left, startVal);
result := crcFromRopeAux(r.right, result);
end
end;
function newCrcFromRopeAux(r: PRope; startVal: TCrc32): TCrc32;
var
stack: TRopeSeq;
it: PRope;
L, i: int;
begin
initStack(stack);
push(stack, r);
result := startVal;
while length(stack) > 0 do begin
it := pop(stack);
while it.data = snil do begin
push(stack, it.right);
it := it.left;
end;
assert(it.data <> snil);
i := strStart;
L := length(it.data)+strStart;
while i < L do begin
result := updateCrc32(it.data[i], result);
inc(i);
end
end
end;
function crcFromRope(r: PRope): TCrc32;
begin
result := newCrcFromRopeAux(r, initCrc32)
end;
function writeRopeIfNotEqual(r: PRope; const filename: string): boolean;
// returns true if overwritten
var
c: TCrc32;
begin
c := crcFromFile(filename);
if c <> crcFromRope(r) then begin
writeRope(r, filename);
result := true
end
else
result := false
end;
initialization
new(N); // init dummy node for splay algorithm
{@ignore}
fillChar(N^, sizeof(N^), 0);
{@emit}
end.
|