/*	iface.c		Paul Albrecht		Aug 1987
 
	Copyright (C) Paul Albrecht 1988.  All rights reserved.

	Last Update:	13 June 1995 (GBM)
	EMACS_MODES:	tabstop=4

	Interface between higher level plt calls and the lower
	level Unix plot utilities.

	Copyright (C) Paul Albrecht 1988.  All rights reserved.
*/

#include	"plt.h"
#include	"plot.h"
#include	"text.h"

static void label0();
static void label90();

static	Ptype	xax_off, yax_off, xax_roff, yax_roff;


plabel( lbl, x, y, just, angle )
char	*lbl, *just;
Ptype	x, y;
double	angle;
{
	if( lbl[0] == 0 )
		return;

	if( fabs(angle) < 1 )
		label0( lbl, x, y, just );
	else	if( fabs(angle-90) < 1 )
			label90( lbl, x, y, just );
	else	err( NO, "Bad angle %g in for plabel(%s,...)", just, lbl );
}


static void label0( lbl, x, y, just )	/* horizontal label */
char	*lbl, *just;
{
short	len;
Boolean	errflag;

	len = strlen( lbl );
	errflag = NO;

	switch( just[0] ) {
		case 'L':				break;
		case 'R':	x -= len*chwd;		break;
		case 'C':	x -= (3*len-1)*chwd/6;	break;
		default:	errflag = YES;		break;
	}
	switch( just[1] ) {
		case 'B':			break;
		case 'T':	y -= chht;	break;
		case 'C':	y -= chht/3;	break;
		default:	errflag = YES;	break;
	}
	if( errflag ) {
		err( NO, "Unknown justification mode `%s'", just );
	}
	else {	move( x, y );
		label( lbl );
	}
}

static void label90( lbl, x, y, just )	/* vertical label */
char	*lbl, *just;
{
static	char	chr[2] = { 0, 0 };
short	len;
Boolean	errflag;

	len = strlen( lbl );
	errflag = NO;

	switch( just[0] ) {
		case 'L':				break;
		case 'R':	y -= (3*len-1)*chht/3;	break;
		case 'C':	y -= (3*len-1)*chht/6;	break;
		default:	errflag = YES;		break;
	}
	switch( just[1] ) {
		case 'T':			break;	
		case 'B':	x -= chwd;	break;
		case 'C':	x -= chwd/3;	break;
		default:	errflag = YES;	break;
	}
	if( errflag ) {
		err( NO, "Unknown justification mode `%s'", just );
	}
	else {	while( len > 0 ) {
			move( x, y );
			chr[0] = lbl[--len];
			label( chr );
			y += chht;
		}
	}
}

/********* Output the tic labels and axis label for the axes *************/

alabel( what, lbl, x, y )
char	*what, *lbl;
Ptype	x, y;
{
short	n;

	if( strcmp(what,"XN") == 0 ) {
		n = (lbl[0] == '-') ? 2*chwd/3 : 0;
		label0( lbl, x-n, y, "CT" );
		xax_off = chht;
	}
	else	if( strcmp(what,"YN") == 0 ) {
			label0( lbl, x-chwd/3, y, "RC" );
			n = chwd*strlen(lbl)+chwd/3;
			if( n > yax_off )
				yax_off = n;
		}
	else	if( strcmp(what,"XNR") == 0 ) {
			n = (lbl[0] == '-') ? 2*chwd/3 : 0;
			label0( lbl, x+n, y+chht/4, "CB" );
			xax_roff = chht;
		}
	else	if( strcmp(what,"YNR") == 0 ) {
			label0( lbl, x+chwd/3, y, "LC" );
			n = chwd*strlen(lbl)+chwd/3;
			if( n > yax_roff )
				yax_roff = n;
		}
	else	if( strcmp(what,"XT") == 0 ) {
			label0( lbl, x, y-xax_off, "CT" );
		}
	else	if( strcmp(what,"YT") == 0 ) {
			label90( lbl, x-yax_off, y, "CB" );
		}
	else	if( strcmp(what,"XTR") == 0 ) {
			label0( lbl, x, y+xax_roff, "CB" );
		}
	else	if( strcmp(what,"YTR") == 0 ) {
			label90( lbl, x+yax_roff, y, "CT" );
		}
	else	err( YES, "Bad call: alabel(%s,...)", what );
}

/****************** Output log tic labels for the axes ********************/

