]> git.lizzy.rs Git - plan9front.git/commitdiff
add qr.c
authorglenda <glenda@PC192-168-2-10>
Wed, 23 Sep 2015 14:36:55 +0000 (16:36 +0200)
committerglenda <glenda@PC192-168-2-10>
Wed, 23 Sep 2015 14:36:55 +0000 (16:36 +0200)
sys/src/cmd/qr.c [new file with mode: 0644]

diff --git a/sys/src/cmd/qr.c b/sys/src/cmd/qr.c
new file mode 100644 (file)
index 0000000..332da7e
--- /dev/null
@@ -0,0 +1,1185 @@
+#include <u.h>
+#include <libc.h>
+#include <ctype.h>
+
+enum {
+       NUM = 1,
+       ANUM = 2,
+       BYTE = 4,
+};
+
+static uchar gpexp[255] =
+{1, 2, 4, 8, 16, 32, 64, 128, 29, 58, 116, 232, 205, 135, 19, 38, 76, 
+152, 45, 90, 180, 117, 234, 201, 143, 3, 6, 12, 24, 48, 96, 192, 157, 
+39, 78, 156, 37, 74, 148, 53, 106, 212, 181, 119, 238, 193, 159, 35, 
+70, 140, 5, 10, 20, 40, 80, 160, 93, 186, 105, 210, 185, 111, 222, 
+161, 95, 190, 97, 194, 153, 47, 94, 188, 101, 202, 137, 15, 30, 60, 
+120, 240, 253, 231, 211, 187, 107, 214, 177, 127, 254, 225, 223, 163, 
+91, 182, 113, 226, 217, 175, 67, 134, 17, 34, 68, 136, 13, 26, 52, 
+104, 208, 189, 103, 206, 129, 31, 62, 124, 248, 237, 199, 147, 59, 
+118, 236, 197, 151, 51, 102, 204, 133, 23, 46, 92, 184, 109, 218, 
+169, 79, 158, 33, 66, 132, 21, 42, 84, 168, 77, 154, 41, 82, 164, 85, 
+170, 73, 146, 57, 114, 228, 213, 183, 115, 230, 209, 191, 99, 198, 
+145, 63, 126, 252, 229, 215, 179, 123, 246, 241, 255, 227, 219, 171, 
+75, 150, 49, 98, 196, 149, 55, 110, 220, 165, 87, 174, 65, 130, 25, 
+50, 100, 200, 141, 7, 14, 28, 56, 112, 224, 221, 167, 83, 166, 81, 
+162, 89, 178, 121, 242, 249, 239, 195, 155, 43, 86, 172, 69, 138, 9, 
+18, 36, 72, 144, 61, 122, 244, 245, 247, 243, 251, 235, 203, 139, 11, 
+22, 44, 88, 176, 125, 250, 233, 207, 131, 27, 54, 108, 216, 173, 71, 
+142};
+
+static uchar gplog[256] = 
+{0, 0, 1, 25, 2, 50, 26, 198, 3, 223, 51, 238, 27, 104, 199, 75, 4, 100, 
+224, 14, 52, 141, 239, 129, 28, 193, 105, 248, 200, 8, 76, 113, 5, 
+138, 101, 47, 225, 36, 15, 33, 53, 147, 142, 218, 240, 18, 130, 69, 
+29, 181, 194, 125, 106, 39, 249, 185, 201, 154, 9, 120, 77, 228, 114, 
+166, 6, 191, 139, 98, 102, 221, 48, 253, 226, 152, 37, 179, 16, 145, 
+34, 136, 54, 208, 148, 206, 143, 150, 219, 189, 241, 210, 19, 92, 
+131, 56, 70, 64, 30, 66, 182, 163, 195, 72, 126, 110, 107, 58, 40, 
+84, 250, 133, 186, 61, 202, 94, 155, 159, 10, 21, 121, 43, 78, 212, 
+229, 172, 115, 243, 167, 87, 7, 112, 192, 247, 140, 128, 99, 13, 103, 
+74, 222, 237, 49, 197, 254, 24, 227, 165, 153, 119, 38, 184, 180, 
+124, 17, 68, 146, 217, 35, 32, 137, 46, 55, 63, 209, 91, 149, 188, 
+207, 205, 144, 135, 151, 178, 220, 252, 190, 97, 242, 86, 211, 171, 
+20, 42, 93, 158, 132, 60, 57, 83, 71, 109, 65, 162, 31, 45, 67, 216, 
+183, 123, 164, 118, 196, 23, 73, 236, 127, 12, 111, 246, 108, 161, 
+59, 82, 41, 157, 85, 170, 251, 96, 134, 177, 187, 204, 62, 90, 203, 
+89, 95, 176, 156, 169, 160, 81, 11, 245, 22, 235, 122, 117, 44, 215, 
+79, 174, 213, 233, 230, 231, 173, 232, 116, 214, 244, 234, 168, 80, 
+88, 175};
+
+static uchar *polys[69];
+static int cap[160];
+static uchar *blocks[160];
+
+static uchar side[40] = {
+       21, 25, 29, 33, 37, 41, 45, 49,
+       53, 57, 61, 65, 69, 73, 77, 81,
+       85, 89, 93, 97, 101, 105, 109, 113,
+       117, 121, 125, 129, 133, 137, 141, 145,
+       149, 153, 157, 161, 165, 169, 173, 177
+};
+
+static uchar align[40][8] = {
+       {0},
+       {6, 18, 0},
+       {6, 22, 0},
+       {6, 26, 0},
+       {6, 30, 0},
+       {6, 34, 0},
+       {6, 22, 38, 0},
+       {6, 24, 42, 0},
+       {6, 26, 46, 0},
+       {6, 28, 50, 0},
+       {6, 30, 54, 0},
+       {6, 32, 58, 0},
+       {6, 34, 62, 0},
+       {6, 26, 46, 66, 0},
+       {6, 26, 48, 70, 0},
+       {6, 26, 50, 74, 0},
+       {6, 30, 54, 78, 0},
+       {6, 30, 56, 82, 0},
+       {6, 30, 58, 86, 0},
+       {6, 34, 62, 90, 0},
+       {6, 28, 50, 72, 94, 0},
+       {6, 26, 50, 74, 98, 0},
+       {6, 30, 54, 78, 102, 0},
+       {6, 28, 54, 80, 106, 0},
+       {6, 32, 58, 84, 110, 0},
+       {6, 30, 58, 86, 114, 0},
+       {6, 34, 62, 90, 118, 0},
+       {6, 26, 50, 74, 98, 122, 0},
+       {6, 30, 54, 78, 102, 126, 0},
+       {6, 26, 52, 78, 104, 130, 0},
+       {6, 30, 56, 82, 108, 134, 0},
+       {6, 34, 60, 86, 112, 138, 0},
+       {6, 30, 58, 86, 114, 142, 0},
+       {6, 34, 62, 90, 118, 146, 0},
+       {6, 30, 54, 78, 102, 126, 150, 0},
+       {6, 24, 50, 76, 102, 128, 154, 0},
+       {6, 28, 54, 80, 106, 132, 158, 0},
+       {6, 32, 58, 84, 110, 136, 162, 0},
+       {6, 26, 54, 82, 110, 138, 166, 0},
+       {6, 30, 58, 86, 114, 142, 170, 0},
+};
+
+static void
+qrinit(void)
+{
+       uchar *p;
+       int in;
+       extern uchar qrrawpolys[];
+       extern uchar qrsizes[];
+       
+       in = 0;
+       for(p = qrrawpolys; in || *p != 0; p++){
+               if(!in)
+                       polys[*p] = p + 1;
+               in = *p;
+       }
+       memset(cap, 0, sizeof(cap));
+       for(p = qrsizes; p[0] != 0 || p[1] != 0; p += 4){
+               cap[p[0]] += p[1] * p[3];
+               if(blocks[p[0]] == nil)
+                       blocks[p[0]] = p;
+       }
+}
+
+static void
+bits(uchar *b, int *elen, int *bctr, ulong data, int n)
+{
+       int l, bc, s;
+       
+       l = *elen;
+       bc = *bctr;
+       
+       s = b[l] << n | data;
+       bc += n;
+       while(bc > 8)
+               b[l++] = s >> (bc -= 8);
+       b[l] = s;
+       if(bc == 8){
+               l++;
+               bc = 0;
+       }
+       
+       *elen = l;
+       *bctr = bc;
+}
+
+static int
+getlbits(int version, int mode)
+{
+       if(version <= 9)
+               switch(mode){
+               case NUM: return 10;
+               case ANUM: return 9;
+               case BYTE: return 8;
+               }
+       else if(version <= 26)
+               switch(mode){
+               case NUM: return 12;
+               case ANUM: return 11;
+               case BYTE: return 16;
+               }
+       else
+               switch(mode){
+               case NUM: return 14;
+               case ANUM: return 13;
+               case BYTE: return 16;
+               }
+       werrstr("unknown mode %#x", mode);
+       return -1;
+}
+
+static char
+atab[] = {
+       ['0'] 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,
+       ['A'] 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,
+       ['a'] 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 
+};
+
+static uchar
+formsel(int version, int level, int mode, uchar *data, int dlen, int *nchar)
+{
+       int bits, lbits, i, n;
+       
+       bits = 8;
+       switch(mode){
+       case NUM:
+               n = 0;
+               for(i = 0; i < dlen; i++)
+                       if(isdigit(data[i]))
+                               n++;
+               *nchar = n;
+               bits += 10 * (n / 3);
+               switch(n%3){
+               case 1: bits += 4; break;
+               case 2: bits += 7; break;
+               }
+               break;
+       case ANUM:
+               n = 0;
+               for(i = 0; i < dlen; i++)
+                       if(atab[data[i]] != 0)
+                               n++;
+               *nchar = n;
+               bits += 11 * (n >> 1) + 6 * (n & 1);
+               break;
+       case BYTE:
+               bits += 8 * dlen;
+               *nchar = dlen;
+               break;
+       }
+       if(version == 0){
+               SET(lbits);
+               for(i = 0; i < 40; i++){
+                       if(i == 0 || i == 10 || i == 27){
+                               lbits = getlbits(i + 1, mode);
+                               if(lbits < 0)
+                                       return -1;
+                       }
+                       if(cap[i << 2 | level] >= (bits + lbits + 7 >> 3))
+                               return i << 2 | level;
+               }
+       }else{
+               lbits = getlbits(version, mode);
+               if(lbits < 0)
+                       return -1;
+               if(cap[version - 1 << 2 | level] >= (bits + lbits + 7 >> 3))
+                       return version - 1 << 2 | level;
+       }
+       werrstr("too many bits");
+       return -1;
+}
+
+static uchar *
+encode(uchar verlev, uchar *data, int dlen, int mode, int nchar)
+{
+       uchar *b;
+       int bc, n, c, nd;
+       int acc;
+       
+       n = 0;
+       c = cap[verlev];
+       b = malloc(c);
+       if(b == nil)
+               return nil;
+       bits(b, &n, &bc, mode, 4);
+       bits(b, &n, &bc, nchar, getlbits((verlev >> 2) + 1, mode));
+       switch(mode){
+       case NUM:
+               nd = 0;
+               acc = 0;
+               for(; dlen-- != 0; data++)
+                       if(isdigit(*data)){
+                               acc = acc * 10 + *data - '0';
+                               if(++nd == 3){
+                                       nd = 0;
+                                       bits(b, &n, &bc, acc, 10);
+                                       acc = 0;
+                               }
+                       }
+               if(nd == 1)
+                       bits(b, &n, &bc, acc, 4);
+               else if(nd == 2)
+                       bits(b, &n, &bc, acc, 7);
+               break;
+       case ANUM:
+               nd = 0;
+               acc = 0;
+               for(; dlen-- != 0; data++)
+                       if(atab[*data] != 0){
+                               acc = acc * 45 + atab[*data] - 1;
+                               if(++nd == 2){
+                                       nd = 0;
+                                       bits(b, &n, &bc, acc, 11);
+                                       acc = 0;
+                               }
+                       }
+               if(nd == 1)
+                       bits(b, &n, &bc, acc, 6);
+               break;
+       case BYTE:
+               while(dlen-- != 0)
+                       bits(b, &n, &bc, *data++, 8);
+               break;
+       default:
+               werrstr("unknown mode %#x", mode);
+               return nil;
+       }
+       bits(b, &n, &bc, 0, 4);
+       if(bc != 0)
+               b[n++] <<= 8 - bc;
+       while(n < c){
+               b[n++] = 0xec;
+               if(n == c)
+                       break;
+               b[n++] = 0x11;
+       }
+       return b;
+}
+
+/* all stored highest coefficient first
+   poly stored as log and highest coefficient implicit 1 */
+static int
+ecc(uchar *data, int dl, uchar *poly, uchar *err, int el, int str)
+{
+       uchar *q;
+       uchar *dp;
+       uchar l;
+       int i, x;
+       
+       q = malloc(el);
+       if(q == nil)
+               return -1;
+       for(i = 0; i < el; i++){
+               q[i] = gpexp[poly[i]];
+               err[i * str] = 0;
+       }
+       for(dp = data + dl - 1; dp != data - 1; dp--){
+               if(*dp != 0){
+                       l = gplog[*dp];
+                       for(i = 0; i < el; i++){
+                               if(q[i] != 0){
+                                       x = gplog[q[i]] + l;
+                                       if(x >= 255)
+                                               x -= 255;
+                                       err[i * str] ^= gpexp[x];
+                               }
+                       }
+               }
+               l = q[0];
+               if(l != 0){
+                       l = gplog[l];
+                       for(i = 0; i < el; i++){
+                               x = l + poly[i];
+                               if(x >= 255)
+                                       x -= 255;
+                               q[i] = (i == el-1 ? 0 : q[i+1]) ^ gpexp[x];
+                       }
+               }else{
+                       for(i = 0; i < el - 1; i++)
+                               q[i] = q[i+1];
+                       q[el - 1] = 0;
+               }
+       }
+       free(q);
+       return 0;
+}
+
+static uchar *
+codewords(uchar verlev, uchar *data, int *ncp)
+{
+       uchar *bp, *p, *dp, *r, *rp;
+       uchar *ep;
+       int ne, i, j, el, nc, nb, ncol;
+       
+       bp = blocks[verlev];
+       el = bp[2] - bp[3]; /* guaranteed constant */
+       ncol = 0;
+       for(p = bp, ne = 0, nc = 0, nb = 0; p[0] == verlev; p += 4){
+               nb += p[1];
+               ne += p[1] * el;
+               nc += p[1] * p[2];
+               ncol = p[3];
+       }
+       r = malloc(nc);
+       if(r == nil)
+               return nil;
+       for(p = bp, dp = data, ep = r + nc - ne; p[0] == verlev; p += 4)
+               for(j = 0; j < p[1]; j++){
+                       if(ecc(dp, p[3], polys[el], ep, el, nb) < 0){
+                               free(r);
+                               return nil;
+                       }
+                       dp += p[3];
+                       ep++;
+               }
+       for(i = 0, rp = r; i < ncol; i++){
+               dp = data + i;
+               for(p = bp; p[0] == verlev; p += 4){
+                       for(j = 0; j < p[1]; j++){
+                               if(p[3] > i)
+                                       *rp++ = *dp;
+                               dp += p[3];
+                       }
+               }
+       }
+       *ncp = nc;
+       return r;
+}
+
+#define max(i,j) ((i)>(j)?(i):(j))
+
+uchar *
+basepat(int verlev)
+{
+       int s, i, j, d;
+       uchar *pat, *p, *ap, *aq;
+       
+       s = side[verlev >> 2];
+       pat = malloc((s + 8) * (s + 8));
+       if(pat == nil)
+               return nil;
+       p = pat;
+       for(j = -4; j < s + 4; j++)
+               for(i = -4; i < s + 4; i++){
+                       if(i < 0 || j < 0 || i >= s || j >= s)
+                               *p++ = 0xff;
+                       else if(i < 8 && j < 8 || i >= s - 8 && j < 8 || i < 8 && j >= s - 8){
+                               d = max(abs(i >= 8 ? i - s + 4 : i-3), abs(j >= 8 ? j - s + 4 : j-3));
+                               *p++ = d == 2 || d == 4 ? 0xff : 0;
+                       }else if(i == 6 || j == 6)
+                               *p++ = (i + j & 1) != 0 ? 0xff : 0;
+                       else if(i <= 8 && j <= 8 || i <= 8 && j >= s - 8 || i >= s - 8 && j <= 8)
+                               *p++ = 0xd0;
+                       else if(verlev >= 24 && (j >= s - 11 && i < 6 || i >= s - 11 && j < 6))
+                               *p++ = 0xc0;
+                       else
+                               *p++ = 0x80;
+               }
+       for(ap = align[verlev >> 2]; *ap != 0; ap++)
+               for(aq = align[verlev >> 2]; *aq != 0; aq++){
+                       if(*ap == 6 && *aq == 6 || *ap >= s - 7 && *aq == 6 || *ap == 6 && *aq >= s - 7)
+                               continue;
+                       for(i = *ap - 2; i <= *ap + 2; i++)
+                               for(j = *aq - 2; j <= *aq + 2; j++)
+                                       pat[i * (s + 8) + j + 4 + (s + 8) * 4] = max(abs(i - *ap), abs(j - *aq)) == 1 ? 0xff : 0;
+               }
+       return pat;
+}
+
+static void
+fill(uchar verlev, uchar *p, uchar *d, int nc)
+{
+       int s, x, y, i, up;
+       uchar *dp, db;
+       
+       s = side[verlev >> 2] + 8;
+       x = s - 5;
+       y = s - 5;
+       up = 1;
+       SET(db);
+       for(dp = d, i = 0;;){
+               if(i == 0){
+                       if(dp - d == nc)
+                               db = 0;
+                       else
+                               db = *dp++;
+               }
+               i = i + 1 & 7;
+               while(p[y * s + x] != 0x80){
+                       if(up){
+                               if((x & 1 ^ (x < 10)) == 0)
+                                       x--;
+                               else if(y >= 4){
+                                       x++;
+                                       y--;
+                               }else{
+                                       x--;
+                                       up = 0;
+                                       if(x == 10)
+                                               x--;
+                               }
+                       }else{
+                               if((x & 1 ^ (x < 10)) == 0)
+                                       x--;
+                               else if(y < s - 4){
+                                       x++;
+                                       y++;
+                               }else{
+                                       x--;
+                                       up = 1;
+                                       if(x == 10)
+                                               x--;
+                               }
+                       }
+                       if(x < 0)
+                               return;
+               }
+               p[y * s + x] = ~((char)db >> 7) ^ 0x0f;
+               db <<= 1;
+       }
+}
+
+static void
+mask(uchar verlev, uchar *p, int m)
+{
+       uchar *q;
+       int s, i, j;
+       
+       s = side[verlev >> 2];
+       for(i = -4, q = p; i < s + 4; i++)
+               for(j = -4; j < s + 4; j++, q++)
+                       if(((*q ^ *q >> 4) & 15) == 15)
+                               switch(m){
+                               case 0: *q ^= (char)((i+j+1) << 7) >> 7; break;
+                               case 1: *q ^= (char)((i+1) << 7) >> 7; break;
+                               case 2: *q ^= (char)(j % 3 - 1) >> 7; break;
+                               case 3: *q ^= (char)((i+j) % 3 - 1) >> 7; break;
+                               case 4: *q ^= (char)(((i/2)+(j/3)+1) << 7) >> 7; break;
+                               case 5: *q ^= (char)((i*j%2)+(i*j%3) - 1) >> 7; break;
+                               case 6: *q ^= (char)((i*j+(i*j%3)+1) << 7) >> 7; break;
+                               case 7: *q ^= (char)((i+j+(i*j%3)+1) << 7) >> 7; break;
+                               }
+}
+
+static int
+evaluate(uchar verlev, uchar *p)
+{
+       int r, i, j, l, s, sc, cnt;
+       uchar q, sh;
+       uchar *pp;
+       enum {
+               N1 = 3,
+               N2 = 3,
+               N3 = 40,
+               N4 = 10,
+       };
+       
+       s = side[verlev >> 2];
+       sc = 0;
+       cnt = 0;
+
+       for(i = 0; i < s; i++){
+               r = 0;
+               l = -1;
+               sh = 0;
+               for(j = 0; j < s; j++){
+                       q = p[i * (s + 8) + j + 4 + (s + 8) * 4] >= 0x80;
+                       cnt += q;
+                       sh = sh << 1 | l;
+                       if(q != l){
+                               if(r >= 5)
+                                       sc += N1 + r;
+                               l = q;
+                               r = 0;
+                       }
+                       if((sh & 0x7f) == 0x5d)
+                               sc += N3;
+               }
+               if(r >= 5)
+                       sc += N1 * r;
+       }
+       sc = abs(cnt - s * s / 2) / (s * s / 20) * N4;
+       
+       for(j = 0; j < s; j++){
+               r = 0;
+               l = -1;
+               sh = 0;
+               for(i = 0; i < s; i++){
+                       q = p[i * (s + 8) + j + 4 + (s + 8) * 4];
+                       sh = sh << 1 | l;
+                       if(q != l){
+                               if(r >= 5)
+                                       sc += N1 + r;
+                               l = q;
+                               r = 0;
+                       }
+                       if((sh & 0x7f) == 0x5d)
+                               sc += N3;
+               }
+               if(r >= 5)
+                       sc += N1 * r;
+       }
+
+       /* approximation */
+       for(j = 1; j < s; j++)
+               for(i = 1; i < s; i++){
+                       pp = p + i * (s + 8) + j + 4 + (s + 8) * 4;
+                       if(pp[0] == pp[-1] && pp[0] == pp[-(s+8)] && pp[0] == pp[-1-(s+8)])
+                               sc += N2;
+               }
+       return sc;
+}
+
+static void
+format(uchar verlev, uchar *p, int m)
+{
+       static ushort format[] = 
+               {21522, 20773, 24188, 23371, 17913, 16590, 20375, 19104, 30660, 
+               29427, 32170, 30877, 26159, 25368, 27713, 26998, 5769, 5054, 7399, 
+               6608, 1890, 597, 3340, 2107, 13663, 12392, 16177, 14854, 9396, 8579, 
+               11994, 11245};
+       int s, i;
+       ushort f;
+       uchar c;
+       
+       s = side[verlev >> 2];
+       f = ~format[verlev << 3 & 0x18 ^ 8 | m];
+       for(i = 0; i < 15; i++){
+               c = (char)(f << 7) >> 7;
+               if(i < 8){
+                       p[(i + (i >= 6)) * (s + 8) + 8 + 4 + (s + 8) * 4] = c;
+                       p[12 * (s + 8) + (s + 3 - i)] = c;
+               }else{
+                       p[(s - 15 + i) * (s + 8) + 8 + 4 + (s + 8) * 4] = c;
+                       p[12 * (s + 8) + (18 - i + (i == 8))] = c;
+               }
+               f >>= 1;
+       }
+       p[(s - 8) * (s + 8) + 8 + 4 + (s + 8) * 4] = 0;
+}
+
+static void
+version(uchar verlev, uchar *p)
+{
+       int s, i, f, q;
+       uchar c;
+       static int version[] = {
+               0x07C94, 0x085BC, 0x09A99, 0x0A4D3, 0x0BBF6, 0x0C762, 0x0D847,
+               0x0E60D, 0x0F928, 0x10B78, 0x1145D, 0x12A17, 0x13532, 0x149A6,
+               0x15683, 0x168C9, 0x177EC, 0x18EC4, 0x191E1, 0x1AFAB, 0x1B08E,
+               0x1CC1A, 0x1D33F, 0x1ED75, 0x1F250, 0x209D5, 0x216F0, 0x228BA,
+               0x2379F, 0x24B0B, 0x2542E, 0x26A64, 0x27541, 0x28C69
+       };
+       
+       if(verlev < 24)
+               return;
+       s = side[verlev >> 2] + 8;
+       f = ~version[(verlev >> 2) - 6];
+       for(i = 0; i < 18; i++){
+               c = (char)(f << 7) >> 7;
+               f >>= 1;
+               q = (s - 15 + i % 3) * s + (4 + i / 3);
+               p[q] = c;
+               p[(q * s) % (s * s - 1)] = c;
+       }
+}
+
+static void
+finalize(uchar verlev, uchar *p)
+{
+       int s;
+       
+       s = side[verlev >> 2] + 8;
+       s *= s;
+       while(s--)
+               *((char *)p++) >>= 4;
+}
+
+uchar *
+qrcode(int ver, int level, int mode, uchar *data, int dlen, int *sidep)
+{
+       int s, nc, m, mv, mbest, mval, nchar;
+       uchar *d, *w, *b;
+       uchar verlev;
+
+       qrinit();
+       verlev = formsel(ver, level, mode, data, dlen, &nchar);
+       if(verlev == 0xff)
+               return nil;
+       s = side[verlev >> 2];
+       d = encode(verlev, data, dlen, mode, nchar);
+       if(d == nil)
+               return nil;
+       w = codewords(verlev, d, &nc);
+       free(d);
+       if(w == nil)
+               return nil;
+       b = basepat(verlev);
+       if(b == nil){
+               free(w);
+               return nil;
+       }
+       fill(verlev, b, w, nc);
+       free(w);
+       SET(mbest); SET(mval);
+       for(m = 0; m < 8; m++){
+               mask(verlev, b, m);
+               mv = evaluate(verlev, b);
+               if(m == 0 || mv < mval){
+                       mbest = m;
+                       mval = mv;
+               }
+               if(m != 7 || mbest != m)
+                       mask(verlev, b, m);
+       }
+       if(mbest != 7)
+               mask(verlev, b, mbest);
+       format(verlev, b, mbest);
+       version(verlev, b);
+       finalize(verlev, b);
+       *sidep = s + 8;
+       return b;
+}
+
+void
+usage(void)
+{
+       fprint(2, "usage: %s [-LMQHn] [-v ver]\n", argv0);
+       exits("usage");
+}
+
+void
+main(int argc, char **argv)
+{
+       uchar buf[8192];
+       int rc;
+       int ver, lev, mode, s;
+       char *c;
+       uchar *b;
+       
+       ver = 0;
+       lev = 1;
+       mode = BYTE;
+       ARGBEGIN {
+       case 'L': lev = 0; break;
+       case 'M': lev = 1; break;
+       case 'Q': lev = 2; break;
+       case 'H': lev = 3; break;
+       case 'v':
+               ver = strtol(EARGF(usage()), &c, 0);
+               if(*c != 0 || ver < 0 || ver > 40)
+                       usage();
+               break;
+       case 'n': mode = NUM; break;
+       case 'a': mode = ANUM; break;
+       default: usage();
+       } ARGEND;
+       
+       if(argc != 0)
+               usage();
+       
+       rc = readn(0, buf, sizeof(buf));
+       if(rc < 0)
+               sysfatal("read: %r");
+       b = qrcode(ver, lev, mode, buf, rc, &s);
+       if(b == nil)
+               sysfatal("%r");
+       print("%11s %11d %11d %11d %11d ", "k8", 0, 0, s, s);
+       write(1, b, s * s);
+}
+
+uchar qrrawpolys[] = {
+       7, 87, 229, 146, 149, 238, 102, 21,
+       0, 10, 251, 67, 46, 61, 118, 70,
+       64, 94, 32, 45, 0, 13, 74, 152,
+       176, 100, 86, 100, 106, 104, 130, 218,
+       206, 140, 78, 0, 15, 8, 183, 61,
+       91, 202, 37, 51, 58, 58, 237, 140,
+       124, 5, 99, 105, 0, 16, 120, 104,
+       107, 109, 102, 161, 76, 3, 91, 191,
+       147, 169, 182, 194, 225, 120, 0, 17,
+       43, 139, 206, 78, 43, 239, 123, 206,
+       214, 147, 24, 99, 150, 39, 243, 163,
+       136, 0, 18, 215, 234, 158, 94, 184,
+       97, 118, 170, 79, 187, 152, 148, 252,
+       179, 5, 98, 96, 153, 0, 20, 17,
+       60, 79, 50, 61, 163, 26, 187, 202,
+       180, 221, 225, 83, 239, 156, 164, 212,
+       212, 188, 190, 0, 22, 210, 171, 247,
+       242, 93, 230, 14, 109, 221, 53, 200,
+       74, 8, 172, 98, 80, 219, 134, 160,
+       105, 165, 231, 0, 24, 229, 121, 135,
+       48, 211, 117, 251, 126, 159, 180, 169,
+       152, 192, 226, 228, 218, 111, 117, 232,
+       87, 96, 227, 21, 0, 26, 173, 125,
+       158, 2, 103, 182, 118, 17, 145, 201,
+       111, 28, 165, 53, 161, 21, 245, 142,
+       13, 102, 48, 227, 153, 145, 218, 70,
+       0, 28, 168, 223, 200, 104, 224, 234,
+       108, 180, 110, 190, 195, 147, 205, 27,
+       232, 201, 21, 43, 245, 87, 42, 195,
+       212, 119, 242, 37, 9, 123, 0, 30,
+       41, 173, 145, 152, 216, 31, 179, 182,
+       50, 48, 110, 86, 239, 96, 222, 125,
+       42, 173, 226, 193, 224, 130, 156, 37,
+       251, 216, 238, 40, 192, 180, 0, 32,
+       10, 6, 106, 190, 249, 167, 4, 67,
+       209, 138, 138, 32, 242, 123, 89, 27,
+       120, 185, 80, 156, 38, 69, 171, 60,
+       28, 222, 80, 52, 254, 185, 220, 241,
+       0, 34, 111, 77, 146, 94, 26, 21,
+       108, 19, 105, 94, 113, 193, 86, 140,
+       163, 125, 58, 158, 229, 239, 218, 103,
+       56, 70, 114, 61, 183, 129, 167, 13,
+       98, 62, 129, 51, 0, 36, 200, 183,
+       98, 16, 172, 31, 246, 234, 60, 152,
+       115, 167, 152, 113, 248, 238, 107, 18,
+       63, 218, 37, 87, 210, 105, 177, 120,
+       74, 121, 196, 117, 251, 113, 233, 30,
+       120, 0, 40, 59, 116, 79, 161, 252,
+       98, 128, 205, 128, 161, 247, 57, 163,
+       56, 235, 106, 53, 26, 187, 174, 226,
+       104, 170, 7, 175, 35, 181, 114, 88,
+       41, 47, 163, 125, 134, 72, 20, 232,
+       53, 35, 15, 0, 42, 250, 103, 221,
+       230, 25, 18, 137, 231, 3, 58, 242,
+       221, 191, 110, 84, 230, 8, 188, 106,
+       96, 147, 15, 131, 139, 34, 101, 223,
+       39, 101, 213, 199, 237, 254, 201, 123,
+       171, 162, 194, 117, 50, 96, 0, 44,
+       190, 7, 61, 121, 71, 246, 69, 55,
+       168, 188, 89, 243, 191, 25, 72, 123,
+       9, 145, 14, 247, 238, 44, 78, 143,
+       62, 224, 126, 118, 114, 68, 163, 52,
+       194, 217, 147, 204, 169, 37, 130, 113,
+       102, 73, 181, 0, 46, 112, 94, 88,
+       112, 253, 224, 202, 115, 187, 99, 89,
+       5, 54, 113, 129, 44, 58, 16, 135,
+       216, 169, 211, 36, 4, 96, 60, 241,
+       73, 104, 234, 8, 249, 245, 119, 174,
+       52, 25, 157, 224, 43, 202, 223, 19,
+       82, 15, 0, 48, 228, 25, 196, 130,
+       211, 146, 60, 24, 251, 90, 39, 102,
+       240, 61, 178, 63, 46, 123, 115, 18,
+       221, 111, 135, 160, 182, 205, 107, 206,
+       95, 150, 120, 184, 91, 21, 247, 156,
+       140, 238, 191, 11, 94, 227, 84, 50,
+       163, 39, 34, 108, 0, 50, 232, 125,
+       157, 161, 164, 9, 118, 46, 209, 99,
+       203, 193, 35, 3, 209, 111, 195, 242,
+       203, 225, 46, 13, 32, 160, 126, 209,
+       130, 160, 242, 215, 242, 75, 77, 42,
+       189, 32, 113, 65, 124, 69, 228, 114,
+       235, 175, 124, 170, 215, 232, 133, 205,
+       0, 52, 116, 50, 86, 186, 50, 220,
+       251, 89, 192, 46, 86, 127, 124, 19,
+       184, 233, 151, 215, 22, 14, 59, 145,
+       37, 242, 203, 134, 254, 89, 190, 94,
+       59, 65, 124, 113, 100, 233, 235, 121,
+       22, 76, 86, 97, 39, 242, 200, 220,
+       101, 33, 239, 254, 116, 51, 0, 54,
+       183, 26, 201, 87, 210, 221, 113, 21,
+       46, 65, 45, 50, 238, 184, 249, 225,
+       102, 58, 209, 218, 109, 165, 26, 95,
+       184, 192, 52, 245, 35, 254, 238, 175,
+       172, 79, 123, 25, 122, 43, 120, 108,
+       215, 80, 128, 201, 235, 8, 153, 59,
+       101, 31, 198, 76, 31, 156, 0, 56,
+       106, 120, 107, 157, 164, 216, 112, 116,
+       2, 91, 248, 163, 36, 201, 202, 229,
+       6, 144, 254, 155, 135, 208, 170, 209,
+       12, 139, 127, 142, 182, 249, 177, 174,
+       190, 28, 10, 85, 239, 184, 101, 124,
+       152, 206, 96, 23, 163, 61, 27, 196,
+       247, 151, 154, 202, 207, 20, 61, 10,
+       0, 58, 82, 116, 26, 247, 66, 27,
+       62, 107, 252, 182, 200, 185, 235, 55,
+       251, 242, 210, 144, 154, 237, 176, 141,
+       192, 248, 152, 249, 206, 85, 253, 142,
+       65, 165, 125, 23, 24, 30, 122, 240,
+       214, 6, 129, 218, 29, 145, 127, 134,
+       206, 245, 117, 29, 41, 63, 159, 142,
+       233, 125, 148, 123, 0, 60, 107, 140,
+       26, 12, 9, 141, 243, 197, 226, 197,
+       219, 45, 211, 101, 219, 120, 28, 181,
+       127, 6, 100, 247, 2, 205, 198, 57,
+       115, 219, 101, 109, 160, 82, 37, 38,
+       238, 49, 160, 209, 121, 86, 11, 124,
+       30, 181, 84, 25, 194, 87, 65, 102,
+       190, 220, 70, 27, 209, 16, 89, 7,
+       33, 240, 0, 62, 65, 202, 113, 98,
+       71, 223, 248, 118, 214, 94, 122, 37,
+       23, 2, 228, 58, 121, 7, 105, 135,
+       78, 243, 118, 70, 76, 223, 89, 72,
+       50, 70, 111, 194, 17, 212, 126, 181,
+       35, 221, 117, 235, 11, 229, 149, 147,
+       123, 213, 40, 115, 6, 200, 100, 26,
+       246, 182, 218, 127, 215, 36, 186, 110,
+       106, 0, 64, 45, 51, 175, 9, 7,
+       158, 159, 49, 68, 119, 92, 123, 177,
+       204, 187, 254, 200, 78, 141, 149, 119,
+       26, 127, 53, 160, 93, 199, 212, 29,
+       24, 145, 156, 208, 150, 218, 209, 4,
+       216, 91, 47, 184, 146, 47, 140, 195,
+       195, 125, 242, 238, 63, 99, 108, 140,
+       230, 242, 31, 204, 11, 178, 243, 217,
+       156, 213, 231, 0, 66, 5, 118, 222,
+       180, 136, 136, 162, 51, 46, 117, 13,
+       215, 81, 17, 139, 247, 197, 171, 95,
+       173, 65, 137, 178, 68, 111, 95, 101,
+       41, 72, 214, 169, 197, 95, 7, 44,
+       154, 77, 111, 236, 40, 121, 143, 63,
+       87, 80, 253, 240, 126, 217, 77, 34,
+       232, 106, 50, 168, 82, 76, 146, 67,
+       106, 171, 25, 132, 93, 45, 105, 0,
+       68, 247, 159, 223, 33, 224, 93, 77,
+       70, 90, 160, 32, 254, 43, 150, 84,
+       101, 190, 205, 133, 52, 60, 202, 165,
+       220, 203, 151, 93, 84, 15, 84, 253,
+       173, 160, 89, 227, 52, 199, 97, 95,
+       231, 52, 177, 41, 125, 137, 241, 166,
+       225, 118, 2, 54, 32, 82, 215, 175,
+       198, 43, 238, 235, 27, 101, 184, 127,
+       3, 5, 8, 163, 238, 0, 0
+};
+
+/* (version - 1) * 4 + level, number of blocks, number of bytes in block, number of data bytes */
+uchar qrsizes[] = {
+       0, 1, 26, 19,
+       1, 1, 26, 16,
+       2, 1, 26, 13,
+       3, 1, 26, 9,
+       4, 1, 44, 34,
+       5, 1, 44, 28,
+       6, 1, 44, 22,
+       7, 1, 44, 16,
+       8, 1, 70, 55,
+       9, 1, 70, 44,
+       10, 2, 35, 17,
+       11, 2, 35, 13,
+       12, 1, 100, 80,
+       13, 2, 50, 32,
+       14, 2, 50, 24,
+       15, 4, 25, 9,
+       16, 1, 134, 108,
+       17, 2, 67, 43,
+       18, 2, 33, 15,
+       18, 2, 34, 16,
+       19, 2, 33, 11,
+       19, 2, 34, 12,
+       20, 2, 86, 68,
+       21, 4, 43, 27,
+       22, 4, 43, 19,
+       23, 4, 43, 15,
+       24, 2, 98, 78,
+       25, 4, 49, 31,
+       26, 2, 32, 14,
+       26, 4, 33, 15,
+       27, 4, 39, 13,
+       27, 1, 40, 14,
+       28, 2, 121, 97,
+       29, 2, 60, 38,
+       29, 2, 61, 39,
+       30, 4, 40, 18,
+       30, 2, 41, 19,
+       31, 4, 40, 14,
+       31, 2, 41, 15,
+       32, 2, 146, 116,
+       33, 3, 58, 36,
+       33, 2, 59, 37,
+       34, 4, 36, 16,
+       34, 4, 37, 17,
+       35, 4, 36, 12,
+       35, 4, 37, 13,
+       36, 2, 86, 68,
+       36, 2, 87, 69,
+       37, 4, 69, 43,
+       37, 1, 70, 44,
+       38, 6, 43, 19,
+       38, 2, 44, 20,
+       39, 6, 43, 15,
+       39, 2, 44, 16,
+       40, 4, 101, 81,
+       41, 1, 80, 50,
+       41, 4, 81, 51,
+       42, 4, 50, 22,
+       42, 4, 51, 23,
+       43, 3, 36, 12,
+       43, 8, 37, 13,
+       44, 2, 116, 92,
+       44, 2, 117, 93,
+       45, 6, 58, 36,
+       45, 2, 59, 37,
+       46, 4, 46, 20,
+       46, 6, 47, 21,
+       47, 7, 42, 14,
+       47, 4, 43, 15,
+       48, 4, 133, 107,
+       49, 8, 59, 37,
+       49, 1, 60, 38,
+       50, 8, 44, 20,
+       50, 4, 45, 21,
+       51, 12, 33, 11,
+       51, 4, 34, 12,
+       52, 3, 145, 115,
+       52, 1, 146, 116,
+       53, 4, 64, 40,
+       53, 5, 65, 41,
+       54, 11, 36, 16,
+       54, 5, 37, 17,
+       55, 11, 36, 12,
+       55, 5, 37, 13,
+       56, 5, 109, 87,
+       56, 1, 110, 88,
+       57, 5, 65, 41,
+       57, 5, 66, 42,
+       58, 5, 54, 24,
+       58, 7, 55, 25,
+       59, 11, 36, 12,
+       59, 7, 37, 13,
+       60, 5, 122, 98,
+       60, 1, 123, 99,
+       61, 7, 73, 45,
+       61, 3, 74, 46,
+       62, 15, 43, 19,
+       62, 2, 44, 20,
+       63, 3, 45, 15,
+       63, 13, 46, 16,
+       64, 1, 135, 107,
+       64, 5, 136, 108,
+       65, 10, 74, 46,
+       65, 1, 75, 47,
+       66, 1, 50, 22,
+       66, 15, 51, 23,
+       67, 2, 42, 14,
+       67, 17, 43, 15,
+       68, 5, 150, 120,
+       68, 1, 151, 121,
+       69, 9, 69, 43,
+       69, 4, 70, 44,
+       70, 17, 50, 22,
+       70, 1, 51, 23,
+       71, 2, 42, 14,
+       71, 19, 43, 15,
+       72, 3, 141, 113,
+       72, 4, 142, 114,
+       73, 3, 70, 44,
+       73, 11, 71, 45,
+       74, 17, 47, 21,
+       74, 4, 48, 22,
+       75, 9, 39, 13,
+       75, 16, 40, 14,
+       76, 3, 135, 107,
+       76, 5, 136, 108,
+       77, 3, 67, 41,
+       77, 13, 68, 42,
+       78, 15, 54, 24,
+       78, 5, 55, 25,
+       79, 15, 43, 15,
+       79, 10, 44, 16,
+       80, 4, 144, 116,
+       80, 4, 145, 117,
+       81, 17, 68, 42,
+       82, 17, 50, 22,
+       82, 6, 51, 23,
+       83, 19, 46, 16,
+       83, 6, 47, 17,
+       84, 2, 139, 111,
+       84, 7, 140, 112,
+       85, 17, 74, 46,
+       86, 7, 54, 24,
+       86, 16, 55, 25,
+       87, 34, 37, 13,
+       88, 4, 151, 121,
+       88, 5, 152, 122,
+       89, 4, 75, 47,
+       89, 14, 76, 48,
+       90, 11, 54, 24,
+       90, 14, 55, 25,
+       91, 16, 45, 15,
+       91, 14, 46, 16,
+       92, 6, 147, 117,
+       92, 4, 148, 118,
+       93, 6, 73, 45,
+       93, 14, 74, 46,
+       94, 11, 54, 24,
+       94, 16, 55, 25,
+       95, 30, 46, 16,
+       95, 2, 47, 17,
+       96, 8, 132, 106,
+       96, 4, 133, 107,
+       97, 8, 75, 47,
+       97, 13, 76, 48,
+       98, 7, 54, 24,
+       98, 22, 55, 25,
+       99, 22, 45, 15,
+       99, 13, 46, 16,
+       100, 10, 142, 114,
+       100, 2, 143, 115,
+       101, 19, 74, 46,
+       101, 4, 75, 47,
+       102, 28, 50, 22,
+       102, 6, 51, 23,
+       103, 33, 46, 16,
+       103, 4, 47, 17,
+       104, 8, 152, 122,
+       104, 4, 153, 123,
+       105, 22, 73, 45,
+       105, 3, 74, 46,
+       106, 8, 53, 23,
+       106, 26, 54, 24,
+       107, 12, 45, 15,
+       107, 28, 46, 16,
+       108, 3, 147, 117,
+       108, 10, 148, 118,
+       109, 3, 73, 45,
+       109, 23, 74, 46,
+       110, 4, 54, 24,
+       110, 31, 55, 25,
+       111, 11, 45, 15,
+       111, 31, 46, 16,
+       112, 7, 146, 116,
+       112, 7, 147, 117,
+       113, 21, 73, 45,
+       113, 7, 74, 46,
+       114, 1, 53, 23,
+       114, 37, 54, 24,
+       115, 19, 45, 15,
+       115, 26, 46, 16,
+       116, 5, 145, 115,
+       116, 10, 146, 116,
+       117, 19, 75, 47,
+       117, 10, 76, 48,
+       118, 15, 54, 24,
+       118, 25, 55, 25,
+       119, 23, 45, 15,
+       119, 25, 46, 16,
+       120, 13, 145, 115,
+       120, 3, 146, 116,
+       121, 2, 74, 46,
+       121, 29, 75, 47,
+       122, 42, 54, 24,
+       122, 1, 55, 25,
+       123, 23, 45, 15,
+       123, 28, 46, 16,
+       124, 17, 145, 115,
+       125, 10, 74, 46,
+       125, 23, 75, 47,
+       126, 10, 54, 24,
+       126, 35, 55, 25,
+       127, 19, 45, 15,
+       127, 35, 46, 16,
+       128, 17, 145, 115,
+       128, 1, 146, 116,
+       129, 14, 74, 46,
+       129, 21, 75, 47,
+       130, 29, 54, 24,
+       130, 19, 55, 25,
+       131, 11, 45, 15,
+       131, 46, 46, 16,
+       132, 13, 145, 115,
+       132, 6, 146, 116,
+       133, 14, 74, 46,
+       133, 23, 75, 47,
+       134, 44, 54, 24,
+       134, 7, 55, 25,
+       135, 59, 46, 16,
+       135, 1, 47, 17,
+       136, 12, 151, 121,
+       136, 7, 152, 122,
+       137, 12, 75, 47,
+       137, 26, 76, 48,
+       138, 39, 54, 24,
+       138, 14, 55, 25,
+       139, 22, 45, 15,
+       139, 41, 46, 16,
+       140, 6, 151, 121,
+       140, 14, 152, 122,
+       141, 6, 75, 47,
+       141, 34, 76, 48,
+       142, 46, 54, 24,
+       142, 10, 55, 25,
+       143, 2, 45, 15,
+       143, 64, 46, 16,
+       144, 17, 152, 122,
+       144, 4, 153, 123,
+       145, 29, 74, 46,
+       145, 14, 75, 47,
+       146, 49, 54, 24,
+       146, 10, 55, 25,
+       147, 24, 45, 15,
+       147, 46, 46, 16,
+       148, 4, 152, 122,
+       148, 18, 153, 123,
+       149, 13, 74, 46,
+       149, 32, 75, 47,
+       150, 48, 54, 24,
+       150, 14, 55, 25,
+       151, 42, 45, 15,
+       151, 32, 46, 16,
+       152, 20, 147, 117,
+       152, 4, 148, 118,
+       153, 40, 75, 47,
+       153, 7, 76, 48,
+       154, 43, 54, 24,
+       154, 22, 55, 25,
+       155, 10, 45, 15,
+       155, 67, 46, 16,
+       156, 19, 148, 118,
+       156, 6, 149, 119,
+       157, 18, 75, 47,
+       157, 31, 76, 48,
+       158, 34, 54, 24,
+       158, 34, 55, 25,
+       159, 20, 45, 15,
+       159, 61, 46, 16,
+       0, 0, 0, 0,
+};