/*	dispn.c	associated subroutines for matching sequences */
/* 	modified for Tek4014 terminal, not 4027 */

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

#define XTERNAL
#include "upam.gbl"

#define TRUE 1
#define FALSE 0
extern FILE *outfd;

extern int iscore, gscore;
extern char name0[], ltitle[];
/*extern int colflg;*/

int linbw[]={0,4,2,1,3};
int lincol[]={0,1,2,3,4};
char linchar[]={96,'a','b','c','d','e','f','g'};
int *linarr;
int nlinarr=5;

extern char lvstr[];

int plinval[3]={200,100,50};
int dlinval[3]={400,200,100};
int *linval;

extern int dnaseq;

extern int smin0,smin1;
extern long loffset;
extern int revflg;
extern int pmirror;

#define DIAG 1
#define INS0 2
#define INS1 4

long pminx, pmaxx, pminy, pmaxy;
int max_x=1000, max_y=770;

double fxscal, fyscal, fxoff, fyoff;
#define SX(x) (int)((double)(x)*fxscal+fxoff)
#define SY(y) (int)((double)(y)*fyscal+fyoff)

discons(seqc0, seqc1, nc)
     char *seqc0, *seqc1;
     int nc;
{
  long x0, x1, y0, y1;
  int direct, ii;

  x1 = x0 = smin1 + loffset;
  if (!revflg) y1 = y0 = smin0;
  else {
    if (!pmirror) y1 = y0 = smin0 + nc;
    else {y1 = y0 = smin0;
	  x1 = x0 = smin1 + nc;
	}
  }

  direct = DIAG;
	
  move(SX(x0),SY(y0));

  for (ii=0; ii<nc; ii++) {
    if (seqc0[ii]==' ' || seqc1[ii]==' ') continue;
    if (seqc0[ii]!='-'&&seqc1[ii]!='-') {
      if (direct!=DIAG) {
	draw(SX(x1),SY(y1));
	x0 = x1; y0 = y1;
	direct = DIAG;
      }
      if (!revflg) {x1++; y1++;}
      else {
	if (!pmirror) {x1++; y1--;}
	else {x1--; y1++;}
      }
    }
    else if (seqc0[ii]=='-') {
      if (direct!=INS0) {
	draw(SX(x1),SY(y1));
	x0 = x1; y0 = y1;
	direct = INS0;
      }
      if (revflg && pmirror) x1--;
      else x1++;
    }
    else if (seqc1[ii]=='-') {
      if (direct!=INS1) {
	draw(SX(x1),SY(y1));
	x0 = x1; y0 = y1;
	direct = INS1;
      }
      if (revflg && !pmirror) y1--;
      else y1++;
    }
  }

  draw(SX(x1),SY(y1));
  if (y1 > pmaxy) printf("\r* n0 * %3ld %3ld\n",y1,pmaxy);
  if (y1 < pminy) printf("\r* n0 * %3ld %3ld\n",y1,pminy);
  if (x1 > pmaxx) printf("\r* n1 * %3ld %3ld\n",x1,pmaxx);
}

disgraph() {};

aancpy(to,from,count)
	char *to, *from;
	int count;
{
	char *tp;
	tp=to;
	while (count--&& *from>=0) {
		if (*from<nsq) *tp++ = sq[*(from++)];
		else *tp++ = *from++;
		}
	*tp=0;
	}

iidex(str, chr)
	char *str, chr;
{
	int i;
	for (i=0; str[i]; i++) if (str[i]==chr) return i;
	return (-1);
	}

toupper(c)
	char c;
{
	if (c>='a' && c<='z') return (c-'a'+'A');
	return c;
	}

min(arg1, arg2)
	int arg1, arg2;
{
	return (arg1<=arg2) ? arg1 : arg2;
	}

openplt(n0, n1)
	long n0, n1;
{
	char *getenv(), *sptr;

	if (strlen(lvstr)>0) {
		sscanf(lvstr,"%d %d %d",&plinval[0],&plinval[1],&plinval[2]);
		linval = plinval;
		}
	else if ((sptr=getenv("LINEVAL"))!=NULL && strlen(sptr)>0) {
		sscanf(sptr,"%d %d %d",&plinval[0],&plinval[1],&plinval[2]);
		linval = plinval;
		}
	else {
		if (dnaseq==1) linval=dlinval;
		else linval=plinval;
		}

	printf("\033\014\n");
	printf("\035");		/* send ESC FF GS */

/*	if (colflg) linarr=lincol; else linarr = linbw; */
	linarr = linbw;

	pmaxx = n1;
	pmaxy = n0;

	fxscal = (double)(max_x-1)/(double)(n1);
	fyscal = (double)(max_y-1)/(double)(n0);

	if (fxscal > fyscal) fxscal = fyscal; else fyscal = fxscal;

	if (fyscal * n0 < (double)max_y/5.0) 
		fyscal = (double)(max_y-1)/((double)(n0)*5.0);

	fxscal *= 0.9; fxoff = (double)(max_x-1)/11.0;
	fyscal *= 0.9; fyoff = (double)(max_y-1)/11.0;

	linetype(0);
	move(SX(0),SY(0));
	draw(SX(0),SY(n0));
	draw(SX(n1),SY(n0));
	draw(SX(n1),SY(0));
	draw(SX(0),SY(0));
	xaxis(n1);
	yaxis(n0);
	legend();
	}
	
