home *** CD-ROM | disk | FTP | other *** search
/ Meeting Pearls 3 / Meeting_Pearls_III.iso / Pearls / texmf / source / dvips / repack.c < prev    next >
C/C++ Source or Header  |  1994-03-29  |  14KB  |  464 lines
  1. /*
  2.  *   Compressed TeX fonts in PostScript.
  3.  *   By Radical Eye Software.
  4.  *   (Slight mods by Don Knuth in December 89.)
  5.  */
  6. #include "dvips.h" /* The copyright notice in that file is included too! */
  7. #ifdef DEBUG
  8. extern integer debug_flag;
  9. #endif /* DEBUG */
  10.  
  11. /*   Given a raster that has been unpacked from PK format,
  12.  *   we compress it by another scheme that is suitable for
  13.  *   a PostScript-oriented unpacking. We write instructions for
  14.  *   a little interpreter whose one-byte instructions have the form
  15.  *   18*opcode+parameter. The interpreter forms each row based
  16.  *   on simple transformations to the previous row. */
  17.  
  18. #define MAXOUT (18)
  19. #define CMD(n) (MAXOUT*(n))
  20. #define ADVXCHG1 CMD(0)
  21. #define ADVXCHG1END CMD(1)
  22. #define CHGX CMD(2)-1
  23. #define CHGXEND CMD(3)-1
  24. #define ADVXLSH CMD(4)
  25. #define ADVXLSHEND CMD(5)
  26. #define ADVXRSH CMD(6)
  27. #define ADVXRSHEND CMD(7)
  28. #define ADVX CMD(8)-1
  29. #define REPX CMD(9)-1
  30. #define SETX CMD(10)-1
  31. #define CLRX CMD(11)-1
  32. #define ADVXCHG2 CMD(12)
  33. #define ADVXCHG2END CMD(13)
  34. #define END CMD(14)
  35.  
  36. #ifdef AMIGA
  37. #include "repack_protos.h"
  38. #include "dvips_protos.h"
  39. #include "unpack_protos.h"
  40. static void addts(register unsigned char);
  41. #else
  42. extern void error() ;
  43. extern long getlong() ;
  44. extern long unpack() ;
  45. #endif
  46.  
  47. static int rowlength = -1 ;
  48. static unsigned char *specdata ;
  49. static long tslen = 0 ;
  50. static unsigned char *tempstore, *tsp, *tsend ;
  51.  
  52. void putlong(a, i)
  53. register char *a ;
  54. long i ;
  55. {
  56.    a[0] = i >> 24 ;
  57.    a[1] = i >> 16 ;
  58.    a[2] = i >> 8 ;
  59.    a[3] = i ;
  60. }
  61.  
  62. long getlong(a)
  63. register unsigned char *a ;
  64. {
  65.    return ((((((a[0] << 8L) + a[1]) << 8L) + a[2]) << 8L) + a[3]) ;
  66. }
  67.  
  68. /* First, a routine that appends one byte to the compressed output. */
  69.  
  70. #define addtse(n) {lcm=tsp-tempstore;addts(n);} /* mark END option position */
  71.  
  72. static void addts(what)
  73. register unsigned char what ;
  74. {
  75.    register unsigned char *p, *q ;
  76.  
  77.    if (tsp >= tsend) {
  78.       if (tempstore == NULL) {
  79.          tslen = 2020 ;
  80.          tempstore = (unsigned char *)mymalloc((integer)tslen) ;
  81.          tsp = tempstore ;
  82.       } else {
  83.          tslen = 2 * tslen ;
  84.          tsp = (unsigned char *)mymalloc((integer)tslen) ;
  85.          for (p=tempstore, q=tsp; p<tsend; p++, q++)
  86.             *q = *p ;
  87.          free((char *)tempstore) ;
  88.          tempstore = tsp ;
  89.          tsp = q ;
  90.       }
  91.       tsend = tempstore + tslen ;
  92.    }
  93.    *tsp++ = what ;
  94. }
  95.  
  96. /* Next, a routine that discovers how to do the compression. */
  97.  
  98. #define rsh(a,b) ( ((a)==0) ? ((b)==128) : ( ((a)==255) ? ((b)==127) :\
  99.                                     ((b)==(((a)>>1)|((a)&128))) ))
  100. #define lsh(a,b) ( ((a)==0) ? ((b)==1) : ( ((a)==255) ? ((b)==254) :\
  101.                                     ((b)==((((a)<<1)&255)|((a)&1))) ))
  102. #define DIFFERENT (1)
  103. #define LSHPOSSIB (2)
  104. #define RSHPOSSIB (4)
  105. #define BLKPOSSIB (8)
  106. #define WHTPOSSIB (16)
  107. #define ENDROW (32)
  108. #define NOPOSSIB(n) ((n&(LSHPOSSIB|RSHPOSSIB|BLKPOSSIB|WHTPOSSIB))==0)
  109. #define NOSHIFT(n) ((n&(LSHPOSSIB|RSHPOSSIB))==0)
  110. /*
  111.  *   Our input bytes are packed to the 16-bit word.  On output,
  112.  *   they're packed to bytes. Thus, we may have to skip a byte at
  113.  *   the end of each row.
  114.  */
  115. void dochar(from, width, height)
  116. unsigned char *from ;
  117. short width, height ; /* in bytes */
  118. {
  119.    register int i ;
  120.    register unsigned char *f, *t, *d ;
  121.    register unsigned char *e ;
  122.    register int accum ;
  123.    int j, k, kk ;
  124.    int diffrow ;
  125.    int repeatcount ;
  126.    int lit, pos, cmd = 0 ;
  127.    long lcm ;
  128.    int widthc ;
  129.  
  130.    widthc = width + (width & 1) ; /* halfword correction */
  131.    lcm = -1 ;
  132.    if (widthc > rowlength) {
  133.       if (rowlength > 0)
  134.          free((char *)specdata) ;
  135.       rowlength = widthc + 30 ;
  136.       specdata = (unsigned char *)mymalloc((integer)(rowlength + 15)) ;
  137.    }
  138.    for (i= -15, t=specdata; i<=widthc; i++, t++)
  139.       *t = 0 ;
  140.    repeatcount = 0 ;
  141.    f = specdata + 2 ;
  142.    for (j=0; j<height; j++, f = from, from += widthc) {
  143.       diffrow = 0 ;
  144.       for (i=0, t=from, d=specdata; i<width; i++, f++, t++, d++) {
  145.          if (*f == *t) {
  146.             if (*t == 0)
  147.                *d = WHTPOSSIB ;
  148.             else if (*t == 255)
  149.                *d = BLKPOSSIB ;
  150.             else
  151.                *d = 0 ;
  152.          } else {
  153.             accum = DIFFERENT ;
  154.             if (rsh(*f, *t))
  155.                accum |= RSHPOSSIB ;
  156.             else if (lsh(*f, *t))
  157.                accum |= LSHPOSSIB ;
  158.             if (*t == 0)
  159.                accum |= WHTPOSSIB ;
  160.             else if (*t == 255)
  161.                accum |= BLKPOSSIB ;
  162.             *d = accum ;
  163.             diffrow++ ;
  164.          }
  165.       } /* end 'for i' */
  166.       *d = ENDROW ;
  167.       if (diffrow == 0) {
  168.          repeatcount++ ;
  169.       } else {
  170.          if (repeatcount) {
  171.             while (repeatcount > MAXOUT) {
  172.                addts((unsigned char)(REPX+MAXOUT)) ;
  173.                repeatcount -= MAXOUT ;
  174.             }
  175.             addts((unsigned char)(REPX+repeatcount)) ;
  176.             repeatcount = 0 ;
  177.          }
  178.          pos = 0 ;
  179.          for (i=0, d=specdata, f=t-width; i<width;) {
  180.             if ((*d & DIFFERENT) == 0) {
  181.                i++ ;
  182.                d++ ;
  183.                f++ ;
  184.             } else {
  185.                accum = 0 ;
  186.                if (pos != i)
  187.                   lit = NOSHIFT(*d) ;
  188.                else /* N.B.: 'lit' does not imply literate programming here */
  189.                   lit = NOPOSSIB(*d) ;
  190.                for (e=d; ;e++) {
  191.                   if (NOPOSSIB(*e))
  192.                      lit = 1 ;
  193.                   if ((*e & DIFFERENT) == 0)
  194.                      break ;
  195.                   if ((*e & WHTPOSSIB) &&
  196.                       (e[1] & WHTPOSSIB)) {
  197.                      while (*e & WHTPOSSIB) {
  198.                         e++ ;
  199.                         accum++ ;
  200.                      }
  201.                      cmd = CLRX ;
  202.                      e -= accum ;
  203.                      break ;
  204.                   } else if ((*e & BLKPOSSIB) &&
  205.                       (e[1] & BLKPOSSIB)) {
  206.                      while (*e & BLKPOSSIB) {
  207.                         e++ ;
  208.                         accum++ ;
  209.                      }
  210.                      cmd = SETX ;
  211.                      e -= accum ;
  212.                      break ;
  213.                   }
  214.                } /* end 'for e'; d pts to first bad byte, e to next good one */
  215.                while (i - pos > MAXOUT) {
  216.                   addts((unsigned char)(ADVX+MAXOUT)) ;
  217.                   pos += MAXOUT ;
  218.                }
  219.                if (0 != (k = (e - d))) {
  220.                   if (lit) {
  221.                      if (k > 2) {
  222.                         if (i > pos) {
  223.                            addts((unsigned char)(ADVX + i - pos)) ;
  224.                            pos = i ;
  225.                         }
  226.                         while (k > MAXOUT) {
  227.                            addts((unsigned char)(CHGX + MAXOUT)) ;
  228.                            for (kk=0; kk<MAXOUT; kk++)
  229.                               addts((unsigned char)(*f++)) ;
  230.                            d += MAXOUT ;
  231.                            pos += MAXOUT ;
  232.                            i += MAXOUT ;
  233.                            k -= MAXOUT ;
  234.                         }
  235.                         addtse((unsigned char)(CHGX + k)) ;
  236.                         pos += k ;
  237.                         for (; d<e; d++, i++, f++)
  238.                            addts((unsigned char)(*f)) ;
  239.                      } else {
  240.                         if (k == 1) {
  241.                            if (i == pos+MAXOUT) {
  242.                               addts((unsigned char)(ADVX + MAXOUT)) ;
  243.                               pos = i ;
  244.                            }
  245.                            addtse((unsigned char)(ADVXCHG1 + i - pos)) ;
  246.                            addts((unsigned char)(*f)) ;
  247.                            i++ ;
  248.                            pos = i ;
  249.                            d++ ;
  250.                            f++ ;
  251.                         } else {
  252.                            if (i == pos+MAXOUT) {
  253.                               addts((unsigned char)(ADVX + MAXOUT)) ;
  254.                               pos = i ;
  255.                            }
  256.                            addtse((unsigned char)(ADVXCHG2 + i - pos)) ;
  257.                            addts((unsigned char)(*f)) ;
  258.                            addts((unsigned char)(f[1])) ;
  259.                            i += 2 ;
  260.                            pos = i ;
  261.                            d += 2 ;
  262.                            f += 2 ;
  263.                         }
  264.                      }
  265.                   } else {
  266.                      while (e > d) {
  267.                         if (*d & LSHPOSSIB) {
  268.                            if (i == pos+MAXOUT) {
  269.                               addts((unsigned char)(ADVX + MAXOUT)) ;
  270.                               pos = i ;
  271.                            }
  272.                            addtse((unsigned char)(ADVXLSH + i - pos)) ;
  273.                         } else if (*d & RSHPOSSIB) {
  274.                            if (i == pos+MAXOUT) {
  275.                               addts((unsigned char)(ADVX + MAXOUT)) ;
  276.                               pos = i ;
  277.                            }
  278.                            addtse((unsigned char)(ADVXRSH + i - pos)) ;
  279.                         } else if (*d & WHTPOSSIB) { /* i==pos */
  280.                            addts((unsigned char)(CLRX + 1)) ; lcm = -1 ;
  281.                         } else if (*d & BLKPOSSIB) { /* i==pos */
  282.                            addts((unsigned char)(SETX + 1)) ; lcm = -1 ;
  283.                         } else
  284.                            error("! bug") ; /* why wasn't lit true? */
  285.                         d++ ;
  286.                         f++ ;
  287.                         i++ ;
  288.                         pos = i ;
  289.                      } /* end 'while e>d' */
  290.                   }
  291.                } /* end 'if e>d' */
  292.                if (accum > 0) {
  293.                   if (i > pos) {
  294.                      addts((unsigned char)(ADVX + i - pos)) ;
  295.                      pos = i ;
  296.                   }
  297.                   i += accum ;
  298.                   d += accum ;
  299.                   f += accum ;
  300.                   while (accum > MAXOUT) {
  301.                      addts((unsigned char)(cmd + MAXOUT)) ;
  302.                      accum -= MAXOUT ;
  303.                      pos += MAXOUT ;
  304.                   }
  305.                   lcm = -1 ;
  306.                   addts((unsigned char)(cmd + accum)) ;
  307.                   pos += accum ;
  308.                }
  309.             } /* end 'else DIFFERENT' */
  310.          } /* end 'for i' */
  311.          if (lcm != -1) {
  312.             tempstore[lcm] += MAXOUT ;  /* append END */
  313.             lcm = -1 ;
  314.          } else {
  315.             addts((unsigned char)(END)) ;
  316.          }
  317.       } /* end 'else rows different' */
  318. #ifdef DEBUG
  319.       if (d != specdata + width)
  320.          error("! internal inconsistency in repack") ;
  321. #endif
  322.    } /* end 'for j' */
  323.    if (repeatcount) {
  324.       while (repeatcount > MAXOUT) {
  325.          addts((unsigned char)(REPX+MAXOUT)) ;
  326.          repeatcount -= MAXOUT ;
  327.       }
  328.       addts((unsigned char)(REPX+repeatcount)) ;
  329.       repeatcount = 0 ;
  330.    }
  331. }
  332.  
  333. extern long bytesleft ;
  334. extern quarterword *raster ;
  335. long mbytesleft ;
  336. quarterword *mraster ;
  337.  
  338. char *
  339. makecopy(what, len, p)
  340. register unsigned char *what ;
  341. register long len ;
  342. register unsigned char *p ;
  343. {
  344.    register unsigned char *q ;
  345.  
  346.    if (p == NULL) {
  347.       if (len > MINCHUNK)
  348.          p = (unsigned char *)mymalloc((integer)len) ;
  349.       else {
  350.          if (bytesleft < len) {
  351.             raster = (quarterword *)mymalloc((integer)RASTERCHUNK) ;
  352.             bytesleft = RASTERCHUNK ;
  353.          }
  354.          p = (unsigned char *)raster ;
  355.          bytesleft -= len ;
  356.          raster += len ;
  357.       }
  358.    }
  359.    q = p ;
  360.    while (len > 0) {
  361.       *p++ = *what++ ;
  362.       len -- ;
  363.    }
  364.    return ((char *)q) ;
  365. }
  366.  
  367. /* Now the main routine, which is called when a character is used
  368.    for the very first time. */
  369. void repack(cp)
  370. register chardesctype *cp ;
  371. {
  372.    register long i, j ;
  373.    register unsigned char *p ;
  374.    register int width, height ;
  375.    register int wwidth ;
  376.    char startbytes ;
  377.    int smallchar ;
  378.  
  379.    p = cp->packptr ;
  380.    if (p == NULL)
  381.       error("! no raster?") ;
  382.    tsp = tempstore ;
  383.    tsend = tempstore + tslen ;
  384.    addts(*p) ; /* copy the PK flag byte */
  385.    if (*p & 4) {
  386.       if ((*p & 7) == 7) {
  387.          startbytes = 17 ;
  388.          width = getlong(p + 1) ;
  389.          height = getlong(p + 5) ;
  390.          for (i=0; i<12; i++)
  391.             addts(*++p) ;
  392.       } else {
  393.          startbytes = 9 ;
  394.          width = p[1] * 256 + p[2] ;
  395.          height = p[3] * 256 + p[4] ;
  396.          addts(*++p) ;
  397.          addts(*++p) ;
  398.          addts(*++p) ;
  399.          addts(*++p) ;
  400.       }
  401.    } else {
  402.       startbytes = 5 ;
  403.       width = p[1] ;
  404.       height = p[2] ;
  405.    }
  406.    addts(*++p) ;
  407.    addts(*++p) ;
  408.    addts(*++p) ;
  409.    addts(*++p) ;
  410.    p++ ; /* OK, p now points to beginning of the nibbles */
  411.    addts((unsigned char)0) ;
  412.    addts((unsigned char)0) ;
  413.    addts((unsigned char)0) ;
  414.    addts((unsigned char)0) ; /* leave room to stick in a new length field */
  415.    wwidth = (width + 15) / 16 ;
  416.    i = 2 * height * (long)wwidth ;
  417.    if (i <= 0)
  418.       i = 2 ;
  419.    if ((cp->flags & BIGCHAR) == 0)
  420.       smallchar = 5 ;
  421.    else
  422.       smallchar = 0 ;
  423.    i += smallchar ;
  424.    if (mbytesleft < i) {
  425.       if (mbytesleft >= RASTERCHUNK)
  426.          (void) free((char *) mraster) ;
  427.       if (RASTERCHUNK > i) {
  428.          mraster = (quarterword *)mymalloc((integer)RASTERCHUNK) ;
  429.          mbytesleft = RASTERCHUNK ;
  430.       } else {
  431.          i += i / 4 ;
  432.          mraster = (quarterword *)mymalloc((integer)(i + 3)) ;
  433.          mbytesleft = i ;
  434.       }
  435.    }
  436.    while (i > 0)
  437.       mraster[--i] = 0 ;
  438.    i = startbytes + unpack(p, (halfword *)mraster, (halfword)width,
  439.                 (halfword)height,  *(cp->packptr)) ;
  440.    dochar(mraster, (width + 7) >> 3, height) ;
  441.    if (smallchar) {
  442.       addts((unsigned char)0) ;
  443.       addts((unsigned char)0) ;
  444.       addts((unsigned char)0) ;
  445.       addts((unsigned char)0) ;
  446.       addts((unsigned char)0) ;
  447.    }
  448.    j = tsp - tempstore ;
  449. #ifdef DEBUG
  450.    if (dd(D_COMPRESS))
  451.         (void)fprintf(stderr,"PK %ld bytes, unpacked %ld, compressed %ld\n",
  452.        i-(long)startbytes, (long)((width+7L)/8)*height, j-(long)startbytes-4) ;
  453. #endif /* DEBUG */
  454.    if ( i < j ) {
  455.       if (i > MINCHUNK)
  456.          free(cp->packptr) ;
  457.       cp->packptr =
  458.               (unsigned char *)makecopy(tempstore, j, (unsigned char *)0) ;
  459.    } else
  460.       makecopy(tempstore, j, (unsigned char *)cp->packptr) ;
  461.    putlong((char *)(cp->packptr+startbytes), j-startbytes-4-smallchar) ;
  462.    cp->flags |= REPACKED ;
  463. }
  464.