elabel( what, base, exponent, x, y )
char	*what, *base, *exponent;
Ptype	x, y;
{
Ptype	off;

	off = chwd * (strlen(base) + strlen(exponent));

	if( strcmp(what,"XE") == 0 ) {
		x += chwd/3;
		label0( exponent,  x+off/2, y, "RT" );
		label0( base, x-off/2, y-chht/2, "LT" );
		xax_off = 3*chht/2;
	}
	else 	if( strcmp(what,"YE") == 0 ) {
			x -= chwd/3;
			label0( exponent, x, y+(chht+3)/6, "RB" );
			label0( base, x-off, y-chht/3, "LB" );
			if( off > yax_off )
				yax_off = off;
		}
	else	if( strcmp(what,"XER") == 0 ) {
			x += chwd/3;
			label0( exponent,  x+off/2, y+3*chht/4, "RB" );
			label0( base, x-off/2, y+chht/4, "LB" );
			xax_roff = 3*chht/2;
		}
	else 	if( strcmp(what,"YER") == 0 ) {
			x += chwd/2;
			label0( exponent, x+off, y+(chht+3)/6, "RB" );
			label0( base, x, y-chht/3, "LB" );
			if( off > yax_roff )
				yax_roff = off;
		}
	else	err( YES, "Bad call: elabel( %s,...)", what );
}

/****** Scatterplot characters or symbols with or without error bars *******/

static	char	circle[] = {100, 0, 89, 46, 57, 82, 12, 99, -34, 94, -74, 66,
 -96, 24, -96, -23, -74, -65, -34, -93, 12, -98, 57, -81, 89, -45, 100, 0};
static	char	diamond[] = {0, -15, 13, 0, 0, 15, -10, 0, 0, -12 };
static	char	triangle[] = {0, 1, -1, -1, 1, -1, 0, 1};
static	char	utriangle[] = {0, -1, -1, 1, 1, 1, 0, -1};
static	char	square[] = {-1, -1, -1, 1, 1, 1, 1, -1, -1, -1};

static	struct	syminf {
	char	*pts;	/* x,y points specifying the symbol contour */
	short	npts;
	double	scl;	/* factor for converting *pts to inches */
	}	*sym, syms[] = {
		circle, sizeof(circle), 0.11,
		utriangle, sizeof(utriangle), 12,
		diamond, sizeof(diamond), 1,
		square, sizeof(square), 10,
		triangle, sizeof(triangle), 12 };


scatterplot( plt, x, y, yneg, ypos )
PltPtr	plt;
Ptype	x, y, yneg, ypos;
{
static	char	lbl[2] = {0, 0};
double	xscl, yscl;
Ptype	n, xd, yd, xdmin, xdmax, ydmin, ydmax;

	if( plt->symbol ) {
		sym = &syms[(plt->pc-'0') % (sizeof(syms)/sizeof(syms[0]))];
		xdmin = ydmin =  30000;
		xdmax = ydmax = -30000;
		xscl = xinch*sym->scl/250;
		yscl = yinch*sym->scl/250;

		for( n=0;  n < sym->npts;  n+=2 ) {
			xd = sym->pts[n];
			yd = sym->pts[n+1];
			xd = xscl*xd + (xd < 0 ? -0.5 : 0.5);
			yd = yscl*yd + (yd < 0 ? -0.5 : 0.5);
			if( n == 0 )
				move( x+xd, y+yd );
			else
				cont( x+xd, y+yd );
			if( xd < xdmin ) xdmin = xd;
			if( xd > xdmax ) xdmax = xd;
			if( yd < ydmin ) ydmin = yd;
			if( yd > ydmax ) ydmax = yd;
		}
	}
	else {	if( plt->pc != '.' ) {
			lbl[0] = plt->pc;
			label0( lbl, x, y, "CC" );
		}
		else	point( x, y );

		xdmax = chwd/2;
		xdmin = -xdmax;
		ydmax = chht/3;
		ydmin = -ydmax;
	}

	if( ypos-y >= ydmax ) {
		move( x+xdmin, ypos );
		cont( x+xdmax, ypos );
		move( x, ypos );
		cont( x, y+ydmax );
	}
	if( yneg-y <= ydmin ) {
		move( x+xdmin, yneg );
		cont( x+xdmax, yneg );
		move( x, yneg );
		cont( x, y+ydmin );
	}
}


PolyDef( x, y, what )
Ptype	x, y;
Const	what;
{
static	int	x0, y0;

	switch( what ) {
		case CMOVE:
			move( x, y );
			x0 = x;
			y0 = y;
			break;
		case CCONT:
			if( oldGrayLevel != PS_WHITE )
				cont( x, y );
			break;
		case CBBCONT:
			cont( x, y );
			break;
		case CFILL:
		case CSTROKE:
		case CFILLI:
			if( oldGrayLevel != PS_WHITE ) {
				cont( x, y );
				cont( x0, y0 );
			}
			break;
		case COSTROKE:
			break;
		case CBBFILL:
		case CBBFILLI:
			cont( x, y );
			cont( x0, y0 );
			break;
		default:
			err( YES, "Bad PolyDef() code %d", what );
	}
}

/******* Give the x and y dimension of a string in device units ******/

strsize( str, xsize, ysize, angle )
char	*str;
Ptype	*xsize, *ysize;
double	angle;
{
	if( fabs(angle) < 1  ) {
		*xsize = chwd * strlen(str);
		*ysize = chht;
	}
	else	if( fabs(angle-90) < 1 ) {
		*xsize = chwd;
		*ysize = chht * strlen(str);
	}
	else	err( NO, "Bad angle %lg in for strsize(%s,...)", angle, str );
}