drawdiag(n0,n1)
	long n0, n1;
{
	linetype(0);
	move(SX(0),SY(0));
	draw(SX(n0),SY(n1));
}

int tarr[] = {10,20,50,100,200,500,1000,2000,5000};
int ntarr = sizeof(tarr);
xaxis(n)
     long n;
{
	int i, jm, tick;
	long js;
	char numstr[20];

	tick = 10;

	for (i=0; i<ntarr; i++) if ((jm = n/tarr[i])<21) goto found;
	jm = n/5000l; i=ntarr-1;
found:	js = tarr[i];

	for (i=1; i<=jm; i++) {
		move(SX((long)i*js),SY(0));
		draw(SX((long)i*js),SY(0)-tick);
		}
	sprintf(numstr,"%ld",js);
	move(SX(js)-6*strlen(numstr),SY(0)-tick-18);
	drawstr(numstr);
	sprintf(numstr,"%ld",jm*js);
	move(SX((long)jm*js)-6*strlen(numstr),SY(0)-tick-18);
	drawstr(numstr);

	move(SX(n/2)-6*strlen(ltitle),SY(0)-tick-40);
	drawstr(ltitle);
	}
		
yaxis(n)
     long n;
{
	int i, jm, tick;
	long js;
	char numstr[20];
	
	tick = 10;

	for (i=0; i<ntarr; i++) if ((jm = n/tarr[i])<21) goto found;
	jm = n/5000l; i=ntarr-1;
found:	js = (long)tarr[i];

	for (i=1; i<=jm; i++) {
		move(SX(0),SY((long)i*js));
		draw(SX(0)-tick,SY((long)i*js));
		}
	sprintf(numstr,"%ld",js);
	move(SX(0)-tick-14*strlen(numstr),SY(js)-7);
	drawstr(numstr);
	sprintf(numstr,"%ld",(long)jm*js);
	move(SX(0)-tick-14*strlen(numstr),SY((long)jm*js)-7);
	drawstr(numstr);
	}

legend()
{
	int i, del;
	int ixp, iyp;
	char numstr[10];

	del = 10;
	for (i=0; i<4; i++) {
		linetype(i);
		ixp = max_x-11*del;
		iyp = max_y*(4-i)/5;
		move(ixp,iyp);
		draw(ixp+5*del,iyp);
		move(ixp+6*del,iyp);
		if (i==3) sprintf(numstr,"<%3d",linval[2]);
		else sprintf(numstr,">%3d",linval[i]);
		drawstr(numstr);
		}
	}

linetype(type)
	int type;
{
/*	if (colflg) printf("\r!COL C%1d",type);	*/
	printf("\033%c",linchar[linarr[type]]);
		}

closeplt()
{
	move(0,0);
	putchar('\r');
	}

opnline(x,y,s,percent,nc)
     long x, y;
     int s, nc;
     float percent;
{
  if (s>linval[0]) linetype(0);
  else if (s>linval[1]) linetype(1);
  else if (s>linval[2]) linetype(2);
  else linetype(3);
}

clsline(x,y,s)
     long x, y;
     int s;
{
}

move(x,y)
	unsigned int x, y;
{
	fputc(0x1d,stdout);
	fputc(((y&0x3e0)>>5)+0x20,stdout);
	fputc((y&0x1f)+0x60,stdout);
	fputc(((x&0x3e0)>>5)+0x20,stdout);
	fputc((x&0x1f)+0x40,stdout);
	}

draw(x,y)
	unsigned int x, y;
{
	fputc(((y&0x3e0)>>5)+0x20,stdout);
	fputc((y&0x1f)+0x60,stdout);
	fputc(((x&0x3e0)>>5)+0x20,stdout);
	fputc((x&0x1f)+0x40,stdout);
	}

drawstr(str)
	char *str;
{
 	fputc(0x1f,stdout);
	fputs(str,stdout);
	fputc(0x1d,stdout);
	}

#ifdef VMS
memset(str, c, cnt)
	char *str; int cnt; char c;
{
	while (cnt--) *str++ = c;
	}
#endif