wok-tiny annotate memtest/stuff/unlzma.S @ rev 170

memtest,plop: tune lzma compression (again)
author Pascal Bellard <pascal.bellard@slitaz.org>
date Fri May 07 15:00:24 2021 +0000 (2021-05-07)
parents a88405d94c1e
children eb617e43dc08
rev   line source
pascal@168 1 // #define RC_NORMALIZE if (Range < kTopValue) { Range <<= 8; Code = (Code << 8) | RC_READ_BYTE; }
pascal@168 2 //
pascal@168 3 // #define IfBit0(p) RC_NORMALIZE; bound = (Range >> kNumBitModelTotalBits) * *(p); if (Code < bound)
pascal@168 4 // #define UpdateBit0(p) Range = bound; *(p) += (kBitModelTotal - *(p)) >> kNumMoveBits;
pascal@168 5 // #define UpdateBit1(p) Range -= bound; Code -= bound; *(p) -= (*(p)) >> kNumMoveBits;
pascal@168 6 //
pascal@168 7 //#define RC_GET_BIT2(p, mi, A0, A1) IfBit0(p) \
pascal@168 8 // { UpdateBit0(p); mi <<= 1; A0; } else \
pascal@168 9 // { UpdateBit1(p); mi = (mi + mi) + 1; A1; }
pascal@168 10 //
pascal@168 11 // #define RC_GET_BIT(p, mi) RC_GET_BIT2(p, mi, ; , ;)
pascal@168 12 //
pascal@168 13 // #define RangeDecoderBitTreeDecode(probs, numLevels, res) \
pascal@168 14 // { int i = numLevels; res = 1; \
pascal@168 15 // do { CProb *p = probs + res; RC_GET_BIT(p, res) } while(--i != 0); \
pascal@168 16 // res -= (1 << numLevels); }
pascal@168 17 /*
pascal@168 18 * Compression with : lzma e src dst -eos -pb2 -lp0 -lc3
pascal@168 19 */
pascal@168 20
pascal@170 21 #define PROP_PB 3
pascal@168 22 #define PROP_LP 0
pascal@168 23 #define PROP_LC 3
pascal@168 24 #define PROPS (PROP_LC+(PROP_LP*9)+(PROP_PB*45))
pascal@168 25
pascal@168 26 // static const Byte *Buffer;
pascal@168 27 // static UInt32 bound, Code, Range;
pascal@168 28
pascal@168 29 /*
pascal@168 30 * Buffer register DS:SI
pascal@168 31 * all var based ws=ss:bp
pascal@168 32 */
pascal@168 33
pascal@168 34 rep0 = -4 // long
pascal@168 35 rep1 = rep0-4 // long
pascal@168 36 rep2 = rep0-8 // long
pascal@168 37 rep3 = rep0-12 // long
pascal@168 38 state = -17 // byte, 0..11
pascal@168 39 posState = state-1 // byte, 0..15
pascal@168 40 posState2 = posState-1 // byte, 0..15
pascal@168 41 scratched = rep0-16 // byte = 1
pascal@168 42 Code = -24 // long
pascal@168 43 outStream = -28 // long
pascal@168 44 nowPos = outStream // long
pascal@168 45 Range = Code-8 // long
pascal@168 46 #define LOCALS 32
pascal@168 47
pascal@168 48 // int LzmaDecode(CLzmaDecoderState *vs,
pascal@168 49 // const unsigned char *inStream,
pascal@168 50 // unsigned char *outStream)
pascal@168 51 // {
pascal@168 52 // CProb *p = vs->Probs;
pascal@168 53 // SizeT nowPos = 0;
pascal@168 54 // #define posStateMask = (1 << (vs->Properties.pb)) - 1;
pascal@168 55 // #define literalPosMask = (1 << (vs->Properties.lp)) - 1;
pascal@168 56 // int lc = vs->Properties.lc, state = 0, len = 0;
pascal@168 57 // UInt32 rep0 = 1, rep1 = 1, rep2 = 1, rep3 = 1;
pascal@168 58 //
pascal@168 59 // {
pascal@168 60 // UInt32 i, numProbs = Literal /*1846*/
pascal@168 61 // + ((UInt32)LZMA_LIT_SIZE /*768*/ << (lc + vs->Properties.lp));
pascal@168 62 // for (i = 0; i < numProbs; i++) p[i] = kBitModelTotal /*2048*/ >> 1;
pascal@168 63
pascal@168 64 #define WS (1846+(768<<(PROP_LC+PROP_LP)))
pascal@168 65 #if (WS+WS+LOCALS) >= 65000
pascal@168 66 /* MAX WS = (1846+(768<<(8+4))) > 3MB! */
pascal@168 67 #error invalid (lc,lp,pb) : out of memory
pascal@168 68 #endif
pascal@168 69
pascal@168 70 ws1 = WS
pascal@168 71 ws2 = ws1*2
pascal@168 72 ws = ws2+LOCALS+15
pascal@168 73
pascal@168 74 #ifndef FLAT32
pascal@168 75 #define AX %ax
pascal@168 76 #define BX %bx
pascal@168 77 #define CX %cx
pascal@168 78 #define DX %dx
pascal@168 79 #define SI %si
pascal@168 80 #define DI %di
pascal@168 81 #define BP %bp
pascal@168 82 #define SP %sp
pascal@168 83 #define CWD cwd
pascal@168 84 #else
pascal@168 85 #define AX %eax
pascal@168 86 #define BX %ebx
pascal@168 87 #define CX %ecx
pascal@168 88 #define DX %edx
pascal@168 89 #define SI %esi
pascal@168 90 #define DI %edi
pascal@168 91 #define BP %ebp
pascal@168 92 #define SP %esp
pascal@168 93 #define CWD cdq
pascal@168 94 #endif
pascal@168 95 /*
pascal@168 96 * LzmaDecode:
pascal@168 97 #ifndef FLAT32
pascal@168 98 * input ds:si=inStream, es:di=outStream
pascal@168 99 * output outStream[], ds:si, es:di
pascal@168 100 .code 16
pascal@168 101 #else
pascal@168 102 * input esi=inStream, edi=outStream
pascal@168 103 * output outStream[], esi, edi
pascal@168 104 .code 32
pascal@168 105 #endif
pascal@168 106 */
pascal@168 107
pascal@168 108 mov $ws1, CX
pascal@168 109 #ifdef ONLY8086
pascal@168 110 movw $2048/2, %ax
pascal@168 111 lzd1:
pascal@168 112 pushw %ax
pascal@168 113 #else
pascal@168 114 lzd1:
pascal@168 115 pushw $2048/2
pascal@168 116 #endif
pascal@168 117 loop lzd1
pascal@168 118 mov SP, BP
pascal@168 119 movb $((LOCALS+3)/4)*2, %cl
pascal@168 120 #ifdef ONLY8086
pascal@168 121 movw $1, %ax
pascal@168 122 cwd
pascal@168 123 initlocals:
pascal@168 124 pushw %dx
pascal@168 125 pushw %ax
pascal@168 126 #else
pascal@168 127 initlocals:
pascal@168 128 pushl $1
pascal@168 129 #endif
pascal@168 130 loop initlocals
pascal@168 131
pascal@168 132 #if !defined(FLAT32) && !defined(FLAT16OUT)
pascal@168 133 movb $4, %cl
pascal@168 134 movw %es, %bx
pascal@168 135 shrw %cl, %bx
pascal@168 136 movw %es, %dx
pascal@168 137 shlw %cl, %dx
pascal@168 138 addw %dx, %di
pascal@168 139 movw %di, outStream(%bp)
pascal@168 140 adcb %bh, outStream+2(%bp)
pascal@168 141 incw %cx
pascal@168 142 #else
pascal@168 143 movb $5, %cl
pascal@168 144 mov DI, outStream(BP)
pascal@168 145 #endif
pascal@168 146
pascal@168 147 // Byte previousByte = 0;
pascal@168 148 xor BX, BX
pascal@168 149
pascal@168 150 // #define RC_INIT(buffer)
pascal@168 151 // Buffer = buffer; Code = 0; Range = 0xFFFFFFFF;
pascal@168 152 // { int i; for(i=0; i<5; i++) { Code = (Code<<8) | RC_READ_BYTE; }}
pascal@168 153 // }
pascal@168 154 // RC_INIT(inStream);
pascal@168 155
pascal@168 156 #ifndef NO_LZMA_HEADER
pascal@168 157 #ifdef CHECK_LZMA_HEADER
pascal@168 158 cmp.w $0x5A4C, (SI) // lzip header ('LZIP' version:1 dicobits:1)
pascal@168 159 je lzip_header
pascal@168 160 cmp.w $0x5D, (SI) // lzma header (0x5D dicosz:4 orgsz:8)
pascal@168 161 jne no_header
pascal@168 162 add $13-6, SI // skip lzma header
pascal@168 163 lzip_header:
pascal@168 164 add $6, SI // skip lzip header
pascal@168 165 no_header:
pascal@168 166 #else
pascal@168 167 add $13, SI // skip lzma header (0x5D dicosz:4 orgsz:8)
pascal@168 168 #endif
pascal@168 169 #endif
pascal@168 170 setrep:
pascal@168 171 call RC_LOAD_BYTE
pascal@168 172 decb Range(BP)
pascal@168 173 loop setrep
pascal@168 174
pascal@168 175 lzdmainlp:
pascal@168 176 // while(1) {
pascal@168 177 // CProb *prob;
pascal@168 178 // int posState = (int)((nowPos) & posStateMask);
pascal@168 179 //
pascal@168 180 // prob = p + IsMatch /*0*/ + (state << kNumPosBitsMax /*4*/) + posState;
pascal@168 181 // if (Bit0(prob)) { /* char */
pascal@168 182
pascal@168 183 xor DX, DX
pascal@168 184 call Bit1state // Bit1(dx + (state << kNumPosBitsMax /*4*/) + posState)
pascal@168 185 mov $state, DI
pascal@168 186 jc lzdstring
pascal@168 187
pascal@168 188 // prob = p + Literal /*1846*/ + (LZMA_LIT_SIZE /*768*/ *
pascal@168 189 // ((((nowPos) & literalPosMask) << lc) + (previousByte >> (8 - lc))));
pascal@168 190
pascal@168 191 #if PROP_LC != 0
pascal@168 192 # ifdef ONLY8086
pascal@168 193 movb $8-PROP_LC, %cl
pascal@168 194 shrb %cl, %bl
pascal@168 195 # else
pascal@168 196 shrb $8-PROP_LC, %bl
pascal@168 197 # endif
pascal@170 198 #else
pascal@170 199 xor %bx,%bx
pascal@168 200 #endif
pascal@168 201
pascal@168 202 #if PROP_LP != 0
pascal@168 203 movb posState2(BP), %dl
pascal@170 204 # if PROP_LC != 0
pascal@170 205 # ifdef ONLY8086
pascal@168 206 movb $PROP_LC, %cl
pascal@168 207 shl %cl, DX
pascal@170 208 # else
pascal@168 209 shl $PROP_LC, DX
pascal@170 210 # endif
pascal@170 211 movb $0, %bh
pascal@168 212 # endif
pascal@168 213 add BX, DX
pascal@168 214 #endif
pascal@168 215
pascal@168 216 movb $3, %ah
pascal@168 217 mul BX // dx = 3*bh
pascal@168 218 add $1846, AX
pascal@168 219
pascal@168 220 // int symbol = 1;
pascal@168 221
pascal@168 222 CWD
pascal@168 223 inc DX // symbol = 1
pascal@168 224 xchg AX, CX // save prob
pascal@168 225
pascal@168 226 // if (state >= kNumLitStates /*7*/) { /* previous was string */
pascal@168 227 // if (state < 4) state = 0;
pascal@168 228
pascal@168 229 lzd6z:
pascal@168 230 subb $3, (BP, DI)
pascal@168 231
pascal@168 232 // if (state < 4) state = 0;
pascal@168 233
pascal@168 234 jnc lzd6
pascal@168 235 movb %dh, (BP, DI) // %dh = 0
pascal@168 236
pascal@168 237 lzd6:
pascal@168 238 // else if (state < 10) state -= 3;
pascal@168 239
pascal@168 240 cmpb $10-3, (BP, DI)
pascal@168 241
pascal@168 242 // else state -= 6;
pascal@168 243
pascal@168 244 jnb lzd6z
pascal@168 245 cmpb $7-3-1, (BP, DI)
pascal@168 246 jbe lzd3
pascal@168 247
pascal@168 248 // int matchByte = outStream[nowPos - rep0];
pascal@168 249
pascal@168 250 call DicoRep02ESDI // %bl = outStream[nowPos - rep0];
pascal@168 251
pascal@168 252 // do {
pascal@168 253 // int bit;
pascal@168 254 // CProb *probLit;
pascal@168 255 // matchByte <<= 1; bit = (matchByte & 0x100);
pascal@168 256
pascal@168 257 movb $1, %bh
pascal@168 258 lzd4:
pascal@168 259 shlb $1, %bl // matchByte <<= 1
pascal@168 260 sbb DI, DI // save bit=C
pascal@168 261
pascal@168 262 // probLit = prob + 0x100 + bit + symbol;
pascal@168 263
pascal@168 264 mov CX, AX // restore prob
pascal@168 265 adcb %bh, %ah // + bit + 0x100
pascal@168 266
pascal@168 267 // RC_GET_BIT2(probLit, symbol, if (bit) break, if (!bit) break)
pascal@168 268
pascal@168 269 call Bit1axdx // C,%ax = Bit1(prob+%ax)
pascal@168 270 rclb $1, %dl // symbol <<= 1; symbol |= C
pascal@168 271 jc lzd5 // if symbol >= 0x100
pascal@168 272 cmp DI, AX
pascal@168 273 jz lzd4 // if bit == Bit1(prob+%ax)
pascal@168 274
pascal@168 275 // } while (symbol < 0x100);
pascal@168 276 // }
pascal@168 277 lzd3:
pascal@168 278 // while (symbol < 0x100) {
pascal@168 279 // CProb *probLit = prob + symbol;
pascal@168 280 // RC_GET_BIT(probLit, symbol)
pascal@168 281 // }
pascal@168 282
pascal@168 283 xor BX, BX
pascal@168 284 jmp lzd4
pascal@168 285 lzd5:
pascal@168 286
pascal@168 287 // outStream[nowPos++] = previousByte = (Byte)symbol;
pascal@168 288
pascal@168 289 xchg AX, DX
pascal@168 290 call outchar // %bl = outStream[nowPos++] = %al;
pascal@168 291 jmp lzdmainlp
pascal@168 292
pascal@168 293 // }
pascal@168 294
pascal@168 295 lzdstring:
pascal@168 296 mov $1, CX
pascal@168 297
pascal@168 298 // else { /* string */
pascal@168 299 // prob = p + IsRep /*192*/ + state;
pascal@168 300
pascal@168 301 movb $192, %dl
pascal@168 302 addb (BP, DI), %dl
pascal@168 303 mov $rep0, DI
pascal@168 304
pascal@168 305 // if (Bit0(prob)) {
pascal@168 306
pascal@168 307 call Bit1dx // Bit1(prob)
pascal@168 308 jc lzd8
pascal@168 309
pascal@168 310 // rep3 = rep2; rep2 = rep1; rep1 = rep0;
pascal@168 311 // state = (state < kNumLitStates /*7*/) ? 0 : 3;
pascal@168 312
pascal@168 313 stc
pascal@168 314
pascal@168 315 // prob = p + LenCoder /*818*/;
pascal@168 316
pascal@168 317 mov $818, DX
pascal@168 318
pascal@168 319 // }
pascal@168 320
pascal@168 321 jmp lzd11a
pascal@168 322
pascal@168 323 // else {
pascal@168 324 lzd8:
pascal@168 325 // prob += kNumStates /*12*/;
pascal@168 326 // if (Bit0(prob)) {
pascal@168 327 call Bit1dx12 // prob += 12; Bit1(prob)
pascal@168 328 jc lzd11
pascal@168 329 // prob = p + IsRep0Long /*240*/ + (state << kNumPosBitsMax /*4*/)
pascal@168 330 // + posState;
pascal@168 331 movb $240, %dl // dh=0
pascal@168 332
pascal@168 333 // if (Bit0(prob)) {
pascal@168 334
pascal@168 335 call Bit1state // Bit1(dx + (state << kNumPosBitsMax /*4*/) + posState)
pascal@168 336 jc lzd12
pascal@168 337
pascal@168 338 // // if (nowPos == 0) return LZMA_RESULT_DATA_ERROR;
pascal@168 339 // state = (state < kNumLitStates /*7*/) ? 9 : 11;
pascal@168 340
pascal@168 341 movb $9, %dl
pascal@168 342
pascal@168 343 // len++; goto string;
pascal@168 344 jmp lzd13string // ax = 0
pascal@168 345 // }
pascal@168 346 // }
pascal@168 347 // else {
pascal@168 348 lzd11:
pascal@168 349 // UInt32 distance = rep1;
pascal@168 350 // prob += kNumStates /*12*/;
pascal@168 351 // if (!Bit0(prob)) {
pascal@168 352
pascal@168 353 call Bit1dx12 // prob += 12; Bit1(prob)
pascal@168 354 jnc lzd11z
pascal@168 355
pascal@168 356 // prob += kNumStates /*12*/;
pascal@168 357 // if (Bit0(prob)) distance = rep2;
pascal@168 358
pascal@168 359 call Bit1dx12 // prob += 12; Bit1(prob)
pascal@168 360 lzd11a:
pascal@168 361 adcb %cl, %cl
pascal@168 362
pascal@168 363 // else { distance = rep3; rep3 = rep2; }
pascal@168 364 // rep2 = rep1;
pascal@168 365 // }
pascal@168 366 // rep1 = rep0; rep0 = distance;
pascal@168 367
pascal@168 368 lzd11z:
pascal@168 369 # ifdef ONLY8086
pascal@168 370 shl $1, CX
pascal@168 371 shl $1, CX // 8->32 bits
pascal@168 372 sub CX, DI // &rep[cx]
pascal@168 373 movw (BP, DI), %ax
pascal@168 374 pushw 2(BP, DI)
pascal@168 375 rotreplp:
pascal@168 376 movb 4(BP, DI), %bl
pascal@168 377 movb %bl, (BP, DI)
pascal@168 378 inc DI
pascal@168 379 loop rotreplp
pascal@168 380 popw %bx
pascal@168 381 testb %dh, %dh
pascal@168 382 jnz lzd10
pascal@168 383 movw %ax, (BP, DI)
pascal@168 384 movw %bx, 2(BP, DI)
pascal@168 385 # else
pascal@168 386 shl $2, CX // 8->32 bits
pascal@168 387 sub CX, DI // &rep[cx]
pascal@168 388 movl (BP, DI), %eax
pascal@168 389 rotreplp:
pascal@168 390 movb 4(BP, DI), %bl
pascal@168 391 movb %bl, (BP, DI)
pascal@168 392 inc DI
pascal@168 393 loop rotreplp
pascal@168 394 testb %dh, %dh
pascal@168 395 jnz lzd10
pascal@168 396 movl %eax, (BP, DI)
pascal@168 397 # endif
pascal@168 398
pascal@168 399 // }
pascal@168 400 lzd12:
pascal@168 401 // state = (state < kNumLitStates /*7*/) ? 8 : 11;
pascal@168 402
pascal@168 403 movb $0x08, %cl
pascal@168 404
pascal@168 405 // prob = p + RepLenCoder /*1332*/;
pascal@168 406
pascal@168 407 mov $1332, DX
pascal@168 408
pascal@168 409 // }
pascal@168 410 lzd10:
pascal@168 411 push CX // CX = 0
pascal@168 412
pascal@168 413 // { /* get len */
pascal@168 414 // int numBits, offset;
pascal@168 415 // CProb *probLen = prob + LenChoice /*0*/;
pascal@168 416 // numBits = kLenNumLowBits /*3*/;
pascal@168 417
pascal@168 418 movb $8, %cl // numBits : 3,3,8
pascal@168 419
pascal@168 420 // if (Bit0(probLen)) {
pascal@168 421
pascal@168 422 call Bit1dx // Bit1(prob)
pascal@168 423 xchg AX, BX
pascal@168 424 inc DX
pascal@168 425 jnc lzd15 // bx=0
pascal@168 426
pascal@168 427 // probLen = prob + LenLow/*2*/ + (posState << kLenNumLowBits/*3*/);
pascal@168 428 // offset = 0;
pascal@168 429 // }
pascal@168 430 // else {
pascal@168 431 // probLen = prob + LenChoice2 /*1*/;
pascal@168 432
pascal@168 433 call Bit1dx // Bit1(prob)
pascal@168 434 add AX, BX
pascal@168 435
pascal@168 436 #if PROP_PB != 0
pascal@168 437 inc AX // ah=0
pascal@168 438 #endif
pascal@168 439 jc lzd16 // %ax=0, %bx=-2
pascal@168 440 lzd15:
pascal@168 441 #if PROP_PB != 0
pascal@168 442 movb $8, %al
pascal@168 443 mulb posState(BP)
pascal@168 444 #endif
pascal@168 445
pascal@168 446 // if (Bit0(probLen)) {
pascal@168 447 // probLen = prob + LenMid/*130*/ + (posState << kLenNumMidBits/*3*/);
pascal@168 448
pascal@168 449 movb $3, %cl // numBits : 3,3,8
pascal@168 450 lzd16:
pascal@168 451 #if PROP_PB != 0
pascal@168 452 add $2-128-1, AX // probLen : 2,130,258
pascal@168 453 #else
pascal@168 454 mov $2-128-1, AX // probLen : 2,130,258
pascal@168 455 #endif
pascal@168 456 add DX, AX
pascal@168 457 mov $-8+1, DX // offset : 0,8,16
pascal@168 458 lzdargslp:
pascal@168 459 add $8, DX
pascal@168 460 add $128, AX
pascal@168 461 inc BX
pascal@168 462 jle lzdargslp // leave with bx=1
pascal@168 463
pascal@168 464 // offset = kLenNumLowSymbols /*8*/;
pascal@168 465 // //numBits = kLenNumMidBits /*3*/;
pascal@168 466 // }
pascal@168 467 // else {
pascal@168 468 // probLen = prob + LenHigh /*258*/;
pascal@168 469 // offset = kLenNumLowSymbols /*8*/ + kLenNumMidSymbols /*8*/;
pascal@168 470 // numBits = kLenNumHighBits /*8*/;
pascal@168 471 // }
pascal@168 472 // }
pascal@168 473 // RangeDecoderBitTreeDecode(probLen, numBits, len); len += offset;
pascal@168 474
pascal@168 475 push DX
pascal@168 476 call RangeDecoder // %ax=probs, %cx=numLevels, %ax=res
pascal@168 477 pop DX
pascal@168 478 add DX, AX // offset
pascal@168 479 pop DX // 0
pascal@168 480 lzd13string:
pascal@168 481 push AX
pascal@168 482
pascal@168 483 // state = (state < kNumLitStates /*7*/) ? dl : dl|3;
pascal@168 484
pascal@168 485 movb $7, %cl
pascal@168 486 cmpb %cl, state(BP)
pascal@168 487 jb new_state
pascal@168 488 orb $3, %dl
pascal@168 489 new_state:
pascal@168 490 movb %dl, state(BP)
pascal@168 491
pascal@168 492 // } /* get len */
pascal@168 493 // if (state < 4) {
pascal@168 494
pascal@168 495 cmpb $4-1, %dl
pascal@168 496 ja lzd19
pascal@168 497
pascal@168 498 // int posSlot;
pascal@168 499 // state += kNumLitStates /*7*/;
pascal@168 500
pascal@168 501 addb %cl, state(BP)
pascal@168 502
pascal@168 503 // prob = p + PosSlot /*432*/ + (((len < kNumLenToPosStates /*4*/) ?
pascal@168 504 // len : kNumLenToPosStates - 1) << kNumPosSlotBits /*6*/);
pascal@168 505
pascal@168 506 cmp $4+1, AX
pascal@168 507 jb lzd21
pascal@168 508 mov $3+1, AX
pascal@168 509
pascal@168 510 lzd21:
pascal@168 511
pascal@168 512 dec CX // cx = 6
pascal@168 513 shl %cl, AX
pascal@168 514 add $432-64, AX
pascal@168 515
pascal@168 516 // RangeDecoderBitTreeDecode(prob, kNumPosSlotBits /*6*/, posSlot);
pascal@168 517
pascal@168 518 call RangeDecoder // %ax=probs, %cx=numLevels, %ax=res
pascal@168 519
pascal@168 520 // if (posSlot >= kStartPosModelIndex /*4*/) {
pascal@168 521 // int numDirectBits = ((posSlot >> 1) - 1);
pascal@168 522
pascal@168 523 #ifndef FLAT32
pascal@168 524 movw %cx, 2(%bp, %di) // %cx = 0
pascal@168 525 #endif
pascal@168 526 mov AX, (BP, DI)
pascal@168 527 mov AX, CX
pascal@168 528 shrw $1, CX
pascal@168 529 dec CX
pascal@168 530 cmpb $4, %al
pascal@168 531 jb lzd22
pascal@168 532
pascal@168 533 // rep0 = (2 | ((UInt32)posSlot & 1));
pascal@168 534
pascal@168 535 andb %bl, (BP, DI) // %bx=1
pascal@168 536 orb $2, (BP, DI)
pascal@168 537
pascal@168 538 // if (posSlot < kEndPosModelIndex /*14*/) {
pascal@168 539
pascal@168 540 cmpb $14, %al
pascal@168 541 jnb lzd23
pascal@168 542
pascal@168 543 // rep0 <<= numDirectBits;
pascal@168 544
pascal@168 545 neg AX
pascal@168 546 # ifdef ONLY8086
pascal@168 547 pushw %cx
pascal@168 548 movb $0, %ch
pascal@168 549 shllrep0:
pascal@168 550 shlw $1, (BP, DI)
pascal@168 551 rclw $1, 2(BP, DI)
pascal@168 552 loop shllrep0
pascal@168 553 popw %cx
pascal@168 554 # else
pascal@168 555 shll %cl, (BP, DI)
pascal@168 556 # endif
pascal@168 557 add (BP, DI), AX
pascal@168 558
pascal@168 559 // prob = p + SpecPos /*688*/ + rep0 - posSlot - 1;
pascal@168 560
pascal@168 561 add $687, AX
pascal@168 562 jmp lzd24
pascal@168 563
pascal@168 564 // }
pascal@168 565 // else {
pascal@168 566 lzd23:
pascal@168 567 // numDirectBits -= kNumAlignBits /*4*/;
pascal@168 568 // do {
pascal@168 569 // RC_NORMALIZE; Range >>= 1; rep0 <<= 1;
pascal@168 570 // if (Code >= Range) { Code -= Range; rep0 |= 1; }
pascal@168 571
pascal@168 572 lzd23z:
pascal@168 573 call RC_NORMALIZE
pascal@168 574 # ifdef ONLY8086
pascal@168 575 pushw %dx
pascal@168 576 shrw $1, Range+2(BP)
pascal@168 577 rcrw $1, Range(BP)
pascal@168 578 movw Range(BP), %ax
pascal@168 579 movw Range+2(BP), %dx
pascal@168 580 cmpw Code+2(BP), %dx
pascal@168 581 ja lzd25
pascal@168 582 jb lzd25x
pascal@168 583 cmpw Code(BP), %ax
pascal@168 584 ja lzd25
pascal@168 585 lzd25x:
pascal@168 586 subw %ax, Code(BP)
pascal@168 587 sbbw %dx, Code+2(BP)
pascal@168 588 stc
pascal@168 589 lzd25:
pascal@168 590 popw %dx
pascal@168 591 rclw $1, (BP, DI)
pascal@168 592 rclw $1, 2(BP, DI)
pascal@168 593 # else
pascal@168 594 shrl $1, Range(BP)
pascal@168 595 movl Range(BP), %eax
pascal@168 596 cmpl Code(BP), %eax
pascal@168 597 ja lzd25
pascal@168 598 subl %eax, Code(BP)
pascal@168 599 stc
pascal@168 600 lzd25:
pascal@168 601 rcll $1, (BP, DI)
pascal@168 602 # endif
pascal@168 603
pascal@168 604 // } while (--numDirectBits != 0);
pascal@168 605
pascal@168 606 cmpb $4+1, %cl
pascal@168 607 loopne lzd23z
pascal@168 608
pascal@168 609 // prob = p + Align /* 802 */; numDirectBits = kNumAlignBits /*4*/;
pascal@168 610 // rep0 <<= numDirectBits;
pascal@168 611
pascal@168 612 # ifdef ONLY8086
pascal@168 613 pushw %cx
pascal@168 614 movb $0, %ch
pascal@168 615 shlrep0:
pascal@168 616 shlw $1, (BP, DI)
pascal@168 617 rclw $1, 2(BP, DI)
pascal@168 618 loop shlrep0
pascal@168 619 popw %cx
pascal@168 620 # else
pascal@168 621 shll %cl, (BP, DI)
pascal@168 622 # endif
pascal@168 623 mov $802, AX
pascal@168 624 // }
pascal@168 625
pascal@168 626 lzd24:
pascal@168 627 call RangeDecoder // %ax=probs, %cx=numLevels, %ax=res
pascal@168 628
pascal@168 629 // {
pascal@168 630 // int i = 1, mi = 1;
pascal@168 631 // do {
pascal@168 632 // CProb *prob3 = prob + mi;
pascal@168 633 // RC_GET_BIT2(prob3, mi, ; , rep0 |= i);
pascal@168 634
pascal@168 635 orb %dh, (BP, DI) // update rep0 with DirectBits
pascal@168 636
pascal@168 637 // i <<= 1;
pascal@168 638 // } while(--numDirectBits != 0);
pascal@168 639 // }
pascal@168 640 // } else rep0 = posSlot;
pascal@168 641 lzd22:
pascal@168 642 // if (++rep0 == (UInt32)(0)) break; /* EOF */
pascal@168 643
pascal@168 644 # ifdef ONLY8086
pascal@168 645 incw (BP, DI)
pascal@168 646 jnz lzd19
pascal@168 647 incw 2(BP, DI)
pascal@168 648 # else
pascal@168 649 incl (BP, DI)
pascal@168 650 # endif
pascal@168 651
pascal@168 652 lzd19:
pascal@168 653 pop CX
pascal@168 654 jz lzdone
pascal@168 655
pascal@168 656 // }
pascal@168 657 // len += kMatchMinLen;/*2*/
pascal@168 658
pascal@168 659 inc CX
pascal@168 660
pascal@168 661 // string: // if (rep0 > nowPos) return LZMA_RESULT_DATA_ERROR;
pascal@168 662 // do {
pascal@168 663 lzd13z:
pascal@168 664 // previousByte = outStream[nowPos - rep0];
pascal@168 665 // outStream[nowPos++] = previousByte;
pascal@168 666
pascal@168 667 call outcharDico // %bl = outStream[nowPos++] = outStream[nowPos - rep0]
pascal@168 668
pascal@168 669 // } while(--len != 0);
pascal@168 670
pascal@168 671 loop lzd13z
pascal@168 672
pascal@168 673 // } /* char/string */
pascal@168 674 // }
pascal@168 675
pascal@168 676 jmp lzdmainlp
pascal@168 677
pascal@168 678 lzdone:
pascal@168 679 // //RC_NORMALIZE;
pascal@168 680 // //*inSizeProcessed = (SizeT)(Buffer - inStream); *outSizeProcessed = nowPos;
pascal@168 681 // return LZMA_RESULT_OK;
pascal@168 682 call Dico2ESDI // set es & di (rep0 = 0)
pascal@168 683 lea ws2(BP), SP // dealloc
pascal@168 684 ret
pascal@168 685 // }
pascal@168 686
pascal@168 687 // al = outStream[nowPos - rep0];
pascal@168 688
pascal@168 689 /*
pascal@168 690 * output es:di, al
pascal@168 691 * scratch bh, cl, flags
pascal@168 692 */
pascal@168 693
pascal@168 694 DicoRep02ESDI:
pascal@168 695 stc
pascal@168 696
pascal@168 697 // bl = outStream[nowPos];
pascal@168 698
pascal@168 699 /*
pascal@168 700 * output es:di, bl
pascal@168 701 * scratch bh, cl, flags
pascal@168 702 */
pascal@168 703
pascal@168 704 Dico2ESDI:
pascal@168 705 #if !defined(FLAT32) && !defined(FLAT16OUT)
pascal@168 706 # ifdef ONLY8086
pascal@168 707 pushw %ax
pascal@168 708 movw nowPos(%bp), %bx
pascal@168 709 movw nowPos+2(%bp), %ax
pascal@168 710 jnc Dico2ESDIz
pascal@168 711 subw rep0(%bp), %bx
pascal@168 712 sbbw rep0+2(%bp), %ax
pascal@168 713 Dico2ESDIz:
pascal@168 714 movw $0xF, %di
pascal@168 715 andw %bx, %di
pascal@168 716 pushw %cx
pascal@168 717 movb $4, %cl
pascal@168 718 shrw %cl, %bx
pascal@168 719 shlw %cl, %ax
pascal@168 720 popw %cx
pascal@168 721 addb %al, %bh
pascal@168 722 popw %ax
pascal@168 723 # else
pascal@168 724 movl nowPos(%bp), %ebx
pascal@168 725 jnc Dico2ESDIz
pascal@168 726 subl rep0(%bp), %ebx
pascal@168 727 Dico2ESDIz:
pascal@168 728 movw %bx, %di
pascal@168 729 xorw %bx, %bx
pascal@168 730 shrl $4, %ebx
pascal@168 731 # endif
pascal@168 732 movw %bx, %es
pascal@168 733 #else
pascal@168 734 mov nowPos(BP), DI
pascal@168 735 jnc Dico2ESDIz
pascal@168 736 sub rep0(BP), DI
pascal@168 737 Dico2ESDIz:
pascal@168 738 #endif
pascal@168 739 #ifdef FLAT32
pascal@168 740 movb (DI), %bl
pascal@168 741 #else
pascal@168 742 movb %es:(%di), %bl
pascal@168 743 #endif
pascal@168 744 ret
pascal@168 745
pascal@168 746 outcharDico:
pascal@168 747
pascal@168 748 // bl = outStream[nowPos++] = outStream[nowPos - rep0]
pascal@168 749
pascal@168 750 /*
pascal@168 751 * output es:di, bl
pascal@168 752 * update nowPos
pascal@168 753 * scratch ax, dx, bh, cl, flags
pascal@168 754 */
pascal@168 755
pascal@168 756 call DicoRep02ESDI // %bl = outStream[nowPos - rep0]
pascal@168 757 xchg AX, BX
pascal@168 758 outchar:
pascal@168 759
pascal@168 760 // bl = outStream[nowPos++] = previousByte = al;
pascal@168 761
pascal@168 762 /*
pascal@168 763 * output bl
pascal@168 764 * update nowPos
pascal@168 765 * scratch ax, dx, bh, di, cl, flags
pascal@168 766 */
pascal@168 767
pascal@168 768 clc
pascal@168 769 call Dico2ESDI
pascal@168 770 stosb
pascal@168 771 xchg AX, BX // previous byte
pascal@168 772
pascal@168 773 // int posState = (int)((nowPos) & posStateMask);
pascal@168 774
pascal@168 775 #if PROP_PB != 0 && PROP_LP != 0
pascal@168 776 addw $0x0101, posState2(BP)
pascal@168 777 andb $(((1 << PROP_PB) -1)<<8)+((1 << PROP_LP) -1), posState2(BP)
pascal@168 778 #else
pascal@168 779 # if PROP_PB != 0
pascal@168 780 incb posState(BP)
pascal@168 781 andb $((1 << PROP_PB) -1), posState(BP)
pascal@168 782 # endif
pascal@168 783 # if PROP_LP != 0
pascal@168 784 incb posState2(BP)
pascal@168 785 andb $((1 << PROP_LP) -1), posState2(BP)
pascal@168 786 # endif
pascal@168 787 #endif
pascal@168 788 #ifdef ONLY8086
pascal@168 789 incw nowPos(BP)
pascal@168 790 jnz incnowPosDone
pascal@168 791 incw nowPos+2(BP)
pascal@168 792 incnowPosDone:
pascal@168 793 #else
pascal@168 794 incl nowPos(BP)
pascal@168 795 #endif
pascal@168 796 ret
pascal@168 797
pascal@168 798 //
pascal@168 799 // #define RC_NORMALIZE if (Range < kTopValue)
pascal@168 800 // { Range <<= 8; Code = (Code << 8) | RC_READ_BYTE; }
pascal@168 801
pascal@168 802 /*
pascal@168 803 * update Range, Code, ds:si
pascal@168 804 * scratch flags
pascal@168 805 */
pascal@168 806
pascal@168 807 RC_NORMALIZE:
pascal@168 808 cmpb $0, Range+3(BP)
pascal@168 809 jne RC_NORMALIZE_1
pascal@168 810 RC_LOAD_BYTE:
pascal@168 811 push AX
pascal@168 812 #ifdef ONLY8086
pascal@168 813 movw Range+1(BP), %ax
pascal@168 814 movw %ax, Range+2(BP)
pascal@168 815 movw Code+1(BP), %ax
pascal@168 816 movw %ax, Code+2(BP)
pascal@168 817 xorw %ax, %ax
pascal@168 818 movb Range(BP), %ah
pascal@168 819 movw %ax, Range(BP)
pascal@168 820 movb Code(BP), %ah
pascal@168 821 movw %ax, Code(BP)
pascal@168 822 #else
pascal@168 823 shll $8, Range(BP)
pascal@168 824 shll $8, Code(BP)
pascal@168 825 #endif
pascal@168 826 #if !defined(FLAT16) && !defined(FLAT32)
pascal@168 827 testw %si, %si
pascal@168 828 jns RC_READ_BYTE
pascal@168 829 movw %ds, %ax
pascal@168 830 incw %ax
pascal@168 831 movw %ax, %ds
pascal@168 832 addw $-16, %si
pascal@168 833 RC_READ_BYTE:
pascal@168 834 #endif
pascal@168 835 lodsb
pascal@168 836 movb %al, Code(BP)
pascal@168 837 pop AX
pascal@168 838 RC_NORMALIZE_1:
pascal@168 839 ret
pascal@168 840
pascal@168 841 // Bit1(dx + (state << kNumPosBitsMax /*4*/) + posState)
pascal@168 842
pascal@168 843 Bit1state:
pascal@168 844 movb $16, %al
pascal@168 845 mulb state(BP)
pascal@168 846 # if PROP_PB != 0
pascal@168 847 addb posState(BP), %al
pascal@168 848 # endif
pascal@168 849 Bit1axdx:
pascal@168 850 add DX, AX
pascal@168 851 jmp Bit1
pascal@168 852
pascal@168 853 // prob += 12; Bit1(prob)
pascal@168 854
pascal@168 855 Bit1dx12:
pascal@168 856 add $12, DX
pascal@168 857 Bit1dx:
pascal@168 858 mov DX, AX
pascal@168 859
pascal@168 860 // static int Bit1(CProb *p)
pascal@168 861
pascal@168 862 Bit1:
pascal@168 863 /*
pascal@168 864 * input ax=p
pascal@168 865 * output C, ax
pascal@168 866 * update bound, Range, Code, ds:si
pascal@168 867 * scratch flags
pascal@168 868 */
pascal@168 869
pascal@168 870 // {
pascal@168 871 // RC_NORMALIZE;
pascal@168 872
pascal@168 873 call RC_NORMALIZE // kill %ax, update %si
pascal@168 874
pascal@168 875 #ifdef ONLY8086
pascal@168 876 pushw %ax
pascal@168 877 pushw %cx
pascal@168 878 pushw %dx
pascal@168 879 pushw %di
pascal@168 880 #else
pascal@168 881 pushal
pascal@168 882 #endif
pascal@168 883
pascal@168 884 xchg AX, DI
pascal@168 885 add DI, DI // short *
pascal@168 886
pascal@168 887
pascal@168 888 // bound = (Range>>kNumBitModelTotalBits /*11*/) * *(p);
pascal@168 889
pascal@168 890 #ifdef ONLY8086
pascal@168 891 movw Range(BP), %dx
pascal@168 892 movw Range+2(BP), %ax
pascal@168 893 movw $11, %cx
pascal@168 894 shr11lp:
pascal@168 895 shrw $1, %ax
pascal@168 896 rcrw $1, %dx
pascal@168 897 loop shr11lp
pascal@168 898 movw %dx, %cx
pascal@168 899 mulw (BP, DI)
pascal@168 900 xchgw %ax, %cx
pascal@168 901 mulw (BP, DI)
pascal@168 902 addw %cx, %dx
pascal@168 903 #else
pascal@168 904 movl Range(BP), %eax
pascal@168 905 shrl $11, %eax
pascal@168 906 movzwl (BP, DI), %edx
pascal@168 907 mull %edx
pascal@168 908 #endif
pascal@168 909
pascal@168 910 // if (Code < bound) {
pascal@168 911
pascal@168 912 #ifdef ONLY8086
pascal@168 913 cmpw Code+2(BP), %dx
pascal@168 914 jb Bit1_1
pascal@168 915 ja Bit1_1x
pascal@168 916 cmpw Code(BP), %ax
pascal@168 917 jbe Bit1_1
pascal@168 918 Bit1_1x:
pascal@168 919
pascal@168 920 // Range = bound;
pascal@168 921
pascal@168 922 movw %ax, Range(BP)
pascal@168 923 movw %dx, Range+2(BP)
pascal@168 924 #else
pascal@168 925 cmpl Code(BP), %eax
pascal@168 926 jbe Bit1_1
pascal@168 927
pascal@168 928 // Range = bound;
pascal@168 929
pascal@168 930 movl %eax, Range(BP)
pascal@168 931 #endif
pascal@168 932
pascal@168 933 // *(p) += (kBitModelTotal /*2048*/ - *(p)) >> kNumMoveBits /*5*/;
pascal@168 934
pascal@168 935 movw $2048, %ax
pascal@168 936
pascal@168 937 // return 0;
pascal@168 938
pascal@168 939 jmp Bit1_2
pascal@168 940
pascal@168 941 // }
pascal@168 942 // else {
pascal@168 943
pascal@168 944 Bit1_1:
pascal@168 945
pascal@168 946 // Range -= bound; Code -= bound;
pascal@168 947
pascal@168 948 #ifdef ONLY8086
pascal@168 949 subw %ax, Range(BP)
pascal@168 950 sbbw %dx, Range+2(BP)
pascal@168 951 subw %ax, Code(BP)
pascal@168 952 sbbw %dx, Code+2(BP)
pascal@168 953 #else
pascal@168 954 subl %eax, Range(BP)
pascal@168 955 subl %eax, Code(BP)
pascal@168 956 #endif
pascal@168 957
pascal@168 958 // *(p) -= (*(p)) >> kNumMoveBits /*5*/;
pascal@168 959
pascal@168 960 movw $31, %ax
pascal@168 961
pascal@168 962 // return 1;
pascal@168 963
pascal@168 964 stc
pascal@168 965 Bit1_2:
pascal@168 966 pushf
pascal@168 967 subw (BP, DI), %ax
pascal@168 968 #ifdef ONLY8086
pascal@168 969 movb $5, %cl
pascal@168 970 sarw %cl, %ax
pascal@168 971 #else
pascal@168 972 sarw $5, %ax
pascal@168 973 #endif
pascal@168 974 addw %ax, (BP, DI)
pascal@168 975 popf
pascal@168 976 #ifdef ONLY8086
pascal@168 977 popw %di
pascal@168 978 popw %dx
pascal@168 979 popw %cx
pascal@168 980 popw %ax
pascal@168 981 #else
pascal@168 982 popal
pascal@168 983 #endif
pascal@168 984 sbb AX, AX
pascal@168 985
pascal@168 986 // }
pascal@168 987 // }
pascal@168 988
pascal@168 989 ret
pascal@168 990
pascal@168 991 RangeDecoder:
pascal@168 992
pascal@168 993 /*
pascal@168 994 * input ax=probs cx=numLevels (< 8) bx=1
pascal@168 995 * output ax=res (backward), dh (forward)
pascal@168 996 * update bound, Range, Code, ds:si
pascal@168 997 * scratch flags, cx=0, dl
pascal@168 998 */
pascal@168 999
pascal@168 1000 push BX
pascal@168 1001
pascal@168 1002 // { int i = numLevels; res = 1;
pascal@168 1003 mov BX, DX // res = 1
pascal@168 1004
pascal@168 1005 // do { CProb *p = probs + res; RC_GET_BIT(p, res) } while(--i != 0);
pascal@168 1006
pascal@168 1007 RangeDecoder_1:
pascal@168 1008 push AX
pascal@168 1009 call Bit1axdx // C,%ax = Bit1(prob+%ax)
pascal@168 1010 rclb $1, %dl // res <<= 1; res |= C
pascal@168 1011 andb %bl, %al // current bit
pascal@168 1012 orb %al, %bh // store in bh
pascal@168 1013 shlb $1, %bl // update max
pascal@168 1014 pop AX
pascal@168 1015 loop RangeDecoder_1
pascal@168 1016
pascal@168 1017 // res -= (1 << numLevels); }
pascal@168 1018
pascal@168 1019 xchg AX, BX // move bh to dh
pascal@168 1020 xchg AX, DX // and dl to al
pascal@168 1021 sub %dl, %al // sub max
pascal@168 1022 pop BX
pascal@168 1023 ret