LAa2b.c

#include <stdlib.h>
#include <stdio.h>

#include "DB.h"
#include "align.h"

int main(int argc, char *argv[])
{ char    code, which;
  int64   total;
  int     aread, bread;
  char    orient, chain;
  int     alen, blen;
  int     ab, ae, bb, be;
  int     diffs;
  int     len;
  int     tspace, small;
  uint8  *tbuffer = NULL;
  uint16 *sbuffer = NULL;

  (void) argv;

  //  Process arguments

  if (argc > 1)
    { fprintf(stderr,"Usage: LAa2b <(ascii) >(binary)\n");
      exit (1);
    }

  while (scanf(" %c",&code) == 1)       //  Header lines
    if (code == '@' || code == '+' || code == '%')
      { scanf(" %c %lld",&which,&total);
        fwrite(&code,sizeof(char),1,stdout);
        fwrite(&which,sizeof(char),1,stdout);
        fwrite(&total,sizeof(int64),1,stdout);
        if (code == '@')
          { tbuffer = (uint8 *) malloc(2*total*sizeof(uint16));
            sbuffer = (uint16 *) tbuffer;
          }
      }
    else
      { ungetc(code,stdin);
        break;
      }
  small = 0;
  if (tbuffer != NULL)
    { if (code != 'X')
        { fprintf(stderr,"LAa2b: .las dump has traces but no X-line\n");
          exit (1);
        }
      scanf(" X %d",&tspace);
      small = (tspace <= TRACE_XOVR && tspace != 0);
      fwrite(&code,sizeof(char),1,stdout);
      fwrite(&tspace,sizeof(int),1,stdout);
    }

  while (scanf(" %c",&code) == 1)       //  For each data line do
    { fwrite(&code,sizeof(char),1,stdout);
      switch (code)
      { case 'P':                         //  Alignment pair
          scanf(" %d %d %c %c",&aread,&bread,&orient,&chain);
          fwrite(&aread,sizeof(int),1,stdout);
          fwrite(&bread,sizeof(int),1,stdout);
          fwrite(&orient,sizeof(char),1,stdout);
          fwrite(&chain,sizeof(char),1,stdout);
          break;
        case 'L':                         //  Read lengths
          scanf(" %d %d",&alen,&blen);
          fwrite(&len,sizeof(int),1,stdout);
          fwrite(&blen,sizeof(int),1,stdout);
          break;
        case 'C':                         //  Coordinate intervals
          scanf(" %d %d %d %d",&ab,&ae,&bb,&be);
          fwrite(&ab,sizeof(int),1,stdout);
          fwrite(&ae,sizeof(int),1,stdout);
          fwrite(&bb,sizeof(int),1,stdout);
          fwrite(&be,sizeof(int),1,stdout);
          break;
        case 'D':                         //  Differences
          scanf(" %d",&diffs);
          fwrite(&diffs,sizeof(int),1,stdout);
          break;
        case 'T':                         //  Mask
          if (tbuffer == NULL)
            { fprintf(stderr,"LAa2b: .las dump has traces but no @ T-line\n");
              exit (1);
            }
          scanf(" %d",&len);
          fwrite(&len,sizeof(int),1,stdout);
          len *= 2;
          if (small)
            { for (int i = 0; i < len; i += 2)
                scanf(" %hhd %hhd",tbuffer+i,tbuffer+(i+1));
              fwrite(tbuffer,sizeof(uint8),len,stdout);
            }
          else
            { for (int i = 0; i < len; i += 2)
                scanf(" %hd %hd",sbuffer+i,sbuffer+(i+1));
              fwrite(sbuffer,sizeof(uint16),len,stdout);
            }
      }
    }

  exit (0);
}