diff options
Diffstat (limited to 'research/experiments/044-specs')
| -rw-r--r-- | research/experiments/044-specs/2026-10-18-stage2c.jsonl | 6 | ||||
| -rw-r--r-- | research/experiments/044-specs/json_rec.py | 159 | ||||
| -rw-r--r-- | research/experiments/044-specs/json_rec_time.c | 119 |
3 files changed, 284 insertions, 0 deletions
diff --git a/research/experiments/044-specs/2026-10-18-stage2c.jsonl b/research/experiments/044-specs/2026-10-18-stage2c.jsonl new file mode 100644 index 0000000..99972d6 --- /dev/null +++ b/research/experiments/044-specs/2026-10-18-stage2c.jsonl @@ -0,0 +1,6 @@ +{"size":103781,"jit":"0","batch":"","c_ns_per_byte":13.050,"pure_ns_per_byte":599.753,"ratio":46.0,"instructions":25733137,"instructions_per_byte":248.0} +{"size":103781,"python_ns_per_byte":367.8,"python_json_module_ns_per_byte":54.0} +{"size":103781,"jit":"0","batch":"","c_ns_per_byte":13.814,"pure_ns_per_byte":600.877,"ratio":43.5,"instructions":25733137,"instructions_per_byte":248.0} +{"size":103781,"python_ns_per_byte":365.5,"python_json_module_ns_per_byte":53.8} +{"size":103781,"jit":"0","batch":"","c_ns_per_byte":13.141,"pure_ns_per_byte":598.454,"ratio":45.5,"instructions":25733137,"instructions_per_byte":248.0} +{"size":103781,"python_ns_per_byte":369.3,"python_json_module_ns_per_byte":54.8} diff --git a/research/experiments/044-specs/json_rec.py b/research/experiments/044-specs/json_rec.py new file mode 100644 index 0000000..c7951e5 --- /dev/null +++ b/research/experiments/044-specs/json_rec.py @@ -0,0 +1,159 @@ +# A JSON text read into a tree by recursive descent and walked depth first for a digest: the design of aup's +# examples/json_tree.aup in 027's stage 2c, in plain Python (the reference 140 counts its lines against; the tree is +# Python's own values, an object kept as pairs). A text that is no JSON raises (ValueError, or IndexError past its +# end), as the event fails. +import re + + +class Pairs(list): + pass + + +def space(b): + return b in (0x20, 0x09, 0x0a, 0x0d) + + +def digit(b): + return 0x30 <= b <= 0x39 + + +def hexdigit(b): + return digit(b) or 0x61 <= b <= 0x66 or 0x41 <= b <= 0x46 + + +def require(ok): + if not ok: + raise ValueError + + +def code(d, i): + return int(d[i:i + 4], 16) + + +def escape(d, i): + return d[i] == 0x5c and d[i + 1] == 0x75 and all(hexdigit(d[k]) for k in range(i + 2, i + 6)) + + +def skip(d, i): + while i < len(d) and space(d[i]): + i += 1 + return i + + +def string(d, i): + require(d[i] == 0x22) + s, j = bytearray(), i + 1 + while d[j] != 0x22: + start = j + while d[j] != 0x22 and d[j] != 0x5c and d[j] >= 0x20: + j += 1 + s += d[start:j] + require(d[j] >= 0x20) + if d[j] == 0x5c: + e = d[j + 1] + if e == 0x75: + require(escape(d, j)) + u = code(d, j + 2) + j += 6 + if 0xd800 <= u < 0xdc00 and j + 5 < len(d) and escape(d, j) and 0xdc00 <= code(d, j + 2) < 0xe000: + u = 0x10000 + ((u - 0xd800) << 10 | code(d, j + 2) - 0xdc00) + j += 6 + s += chr(u).encode() # a lone surrogate raises + else: + require(e in b'"\\/bfnrt') + s += {0x62: b'\b', 0x66: b'\f', 0x6e: b'\n', 0x72: b'\r', 0x74: b'\t'}.get(e, bytes([e])) + j += 2 + return s.decode(), j + 1 + + +def number(t): + require(re.fullmatch(r'-?(0|[1-9]\d*)(\.\d+)?([eE][+-]?\d+)?', t)) + x = float(t) + require(abs(x) != float('inf')) + return x + + +def value(d, i): + j = skip(d, i) + b = d[j] + if b in (0x5b, 0x7b): + items = [] if b == 0x5b else Pairs() + j = skip(d, j + 1) + if d[j] != b + 2: + while True: + if b == 0x5b: + v, k = value(d, j) + items.append(v) + else: + name, k = string(d, j) + k = skip(d, k) + require(d[k] == 0x3a) + v, k = value(d, k + 1) + items.append((name, v)) + j = skip(d, k) + if d[j] != 0x2c: + break + j = skip(d, j + 1) + require(d[j] == b + 2) + return items, j + 1 + if b == 0x22: + return string(d, j) + if b == 0x2d or digit(b): + start, j = j, j + 1 + while j < len(d) and (digit(d[j]) or d[j] in b'.eE+-'): + j += 1 + return number(d[start:j].decode()), j + for word, v in ((b'true', True), (b'false', False), (b'null', None)): + if d[j:j + len(word)] == word: + return v, j + len(word) + raise ValueError + + +def walk(x, d): + if x is None: d['nulls'] += 1 + elif x is True: d['trues'] += 1 + elif x is False: d['falses'] += 1 + elif isinstance(x, float): d['numbers'] += 1; d['sum'] += x + elif isinstance(x, str): d['strings'] += 1; d['text_bytes'] += len(x.encode()) + elif isinstance(x, Pairs): + d['objects'] += 1 + for k, v in x: + d['key_bytes'] += len(k.encode()) + walk(v, d) + else: + d['arrays'] += 1 + for y in x: + walk(y, d) + + +def check(data): + data.decode() # UTF-8, as text(check.data) is + root, end = value(data, 0) + require(skip(data, end) == len(data)) + d = dict(nulls=0, trues=0, falses=0, numbers=0, strings=0, arrays=0, objects=0, sum=0.0, text_bytes=0, key_bytes=0) + walk(root, d) + return d + + +if __name__ == '__main__': + import json, os, sys, time + sys.setrecursionlimit(100000) + suite = sys.argv[1] + for name in sorted(os.listdir(suite)): + try: + got = check(open(os.path.join(suite, name), 'rb').read()) + except (ValueError, IndexError, RecursionError): + got = None + assert not name.startswith('y_') or got is not None, name + assert not name.startswith('n_') or got is None, name + big = b'[' + b','.join(b'{"key%d": [1.5e3, true, null, "text \\u00e9 %d"]}' % (k, k) for k in range(2000)) + b']' + runs, loads = [], [] + for _ in range(7): + start = time.perf_counter() + check(big) + runs.append((time.perf_counter() - start) * 1e9 / len(big)) + start = time.perf_counter() + d = dict(nulls=0, trues=0, falses=0, numbers=0, strings=0, arrays=0, objects=0, sum=0.0, text_bytes=0, key_bytes=0) + walk(json.loads(big, object_pairs_hook=Pairs, parse_int=float), d) + loads.append((time.perf_counter() - start) * 1e9 / len(big)) + print('{"size":%d,"python_ns_per_byte":%.1f,"python_json_module_ns_per_byte":%.1f}' % (len(big), sorted(runs)[3], sorted(loads)[3])) diff --git a/research/experiments/044-specs/json_rec_time.c b/research/experiments/044-specs/json_rec_time.c new file mode 100644 index 0000000..03c02c4 --- /dev/null +++ b/research/experiments/044-specs/json_rec_time.c @@ -0,0 +1,119 @@ +/* Research 140, design 027 stage 2c: aup's examples/json_tree.aup (run by + * Pure: a JSON text read into a tree by recursive descent, then a digest of + * the tree, depth first) against C: + * cJSON (the system's libcjson.so.1, its few declarations written here) + * reading the same document into its tree, and the same digest. Both digests + * must be equal; then each is timed (the median of seven runs of rounds). + * + * cc -O2 -std=c17 -I../../../core json_rec_time.c ../../../build/libpure.a -l:libcjson.so.1 -lm -pthread -o json_rec_time + * ./json_rec_time json_tree.pmod ROUNDS + * + * The document is json_time.c's: 2,000 objects of a key and an array of a + * number, true, null and a string with an escape. One JSON line. */ +#define _POSIX_C_SOURCE 200809L +#include "module.h" +#include <stdio.h> +#include <stdlib.h> +#include <string.h> +#include <time.h> + +typedef struct cJSON { + struct cJSON *next, *prev, *child; + int type; + char *valuestring; + int valueint; + double valuedouble; + char *string; +} cJSON; +cJSON *cJSON_Parse(const char *value); +void cJSON_Delete(cJSON *item); +enum { FALSE_ = 1, TRUE_ = 2, NULL_ = 4, NUMBER = 8, STRING = 16, ARRAY = 32, OBJECT = 64 }; + +/* nulls, trues, falses, numbers, strings, arrays, objects, sum, text bytes, key bytes: depth first, as json_tree.aup. */ +typedef struct { double count[7], sum, text_bytes, key_bytes; } Digest; +static void walk(const cJSON *x, Digest *d) { + switch (x->type & 0xFF) { + case NULL_: d->count[0]++; break; + case TRUE_: d->count[1]++; break; + case FALSE_: d->count[2]++; break; + case NUMBER: d->count[3]++, d->sum += x->valuedouble; break; + case STRING: d->count[4]++, d->text_bytes += (double)strlen(x->valuestring); break; + case ARRAY: d->count[5]++; for (const cJSON *y = x->child; y; y = y->next) walk(y, d); break; + default: d->count[6]++; for (const cJSON *y = x->child; y; y = y->next) d->key_bytes += (double)strlen(y->string), walk(y, d); + } +} +static Digest digest(const char *text, cJSON **queue) { + (void)queue; + Digest d = {{0}, 0, 0, 0}; + cJSON *root = cJSON_Parse(text); + walk(root, &d); + cJSON_Delete(root); + return d; +} + +static double now_ns(void) { + struct timespec t; + clock_gettime(CLOCK_MONOTONIC, &t); + return (double)t.tv_sec * 1e9 + (double)t.tv_nsec; +} + +static int by_value(const void *a, const void *b) { + double x = *(const double *)a, y = *(const double *)b; + return (x > y) - (x < y); +} + +int main(int argc, char **argv) { + if (argc != 3) { fprintf(stderr, "usage: %s MODULE ROUNDS\n", argv[0]); return 2; } + int rounds = atoi(argv[2]); + static char doc[1 << 18]; + static cJSON *queue[1 << 16]; + size_t n = 0; + doc[n++] = '['; + for (int k = 0; k < 2000; ++k) + n += (size_t)sprintf(doc + n, "%s{\"key%d\": [1.5e3, true, null, \"text \\u00e9 %d\"]}", k ? "," : "", k, k); + doc[n++] = ']'; + doc[n] = 0; + PureModule module; + char error[160]; + if (pure_module_read(argv[1], &module, error, sizeof(error)) != PURE_OK) { fprintf(stderr, "%s\n", error); return 1; } + PureApp app; + if (pure_app_start(&app, &module, (PureData){0}) != PURE_OK) return 1; + int check = pure_module_input(&module, "check"); + if (pure_app_message(&app, (uint32_t)check, NULL, -1, (PureSpan){(const unsigned char *)doc, n}) != PURE_OK) { + fprintf(stderr, "event failed\n"); + return 1; + } + uint64_t work = app.work; + Digest c = digest(doc, queue); + /* The states after valid: nulls, trues, falses, numbers, strings, arrays, objects, sum, text_bytes, key_bytes. */ + const char *names[] = {"nulls", "trues", "falses", "numbers", "strings", "arrays", "objects", "sum", "text_bytes", "key_bytes"}; + double want[] = {c.count[0], c.count[1], c.count[2], c.count[3], c.count[4], c.count[5], c.count[6], c.sum, c.text_bytes, c.key_bytes}; + for (unsigned k = 0; k < 10; ++k) { + int state = -1; + for (uint32_t s = 0; s < module.states; ++s) if (!strcmp(module.state_names[s], names[k])) state = (int)s; + if (state < 0 || app.values[state].num != want[k]) { + fprintf(stderr, "%s: %g against %g\n", names[k], state < 0 ? -1 : app.values[state].num, want[k]); + return 1; + } + } + double c_runs[7], pure_runs[7]; + volatile double sink = 0; + for (int r = 0; r < 7; ++r) { + double start = now_ns(); + for (int i = 0; i < rounds; ++i) sink += digest(doc, queue).sum; + c_runs[r] = (now_ns() - start) / rounds; + start = now_ns(); + for (int i = 0; i < rounds; ++i) + if (pure_app_message(&app, (uint32_t)check, NULL, -1, (PureSpan){(const unsigned char *)doc, n}) != PURE_OK) return 1; + pure_runs[r] = (now_ns() - start) / rounds; + } + qsort(c_runs, 7, sizeof(double), by_value); + qsort(pure_runs, 7, sizeof(double), by_value); + const char *jit = getenv("PURE_JIT"), *batch = getenv("PURE_BATCH"); + printf("{\"size\":%zu,\"jit\":\"%s\",\"batch\":\"%s\",\"c_ns_per_byte\":%.3f,\"pure_ns_per_byte\":%.3f,\"ratio\":%.1f," + "\"instructions\":%llu,\"instructions_per_byte\":%.1f}\n", + n, jit ? jit : "", batch ? batch : "", c_runs[3] / (double)n, pure_runs[3] / (double)n, pure_runs[3] / c_runs[3], + (unsigned long long)work, (double)work / (double)n); + pure_app_stop(&app); + return 0; +} |
