The R Project SVN – R – Diff – /branches/djm-source/src/main/plot.c


/*
 *  R : A Computer Language for Statistical Data Analysis
 *  Copyright (C) 1995, 1996  Robert Gentleman and Ross Ihaka
 *  Copyright (C) 1997--2001  Robert Gentleman, Ross Ihaka and the
 *			      R Development Core Team
 *  Copyright (C) 2002--2004  The R Foundation
 *
 *  This program is free software; you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation; either version 2 of the License, or
 *  (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  A copy of the GNU General Public License is available via WWW at
 *  http://www.gnu.org/copyleft/gpl.html.  You can also obtain it by
 *  writing to the Free Software Foundation, Inc., 59 Temple Place,
 *  Suite 330, Boston, MA  02111-1307  USA.
 */
 
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
 
#include <Defn.h>
#include <Rmath.h>
#include <Graphics.h>
#include <Rdevices.h>
#include <Print.h>
 
#ifndef HAVE_HYPOT
# define hypot pythag
#endif
 
 
void NewFrameConfirm(void)
{
    unsigned char buf[16];
    R_ReadConsole("Hit <Return> to see next plot: ", buf, 16, 0);
}
 
	/* Remember: +1 and/or -1 because C arrays are */
	/* zero-based and R-vectors are one-based. */
 
#define checkArity_length					\
    checkArity(op, args);					\
    if(!LENGTH(CAR(args)))					\
	errorcall(call, "argument must have positive length")
 
 
SEXP do_devcontrol(SEXP call, SEXP op, SEXP args, SEXP env)
{
    int listFlag;
 
    checkArity(op, args);
    listFlag = asLogical(CAR(args));
    if(listFlag == NA_LOGICAL) errorcall(call, "invalid argument");
    if(listFlag)
	enableDisplayList(CurrentDevice());
    else
	inhibitDisplayList(CurrentDevice());
    return ScalarLogical(listFlag);
}
 
SEXP do_devcopy(SEXP call, SEXP op, SEXP args, SEXP env)
{
    checkArity_length;
    GEcopyDisplayList(INTEGER(CAR(args))[0] - 1);
    return R_NilValue;
}
 
SEXP do_devcur(SEXP call, SEXP op, SEXP args, SEXP env)
{
    SEXP cd = allocVector(INTSXP, 1);
    checkArity(op, args);
    INTEGER(cd)[0] = curDevice() + 1;
    return cd;
}
 
SEXP do_devnext(SEXP call, SEXP op, SEXP args, SEXP env)
{
    SEXP nd = allocVector(INTSXP, 1);
    checkArity_length;
    INTEGER(nd)[0] = nextDevice(INTEGER(CAR(args))[0] - 1) + 1;
    return nd;
}
 
SEXP do_devprev(SEXP call, SEXP op, SEXP args, SEXP env)
{
    SEXP pd = allocVector(INTSXP, 1);
    checkArity_length;
    INTEGER(pd)[0] = prevDevice(INTEGER(CAR(args))[0] - 1) + 1;
    return pd;
}
 
SEXP do_devset(SEXP call, SEXP op, SEXP args, SEXP env)
{
    int devNum = INTEGER(CAR(args))[0] - 1;
    SEXP sd = allocVector(INTSXP, 1);
    checkArity(op, args);
    INTEGER(sd)[0] = selectDevice(devNum) + 1;
    return sd;
}
 
SEXP do_devoff(SEXP call, SEXP op, SEXP args, SEXP env)
{
    checkArity_length;
    killDevice(INTEGER(CAR(args))[0] - 1);
    return R_NilValue;
}
 
 
/*  P A R A M E T E R	 U T I L I T I E S  */
 
 
/* ProcessInLinePars handles inline par specifications in graphics functions.
 * It does this by calling Specify2() from ./par.c */
 
void ProcessInlinePars(SEXP s, DevDesc *dd, SEXP call)
{
    if (isList(s)) {
	while (s != R_NilValue) {
	    if (isList(CAR(s)))
		ProcessInlinePars(CAR(s), dd, call);
	    else if (TAG(s) != R_NilValue)
	    Specify2(CHAR(PRINTNAME(TAG(s))), CAR(s), dd, call);
	    s = CDR(s);
	}
    }
}
 
/*
 * Extract specified par from list of inline pars
 */
SEXP getInlinePar(SEXP s, char *name)
{
    SEXP result = R_NilValue;
    int found = 0;
    if (isList(s) && !found) {
	while (s != R_NilValue) {
	    if (isList(CAR(s))) {
		result = getInlinePar(CAR(s), name);
		if (result)
		    found = 1;
	    } else
		if (TAG(s) != R_NilValue)
		    if (!strcmp(CHAR(PRINTNAME(TAG(s))), name)) {
			result = CAR(s);
			found = 1;
		    }
	    s = CDR(s);
	}
    }
    return result;
}
 
SEXP FixupPch(SEXP pch, int dflt)
{
    int i, n;
    SEXP ans = R_NilValue;/* -Wall*/
 
    n = length(pch);
    if (n == 0) {
	ans = allocVector(INTSXP, 1);
	INTEGER(ans)[0] = dflt;
    }
    else if (isList(pch)) {
	ans = allocVector(INTSXP, n);
	for (i = 0; pch != R_NilValue;	pch = CDR(pch))
	    INTEGER(ans)[i++] = asInteger(CAR(pch));
    }
    else if (isInteger(pch)) {
	ans = allocVector(INTSXP, n);
	for (i = 0; i < n; i++)
	    INTEGER(ans)[i] = INTEGER(pch)[i];
    }
    else if (isReal(pch)) {
	ans = allocVector(INTSXP, n);
	for (i = 0; i < n; i++)
	    INTEGER(ans)[i] = R_FINITE(REAL(pch)[i]) ?
		REAL(pch)[i] : NA_INTEGER;
    }
    else if (isString(pch)) {
	ans = allocVector(INTSXP, n);
	for (i = 0; i < n; i++)
	    INTEGER(ans)[i] = STRING_ELT(pch, i) != NA_STRING ?
		CHAR(STRING_ELT(pch, i))[0] : NA_INTEGER;
    }
    else if (isLogical(pch)) {/* NA, but not TRUE/FALSE */
	ans = allocVector(INTSXP, n);
	for (i = 0; i < n; i++)
	    if(LOGICAL(pch)[i] == NA_LOGICAL)
		INTEGER(ans)[i] = NA_INTEGER;
	    else error("only NA allowed in logical plotting symbol");
    }
    else error("invalid plotting symbol");
    for (i = 0; i < n; i++) {
	if (INTEGER(ans)[i] < 0 && INTEGER(ans)[i] != NA_INTEGER)
	    INTEGER(ans)[i] = dflt;
    }
    return ans;
}
 
SEXP FixupLty(SEXP lty, int dflt)
{
    int i, n;
    SEXP ans;
    n = length(lty);
    if (n == 0) {
	ans = allocVector(INTSXP, 1);
	INTEGER(ans)[0] = dflt;
    }
    else {
	ans = allocVector(INTSXP, n);
	for (i = 0; i < n; i++)
	    INTEGER(ans)[i] = LTYpar(lty, i);
    }
    return ans;
}
 
SEXP FixupLwd(SEXP lwd, double dflt)
{
    int i, n;
    double w;
    SEXP ans = NULL;
 
    n = length(lwd);
    if (n == 0) {
	ans = allocVector(REALSXP, 1);
	REAL(ans)[0] = dflt;
    }
    else {
	PROTECT(lwd = coerceVector(lwd, REALSXP));
	n = length(lwd);
	ans = allocVector(REALSXP, n);
	for (i = 0; i < n; i++) {
	    w = REAL(lwd)[i];
	    if (w < 0) w = NA_REAL;
	    REAL(ans)[i] = w;
 
	}
	UNPROTECT(1);
    }
    return ans;
}
 
SEXP FixupFont(SEXP font, int dflt)
{
    int i, k, n;
    SEXP ans = R_NilValue;/* -Wall*/
    n = length(font);
    if (n == 0) {
	ans = allocVector(INTSXP, 1);
	INTEGER(ans)[0] = dflt;
    }
    else if (isInteger(font) || isLogical(font)) {
	ans = allocVector(INTSXP, n);
	for (i = 0; i < n; i++) {
	    k = INTEGER(font)[i];
#ifndef Win32
	    if (k < 1 || k > 4) k = NA_INTEGER;
#else
	    if (k < 1 || k > 32) k = NA_INTEGER;
#endif
	    INTEGER(ans)[i] = k;
	}
    }
    else if (isReal(font)) {
	ans = allocVector(INTSXP, n);
	for (i = 0; i < n; i++) {
	    k = REAL(font)[i];
#ifndef Win32
	    if (k < 1 || k > 4) k = NA_INTEGER;
#else
	    if (k < 1 || k > 32) k = NA_INTEGER;
#endif
	    INTEGER(ans)[i] = k;
	}
    }
    else error("invalid font specification");
    return ans;
}
 
SEXP FixupCol(SEXP col, unsigned int dflt)
{
    int i, n;
    SEXP ans;
    n = length(col);
    if (n == 0) {
	ans = allocVector(INTSXP, 1);
	INTEGER(ans)[0] = dflt;
    }
    else {
	ans = allocVector(INTSXP, n);
	if (isList(col))
	    for (i = 0; i < n; i++) {
		INTEGER(ans)[i] = RGBpar(CAR(col), 0);
		col = CDR(col);
	    }
	else
	    for (i = 0; i < n; i++)
		INTEGER(ans)[i] = RGBpar(col, i);
    }
    return ans;
}
 
SEXP FixupCex(SEXP cex, double dflt)
{
    SEXP ans;
    int i, n;
    n = length(cex);
    if (n == 0) {
	ans = allocVector(REALSXP, 1);
	if (R_FINITE(dflt) && dflt > 0)
	    REAL(ans)[0] = dflt;
	else
	    REAL(ans)[0] = NA_REAL;
    }
    else {
	double c;
	ans = allocVector(REALSXP, n);
	if (isReal(cex))
	    for (i = 0; i < n; i++) {
		c = REAL(cex)[i];
		if (R_FINITE(c) && c > 0)
		    REAL(ans)[i] = c;
		else
		    REAL(ans)[i] = NA_REAL;
	    }
	else if (isInteger(cex) || isLogical(cex))
 
	    for (i = 0; i < n; i++) {
		c = INTEGER(cex)[i];
		if (c == NA_INTEGER || c <= 0)
		    c = NA_REAL;
		REAL(ans)[i] = c;
	    }
    }
    return ans;
}
 
SEXP FixupVFont(SEXP vfont) {
    SEXP ans = R_NilValue;
    if (!isNull(vfont)) {
	SEXP vf;
	int typeface, fontindex;
	int minindex, maxindex=0;/* -Wall*/
	int i;
	PROTECT(vf = coerceVector(vfont, INTSXP));
	if (length(vf) != 2)
	    error("Invalid vfont value");
	typeface = INTEGER(vf)[0];
	if (typeface < 0 || typeface > 7)
	    error("Invalid vfont value [typeface]");
	/* For each of the typefaces {0..7}, there are several fontindices
	   available; how many depends on the typeface.
	   The possible combinations are "given" in ./g_fontdb.c
	   and also listed in help(Hershey).
	 */
	minindex = 1;
	switch (typeface) {
	case 0: /* serif */
	    maxindex = 7;	    break;
	case 1: /* sans serif */
	case 6: /* serif symbol */
	    maxindex = 4;	    break;
	case 2: /* script */
	    maxindex = 3;	    break;
	case 3: /* gothic english */
	case 4: /* gothic german */
	case 5: /* gothic italian */
	    maxindex = 1;	    break;
	case 7: /* sans serif symbol */
	    maxindex = 2;
	}
	fontindex = INTEGER(vf)[1];
	if (fontindex < minindex || fontindex > maxindex)
	    error("Invalid vfont value [fontindex]");
	ans = allocVector(INTSXP, 2);
	for (i=0; i<2; i++)
	    INTEGER(ans)[i] = INTEGER(vf)[i];
	UNPROTECT(1);
    }
    return ans;
}
 
/* GetTextArg() : extract from call and possibly set text arguments
 *  ("label", col=, cex=, font=, vfont=)
 *
 * Main purpose: Treat things like  title(main = list("This Title", font= 4))
 *
 * Called from	do_title()  [only, currently]
 */
static void
GetTextArg(SEXP call, SEXP spec, SEXP *ptxt,
	   int *pcol, double *pcex, int *pfont, SEXP *pvfont)
{
    int i, n, col, font, colspecd;
    double cex;
    SEXP txt, vfont, nms;
 
    txt	  = R_NilValue;
    vfont = R_NilValue;
    cex	  = NA_REAL;
    col	  = NA_INTEGER;
    colspecd = 0;
    font  = NA_INTEGER;
    PROTECT(txt);
 
    switch (TYPEOF(spec)) {
    case LANGSXP:
    case SYMSXP:
	UNPROTECT(1);
	PROTECT(txt = coerceVector(spec, EXPRSXP));
	break;
    case VECSXP:
	if (length(spec) == 0) {
	    *ptxt = R_NilValue;
	}
	else {
	    nms = getAttrib(spec, R_NamesSymbol);
	    if (nms==R_NilValue){ /* PR#1939 */
	       txt = VECTOR_ELT(spec, 0);
	       if (TYPEOF(txt) == LANGSXP ||TYPEOF(txt) == SYMSXP ) {
		    UNPROTECT(1);
		    PROTECT(txt = coerceVector(txt, EXPRSXP));
	       }
	       else if (!isExpression(txt)) {
		    UNPROTECT(1);
		    PROTECT(txt = coerceVector(txt, STRSXP));
	       }
	    } else {
	       n = length(nms);
	       for (i = 0; i < n; i++) {
		if (!strcmp(CHAR(STRING_ELT(nms, i)), "cex")) {
		    cex = asReal(VECTOR_ELT(spec, i));
		}
		else if (!strcmp(CHAR(STRING_ELT(nms, i)), "col")) {
		    SEXP colsxp = VECTOR_ELT(spec, i);
		    if (!isNAcol(colsxp, 0, LENGTH(colsxp))) {
			col = asInteger(FixupCol(colsxp, NA_INTEGER));
			colspecd = 1;
		    }
		}
		else if (!strcmp(CHAR(STRING_ELT(nms, i)), "font")) {
		    font = asInteger(FixupFont(VECTOR_ELT(spec, i), NA_INTEGER));
		}
		else if (!strcmp(CHAR(STRING_ELT(nms, i)), "vfont")) {
		    vfont = FixupVFont(VECTOR_ELT(spec, i));
		}
		else if (!strcmp(CHAR(STRING_ELT(nms, i)), "")) {
		    txt = VECTOR_ELT(spec, i);
		    if (TYPEOF(txt) == LANGSXP || TYPEOF(txt)==SYMSXP) {
			UNPROTECT(1);
			PROTECT(txt = coerceVector(txt, EXPRSXP));
		    }
		    else if (!isExpression(txt)) {
			UNPROTECT(1);
			PROTECT(txt = coerceVector(txt, STRSXP));
		    }
		}
		else errorcall(call, "invalid graphics parameter");
	       }
	    }
	}
	break;
    case STRSXP:
    case EXPRSXP:
	txt = spec;
	break;
    default:
	txt = coerceVector(spec, STRSXP);
	break;
    }
    UNPROTECT(1);
    if (txt != R_NilValue) {
	*ptxt = txt;
	if (R_FINITE(cex))	 *pcex	 = cex;
	if (colspecd)	         *pcol	 = col;
	if (font != NA_INTEGER)	 *pfont	 = font;
	if (vfont != R_NilValue) *pvfont = vfont;
    }
}/* GetTextArg */
 
 
    /* GRAPHICS FUNCTION ENTRY POINTS */
 
 
SEXP do_plot_new(SEXP call, SEXP op, SEXP args, SEXP env)
{
    /* plot.new() - create a new plot "frame" */
 
    DevDesc *dd;
 
    checkArity(op, args);
 
    dd = GNewPlot(GRecording(call));
 
    Rf_dpptr(dd)->xlog = Rf_gpptr(dd)->xlog = FALSE;
    Rf_dpptr(dd)->ylog = Rf_gpptr(dd)->ylog = FALSE;
 
    GScale(0.0, 1.0, 1, dd);
    GScale(0.0, 1.0, 2, dd);
    GMapWin2Fig(dd);
    GSetState(1, dd);
 
    if (GRecording(call))
	recordGraphicOperation(op, args, dd);
    return R_NilValue;
}
 
 
/*
 *  SYNOPSIS
 *
 *	plot.window(xlim, ylim, log="", asp=NA)
 *
 *  DESCRIPTION
 *
 *	This function sets up the world coordinates for a graphics
 *	window.	 Note that if asp is a finite positive value then
 *	the window is set up so that one data unit in the y direction
 *	is equal in length to one data unit in the x direction divided
 *	by asp.
 *
 *	The special case asp == 1 produces plots where distances
 *	between points are represented accurately on screen.
 *
 *  NOTE
 *
 *	The use of asp can have weird effects when axis is an
 *	interpreted function.  It has to be internal so that the
 *	full computation is captured in the display list.
 */
 
SEXP do_plot_window(SEXP call, SEXP op, SEXP args, SEXP env)
{
    SEXP xlim, ylim, logarg;
    double asp, xmin, xmax, ymin, ymax;
    Rboolean logscale;
    char *p;
    SEXP originalArgs = args;
    DevDesc *dd = CurrentDevice();
 
    if (length(args) < 3)
	errorcall(call, "at least 3 arguments required");
 
    xlim = CAR(args);
    if (!isNumeric(xlim) || LENGTH(xlim) != 2)
	errorcall(call, "invalid xlim");
    args = CDR(args);
 
    ylim = CAR(args);
    if (!isNumeric(ylim) || LENGTH(ylim) != 2)
	errorcall(call, "invalid ylim");
    args = CDR(args);
 
    logscale = FALSE;
    logarg = CAR(args);
    if (!isString(logarg))
	errorcall(call, "\"log=\" specification must be character");
    p = CHAR(STRING_ELT(logarg, 0));
    while (*p) {
	switch (*p) {
	case 'x':
	    Rf_dpptr(dd)->xlog = Rf_gpptr(dd)->xlog = logscale = TRUE;
	    break;
	case 'y':
	    Rf_dpptr(dd)->ylog = Rf_gpptr(dd)->ylog = logscale = TRUE;
	    break;
	default:
	    errorcall(call,"invalid \"log=%s\" specification",p);
	}
	p++;
    }
    args = CDR(args);
 
    asp = (logscale) ? NA_REAL : asReal(CAR(args));;
    args = CDR(args);
 
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
 
    if (isInteger(xlim)) {
	if (INTEGER(xlim)[0] == NA_INTEGER || INTEGER(xlim)[1] == NA_INTEGER)
	    errorcall(call, "NAs not allowed in xlim");
	xmin = INTEGER(xlim)[0];
	xmax = INTEGER(xlim)[1];
    }
    else {
	if (!R_FINITE(REAL(xlim)[0]) || !R_FINITE(REAL(xlim)[1]))
	    errorcall(call, "need finite xlim values");
	xmin = REAL(xlim)[0];
	xmax = REAL(xlim)[1];
    }
    if (isInteger(ylim)) {
	if (INTEGER(ylim)[0] == NA_INTEGER || INTEGER(ylim)[1] == NA_INTEGER)
	    errorcall(call, "NAs not allowed in ylim");
	ymin = INTEGER(ylim)[0];
	ymax = INTEGER(ylim)[1];
    }
    else {
	if (!R_FINITE(REAL(ylim)[0]) || !R_FINITE(REAL(ylim)[1]))
	    errorcall(call, "need finite ylim values");
	ymin = REAL(ylim)[0];
	ymax = REAL(ylim)[1];
    }
    if ((Rf_dpptr(dd)->xlog && (xmin < 0 || xmax < 0)) ||
       (Rf_dpptr(dd)->ylog && (ymin < 0 || ymax < 0)))
	    errorcall(call, "Logarithmic axis must have positive limits");
 
    if (R_FINITE(asp) && asp > 0) {
	double pin1, pin2, scale, xdelta, ydelta, xscale, yscale, xadd, yadd;
	pin1 = GConvertXUnits(1.0, NPC, INCHES, dd);
	pin2 = GConvertYUnits(1.0, NPC, INCHES, dd);
	xdelta = fabs(xmax - xmin) / asp;
	ydelta = fabs(ymax - ymin);
	if(xdelta == 0.0 && ydelta == 0.0) {
	    /* We really do mean zero: small non-zero values work
	       Mimic the behaviour of GScale for the x axis. */
	    xadd = yadd = ((xmin == 0.0) ? 1 : 0.4) * asp;
	    xadd *= asp;
	} else {
	    xscale = pin1 / xdelta;
	    yscale = pin2 / ydelta;
	    scale = (xscale < yscale) ? xscale : yscale;
	    xadd = .5 * (pin1 / scale - xdelta) * asp;
	    yadd = .5 * (pin2 / scale - ydelta);
	}
	GScale(xmin - xadd, xmax + xadd, 1, dd);
	GScale(ymin - yadd, ymax + yadd, 2, dd);
    }
    else {
	GScale(xmin, xmax, 1, dd);
	GScale(ymin, ymax, 2, dd);
    }
    GMapWin2Fig(dd);
    GRestorePars(dd);
    /* NOTE: the operation is only recorded if there was no "error" */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}
 
void GetAxisLimits(double left, double right, double *low, double *high)
{
/*	Called from do_axis()	such as
 *	GetAxisLimits(Rf_gpptr(dd)->usr[0], Rf_gpptr(dd)->usr[1], &low, &high)
 *
 *	Computes  *low < left, right < *high  (even if left=right)
 */
    double eps;
    if (left > right) {/* swap */
	eps = left; left = right; right = eps;
    }
    eps = right - left;
    if (eps == 0.)
	eps = 0.5 * FLT_EPSILON;
    else
	eps *= FLT_EPSILON;
    *low = left - eps;
    *high = right + eps;
}
 
 
/* axis(side, at, labels, ...) */
 
SEXP labelformat(SEXP labels)
{
    /* format(labels): i.e. from numbers to strings */
    SEXP ans = R_NilValue;/* -Wall*/
    int i, n, w, d, e, wi, di, ei;
    char *strp;
    n = length(labels);
    R_print.digits = 7;/* maximally 7 digits -- ``burnt in'';
			  S-PLUS <= 5.x has about 6
			  (but really uses single precision..) */
    switch(TYPEOF(labels)) {
    case LGLSXP:
	PROTECT(ans = allocVector(STRSXP, n));
	for (i = 0; i < n; i++) {
	    strp = EncodeLogical(LOGICAL(labels)[i], 0);
	    SET_STRING_ELT(ans, i, mkChar(strp));
	}
	UNPROTECT(1);
	break;
    case INTSXP:
	PROTECT(ans = allocVector(STRSXP, n));
	for (i = 0; i < n; i++) {
	    strp = EncodeInteger(INTEGER(labels)[i], 0);
	    SET_STRING_ELT(ans, i, mkChar(strp));
	}
	UNPROTECT(1);
	break;
    case REALSXP:
	formatReal(REAL(labels), n, &w, &d, &e, 0);
	PROTECT(ans = allocVector(STRSXP, n));
	for (i = 0; i < n; i++) {
	    strp = EncodeReal(REAL(labels)[i], 0, d, e);
	    SET_STRING_ELT(ans, i, mkChar(strp));
	}
	UNPROTECT(1);
	break;
    case CPLXSXP:
	formatComplex(COMPLEX(labels), n, &w, &d, &e, &wi, &di, &ei, 0);
	PROTECT(ans = allocVector(STRSXP, n));
	for (i = 0; i < n; i++) {
	    strp = EncodeComplex(COMPLEX(labels)[i], 0, d, e, 0, di, ei);
	    SET_STRING_ELT(ans, i, mkChar(strp));
	}
	UNPROTECT(1);
	break;
    case STRSXP:
	PROTECT(ans = allocVector(STRSXP, n));
	for (i = 0; i < n; i++) {
	    SET_STRING_ELT(ans, i, STRING_ELT(labels, i));
	}
	UNPROTECT(1);
	break;
    default:
	error("invalid type for axis labels");
    }
    return ans;
}
 
SEXP CreateAtVector(double *axp, double *usr, int nint, Rboolean logflag)
{
/*	Create an  'at = ...' vector for  axis(.) / do_axis,
 *	i.e., the vector of tick mark locations,
 *	when none has been specified (= default).
 *
 *	axp[0:2] = (x1, x2, nInt), where x1..x2 are the extreme tick marks
 *                 {unless in log case, where nint \in {1,2,3 ; -1,-2,....}
 *                  and the `nint' argument is used.
 
 *	The resulting REAL vector must have length >= 1, ideally >= 2
 */
    SEXP at = R_NilValue;/* -Wall*/
    double umin, umax, dn, rng, small;
    int i, n, ne;
    if (!logflag || axp[2] < 0) { /* --- linear axis --- Only use axp[] arg. */
	n = fabs(axp[2]) + 0.25;/* >= 0 */
	dn = imax2(1, n);
	rng = axp[1] - axp[0];
	small = fabs(rng)/(100.*dn);
	at = allocVector(REALSXP, n + 1);
	for (i = 0; i <= n; i++) {
	    REAL(at)[i] = axp[0] + (i / dn) * rng;
	    if (fabs(REAL(at)[i]) < small)
		REAL(at)[i] = 0;
	}
    }
    else { /* ------ log axis ----- */
	n = (axp[2] + 0.5);
	/* {xy}axp[2] for 'log': GLpretty() [../graphics.c] sets
	   n < 0: very small scale ==> linear axis, above, or
	   n = 1,2,3.  see switch() below */
	umin = usr[0];
	umax = usr[1];
	/* Debugging: When does the following happen... ? */
	if (umin > umax)
	    warning("CreateAtVector \"log\"(from axis()): "
		    "usr[0] = %g > %g = usr[1] !", umin, umax);
	dn = axp[0];
	if (dn < 1e-300)
	    warning("CreateAtVector \"log\"(from axis()): axp[0] = %g !", dn);
 
	/* You get the 3 cases below by
	 *  for (y in 1e-5*c(1,2,8))  plot(y, log = "y")
	 */
	switch(n) {
	case 1: /* large range:	1	 * 10^k */
	    i = floor(log10(axp[1])) - ceil(log10(axp[0])) + 0.25;
	    ne = i / nint + 1;
	    if (ne < 1)
		error("log - axis(), 'at' creation, _LARGE_ range: "
		      "ne = %d <= 0 !!\n"
		      "\t axp[0:1]=(%g,%g) ==> i = %d;	nint = %d",
		      ne, axp[0],axp[1], i, nint);
	    rng = pow(10., (double)ne);/* >= 10 */
	    n = 0;
	    while (dn < umax) {
		n++;
		dn *= rng;
	    }
	    if (!n)
		error("log - axis(), 'at' creation, _LARGE_ range: "
		      "illegal {xy}axp or par; nint=%d\n"
		      "	 axp[0:1]=(%g,%g), usr[0:1]=(%g,%g); i=%d, ni=%d",
		      nint, axp[0],axp[1], umin,umax, i,ne);
	    at = allocVector(REALSXP, n);
	    dn = axp[0];
	    n = 0;
	    while (dn < umax) {
		REAL(at)[n++] = dn;
		dn *= rng;
	    }
	    break;
 
	case 2: /* medium range:  1, 5	  * 10^k */
	    n = 0;
	    if (0.5 * dn >= umin) n++;
	    for (;;) {
		if (dn > umax) break;		n++;
		if (5 * dn > umax) break;	n++;
		dn *= 10;
	    }
	    if (!n)
		error("log - axis(), 'at' creation, _MEDIUM_ range: "
		      "illegal {xy}axp or par;\n"
		      "	 axp[0]= %g, usr[0:1]=(%g,%g)",
		      axp[0], umin,umax);
 
	    at = allocVector(REALSXP, n);
	    dn = axp[0];
	    n = 0;
	    if (0.5 * dn >= umin) REAL(at)[n++] = 0.5 * dn;
	    for (;;) {
		if (dn > umax) break;		REAL(at)[n++] = dn;
		if (5 * dn > umax) break;	REAL(at)[n++] = 5 * dn;
		dn *= 10;
	    }
	    break;
 
	case 3: /* small range:	 1,2,5,10 * 10^k */
	    n = 0;
	    if (0.2 * dn >= umin) n++;
	    if (0.5 * dn >= umin) n++;
	    for (;;) {
		if (dn > umax) break;		n++;
		if (2 * dn > umax) break;	n++;
		if (5 * dn > umax) break;	n++;
		dn *= 10;
	    }
	    if (!n)
		error("log - axis(), 'at' creation, _SMALL_ range: "
		      "illegal {xy}axp or par;\n"
		      "	 axp[0]= %g, usr[0:1]=(%g,%g)",
		      axp[0], umin,umax);
	    at = allocVector(REALSXP, n);
	    dn = axp[0];
	    n = 0;
	    if (0.2 * dn >= umin) REAL(at)[n++] = 0.2 * dn;
	    if (0.5 * dn >= umin) REAL(at)[n++] = 0.5 * dn;
	    for (;;) {
		if (dn > umax) break;		REAL(at)[n++] = dn;
		if (2 * dn > umax) break;	REAL(at)[n++] = 2 * dn;
		if (5 * dn > umax) break;	REAL(at)[n++] = 5 * dn;
		dn *= 10;
	    }
	    break;
	default:
	    error("log - axis(), 'at' creation: ILLEGAL {xy}axp[3] = %g",
		  axp[2]);
	}
    }
    return at;
}
 
SEXP do_axis(SEXP call, SEXP op, SEXP args, SEXP env)
{
    /* axis(side, at, labels, tick, line, pos,
     *	    outer, font, vfont, lty, lwd, col, ...) */
 
    SEXP at, lab, vfont;
    int col, font, lty;
    int i, n, nint = 0, ntmp, side, *ind, outer, lineoff = 0;
    int istart, iend, incr;
    Rboolean dolabels, doticks, logflag = FALSE;
    Rboolean create_at, vectorFonts = FALSE;
    double x, y, temp, tnew, tlast;
    double axp[3], usr[2];
    double gap, labw, low, high, line, pos, lwd;
    double axis_base, axis_tick, axis_lab, axis_low, axis_high;
 
    SEXP originalArgs = args, label;
    DevDesc *dd = CurrentDevice();
 
    /* Arity Check */
    /* This is a builtin function, so it should always have */
    /* the correct arity, but it doesn't hurt to be defensive. */
 
    if (length(args) < 9)
	errorcall(call, "too few arguments");
    GCheckState(dd);
 
    /* Required argument: "side" */
    /* Which side of the plot the axis is to appear on. */
    /* side = 1 | 2 | 3 | 4. */
 
    side = asInteger(CAR(args));
    if (side < 1 || side > 4)
	errorcall(call, "invalid axis number %d", side);
    args = CDR(args);
 
    /* Required argument: "at" */
    /* This gives the tick-label locations. */
    /* Note that these are coerced to the correct type below. */
 
    at = CAR(args); args = CDR(args);
 
    /* Required argument: "labels" */
    /* Labels can be a logical, indicating whether or not */
    /* to label the axis; or it can be a vector of character */
    /* strings or expressions which give the labels explicitly. */
    /* The expressions are used to set mathematical labelling. */
 
    dolabels = TRUE;
    if (isLogical(CAR(args)) && length(CAR(args)) > 0) {
	i = asLogical(CAR(args));
	if (i == 0 || i == NA_LOGICAL)
	    dolabels = FALSE;
	PROTECT(lab = R_NilValue);
    }
    else if (isExpression(CAR(args))) {
	PROTECT(lab = CAR(args));
    }
    else {
	PROTECT(lab = coerceVector(CAR(args), STRSXP));
    }
    args = CDR(args);
 
    /* Required argument: "tick" */
    /* This indicates whether or not ticks and the axis line */
    /* should be plotted: TRUE => show, FALSE => don't show. */
 
    doticks = asLogical(CAR(args));
    doticks = (doticks == NA_LOGICAL) ? TRUE : (Rboolean) doticks;
    args = CDR(args);
 
    /* Optional argument: "line" */
 
    /* Specifies an offset outward from the plot for the axis.
     * The values in the par value "mgp" are interpreted
     * relative to this value. */
    line = asReal(CAR(args));
    if (!R_FINITE(line)) {
	/* Except that here mgp values are not relative to themselves */
	line = Rf_gpptr(dd)->mgp[2];
	lineoff = line;
    }
    args = CDR(args);
 
    /* Optional argument: "pos" */
    /* Specifies a user coordinate at which the axis should be drawn. */
    /* This overrides the value of "line".  Again the "mgp" par values */
    /* are interpreted relative to this value. */
 
    pos = asReal(CAR(args));
    if (!R_FINITE(pos)) pos = NA_REAL; else lineoff = 0;
    args = CDR(args);
 
    /* Optional argument: "outer" */
    /* Should the axis be drawn in the outer margin. */
    /* This only affects the computation of axis_base. */
 
    outer = asLogical(CAR(args));
    if (outer == NA_LOGICAL || outer == 0)
	outer = NPC;
    else
	outer = NIC;
    args = CDR(args);
 
    /* Optional argument: "font" */
    font = asInteger(FixupFont(CAR(args), NA_INTEGER));
    args = CDR(args);
 
    /* Optional argument: "vfont" */
    /* Allows Hershey vector fonts to be used */
    PROTECT(vfont = FixupVFont(CAR(args)));
    if (!isNull(vfont))
	vectorFonts = TRUE;
    args = CDR(args);
 
    /* Optional argument: "lty" */
    lty = asInteger(FixupLty(CAR(args), NA_INTEGER));
    args = CDR(args);
 
    /* Optional argument: "lwd" */
    lwd = asReal(FixupLwd(CAR(args), Rf_gpptr(dd)->lwd));
    args = CDR(args);
 
    /* Optional argument: "col" */
    col = asInteger(FixupCol(CAR(args), Rf_gpptr(dd)->fg));
    args = CDR(args);
 
    /* Retrieve relevant "par" values. */
 
    switch(side) {
    case 1:
    case 3:
	axp[0] = Rf_dpptr(dd)->xaxp[0];
	axp[1] = Rf_dpptr(dd)->xaxp[1];
	axp[2] = Rf_dpptr(dd)->xaxp[2];
	usr[0] = Rf_dpptr(dd)->usr[0];
	usr[1] = Rf_dpptr(dd)->usr[1];
	logflag = Rf_dpptr(dd)->xlog;
	nint = Rf_dpptr(dd)->lab[0];
	break;
    case 2:
    case 4:
	axp[0] = Rf_dpptr(dd)->yaxp[0];
	axp[1] = Rf_dpptr(dd)->yaxp[1];
	axp[2] = Rf_dpptr(dd)->yaxp[2];
	usr[0] = Rf_dpptr(dd)->usr[2];
	usr[1] = Rf_dpptr(dd)->usr[3];
	logflag = Rf_dpptr(dd)->ylog;
	nint = Rf_dpptr(dd)->lab[1];
	break;
    }
 
    /* Determine the tickmark positions.  Note that these may fall */
    /* outside the plot window. We will clip them in the code below. */
 
    create_at = (length(at) == 0);
    if (create_at) {
	PROTECT(at = CreateAtVector(axp, usr, nint, logflag));
    }
    else {
	if (isReal(at)) PROTECT(at = duplicate(at));
	else PROTECT(at = coerceVector(at, REALSXP));
    }
    n = length(at);
 
    /* Check/setup the tick labels.  This can mean using user-specified */
    /* labels, or encoding the "at" positions as strings. */
 
    if (dolabels) {
	if (length(lab) == 0)
	    lab = labelformat(at);
	else if (!isExpression(lab))
	    lab = labelformat(lab);
	if (length(at) != length(lab))
	    errorcall(call, "location and label lengths differ, %d != %d",
		      length(at), length(lab));
    }
    PROTECT(lab);
 
    /* Check there are no NA, Inf or -Inf values for tick positions. */
    /* The code here is long-winded.  Couldn't we just inline things */
    /* below.  Hmmm - we need the min and max of the finite values ... */
 
    ind = (int *) R_alloc(n, sizeof(int));
    for(i = 0; i < n; i++) ind[i] = i;
    rsort_with_index(REAL(at), ind, n);
    ntmp = 0;
    for(i = 0; i < n; i++) {
	if(R_FINITE(REAL(at)[i])) ntmp = i+1;
    }
    n = ntmp;
    if (n == 0)
	errorcall(call, "no locations are finite");
 
    /* Ok, all systems are "GO".  Let's get to it.
     * First we process all the remaining inline par values */
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
 
    /* At this point we know the value of "xaxt" and "yaxt",
     * so we test to see whether the relevant one is "n".
     * If it is, we just bail out at this point. */
 
    if (((side == 1 || side == 3) && Rf_gpptr(dd)->xaxt == 'n') ||
	((side == 2 || side == 4) && Rf_gpptr(dd)->yaxt == 'n')) {
	GRestorePars(dd);
	UNPROTECT(4);
	return R_NilValue;
    }
 
 
    /* no! we do allow an `lty' argument -- will not be used often though
     *	Rf_gpptr(dd)->lty = LTY_SOLID; */
    Rf_gpptr(dd)->lty = lty;
    Rf_gpptr(dd)->lwd = lwd;
 
    /* Override par("xpd") and force clipping to figure region.
     * NOTE: don't override to _reduce_ clipping region */
 
    Rf_gpptr(dd)->xpd = 2;
 
    Rf_gpptr(dd)->adj = 0.5;
    Rf_gpptr(dd)->font = (font == NA_INTEGER)? Rf_gpptr(dd)->fontaxis : font;
    Rf_gpptr(dd)->cex = Rf_gpptr(dd)->cexbase * Rf_gpptr(dd)->cexaxis;
    /* no!   col = Rf_gpptr(dd)->col; */
 
    /* Draw the axis */
    GMode(1, dd);
    switch (side) {
    case 1: /*--- x-axis -- horizontal --- */
    case 3:
	GetAxisLimits(Rf_gpptr(dd)->usr[0], Rf_gpptr(dd)->usr[1], &low, &high);
	axis_low  = GConvertX(fmax2(low, REAL(at)[0]), USER, NFC, dd);
	axis_high = GConvertX(fmin2(high, REAL(at)[n-1]), USER, NFC, dd);
	if (side == 1) {
	    if (R_FINITE(pos))
		axis_base = GConvertY(pos, USER, NFC, dd);
	    else
		axis_base = GConvertY(0.0, outer, NFC, dd)
		    - GConvertYUnits(line, LINES, NFC, dd);
	    if (R_FINITE(Rf_gpptr(dd)->tck)) {
		double len, xu, yu;
		if(Rf_gpptr(dd)->tck > 0.5)
		    len = GConvertYUnits(Rf_gpptr(dd)->tck, NPC, NFC, dd);
		else {
		    xu = GConvertXUnits(Rf_gpptr(dd)->tck, NPC, INCHES, dd);
		    yu = GConvertYUnits(Rf_gpptr(dd)->tck, NPC, INCHES, dd);
		    xu = (fabs(xu) < fabs(yu)) ? xu : yu;
		    len = GConvertYUnits(xu, INCHES, NFC, dd);
		}
		axis_tick = axis_base + len;
 
	    } else
		axis_tick = axis_base +
			GConvertYUnits(Rf_gpptr(dd)->tcl, LINES, NFC, dd);
	}
	else {
	    if (R_FINITE(pos))
		axis_base = GConvertY(pos, USER, NFC, dd);
	    else
		axis_base =  GConvertY(1.0, outer, NFC, dd)
		    + GConvertYUnits(line, LINES, NFC, dd);
	    if (R_FINITE(Rf_gpptr(dd)->tck)) {
		double len, xu, yu;
		if(Rf_gpptr(dd)->tck > 0.5)
		    len = GConvertYUnits(Rf_gpptr(dd)->tck, NPC, NFC, dd);
		else {
		    xu = GConvertXUnits(Rf_gpptr(dd)->tck, NPC, INCHES, dd);
		    yu = GConvertYUnits(Rf_gpptr(dd)->tck, NPC, INCHES, dd);
		    xu = (fabs(xu) < fabs(yu)) ? xu : yu;
		    len = GConvertYUnits(xu, INCHES, NFC, dd);
		}
		axis_tick = axis_base - len;
	    } else
		axis_tick = axis_base -
		    GConvertYUnits(Rf_gpptr(dd)->tcl, LINES, NFC, dd);
	}
	if (doticks) {
	    Rf_gpptr(dd)->col = col;/*was fg */
	    GLine(axis_low, axis_base, axis_high, axis_base, NFC, dd);
	    for (i = 0; i < n; i++) {
		x = REAL(at)[i];
		if (low <= x && x <= high) {
		    x = GConvertX(x, USER, NFC, dd);
		    GLine(x, axis_base, x, axis_tick, NFC, dd);
		}
	    }
	}
	/* Tickmark labels. */
	Rf_gpptr(dd)->col = Rf_gpptr(dd)->colaxis;
	gap = GStrWidth("m", NFC, dd);	/* FIXUP x/y distance */
	tlast = -1.0;
	if (Rf_gpptr(dd)->las == 2 || Rf_gpptr(dd)->las == 3) {
	    Rf_gpptr(dd)->adj = (side == 1) ? 1 : 0;
	}
	else Rf_gpptr(dd)->adj = 0.5;
	if (side == 1) {
	    axis_lab = - axis_base
		+ GConvertYUnits(Rf_gpptr(dd)->mgp[1]-lineoff, LINES, NFC, dd)
		+ GConvertY(0.0, NPC, NFC, dd);
	}
	else { /* side == 3 */
	    axis_lab = axis_base
		+ GConvertYUnits(Rf_gpptr(dd)->mgp[1]-lineoff, LINES, NFC, dd)
		- GConvertY(1.0, NPC, NFC, dd);
	}
	axis_lab = GConvertYUnits(axis_lab, NFC, LINES, dd);
 
	/* The order of processing is important here. */
	/* We must ensure that the labels are drawn left-to-right. */
	/* The logic here is getting way too convoluted. */
	/* This needs a serious rewrite. */
 
	if (Rf_gpptr(dd)->usr[0] > Rf_gpptr(dd)->usr[1]) {
	    istart = n - 1;
	    iend = -1;
	    incr = -1;
	}
	else {
	    istart = 0;
	    iend = n;
	    incr = 1;
	}
	for (i = istart; i != iend; i += incr) {
	    x = REAL(at)[i];
	    if (!R_FINITE(x)) continue;
	    temp = GConvertX(x, USER, NFC, dd);
	    if (dolabels) {
		/* Clip tick labels to user coordinates. */
		if (x > low && x < high) {
		    if (isExpression(lab)) {
			GMMathText(VECTOR_ELT(lab, ind[i]), side,
				   axis_lab, 0, x, Rf_gpptr(dd)->las, dd);
		    }
		    else {
			label = STRING_ELT(lab, ind[i]);
			if(label != NA_STRING) {
			    labw = GStrWidth(CHAR(label), NFC, dd);
			    tnew = temp - 0.5 * labw;
			    /* Check room for perpendicular labels. */
			    if (Rf_gpptr(dd)->las == 2 ||
				Rf_gpptr(dd)->las == 3 ||
				tnew - tlast >= gap) {
				GMtext(CHAR(label), side, axis_lab, 0, x,
				       Rf_gpptr(dd)->las, dd);
				tlast = temp + 0.5 *labw;
			    }
			}
		    }
		}
	    }
	}
	break;
 
    case 2: /*--- y-axis -- vertical --- */
    case 4:
	GetAxisLimits(Rf_gpptr(dd)->usr[2], Rf_gpptr(dd)->usr[3], &low, &high);
	axis_low = GConvertY(fmax2(low, REAL(at)[0]), USER, NFC, dd);
	axis_high = GConvertY(fmin2(high, REAL(at)[n-1]), USER, NFC, dd);
	if (side == 2) {
	    if (R_FINITE(pos))
		axis_base = GConvertX(pos, USER, NFC, dd);
	    else
		axis_base =  GConvertX(0.0, outer, NFC, dd)
		    - GConvertXUnits(line, LINES, NFC, dd);
	    if (R_FINITE(Rf_gpptr(dd)->tck)) {
		double len, xu, yu;
		if(Rf_gpptr(dd)->tck > 0.5)
		    len = GConvertXUnits(Rf_gpptr(dd)->tck, NPC, NFC, dd);
		else {
		    xu = GConvertXUnits(Rf_gpptr(dd)->tck, NPC, INCHES, dd);
		    yu = GConvertYUnits(Rf_gpptr(dd)->tck, NPC, INCHES, dd);
		    xu = (fabs(xu) < fabs(yu)) ? xu : yu;
		    len = GConvertXUnits(xu, INCHES, NFC, dd);
		}
		axis_tick = axis_base + len;
	    } else
		axis_tick = axis_base +
		    GConvertXUnits(Rf_gpptr(dd)->tcl, LINES, NFC, dd);
	}
	else {
	    if (R_FINITE(pos))
		axis_base = GConvertX(pos, USER, NFC, dd);
	    else
		axis_base =  GConvertX(1.0, outer, NFC, dd)
		    + GConvertXUnits(line, LINES, NFC, dd);
	    if (R_FINITE(Rf_gpptr(dd)->tck)) {
		double len, xu, yu;
		if(Rf_gpptr(dd)->tck > 0.5)
		    len = GConvertXUnits(Rf_gpptr(dd)->tck, NPC, NFC, dd);
		else {
		    xu = GConvertXUnits(Rf_gpptr(dd)->tck, NPC, INCHES, dd);
		    yu = GConvertYUnits(Rf_gpptr(dd)->tck, NPC, INCHES, dd);
		    xu = (fabs(xu) < fabs(yu)) ? xu : yu;
		    len = GConvertXUnits(xu, INCHES, NFC, dd);
		}
		axis_tick = axis_base - len;
	    } else
		axis_tick = axis_base -
		    GConvertXUnits(Rf_gpptr(dd)->tcl, LINES, NFC, dd);
	}
	if (doticks) {
	    Rf_gpptr(dd)->col = col;/*was fg */
	    GLine(axis_base, axis_low, axis_base, axis_high, NFC, dd);
	    for (i = 0; i < n; i++) {
		y = REAL(at)[i];
		if (low <= y && y <= high) {
		    y = GConvertY(y, USER, NFC, dd);
		    GLine(axis_base, y, axis_tick, y, NFC, dd);
		}
	    }
	}
	/* Tickmark labels. */
	Rf_gpptr(dd)->col = Rf_gpptr(dd)->colaxis;
	gap = GStrWidth("m", INCHES, dd);
	gap = GConvertYUnits(gap, INCHES, NFC, dd);
	tlast = -1.0;
	if (Rf_gpptr(dd)->las == 1 || Rf_gpptr(dd)->las == 2) {
	    Rf_gpptr(dd)->adj = (side == 2) ? 1 : 0;
	}
	else Rf_gpptr(dd)->adj = 0.5;
	if (side == 2) {
	    axis_lab = - axis_base
		+ GConvertXUnits(Rf_gpptr(dd)->mgp[1]-lineoff, LINES, NFC, dd)
		+ GConvertX(0.0, NPC, NFC, dd);
	}
	else { /* side == 4 */
	    axis_lab = axis_base
		+ GConvertXUnits(Rf_gpptr(dd)->mgp[1]-lineoff, LINES, NFC, dd)
		- GConvertX(1.0, NPC, NFC, dd);
	}
	axis_lab = GConvertXUnits(axis_lab, NFC, LINES, dd);
 
	/* The order of processing is important here. */
	/* We must ensure that the labels are drawn left-to-right. */
	/* The logic here is getting way too convoluted. */
	/* This needs a serious rewrite. */
 
	if (Rf_gpptr(dd)->usr[2] > Rf_gpptr(dd)->usr[3]) {
	    istart = n - 1;
	    iend = -1;
	    incr = -1;
	}
	else {
	    istart = 0;
	    iend = n;
	    incr = 1;
	}
	for (i = istart; i != iend; i += incr) {
	    y = REAL(at)[i];
	    if (!R_FINITE(y)) continue;
	    temp = GConvertY(y, USER, NFC, dd);
	    if (dolabels) {
		/* Clip tick labels to user coordinates. */
		if (y > low && y < high) {
		    if (isExpression(lab)) {
			GMMathText(VECTOR_ELT(lab, ind[i]), side,
				   axis_lab, 0, y, Rf_gpptr(dd)->las, dd);
		    }
		    else {
			label = STRING_ELT(lab, ind[i]);
			if(label != NA_STRING) {
			    labw = GStrWidth(CHAR(label), INCHES, dd);
			    labw = GConvertYUnits(labw, INCHES, NFC, dd);
			    tnew = temp - 0.5 * labw;
			    /* Check room for perpendicular labels. */
			    if (Rf_gpptr(dd)->las == 1 ||
				Rf_gpptr(dd)->las == 2 ||
				tnew - tlast >= gap) {
				GMtext(CHAR(label), side, axis_lab, 0, y,
				       Rf_gpptr(dd)->las, dd);
				tlast = temp + 0.5 *labw;
			    }
			}
		    }
		}
	    }
	}
	break;
    } /* end  switch(side, ..) */
    UNPROTECT(4); /* lab, vfont, at, lab again */
    GMode(0, dd);
    GRestorePars(dd);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}/* do_axis */
 
 
SEXP do_plot_xy(SEXP call, SEXP op, SEXP args, SEXP env)
{
/*	plot.xy(xy, type, pch, lty, col, cex, ...)
 
 *	plot points or lines of various types
 */
    SEXP sxy, sx, sy, pch, cex, col, bg, lty;
    double *x, *y, xold, yold, xx, yy, thiscex;
    int i, n, npch, ncex, ncol, nbg, nlty, type=0, start=0, thispch, thiscol;
 
    SEXP originalArgs = args;
    DevDesc *dd = CurrentDevice();
 
    /* Basic Checks */
    GCheckState(dd);
    if (length(args) < 6)
	errorcall(call, "too few arguments");
 
    /* Required Arguments */
#define PLOT_XY_DEALING(subname)				\
    sx = R_NilValue;		/* -Wall */			\
    sy = R_NilValue;		/* -Wall */			\
    sxy = CAR(args);						\
    if (isNewList(sxy) && length(sxy) >= 2) {			\
	internalTypeCheck(call, sx = VECTOR_ELT(sxy, 0), REALSXP);	\
	internalTypeCheck(call, sy = VECTOR_ELT(sxy, 1), REALSXP);	\
    }								\
    else if (isList(sxy) && length(sxy) >= 2) {			\
	internalTypeCheck(call, sx = CAR(sxy), REALSXP);	\
	internalTypeCheck(call, sy = CADR(sxy), REALSXP);	\
    }								\
    else							\
	errorcall(call, "invalid plotting structure");		\
    if (LENGTH(sx) != LENGTH(sy))				\
	error("x and y lengths differ in " subname "().");	\
    n = LENGTH(sx);						\
    args = CDR(args)
 
    PLOT_XY_DEALING("plot.xy");
 
    if (isNull(CAR(args))) type = 'p';
    else {
	if (isString(CAR(args)) && LENGTH(CAR(args)) == 1 &&
	    LENGTH(pch = STRING_ELT(CAR(args), 0)) >= 1) {
	    if(LENGTH(pch) > 1)
		warningcall(call, "plot type '%s' truncated to first character",
			    CHAR(pch));
	    type = CHAR(pch)[0];
	}
	else errorcall(call, "invalid plot type");
    }
    args = CDR(args);
 
    PROTECT(pch = FixupPch(CAR(args), Rf_gpptr(dd)->pch));	args = CDR(args);
    npch = length(pch);
 
    PROTECT(lty = FixupLty(CAR(args), Rf_gpptr(dd)->lty));	args = CDR(args);
    nlty = length(lty);
 
    /* Default col was NA_INTEGER (0x80000000) which was interpreted
       as zero (black) or "don't draw" depending on line/rect/circle
       situation. Now we set the default to zero and don't plot at all
       if col==NA.
 
       FIXME: bg needs similar change, but that requires changes to
       the specific drivers. */
 
    PROTECT(col = FixupCol(CAR(args), 0)); args = CDR(args);
    ncol = LENGTH(col);
 
    PROTECT(bg = FixupCol(CAR(args), NA_INTEGER));	args = CDR(args);
    nbg = LENGTH(bg);
 
    PROTECT(cex = FixupCex(CAR(args), 1.0));	args = CDR(args);
    ncex = LENGTH(cex);
 
    /* Miscellaneous Graphical Parameters -- e.g., lwd */
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
 
    x = REAL(sx);
    y = REAL(sy);
 
    if (nlty && INTEGER(lty)[0] != NA_INTEGER)
	Rf_gpptr(dd)->lty = INTEGER(lty)[0];
 
    GMode(1, dd);
    /* removed by paul 26/5/99 because all clipping now happens in graphics.c
     * GClip(dd);
     */
 
    switch(type) {
    case 'l':
    case 'o':
	/* lines and overplotted lines and points */
	Rf_gpptr(dd)->col = INTEGER(col)[0];
	xold = NA_REAL;
	yold = NA_REAL;
	for (i = 0; i < n; i++) {
	    xx = x[i];
	    yy = y[i];
	    /* do the conversion now to check for non-finite */
	    GConvert(&xx, &yy, USER, DEVICE, dd);
	    if ((R_FINITE(xx) && R_FINITE(yy)) &&
		!(R_FINITE(xold) && R_FINITE(yold)))
		start = i;
	    else if ((R_FINITE(xold) && R_FINITE(yold)) &&
		     !(R_FINITE(xx) && R_FINITE(yy))) {
		if (i-start > 1)
		    GPolyline(i-start, x+start, y+start,
			      USER, dd);
	    }
	    else if ((R_FINITE(xold) && R_FINITE(yold)) &&
		     (i == n-1))
		GPolyline(n-start, x+start, y+start, USER, dd);
	    xold = xx;
	    yold = yy;
	}
	break;
 
    case 'b':
    case 'c': /* broken lines (with points in between if 'b') */
    {
	double d, f;
	d = GConvertYUnits(0.5, CHARS, INCHES, dd);
	Rf_gpptr(dd)->col = INTEGER(col)[0];
	xold = NA_REAL;
	yold = NA_REAL;
	for (i = 0; i < n; i++) {
	    xx = x[i];
	    yy = y[i];
	    GConvert(&xx, &yy, USER, INCHES, dd);
	    if (R_FINITE(xold) && R_FINITE(yold) &&
		R_FINITE(xx) && R_FINITE(yy)) {
		if ((f = d/hypot(xx-xold, yy-yold)) < 0.5) {
		    GLine(xold + f * (xx - xold),
			  yold + f * (yy - yold),
			  xx + f * (xold - xx),
			  yy + f * (yold - yy),
			  INCHES, dd);
		}
	    }
	    xold = xx;
	    yold = yy;
	}
    }
    break;
 
    case 's': /* step function	I */
    {
	double xtemp[3], ytemp[3];
	Rf_gpptr(dd)->col = INTEGER(col)[0];
	xold = x[0];
	yold = y[0];
	GConvert(&xold, &yold, USER, DEVICE, dd);
	for (i = 1; i < n; i++) {
	    xx = x[i];
	    yy = y[i];
	    GConvert(&xx, &yy, USER, DEVICE, dd);
	    if (R_FINITE(xold) && R_FINITE(yold) &&
		R_FINITE(xx) && R_FINITE(yy)) {
		xtemp[0] = xold; ytemp[0] = yold;
		xtemp[1] = xx;	ytemp[1] = yold;
		xtemp[2] = xx;	ytemp[2] = yy;
		GPolyline(3, xtemp, ytemp, DEVICE, dd);
	    }
	    xold = xx;
	    yold = yy;
	}
    }
    break;
 
    case 'S': /* step function	II */
    {
	double xtemp[3], ytemp[3];
	Rf_gpptr(dd)->col = INTEGER(col)[0];
	xold = x[0];
	yold = y[0];
	GConvert(&xold, &yold, USER, DEVICE, dd);
	for (i = 1; i < n; i++) {
	    xx = x[i];
	    yy = y[i];
	    GConvert(&xx, &yy, USER, DEVICE, dd);
	    if (R_FINITE(xold) && R_FINITE(yold) &&
		R_FINITE(xx) && R_FINITE(yy)) {
		xtemp[0] = xold; ytemp[0] = yold;
		xtemp[1] = xold; ytemp[1] = yy;
		xtemp[2] = xx; ytemp[2] = yy;
		GPolyline(3, xtemp, ytemp, DEVICE, dd);
	    }
	    xold = xx;
	    yold = yy;
	}
    }
    break;
 
    case 'h': /* h[istogram] (bar plot) */
	if (Rf_gpptr(dd)->ylog)
	    yold = Rf_gpptr(dd)->usr[2];/* DBL_MIN fails.. why ???? */
	else
	    yold = 0.0;
	yold = GConvertY(yold, USER, DEVICE, dd);
	for (i = 0; i < n; i++) {
	    xx = x[i];
	    yy = y[i];
	    GConvert(&xx, &yy, USER, DEVICE, dd);
	    if (R_FINITE(xx) && R_FINITE(yy)
		&& !R_TRANSPARENT(thiscol = INTEGER(col)[i % ncol])) {
		Rf_gpptr(dd)->col = thiscol;
		GLine(xx, yold, xx, yy, DEVICE, dd);
	    }
	}
	break;
 
    case 'p':
    case 'n': /* nothing here */
	break;
 
    default:/* OTHERWISE */
	errorcall(call, "invalid plot type '%c'", type);
 
    } /* End {switch(type)} */
 
    if (type == 'p' || type == 'b' || type == 'o') {
	for (i = 0; i < n; i++) {
	    xx = x[i];
	    yy = y[i];
	    GConvert(&xx, &yy, USER, DEVICE, dd);
	    if (R_FINITE(xx) && R_FINITE(yy)) {
		if (R_FINITE(thiscex = REAL(cex)[i % ncex])
		    && (thispch = INTEGER(pch)[i % npch]) != NA_INTEGER
		    && !R_TRANSPARENT(thiscol = INTEGER(col)[i % ncol]))
		{
		    Rf_gpptr(dd)->cex = thiscex * Rf_gpptr(dd)->cexbase;
		    Rf_gpptr(dd)->col = thiscol;
		    Rf_gpptr(dd)->bg = INTEGER(bg)[i % nbg];
		    GSymbol(xx, yy, DEVICE, thispch, dd);
		}
	    }
	}
    }
    GMode(0, dd);
    GRestorePars(dd);
    UNPROTECT(5);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}/* do_plot_xy */
 
/* Checks for ... , x0, y0, x1, y1 ... */
 
static void xypoints(SEXP call, SEXP args, int *n)
{
    int k=0;/* -Wall */
 
    if (!isNumeric(CAR(args)) || (k = LENGTH(CAR(args))) <= 0)
	errorcall(call, "first argument invalid");
    SETCAR(args, coerceVector(CAR(args), REALSXP));
    *n = k;
    args = CDR(args);
 
    if (!isNumeric(CAR(args)) || (k = LENGTH(CAR(args))) <= 0)
	errorcall(call, "second argument invalid");
    SETCAR(args, coerceVector(CAR(args), REALSXP));
    if (k > *n) *n = k;
    args = CDR(args);
 
    if (!isNumeric(CAR(args)) || (k = LENGTH(CAR(args))) <= 0)
	errorcall(call, "third argument invalid");
    SETCAR(args, coerceVector(CAR(args), REALSXP));
    if (k > *n) *n = k;
    args = CDR(args);
 
    if (!isNumeric(CAR(args)) || (k = LENGTH(CAR(args))) <= 0)
	errorcall(call, "fourth argument invalid");
    SETCAR(args, coerceVector(CAR(args), REALSXP));
    if (k > *n) *n = k;
    args = CDR(args);
}
 
 
SEXP do_segments(SEXP call, SEXP op, SEXP args, SEXP env)
{
    /* segments(x0, y0, x1, y1, col, lty, lwd, ...) */
    SEXP sx0, sx1, sy0, sy1, col, lty, lwd;
    double *x0, *x1, *y0, *y1;
    double xx[2], yy[2];
    int nx0, nx1, ny0, ny1, i, n, ncol, nlty, nlwd;
    SEXP originalArgs = args;
    DevDesc *dd = CurrentDevice();
 
    if (length(args) < 4) errorcall(call, "too few arguments");
    GCheckState(dd);
 
    xypoints(call, args, &n);
 
    sx0 = CAR(args); nx0 = length(sx0); args = CDR(args);
    sy0 = CAR(args); ny0 = length(sy0); args = CDR(args);
    sx1 = CAR(args); nx1 = length(sx1); args = CDR(args);
    sy1 = CAR(args); ny1 = length(sy1); args = CDR(args);
 
    PROTECT(col = FixupCol(CAR(args), NA_INTEGER));
    ncol = LENGTH(col); args = CDR(args);
 
    PROTECT(lty = FixupLty(CAR(args), Rf_gpptr(dd)->lty));
    nlty = length(lty); args = CDR(args);
 
    PROTECT(lwd = FixupLwd(CAR(args), Rf_gpptr(dd)->lwd));
    nlwd = length(lwd); args = CDR(args);
 
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
 
    x0 = REAL(sx0);
    y0 = REAL(sy0);
    x1 = REAL(sx1);
    y1 = REAL(sy1);
 
    GMode(1, dd);
    for (i = 0; i < n; i++) {
	xx[0] = x0[i%nx0];
	yy[0] = y0[i%ny0];
	xx[1] = x1[i%nx1];
	yy[1] = y1[i%ny1];
	GConvert(xx, yy, USER, DEVICE, dd);
	GConvert(xx+1, yy+1, USER, DEVICE, dd);
	if (R_FINITE(xx[0]) && R_FINITE(yy[0]) &&
	    R_FINITE(xx[1]) && R_FINITE(yy[1]))
	{
	    Rf_gpptr(dd)->col = INTEGER(col)[i % ncol];
	    /* NA color should be ok */
	    Rf_gpptr(dd)->lty = INTEGER(lty)[i % nlty];
	    Rf_gpptr(dd)->lwd = REAL(lwd)[i % nlwd];
	    GLine(xx[0], yy[0], xx[1], yy[1], DEVICE, dd);
	}
    }
    GMode(0, dd);
    GRestorePars(dd);
 
    UNPROTECT(3);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}
 
 
SEXP do_rect(SEXP call, SEXP op, SEXP args, SEXP env)
{
    /* rect(xl, yb, xr, yt, col, border, lty, lwd, xpd, ...) */
    SEXP sxl, sxr, syb, syt, sxpd, col, lty, lwd, border;
    double *xl, *xr, *yb, *yt, x0, y0, x1, y1;
    int i, n, nxl, nxr, nyb, nyt, ncol, nlty, nlwd, nborder, xpd;
    SEXP originalArgs = args;
    DevDesc *dd = CurrentDevice();
 
    if (length(args) < 4) errorcall(call, "too few arguments");
    GCheckState(dd);
 
    xypoints(call, args, &n);
    sxl = CAR(args); nxl = length(sxl); args = CDR(args);/* x_left */
    syb = CAR(args); nyb = length(syb); args = CDR(args);/* y_bottom */
    sxr = CAR(args); nxr = length(sxr); args = CDR(args);/* x_right */
    syt = CAR(args); nyt = length(syt); args = CDR(args);/* y_top */
 
    PROTECT(col = FixupCol(CAR(args), NA_INTEGER));
    ncol = LENGTH(col);
    args = CDR(args);
 
    PROTECT(border =  FixupCol(CAR(args), Rf_gpptr(dd)->fg));
    nborder = LENGTH(border);
    args = CDR(args);
 
    PROTECT(lty = FixupLty(CAR(args), Rf_gpptr(dd)->lty));
    nlty = length(lty);
    args = CDR(args);
 
    PROTECT(lwd = FixupLwd(CAR(args), Rf_gpptr(dd)->lwd));
    nlwd = length(lwd);
    args = CDR(args);
 
    sxpd = CAR(args);
    if (sxpd != R_NilValue)
	xpd = asInteger(sxpd);
    else
	xpd = Rf_gpptr(dd)->xpd;
    args = CDR(args);
 
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
 
    if (xpd == NA_INTEGER)
	Rf_gpptr(dd)->xpd = 2;
    else
	Rf_gpptr(dd)->xpd = xpd;
 
    xl = REAL(sxl);
    xr = REAL(sxr);
    yb = REAL(syb);
    yt = REAL(syt);
 
    GMode(1, dd);
    for (i = 0; i < n; i++) {
	if (nlty && INTEGER(lty)[i % nlty] != NA_INTEGER)
	    Rf_gpptr(dd)->lty = INTEGER(lty)[i % nlty];
	else
	    Rf_gpptr(dd)->lty = Rf_dpptr(dd)->lty;
	if (nlwd && REAL(lwd)[i % nlwd] != NA_REAL)
	    Rf_gpptr(dd)->lwd = REAL(lwd)[i % nlwd];
	else
	    Rf_gpptr(dd)->lwd = Rf_dpptr(dd)->lwd;
	x0 = xl[i%nxl];
	y0 = yb[i%nyb];
	x1 = xr[i%nxr];
	y1 = yt[i%nyt];
	GConvert(&x0, &y0, USER, DEVICE, dd);
	GConvert(&x1, &y1, USER, DEVICE, dd);
	if (R_FINITE(x0) && R_FINITE(y0) && R_FINITE(x1) && R_FINITE(y1))
	    GRect(x0, y0, x1, y1, DEVICE, INTEGER(col)[i % ncol],
		  INTEGER(border)[i % nborder], dd);
    }
    GMode(0, dd);
 
    GRestorePars(dd);
    UNPROTECT(4);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}
 
 
SEXP do_arrows(SEXP call, SEXP op, SEXP args, SEXP env)
{
    /* arrows(x0, y0, x1, y1, length, angle, code, col, lty, lwd, xpd) */
    SEXP sx0, sx1, sy0, sy1, sxpd, col, rawcol, lty, lwd;
    double *x0, *x1, *y0, *y1;
    double xx0, yy0, xx1, yy1;
    double hlength, angle;
    int code;
    int nx0, nx1, ny0, ny1, i, n, ncol, nlty, nlwd, xpd;
    SEXP originalArgs = args;
    DevDesc *dd = CurrentDevice();
 
    if (length(args) < 4) errorcall(call, "too few arguments");
    GCheckState(dd);
 
    xypoints(call, args, &n);
 
    sx0 = CAR(args); nx0 = length(sx0); args = CDR(args);
    sy0 = CAR(args); ny0 = length(sy0); args = CDR(args);
    sx1 = CAR(args); nx1 = length(sx1); args = CDR(args);
    sy1 = CAR(args); ny1 = length(sy1); args = CDR(args);
 
    hlength = asReal(CAR(args));
    if (!R_FINITE(hlength) || hlength < 0)
	errorcall(call, "invalid head length");
    args = CDR(args);
 
    angle = asReal(CAR(args));
    if (!R_FINITE(angle))
	errorcall(call, "invalid head angle");
    args = CDR(args);
 
    code = asInteger(CAR(args));
    if (code == NA_INTEGER || code < 0 || code > 3)
	errorcall(call, "invalid arrow head specification");
    args = CDR(args);
 
    /*
     * Need raw colours to be able to check for NAs
     * FixupCol converts NAs to fully transparent
     */
    rawcol = CAR(args);
    PROTECT(col = FixupCol(rawcol, NA_INTEGER));
    ncol = LENGTH(col);
    args = CDR(args);
 
    PROTECT(lty = FixupLty(CAR(args), Rf_gpptr(dd)->lty));
    nlty = length(lty);
    args = CDR(args);
 
    PROTECT(lwd = CAR(args));/*need  FixupLwd ?*/
    nlwd = length(lwd);
    if (nlwd == 0)
	errorcall(call, "'lwd' must be numeric of length >=1");
    args = CDR(args);
 
    sxpd = CAR(args);
    if (sxpd != R_NilValue)
	xpd = asInteger(sxpd);
    else
	xpd = Rf_gpptr(dd)->xpd;
    args = CDR(args);
 
    GSavePars(dd);
 
    if (xpd == NA_INTEGER)
	Rf_gpptr(dd)->xpd = 2;
    else
	Rf_gpptr(dd)->xpd = xpd;
 
    x0 = REAL(sx0);
    y0 = REAL(sy0);
    x1 = REAL(sx1);
    y1 = REAL(sy1);
 
    GMode(1, dd);
    for (i = 0; i < n; i++) {
	xx0 = x0[i%nx0];
	yy0 = y0[i%ny0];
	xx1 = x1[i%nx1];
	yy1 = y1[i%ny1];
	GConvert(&xx0, &yy0, USER, DEVICE, dd);
	GConvert(&xx1, &yy1, USER, DEVICE, dd);
	if (R_FINITE(xx0) && R_FINITE(yy0) && R_FINITE(xx1) && R_FINITE(yy1)) {
 
	    if (isNAcol(rawcol, i, ncol))
		Rf_gpptr(dd)->col = Rf_dpptr(dd)->col;
	    else
		Rf_gpptr(dd)->col = INTEGER(col)[i % ncol];
	    if (nlty == 0 || INTEGER(lty)[i % nlty] == NA_INTEGER)
		Rf_gpptr(dd)->lty = Rf_dpptr(dd)->lty;
	    else
		Rf_gpptr(dd)->lty = INTEGER(lty)[i % nlty];
	    Rf_gpptr(dd)->lwd = REAL(lwd)[i % nlwd];
	    GArrow(xx0, yy0, xx1, yy1, DEVICE,
		   hlength, angle, code, dd);
	}
    }
    GMode(0, dd);
    GRestorePars(dd);
 
    UNPROTECT(3);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}
 
 
static void drawPolygon(int n, double *x, double *y,
			int lty, int fill, int border, DevDesc *dd)
{
    if (lty == NA_INTEGER)
	Rf_gpptr(dd)->lty = Rf_dpptr(dd)->lty;
    else
	Rf_gpptr(dd)->lty = lty;
    GPolygon(n, x, y, USER, fill, border, dd);
}
 
SEXP do_polygon(SEXP call, SEXP op, SEXP args, SEXP env)
{
    /* polygon(x, y, col, border, lty, xpd, ...) */
    SEXP sx, sy, col, border, lty, sxpd;
    int nx;
    int ncol, nborder, nlty, xpd, i, start=0;
    int num = 0;
    double *x, *y, xx, yy, xold, yold;
 
    SEXP originalArgs = args;
    DevDesc *dd = CurrentDevice();
 
    GCheckState(dd);
 
    if (length(args) < 2) errorcall(call, "too few arguments");
    /* (x,y) is checked in R via xy.coords() ; no need here : */
    sx = SETCAR(args, coerceVector(CAR(args), REALSXP));  args = CDR(args);
    sy = SETCAR(args, coerceVector(CAR(args), REALSXP));  args = CDR(args);
    nx = LENGTH(sx);
 
    PROTECT(col = FixupCol(CAR(args), NA_INTEGER));	args = CDR(args);
    ncol = LENGTH(col);
 
    PROTECT(border = FixupCol(CAR(args), Rf_gpptr(dd)->fg));	args = CDR(args);
    nborder = LENGTH(border);
 
    PROTECT(lty = FixupLty(CAR(args), Rf_gpptr(dd)->lty));	args = CDR(args);
    nlty = length(lty);
 
    sxpd = CAR(args);
    if (sxpd != R_NilValue)
	xpd = asInteger(sxpd);
    else
	xpd = Rf_gpptr(dd)->xpd;
    args = CDR(args);
 
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
 
    if (xpd == NA_INTEGER)
	Rf_gpptr(dd)->xpd = 2;
    else
	Rf_gpptr(dd)->xpd = xpd;
 
    GMode(1, dd);
 
    x = REAL(sx);
    y = REAL(sy);
    xold = NA_REAL;
    yold = NA_REAL;
    for (i = 0; i < nx; i++) {
	xx = x[i];
	yy = y[i];
	GConvert(&xx, &yy, USER, DEVICE, dd);
	if ((R_FINITE(xx) && R_FINITE(yy)) &&
	    !(R_FINITE(xold) && R_FINITE(yold)))
	    start = i; /* first point of current segment */
	else if ((R_FINITE(xold) && R_FINITE(yold)) &&
		 !(R_FINITE(xx) && R_FINITE(yy))) {
	    if (i-start > 1) {
		drawPolygon(i-start, x+start, y+start,
			    INTEGER(lty)[num%nlty],
			    INTEGER(col)[num%ncol],
			    INTEGER(border)[num%nborder], dd);
		num++;
	    }
	}
	else if ((R_FINITE(xold) && R_FINITE(yold)) && (i == nx-1)) { /* last */
	    drawPolygon(nx-start, x+start, y+start,
			INTEGER(lty)[num%nlty],
			INTEGER(col)[num%ncol],
			INTEGER(border)[num%nborder], dd);
	    num++;
	}
	xold = xx;
	yold = yy;
    }
 
    GMode(0, dd);
 
    GRestorePars(dd);
    UNPROTECT(3);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}
 
SEXP do_text(SEXP call, SEXP op, SEXP args, SEXP env)
{
/* text(xy, labels, adj, pos, offset,
 *	vfont, cex, col, font, xpd, ...)
 */
    SEXP sx, sy, sxy, sxpd, txt, adj, pos, cex, col, rawcol, font, vfont;
    int i, n, npos, ncex, ncol, nfont, ntxt, xpd;
    double adjx = 0, adjy = 0, offset = 0.5;
    double *x, *y;
    double xx, yy;
    Rboolean vectorFonts = FALSE;
    SEXP string, originalArgs = args;
    DevDesc *dd = CurrentDevice();
 
    GCheckState(dd);
 
    if (length(args) < 3) errorcall(call, "too few arguments");
 
    PLOT_XY_DEALING("text");
 
    /* labels */
    txt = CAR(args);
    if (isSymbol(txt) || isLanguage(txt))
	txt = coerceVector(txt, EXPRSXP);
    else if (!isExpression(txt))
	txt = coerceVector(txt, STRSXP);
    PROTECT(txt);
    if (length(txt) <= 0)
	errorcall(call, "zero length 'labels'");
    args = CDR(args);
 
    PROTECT(adj = CAR(args));
    if (isNull(adj) || (isNumeric(adj) && length(adj) == 0)) {
	adjx = Rf_gpptr(dd)->adj;
	adjy = NA_REAL;
    }
    else if (isReal(adj)) {
	if (LENGTH(adj) == 1) {
	    adjx = REAL(adj)[0];
	    adjy = NA_REAL;
	}
	else {
	    adjx = REAL(adj)[0];
	    adjy = REAL(adj)[1];
	}
    }
    else if (isInteger(adj)) {
	if (LENGTH(adj) == 1) {
	    adjx = INTEGER(adj)[0];
	    adjy = NA_REAL;
	}
	else {
	    adjx = INTEGER(adj)[0];
	    adjy = INTEGER(adj)[1];
	}
    }
    else errorcall(call, "invalid adj value");
    args = CDR(args);
 
    PROTECT(pos = coerceVector(CAR(args), INTSXP));
    npos = length(pos);
    for (i = 0; i < npos; i++)
	if (INTEGER(pos)[i] < 1 || INTEGER(pos)[i] > 4)
	    errorcall(call, "invalid pos value");
    args = CDR(args);
 
    offset = GConvertXUnits(asReal(CAR(args)), CHARS, INCHES, dd);
    args = CDR(args);
 
    PROTECT(vfont = FixupVFont(CAR(args)));
    if (!isNull(vfont))
	vectorFonts = TRUE;
    args = CDR(args);
 
    PROTECT(cex = FixupCex(CAR(args), 1.0));
    ncex = LENGTH(cex);
    args = CDR(args);
 
    rawcol = CAR(args);
    PROTECT(col = FixupCol(rawcol, NA_INTEGER));
    ncol = LENGTH(col);
    args = CDR(args);
 
    PROTECT(font = FixupFont(CAR(args), NA_INTEGER));
    nfont = LENGTH(font);
    args = CDR(args);
 
    sxpd = CAR(args); /* xpd: NULL -> par("xpd") */
    if (sxpd != R_NilValue)
	xpd = asInteger(sxpd);
    else
	xpd = Rf_gpptr(dd)->xpd;
    args = CDR(args);
 
    x = REAL(sx);
    y = REAL(sy);
    n = LENGTH(sx);
    ntxt = LENGTH(txt);
 
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
 
    Rf_gpptr(dd)->xpd = (xpd == NA_INTEGER)? 2 : xpd;
 
    GMode(1, dd);
    for (i = 0; i < n; i++) {
	xx = x[i % n];
	yy = y[i % n];
	GConvert(&xx, &yy, USER, INCHES, dd);
	if (R_FINITE(xx) && R_FINITE(yy)) {
	    if (ncol && !isNAcol(rawcol, i, ncol))
		Rf_gpptr(dd)->col = INTEGER(col)[i % ncol];
	    else
		Rf_gpptr(dd)->col = Rf_dpptr(dd)->col;
	    if (ncex && R_FINITE(REAL(cex)[i%ncex]))
		Rf_gpptr(dd)->cex = Rf_gpptr(dd)->cexbase * REAL(cex)[i % ncex];
	    else
		Rf_gpptr(dd)->cex = Rf_gpptr(dd)->cexbase;
	    if (nfont && INTEGER(font)[i % nfont] != NA_INTEGER)
		Rf_gpptr(dd)->font = INTEGER(font)[i % nfont];
	    else
		Rf_gpptr(dd)->font = Rf_dpptr(dd)->font;
	    if (npos > 0) {
		switch(INTEGER(pos)[i % npos]) {
		case 1:
		    yy = yy - offset;
		    adjx = 0.5;
		    adjy = 1 - (0.5 - Rf_gpptr(dd)->yCharOffset);
		    break;
		case 2:
		    xx = xx - offset;
		    adjx = 1;
		    adjy = Rf_gpptr(dd)->yCharOffset;
		    break;
		case 3:
		    yy = yy + offset;
		    adjx = 0.5;
		    adjy = 0;
		    break;
		case 4:
		    xx = xx + offset;
		    adjx = 0;
		    adjy = Rf_gpptr(dd)->yCharOffset;
		    break;
		}
	    }
	    if (vectorFonts) {
		string = STRING_ELT(txt, i % ntxt);
		if(string != NA_STRING)
		    GVText(xx, yy, INCHES, CHAR(string),
			   INTEGER(vfont)[0], INTEGER(vfont)[1],
			   adjx, adjy, Rf_gpptr(dd)->srt, dd);
	    } else if (isExpression(txt)) {
		GMathText(xx, yy, INCHES, VECTOR_ELT(txt, i % ntxt),
			  adjx, adjy, Rf_gpptr(dd)->srt, dd);
	    } else {
		string = STRING_ELT(txt, i % ntxt);
		if(string != NA_STRING)
		    GText(xx, yy, INCHES, CHAR(string),
			  adjx, adjy, Rf_gpptr(dd)->srt, dd);
	    }
	}
    }
    GMode(0, dd);
 
    GRestorePars(dd);
    UNPROTECT(7);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}
 
static double ComputeAdjValue(double adj, int side, int las)
{
    if (!R_FINITE(adj)) {
	switch(las) {
	case 0:/* parallel to axis */
	    adj = 0.5; break;
	case 1:/* horizontal */
	    switch(side) {
	    case 1:
	    case 3: adj = 0.5; break;
	    case 2: adj = 1.0; break;
	    case 4: adj = 0.0; break;
	    }
	    break;
	case 2:/* perpendicular to axis */
	    switch(side) {
	    case 1:
	    case 2: adj = 1.0; break;
	    case 3:
	    case 4: adj = 0.0; break;
	    }
	    break;
	case 3:/* vertical */
	    switch(side) {
	    case 1: adj = 1.0; break;
	    case 3: adj = 0.0; break;
	    case 2:
	    case 4: adj = 0.5; break;
	    }
	    break;
	}
    }
    return adj;
}
 
static double ComputeAtValueFromAdj(double adj, int side, int outer,
				    DevDesc *dd)
{
    double at = 0;		/* -Wall */
    switch(side % 2) {
    case 0:
	at  = outer ? adj : yNPCtoUsr(adj, dd);
	break;
    case 1:
	at = outer ? adj : xNPCtoUsr(adj, dd);
	break;
    }
    return at;
}
 
static double ComputeAtValue(double at, double adj,
			     int side, int las, int outer,
			     DevDesc *dd)
{
    if (!R_FINITE(at)) {
	/* If the text is parallel to the axis, use "adj" for "at"
	 * Otherwise, centre the text
	 */
	switch(las) {
	case 0:/* parallel to axis */
	    at = ComputeAtValueFromAdj(adj, side, outer, dd);
	    break;
	case 1:/* horizontal */
	    switch(side) {
	    case 1:
	    case 3:
		at = ComputeAtValueFromAdj(adj, side, outer, dd);
		break;
	    case 2:
	    case 4:
		at = outer ? 0.5 : yNPCtoUsr(0.5, dd);
		break;
	    }
	    break;
	case 2:/* perpendicular to axis */
	    switch(side) {
	    case 1:
	    case 3:
		at = outer ? 0.5 : xNPCtoUsr(0.5, dd);
		break;
	    case 2:
	    case 4:
		at = outer ? 0.5 : yNPCtoUsr(0.5, dd);
		break;
	    }
	    break;
	case 3:/* vertical */
	    switch(side) {
	    case 1:
	    case 3:
		at = outer ? 0.5 : xNPCtoUsr(0.5, dd);
		break;
	    case 2:
	    case 4:
		at = ComputeAtValueFromAdj(adj, side, outer, dd);
		break;
	    }
	    break;
	}
    }
    return at;
}
 
/* mtext(text,
	 side = 3,
	 line = 0,
	 outer = TRUE,
	 at = NA,
	 adj = NA,
	 cex = NA,
	 col = NA,
	 font = NA,
	 vfont = NULL,
	 ...) */
 
SEXP do_mtext(SEXP call, SEXP op, SEXP args, SEXP env)
{
    SEXP text, side, line, outer, at, adj, cex, col, font, vfont, string;
    SEXP rawcol;
    int ntext, nside, nline, nouter, nat, nadj, ncex, ncol, nfont;
    Rboolean dirtyplot = FALSE, gpnewsave = FALSE, dpnewsave = FALSE;
    Rboolean vectorFonts = FALSE;
    int i, n, fontsave, colsave;
    double cexsave;
    SEXP originalArgs = args;
    DevDesc *dd = CurrentDevice();
 
    GCheckState(dd);
 
    if (length(args) < 9)
	errorcall(call, "too few arguments");
 
    /* Arg1 : text= */
    text = CAR(args);
    if (isSymbol(text) || isLanguage(text))
	text = coerceVector(text, EXPRSXP);
    else if (!isExpression(text))
	text = coerceVector(text, STRSXP);
    PROTECT(text);
    n = ntext = length(text);
    if (ntext <= 0)
	errorcall(call, "zero length \"text\" specified");
    args = CDR(args);
 
    /* Arg2 : side= */
    PROTECT(side = coerceVector(CAR(args), INTSXP));
    nside = length(side);
    if (nside <= 0) errorcall(call, "zero length \"side\" specified");
    if (n < nside) n = nside;
    args = CDR(args);
 
    /* Arg3 : line= */
    PROTECT(line = coerceVector(CAR(args), REALSXP));
    nline = length(line);
    if (nline <= 0) errorcall(call, "zero length \"line\" specified");
    if (n < nline) n = nline;
    args = CDR(args);
 
    /* Arg4 : outer= */
    /* outer == NA => outer <- 0 */
    PROTECT(outer = coerceVector(CAR(args), INTSXP));
    nouter = length(outer);
    if (nouter <= 0) errorcall(call, "zero length \"outer\" specified");
    if (n < nouter) n = nouter;
    args = CDR(args);
 
    /* Arg5 : at= */
    PROTECT(at = coerceVector(CAR(args), REALSXP));
    nat = length(at);
    if (nat <= 0) errorcall(call, "zero length \"at\" specified");
    if (n < nat) n = nat;
    args = CDR(args);
 
    /* Arg6 : adj= */
    PROTECT(adj = coerceVector(CAR(args), REALSXP));
    nadj = length(adj);
    if (nadj <= 0) errorcall(call, "zero length \"adj\" specified");
    if (n < nadj) n = nadj;
    args = CDR(args);
 
    /* Arg7 : cex */
    PROTECT(cex = FixupCex(CAR(args), 1.0));
    ncex = length(cex);
    if (ncex <= 0) errorcall(call, "zero length \"cex\" specified");
    if (n < ncex) n = ncex;
    args = CDR(args);
 
    /* Arg8 : col */
    rawcol = CAR(args);
    PROTECT(col = FixupCol(rawcol, NA_INTEGER));
    ncol = length(col);
    if (ncol <= 0) errorcall(call, "zero length \"col\" specified");
    if (n < ncol) n = ncol;
    args = CDR(args);
 
    /* Arg9 : font */
    PROTECT(font = FixupFont(CAR(args), NA_INTEGER));
    nfont = length(font);
    if (nfont <= 0) errorcall(call, "zero length \"font\" specified");
    if (n < nfont) n = nfont;
    args = CDR(args);
 
    /* Arg10 : vfont */
    PROTECT(vfont = FixupVFont(CAR(args)));
    if (!isNull(vfont))
	vectorFonts = TRUE;
    args = CDR(args);
 
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
 
    /* If we only scribble in the outer margins, */
    /* we don't want to mark the plot as dirty. */
 
    dirtyplot = FALSE;
    gpnewsave = Rf_gpptr(dd)->new;
    dpnewsave = Rf_dpptr(dd)->new;
    cexsave = Rf_gpptr(dd)->cex;
    fontsave = Rf_gpptr(dd)->font;
    colsave = Rf_gpptr(dd)->col;
 
    /* override par("xpd") and force clipping to figure region */
    /* NOTE: don't override to _reduce_ clipping region */
    if (Rf_gpptr(dd)->xpd < 1)
	Rf_gpptr(dd)->xpd = 1;
 
    if (outer) {
	gpnewsave = Rf_gpptr(dd)->new;
	dpnewsave = Rf_dpptr(dd)->new;
	/* override par("xpd") and force clipping to device region */
	Rf_gpptr(dd)->xpd = 2;
    }
    GMode(1, dd);
 
    for (i = 0; i < n; i++) {
	double atval = REAL(at)[i%nat];
	double adjval = REAL(adj)[i%nadj];
	double cexval = REAL(cex)[i%ncex];
	double lineval = REAL(line)[i%nline];
	int outerval = INTEGER(outer)[i%nouter];
	int sideval = INTEGER(side)[i%nside];
	int fontval = INTEGER(font)[i%nfont];
	int colval = INTEGER(col)[i%ncol];
 
	if (outerval == NA_INTEGER) outerval = 0;
	/* Note : we ignore any shrinking produced */
	/* by mfrow / mfcol specs here.	 I.e. don't */
	/* Rf_gpptr(dd)->cexbase. */
	if (R_FINITE(cexval)) Rf_gpptr(dd)->cex = cexval;
	else cexval = cexsave;
	Rf_gpptr(dd)->font = (fontval == NA_INTEGER) ? fontsave : fontval;
	if (isNAcol(rawcol, i, ncol))
	    Rf_gpptr(dd)->col = colsave;
	else
	    Rf_gpptr(dd)->col = colval;
	Rf_gpptr(dd)->adj = ComputeAdjValue(adjval, sideval, Rf_gpptr(dd)->las);
	atval = ComputeAtValue(atval, Rf_gpptr(dd)->adj, sideval, Rf_gpptr(dd)->las,
			       outerval, dd);
 
	if (vectorFonts) {
	    string = STRING_ELT(text, i%ntext);
#ifdef GMV_implemented
	    if(string != NA_STRING)
		GMVText(CHAR(string),
			INTEGER(vfont)[0], INTEGER(vfont)[1],
			sideval, lineval, outerval, atval,
			Rf_gpptr(dd)->las, dd);
#else
	    warningcall(call,"Hershey fonts not yet implemented for mtext()");
	    if(string != NA_STRING)
		GMtext(CHAR(string), sideval, lineval, outerval, atval,
		       Rf_gpptr(dd)->las, dd);
#endif
	}
	else if (isExpression(text))
	    GMMathText(VECTOR_ELT(text, i%ntext),
		       sideval, lineval, outerval, atval, Rf_gpptr(dd)->las, dd);
	else {
	    string = STRING_ELT(text, i%ntext);
	    if(string != NA_STRING)
		GMtext(CHAR(string), sideval, lineval, outerval, atval,
		       Rf_gpptr(dd)->las, dd);
	}
 
	if (outerval == 0) dirtyplot = TRUE;
    }
    GMode(0, dd);
 
    GRestorePars(dd);
    if (!dirtyplot) {
	Rf_gpptr(dd)->new = gpnewsave;
	Rf_dpptr(dd)->new = dpnewsave;
    }
    UNPROTECT(10);
 
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}/* do_mtext */
 
 
SEXP do_title(SEXP call, SEXP op, SEXP args, SEXP env)
{
/* Annotation for plots :
 
   title(main, sub, xlab, ylab,
	 line, outer,
	 ...) */
 
    SEXP Main, xlab, ylab, sub, vfont, string;
    double adj, adjy, cex, offset, line, hpos, vpos, where;
    int col, font, outer;
    int i, n;
    SEXP originalArgs = args;
    DevDesc *dd = CurrentDevice();
 
    GCheckState(dd);
 
    if (length(args) < 6) errorcall(call, "too few arguments");
 
    Main = sub = xlab = ylab = R_NilValue;
 
    if (CAR(args) != R_NilValue && LENGTH(CAR(args)) > 0)
	Main = CAR(args);
    args = CDR(args);
 
    if (CAR(args) != R_NilValue && LENGTH(CAR(args)) > 0)
	sub = CAR(args);
    args = CDR(args);
 
    if (CAR(args) != R_NilValue && LENGTH(CAR(args)) > 0)
	xlab = CAR(args);
    args = CDR(args);
 
    if (CAR(args) != R_NilValue && LENGTH(CAR(args)) > 0)
	ylab = CAR(args);
    args = CDR(args);
 
    line = asReal(CAR(args));
    args = CDR(args);
 
    outer = asLogical(CAR(args));
    if (outer == NA_LOGICAL) outer = 0;
    args = CDR(args);
 
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
 
    /* override par("xpd") and force clipping to figure region */
    /* NOTE: don't override to _reduce_ clipping region */
    if (Rf_gpptr(dd)->xpd < 1)
	Rf_gpptr(dd)->xpd = 1;
    if (outer)
	Rf_gpptr(dd)->xpd = 2;
    adj = Rf_gpptr(dd)->adj;
 
    GMode(1, dd);
    if (Main != R_NilValue) {
	cex = Rf_gpptr(dd)->cexmain;
	col = Rf_gpptr(dd)->colmain;
	font = Rf_gpptr(dd)->fontmain;
	GetTextArg(call, Main, &Main, &col, &cex, &font, &vfont);
	Rf_gpptr(dd)->col = col;
	Rf_gpptr(dd)->cex = Rf_gpptr(dd)->cexbase * cex;
	Rf_gpptr(dd)->font = font;
	if (outer) {
	    if (R_FINITE(line)) {
		vpos = line;
		adjy = 0;
	    }
	    else {
		vpos = 0.5 * Rf_gpptr(dd)->oma[2];
		adjy = 0.5;
	    }
	    hpos = adj;
	    where = OMA3;
	}
	else {
	    if (R_FINITE(line)) {
		vpos = line;
		adjy = 0;
	    }
	    else {
		vpos = 0.5 * Rf_gpptr(dd)->mar[2];
		adjy = 0.5;
	    }
	    hpos = GConvertX(adj, NPC, USER, dd);
	    where = MAR3;
	}
	if (isExpression(Main)) {
	    GMathText(hpos, vpos, where, VECTOR_ELT(Main, 0),
		      adj, 0.5, 0.0, dd);
	}
	else {
	  n = length(Main);
	  offset = 0.5 * (n - 1) + vpos;
	  for (i = 0; i < n; i++) {
		string = STRING_ELT(Main, i);
		if(string != NA_STRING)
		    GText(hpos, offset - i, where, CHAR(string), adj,
			  adjy, 0.0, dd);
	  }
	}
    }
    if (sub != R_NilValue) {
	cex = Rf_gpptr(dd)->cexsub;
	col = Rf_gpptr(dd)->colsub;
	font = Rf_gpptr(dd)->fontsub;
	GetTextArg(call, sub, &sub, &col, &cex, &font, &vfont);
	Rf_gpptr(dd)->col = col;
	Rf_gpptr(dd)->cex = Rf_gpptr(dd)->cexbase * cex;
	Rf_gpptr(dd)->font = font;
	if (R_FINITE(line))
	    vpos = line;
	else
	    vpos = Rf_gpptr(dd)->mgp[0] + 1;
	if (outer) {
	    hpos = adj;
	    where = 1;
	}
	else {
	    hpos = GConvertX(adj, NPC, USER, dd);
	    where = 0;
	}
	if (isExpression(sub))
	    GMMathText(VECTOR_ELT(sub, 0), 1, vpos, where,
		       hpos, 0, dd);
	else {
	    n = length(sub);
	    for (i = 0; i < n; i++) {
		string = STRING_ELT(sub, i);
		if(string != NA_STRING)
		    GMtext(CHAR(string), 1, vpos, where, hpos, 0, dd);
	    }
	}
    }
    if (xlab != R_NilValue) {
	cex = Rf_gpptr(dd)->cexlab;
	col = Rf_gpptr(dd)->collab;
	font = Rf_gpptr(dd)->fontlab;
	GetTextArg(call, xlab, &xlab, &col, &cex, &font, &vfont);
	Rf_gpptr(dd)->cex = Rf_gpptr(dd)->cexbase * cex;
	Rf_gpptr(dd)->col = col;
	Rf_gpptr(dd)->font = font;
	if (R_FINITE(line))
	    vpos = line;
	else
	    vpos = Rf_gpptr(dd)->mgp[0];
	if (outer) {
	    hpos = adj;
	    where = 1;
	}
	else {
	    hpos = GConvertX(adj, NPC, USER, dd);
	    where = 0;
	}
	if (isExpression(xlab))
	    GMMathText(VECTOR_ELT(xlab, 0), 1, vpos, where,
		       hpos, 0, dd);
	else {
	    n = length(xlab);
	    for (i = 0; i < n; i++) {
		string = STRING_ELT(xlab, i);
		if(string != NA_STRING)
		    GMtext(CHAR(string), 1, vpos + i, where, hpos, 0, dd);
	    }
	}
    }
    if (ylab != R_NilValue) {
	cex = Rf_gpptr(dd)->cexlab;
	col = Rf_gpptr(dd)->collab;
	font = Rf_gpptr(dd)->fontlab;
	GetTextArg(call, ylab, &ylab, &col, &cex, &font, &vfont);
	Rf_gpptr(dd)->cex = Rf_gpptr(dd)->cexbase * cex;
	Rf_gpptr(dd)->col = col;
	Rf_gpptr(dd)->font = font;
	if (R_FINITE(line))
	    vpos = line;
	else
	    vpos = Rf_gpptr(dd)->mgp[0];
	if (outer) {
	    hpos = adj;
	    where = 1;
	}
	else {
	    hpos = GConvertY(adj, NPC, USER, dd);
	    where = 0;
	}
	if (isExpression(ylab))
	    GMMathText(VECTOR_ELT(ylab, 0), 2, vpos, where,
		       hpos, 0, dd);
	else {
	    n = length(ylab);
	    for (i = 0; i < n; i++) {
		string = STRING_ELT(ylab, i);
		if(string != NA_STRING)
		    GMtext(CHAR(string), 2, vpos - i, where, hpos, 0, dd);
	    }
	}
    }
    GMode(0, dd);
    GRestorePars(dd);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}
 
 
/*  abline(a, b, h, v, col, lty, lwd, ...)
    draw lines in intercept/slope form.	 */
 
static void getxlimits(double *x, DevDesc *dd) {
    /*
     * xpd = 0 means clip to current plot region
     * xpd = 1 means clip to current figure region
     * xpd = 2 means clip to device region
     */
    switch (Rf_gpptr(dd)->xpd) {
    case 0:
	x[0] = Rf_gpptr(dd)->usr[0];
	x[1] = Rf_gpptr(dd)->usr[1];
	break;
    case 1:
	x[0] = GConvertX(0, NFC, USER, dd);
	x[1] = GConvertX(1, NFC, USER, dd);
	break;
    case 2:
	x[0] = GConvertX(0, NDC, USER, dd);
	x[1] = GConvertX(1, NDC, USER, dd);
	break;
    }
}
 
static void getylimits(double *y, DevDesc *dd) {
    switch (Rf_gpptr(dd)->xpd) {
    case 0:
	y[0] = Rf_gpptr(dd)->usr[2];
	y[1] = Rf_gpptr(dd)->usr[3];
	break;
    case 1:
	y[0] = GConvertY(0, NFC, USER, dd);
	y[1] = GConvertY(1, NFC, USER, dd);
	break;
    case 2:
	y[0] = GConvertY(0, NDC, USER, dd);
	y[1] = GConvertY(1, NDC, USER, dd);
	break;
    }
}
 
SEXP do_abline(SEXP call, SEXP op, SEXP args, SEXP env)
{
    SEXP a, b, h, v, untf, col, lty, lwd;
    int i, ncol, nlines, nlty, nlwd, lstart, lstop;
    double aa, bb, x[2], y[2];
    SEXP originalArgs = args;
    DevDesc *dd = CurrentDevice();
 
    GCheckState(dd);
 
    if (length(args) < 5) errorcall(call, "too few arguments");
 
    if ((a = CAR(args)) != R_NilValue)
	SETCAR(args, a = coerceVector(a, REALSXP));
    args = CDR(args);
 
    if ((b = CAR(args)) != R_NilValue)
	SETCAR(args, b = coerceVector(b, REALSXP));
    args = CDR(args);
 
    if ((h = CAR(args)) != R_NilValue)
	SETCAR(args, h = coerceVector(h, REALSXP));
    args = CDR(args);
 
    if ((v = CAR(args)) != R_NilValue)
	SETCAR(args, v = coerceVector(v, REALSXP));
    args = CDR(args);
 
    if ((untf = CAR(args)) != R_NilValue)
	SETCAR(args, untf = coerceVector(untf, LGLSXP));
    args = CDR(args);
 
 
    PROTECT(col = FixupCol(CAR(args), NA_INTEGER));	args = CDR(args);
    ncol = LENGTH(col);
 
    PROTECT(lty = FixupLty(CAR(args), Rf_gpptr(dd)->lty)); args = CDR(args);
    nlty = length(lty);
 
    PROTECT(lwd = FixupLwd(CAR(args), Rf_gpptr(dd)->lwd)); args = CDR(args);
    nlwd = length(lwd);
 
    GSavePars(dd);
 
    ProcessInlinePars(args, dd, call);
 
    nlines = 0;
 
    if (a != R_NilValue) {
	if (b == R_NilValue) {
	    if (LENGTH(a) != 2)
		errorcall(call, "invalid a=, b= specification");
	    aa = REAL(a)[0];
	    bb = REAL(a)[1];
	}
	else {
	    aa = asReal(a);
	    bb = asReal(b);
	}
	if (!R_FINITE(aa) || !R_FINITE(bb))
	    errorcall(call, "\"a\" and \"b\" must be finite");
	Rf_gpptr(dd)->col = INTEGER(col)[0];
	Rf_gpptr(dd)->lwd = REAL(lwd)[0];
	if (nlty && INTEGER(lty)[0] != NA_INTEGER)
	    Rf_gpptr(dd)->lty = INTEGER(lty)[0];
	else
	    Rf_gpptr(dd)->lty = Rf_dpptr(dd)->lty;
	GMode(1, dd);
	/* FIXME?
	 * Seems like the logic here is just draw from xmin to xmax
	 * and you're guaranteed to draw at least from ymin to ymax
	 * This MAY cause a problem at some stage when the line being
	 * drawn is VERY steep -- and the problem is worse now that
	 * abline will potentially draw to the extents of the device
	 * (when xpd=NA).  NOTE that R's internal clipping protects the
	 * device drivers from stupidly large numbers, BUT there is
	 * still a risk that we could produce a number which is too
	 * big for the computer's brain.
	 * Paul.
	 *
	 * The problem is worse -- you could get NaN, which at least the
	 * X11 device coerces to -2^31 <TSL>
	 */
	getxlimits(x, dd);
	if (R_FINITE(Rf_gpptr(dd)->lwd)) {
	    if (LOGICAL(untf)[0] == 1 && (Rf_gpptr(dd)->xlog || Rf_gpptr(dd)->ylog)) {
		double xx[101], yy[101];
		double xstep = (x[1] - x[0])/100;
		for (i = 0; i < 100; i++) {
		    xx[i] = x[0] + i*xstep;
		    yy[i] = aa + xx[i] * bb;
		}
		xx[100] = x[1];
		yy[100] = aa + x[1] * bb;
 
		/* now get rid of -ve values */
		lstart=0;lstop=100;
		if (Rf_gpptr(dd)->xlog){
			for(;xx[lstart]<=0 && lstart<101;lstart++);
			for(;xx[lstop]<=0 && lstop>0;lstop--);
		}
		if (Rf_gpptr(dd)->ylog){
			for(;yy[lstart]<=0 && lstart<101;lstart++);
			for(;yy[lstop]<=0 && lstop>0;lstop--);
		}
 
 
		GPolyline(lstop-lstart+1, xx+lstart, yy+lstart, USER, dd);
	    }
	    else {
		double x0, x1;
 
		x0 = ( Rf_gpptr(dd)->xlog ) ?	log10(x[0]) : x[0];
		x1 = ( Rf_gpptr(dd)->xlog ) ?	log10(x[1]) : x[1];
 
		y[0] = aa + x0 * bb;
		y[1] = aa + x1 * bb;
 
		if ( Rf_gpptr(dd)->ylog ){
		    y[0] = pow(10.,y[0]);
		    y[1] = pow(10.,y[1]);
		}
 
		GLine(x[0], y[0], x[1], y[1], USER, dd);
	    }
	}
	GMode(0, dd);
	nlines++;
    }
    if (h != R_NilValue) {
	GMode(1, dd);
	for (i = 0; i < LENGTH(h); i++) {
	    Rf_gpptr(dd)->col = INTEGER(col)[nlines % ncol];
	    if (nlty && INTEGER(lty)[nlines % nlty] != NA_INTEGER)
		Rf_gpptr(dd)->lty = INTEGER(lty)[nlines % nlty];
	    else
		Rf_gpptr(dd)->lty = Rf_dpptr(dd)->lty;
	    Rf_gpptr(dd)->lwd = REAL(lwd)[nlines % nlwd];
	    aa = REAL(h)[i];
	    if (R_FINITE(aa) && R_FINITE(Rf_gpptr(dd)->lwd)) {
		getxlimits(x, dd);
		y[0] = aa;
		y[1] = aa;
		GLine(x[0], y[0], x[1], y[1], USER, dd);
	    }
	    nlines++;
	}
	GMode(0, dd);
    }
    if (v != R_NilValue) {
	GMode(1, dd);
	for (i = 0; i < LENGTH(v); i++) {
	    Rf_gpptr(dd)->col = INTEGER(col)[nlines % ncol];
	    if (nlty && INTEGER(lty)[nlines % nlty] != NA_INTEGER)
		Rf_gpptr(dd)->lty = INTEGER(lty)[nlines % nlty];
	    else
		Rf_gpptr(dd)->lty = Rf_dpptr(dd)->lty;
	    Rf_gpptr(dd)->lwd = REAL(lwd)[nlines % nlwd];
	    aa = REAL(v)[i];
	    if (R_FINITE(aa) && R_FINITE(Rf_gpptr(dd)->lwd)) {
		getylimits(y, dd);
		x[0] = aa;
		x[1] = aa;
		GLine(x[0], y[0], x[1], y[1], USER, dd);
	    }
	    nlines++;
	}
	GMode(0, dd);
    }
    UNPROTECT(3);
    GRestorePars(dd);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}
 
SEXP do_box(SEXP call, SEXP op, SEXP args, SEXP env)
{
/*     box(which="plot", lty="solid", ...)
       --- which is coded, 1 = plot, 2 = figure, 3 = inner, 4 = outer.
*/
    int which, col;
    SEXP colsxp, fgsxp;
    SEXP originalArgs = args;
    DevDesc *dd = CurrentDevice();
 
    GCheckState(dd);
    GSavePars(dd);
    which = asInteger(CAR(args)); args = CDR(args);
    if (which < 1 || which > 4)
	errorcall(call, "invalid \"which\" specification");
    /*
     * If specified non-NA col then use that, else ...
     *
     * if specified non-NA fg then use that, else ...
     *
     * else use par("col")
     */
    col= Rf_gpptr(dd)->col;
    ProcessInlinePars(args, dd, call);
    colsxp = getInlinePar(args, "col");
    if (isNAcol(colsxp, 0, 1)) {
	fgsxp = getInlinePar(args, "fg");
	if (isNAcol(fgsxp, 0, 1))
	    Rf_gpptr(dd)->col = col;
	else
	    Rf_gpptr(dd)->col = Rf_gpptr(dd)->fg;
    }
    /* override par("xpd") and force clipping to device region */
    Rf_gpptr(dd)->xpd = 2;
    GMode(1, dd);
    GBox(which, dd);
    GMode(0, dd);
    GRestorePars(dd);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
}
 
static void drawPointsLines(double xp, double yp, double xold, double yold,
			    char type, int first, DevDesc *dd)
{
    if (type == 'p' || type == 'o')
	GSymbol(xp, yp, DEVICE, Rf_gpptr(dd)->pch, dd);
    if ((type == 'l' || type == 'o') && !first)
	GLine(xold, yold, xp, yp, DEVICE, dd);
}
 
SEXP do_locator(SEXP call, SEXP op, SEXP args, SEXP env)
{
    SEXP x, y, nobs, ans, saveans, stype = R_NilValue;
    int i, n, type='p';
    double xp, yp, xold=0, yold=0;
    DevDesc *dd = CurrentDevice();
 
    /* If replaying, just draw the points and lines that were recorded */
    if (call == R_NilValue) {
	x = CAR(args); args = CDR(args);
	y = CAR(args); args = CDR(args);
	nobs = CAR(args); args = CDR(args);
	n = INTEGER(nobs)[0];
	stype = CAR(args); args = CDR(args);
	type = CHAR(STRING_ELT(stype, 0))[0];
	if (type != 'n') {
	    GMode(1, dd);
	    for (i=0; i<n; i++) {
		xp = REAL(x)[i];
		yp = REAL(y)[i];
		GConvert(&xp, &yp, USER, DEVICE, dd);
		drawPointsLines(xp, yp, xold, yold, type, i==0, dd);
		xold = xp;
		yold = yp;
	    }
	    GMode(0, dd);
	}
	return R_NilValue;
    } else {
	GCheckState(dd);
 
	checkArity(op, args);
	n = asInteger(CAR(args));
	if (n <= 0 || n == NA_INTEGER)
	    error("invalid number of points in locator");
	args = CDR(args);
	if (isString(CAR(args)) && LENGTH(CAR(args)) == 1)
	    stype = CAR(args);
	else
	    errorcall(call, "invalid plot type");
	type = CHAR(STRING_ELT(stype, 0))[0];
	PROTECT(x = allocVector(REALSXP, n));
	PROTECT(y = allocVector(REALSXP, n));
	PROTECT(nobs=allocVector(INTSXP,1));
	i = 0;
 
	GMode(2, dd);
	while (i < n) {
	    if (!GLocator(&(REAL(x)[i]), &(REAL(y)[i]), USER, dd))
		break;
	    if (type != 'n') {
		GMode(1, dd);
		xp = REAL(x)[i];
		yp = REAL(y)[i];
		GConvert(&xp, &yp, USER, DEVICE, dd);
		drawPointsLines(xp, yp, xold, yold, type, i==0, dd);
		GMode(2, dd);
		xold = xp; yold = yp;
	    }
	    i += 1;
	}
	GMode(0, dd);
	INTEGER(nobs)[0] = i;
	while (i < n) {
	    REAL(x)[i] = NA_REAL;
	    REAL(y)[i] = NA_REAL;
	    i += 1;
	}
	PROTECT(ans = allocList(3));
	SETCAR(ans, x);
	SETCADR(ans, y);
	SETCADDR(ans, nobs);
	PROTECT(saveans = allocList(4));
	SETCAR(saveans, x);
	SETCADR(saveans, y);
	SETCADDR(saveans, nobs);
	SETCADDDR(saveans, CAR(args));
	/* Record the points and lines that were drawn in the display list */
	recordGraphicOperation(op, saveans, dd);
	UNPROTECT(5);
	return ans;
    }
}
 
#define THRESHOLD	0.25
 
static void drawLabel(double xi, double yi, int pos, double offset, char *l,
		      DevDesc *dd)
{
    switch (pos) {
    case 4:
	xi = xi+offset;
	GText(xi, yi, INCHES, l, 0.0,
	      Rf_gpptr(dd)->yCharOffset, 0.0, dd);
	break;
    case 2:
	xi = xi-offset;
	GText(xi, yi, INCHES, l, 1.0,
	      Rf_gpptr(dd)->yCharOffset, 0.0, dd);
	break;
    case 3:
	yi = yi+offset;
	GText(xi, yi, INCHES, l, 0.5,
	      0.0, 0.0, dd);
	break;
    case 1:
	yi = yi-offset;
	GText(xi, yi, INCHES, l, 0.5,
	      1-(0.5-Rf_gpptr(dd)->yCharOffset),
	      0.0, dd);
    }
}
 
SEXP do_identify(SEXP call, SEXP op, SEXP args, SEXP env)
{
    SEXP ans, x, y, l, ind, pos, Offset, draw, saveans;
    double xi, yi, xp, yp, d, dmin, offset;
    int i, imin, k, n, npts, plot, posi, warn;
    DevDesc *dd = CurrentDevice();
 
    /* If we are replaying the display list, then just redraw the
       labels beside the identified points */
    if (call == R_NilValue) {
	ind = CAR(args); args = CDR(args);
	pos = CAR(args); args = CDR(args);
	x = CAR(args); args = CDR(args);
	y = CAR(args); args = CDR(args);
	Offset = CAR(args); args = CDR(args);
	l = CAR(args); args = CDR(args);
	draw = CAR(args);
	n = length(x);
	for (i=0; i<n; i++) {
	    plot = LOGICAL(ind)[i];
	    if (LOGICAL(draw)[0] && plot) {
		xi = REAL(x)[i];
		yi = REAL(y)[i];
		GConvert(&xi, &yi, USER, INCHES, dd);
		posi = INTEGER(pos)[i];
		offset = GConvertXUnits(asReal(Offset), CHARS, INCHES, dd);
		drawLabel(xi, yi, posi, offset, CHAR(STRING_ELT(l, i)), dd);
	    }
	}
	return R_NilValue;
    }
    else {
	GCheckState(dd);
 
	checkArity(op, args);
	x = CAR(args);
	args = CDR(args); y = CAR(args);
	args = CDR(args); l = CAR(args);
	args = CDR(args); npts = asInteger(CAR(args));
	args = CDR(args); plot = asLogical(CAR(args));
	args = CDR(args); Offset = CAR(args);
	if (npts <= 0 || npts == NA_INTEGER)
	    error("invalid number of points in identify");
	if (!isReal(x) || !isReal(y) || !isString(l) || !isReal(Offset))
	    errorcall(call, "incorrect argument type");
	if (LENGTH(x) != LENGTH(y) || LENGTH(x) != LENGTH(l))
	    errorcall(call, "different argument lengths");
	n = LENGTH(x);
	if (n <= 0) {
	    R_Visible = 0;
	    return NULL;
	}
 
	offset = GConvertXUnits(asReal(Offset), CHARS, INCHES, dd);
	PROTECT(ind = allocVector(LGLSXP, n));
	PROTECT(pos = allocVector(INTSXP, n));
	for (i = 0; i < n; i++)
	    LOGICAL(ind)[i] = 0;
 
	k = 0;
	GMode(2, dd);
	while (k < npts) {
	    if (!GLocator(&xp, &yp, INCHES, dd)) break;
	    dmin = DBL_MAX;
	    imin = -1;
	    for (i = 0; i < n; i++) {
		xi = REAL(x)[i];
		yi = REAL(y)[i];
		GConvert(&xi, &yi, USER, INCHES, dd);
		if (!R_FINITE(xi) || !R_FINITE(yi)) continue;
		d = hypot(xp-xi, yp-yi);
		if (d < dmin) {
		    imin = i;
		    dmin = d;
		}
	    }
	    /* can't use warning because we want to print immediately  */
	    /* might want to handle warn=2? */
	    warn = asInteger(GetOption(install("warn"), R_NilValue));
	    if (dmin > THRESHOLD) {
	        if(warn >= 0)
		    REprintf("warning: no point with %.2f inches\n",
                                        THRESHOLD);
	    }
	    else if (LOGICAL(ind)[imin]) {
	        if(warn >= 0 )
		    REprintf("warning: nearest point already identified\n");
	    }
	    else {
		k++;
		LOGICAL(ind)[imin] = 1;
 
		xi = REAL(x)[imin];
		yi = REAL(y)[imin];
		GConvert(&xi, &yi, USER, INCHES, dd);
		if (fabs(xp-xi) >= fabs(yp-yi)) {
		    if (xp >= xi) {
			INTEGER(pos)[imin] = 4;
		    }
		    else {
			INTEGER(pos)[imin] = 2;
		    }
		}
		else {
		    if (yp >= yi) {
			INTEGER(pos)[imin] = 3;
		    }
		    else {
			INTEGER(pos)[imin] = 1;
		    }
		}
		if (plot)
		    drawLabel(xi, yi, INTEGER(pos)[imin], offset,
			      CHAR(STRING_ELT(l, imin)), dd);
	    }
	}
	GMode(0, dd);
	PROTECT(ans = allocList(2));
	SETCAR(ans, ind);
	SETCADR(ans, pos);
	PROTECT(saveans = allocList(7));
	SETCAR(saveans, ind);
	SETCADR(saveans, pos);
	SETCADDR(saveans, x);
	SETCADDDR(saveans, y);
	SETCAD4R(saveans, Offset);
	SETCAD4R(CDR(saveans), l);
	PROTECT(draw = allocVector(LGLSXP, 1));
	LOGICAL(draw)[0] = plot;
	SETCAD4R(CDDR(saveans), draw);
 
	/* If we are recording, save enough information to be able to
	   redraw the text labels beside identified points */
	if (GRecording(call))
	    recordGraphicOperation(op, saveans, dd);
	UNPROTECT(5);
 
	return ans;
    }
}
 
/* strheight(str, units)  ||  strwidth(str, units) */
 
#define DO_STR_DIM(KIND) 						\
{									\
    SEXP ans, str, ch;							\
    int i, n, units;							\
    double cex, cexsave;						\
    DevDesc *dd = CurrentDevice();					\
									\
    checkArity(op, args);						\
    GCheckState(dd);							\
									\
    str = CAR(args);							\
    if (isSymbol(str) || isLanguage(str))				\
	str = coerceVector(str, EXPRSXP);				\
    else if (!isExpression(str))					\
	str = coerceVector(str, STRSXP);				\
    PROTECT(str);							\
    args = CDR(args);							\
									\
    if ((units = asInteger(CAR(args))) == NA_INTEGER || units < 0)	\
	errorcall(call, "invalid units");				\
    args = CDR(args);							\
									\
    if (isNull(CAR(args)))						\
	cex = Rf_gpptr(dd)->cex;					\
    else if (!R_FINITE(cex = asReal(CAR(args))) || cex <= 0.0)		\
	errorcall(call, "invalid cex value");				\
									\
    n = LENGTH(str);							\
    PROTECT(ans = allocVector(REALSXP, n));				\
    cexsave = Rf_gpptr(dd)->cex;					\
    Rf_gpptr(dd)->cex = cex * Rf_gpptr(dd)->cexbase;			\
    for (i = 0; i < n; i++)						\
	if (isExpression(str))						\
	    REAL(ans)[i] = GExpression ## KIND(VECTOR_ELT(str, i),	\
					     GMapUnits(units), dd);	\
	else {								\
	    ch = STRING_ELT(str, i);					\
	    REAL(ans)[i] = (ch == NA_STRING) ? 0.0 :			\
		GStr ## KIND(CHAR(ch), GMapUnits(units), dd);		\
	}								\
    Rf_gpptr(dd)->cex = cexsave;					\
    UNPROTECT(2);							\
    return ans;								\
}
 
SEXP do_strheight(SEXP call, SEXP op, SEXP args, SEXP env)
DO_STR_DIM(Height)
 
SEXP do_strwidth (SEXP call, SEXP op, SEXP args, SEXP env)
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
DO_STR_DIM(Width)
 
 
 
 
 
#undef DO_STR_DIM
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
static int *dnd_lptr;
static int *dnd_rptr;
static double *dnd_hght;
static double *dnd_xpos;
static double dnd_hang;
static double dnd_offset;
static SEXP *dnd_llabels;
 
 
static void drawdend(int node, double *x, double *y, DevDesc *dd)
{
/* Recursive function for 'hclust' dendrogram drawing:
 * Do left + Do right + Do myself
 * "do" : 1) label leafs (if there are) and __
 *        2) find coordinates to draw the  | |
 *        3) return (*x,*y) of "my anchor"
 */
    double xl, xr, yl, yr;
    double xx[4], yy[4];
    int k;
 
    *y = dnd_hght[node-1];
    /* left part  */
    k = dnd_lptr[node-1];
    if (k > 0) drawdend(k, &xl, &yl, dd);
    else {
	xl = dnd_xpos[-k-1];
	yl = (dnd_hang >= 0) ? *y - dnd_hang : 0;
 
	if(dnd_llabels[-k-1] != NA_STRING)
	    GText(xl, yl-dnd_offset, USER, CHAR(dnd_llabels[-k-1]),
		  1.0, 0.3, 90.0, dd);
    }
    /* right part */
    k = dnd_rptr[node-1];
    if (k > 0) drawdend(k, &xr, &yr, dd);
    else {
	xr = dnd_xpos[-k-1];
	yr = (dnd_hang >= 0) ? *y - dnd_hang : 0;
 
	if(dnd_llabels[-k-1] != NA_STRING)
	    GText(xr, yr-dnd_offset, USER, CHAR(dnd_llabels[-k-1]),
		  1.0, 0.3, 90.0, dd);
    }
    xx[0] = xl; yy[0] = yl;
    xx[1] = xl; yy[1] = *y;
    xx[2] = xr; yy[2] = *y;
    xx[3] = xr; yy[3] = yr;
    GPolyline(4, xx, yy, USER, dd);
    *x = 0.5 * (xl + xr);
}
 
 
SEXP do_dend(SEXP call, SEXP op, SEXP args, SEXP env)
{
    double x, y;
    int n;
 
    SEXP originalArgs;
    DevDesc *dd;
 
    dd = CurrentDevice();
    GCheckState(dd);
 
    originalArgs = args;
    if (length(args) < 6)
	errorcall(call, "too few arguments");
 
    /* n */
    n = asInteger(CAR(args));
    if (n == NA_INTEGER || n < 2)
	goto badargs;
    args = CDR(args);
 
    /* merge */
    if (TYPEOF(CAR(args)) != INTSXP || length(CAR(args)) != 2*n)
	goto badargs;
    dnd_lptr = &(INTEGER(CAR(args))[0]);
    dnd_rptr = &(INTEGER(CAR(args))[n]);
    args = CDR(args);
 
    /* height */
    if (TYPEOF(CAR(args)) != REALSXP || length(CAR(args)) != n)
	goto badargs;
    dnd_hght = REAL(CAR(args));
    args = CDR(args);
 
    /* ord = order(x$order) */
    if (length(CAR(args)) != n+1)
	goto badargs;
    dnd_xpos = REAL(coerceVector(CAR(args),REALSXP));
    args = CDR(args);
 
    /* hang */
    dnd_hang = asReal(CAR(args));
    if (!R_FINITE(dnd_hang))
	goto badargs;
    dnd_hang = dnd_hang * (dnd_hght[n-1] - dnd_hght[0]);
    args = CDR(args);
 
    /* labels */
    if (TYPEOF(CAR(args)) != STRSXP || length(CAR(args)) != n+1)
	goto badargs;
    dnd_llabels = STRING_PTR(CAR(args));
    args = CDR(args);
 
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
    Rf_gpptr(dd)->cex = Rf_gpptr(dd)->cexbase * Rf_gpptr(dd)->cex;
    dnd_offset = GConvertYUnits(GStrWidth("m", INCHES, dd), INCHES, USER, dd);
 
    /* override par("xpd") and force clipping to figure region */
    /* NOTE: don't override to _reduce_ clipping region */
    if (Rf_gpptr(dd)->xpd < 1)
	Rf_gpptr(dd)->xpd = 1;
 
    GMode(1, dd);
    drawdend(n, &x, &y, dd);
    GMode(0, dd);
    GRestorePars(dd);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    return R_NilValue;
 
  badargs:
    error("invalid dendrogram input");
    return R_NilValue;/* never used; to keep -Wall happy */
}
 
SEXP do_dendwindow(SEXP call, SEXP op, SEXP args, SEXP env)
{
    int i, imax, n;
    double pin, *ll, tmp, yval, *y, ymin, ymax, yrange, m;
    SEXP originalArgs, merge, height, llabels, str;
    char *vmax;
    DevDesc *dd;
 
    dd = CurrentDevice();
    GCheckState(dd);
    originalArgs = args;
    if (length(args) < 5)
	errorcall(call, "too few arguments");
    n = asInteger(CAR(args));
    if (n == NA_INTEGER || n < 2)
	goto badargs;
    args = CDR(args);
    if (TYPEOF(CAR(args)) != INTSXP || length(CAR(args)) != 2 * n)
	goto badargs;
    merge = CAR(args);
 
    args = CDR(args);
    if (TYPEOF(CAR(args)) != REALSXP || length(CAR(args)) != n)
	goto badargs;
    height = CAR(args);
 
 
 
 
    args = CDR(args);
    dnd_hang = asReal(CAR(args));
    if (!R_FINITE(dnd_hang))
	goto badargs;
 
    args = CDR(args);
    if (TYPEOF(CAR(args)) != STRSXP || length(CAR(args)) != n + 1)
	goto badargs;
 
    llabels = CAR(args);
 
    args = CDR(args);
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
    Rf_gpptr(dd)->cex = Rf_gpptr(dd)->cexbase * Rf_gpptr(dd)->cex;
    dnd_offset = GStrWidth("m", INCHES, dd);
    vmax = vmaxget();
    y =  (double*)R_alloc(n, sizeof(double));
    ll = (double*)R_alloc(n, sizeof(double));
    dnd_lptr = &(INTEGER(merge)[0]);
    dnd_rptr = &(INTEGER(merge)[n]);
    ymax = ymin = REAL(height)[0];
    for (i = 1; i < n; i++) {
	m = REAL(height)[i];
	if (m > ymax)
	    ymax = m;
	else if (m < ymin)
	    ymin = m;
    }
    pin = Rf_gpptr(dd)->pin[1];
    for (i = 0; i < n; i++) {
	str = STRING_ELT(llabels, i);
	ll[i] = (str == NA_STRING) ? 0.0 :
	    GStrWidth(CHAR(str), INCHES, dd) + dnd_offset;
    }
 
    imax = -1; yval = -DBL_MAX;
    if (dnd_hang >= 0) {
	ymin = ymax - (1 + dnd_hang) * (ymax - ymin);
	yrange = ymax - ymin;
	/* determine leaf heights */
	for (i = 0; i < n; i++) {
	    if (dnd_lptr[i] < 0)
		y[-dnd_lptr[i] - 1] = REAL(height)[i];
	    if (dnd_rptr[i] < 0)
		y[-dnd_rptr[i] - 1] = REAL(height)[i];
	}
	/* determine the most extreme label depth */
	/* assuming that we are using the full plot */
	/* window for the tree itself */
 
 
	for (i = 0; i < n; i++) {
	    tmp = ((ymax - y[i]) / yrange) * pin + ll[i];
	    if (tmp > yval) {
		yval = tmp;
		imax = i;
	    }
	}
    }
    else {
 
	yrange = ymax;
 
 
	for (i = 0; i < n; i++) {
	    tmp = pin + ll[i];
	    if (tmp > yval) {
		yval = tmp;
		imax = i;
	    }
	}
    }
    /* now determine how much to scale */
    ymin = ymax - (pin/(pin - ll[imax])) * yrange;
    GScale(1.0, n+1.0, 1 /* x */, dd);
    GScale(ymin, ymax, 2 /* y */, dd);
    GMapWin2Fig(dd);
    GRestorePars(dd);
    /* NOTE: only record operation if no "error"  */
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    vmaxset(vmax);
    return R_NilValue;
  badargs:
    error("invalid dendrogram input");
    return R_NilValue;/* never used; to keep -Wall happy */
}
 
 
SEXP do_erase(SEXP call, SEXP op, SEXP args, SEXP env)
{
    SEXP col;
    int ncol;
    DevDesc *dd = CurrentDevice();
    checkArity(op, args);
    PROTECT(col = FixupCol(CAR(args), NA_INTEGER));
    ncol = LENGTH(col);
    GSavePars(dd);
    GMode(1, dd);
    GRect(0.0, 0.0, 1.0, 1.0, NDC, INTEGER(col)[0], NA_INTEGER, dd);
    GMode(0, dd);
    GRestorePars(dd);
    UNPROTECT(1);
    return R_NilValue;
}
 
 
SEXP do_getSnapshot(SEXP call, SEXP op, SEXP args, SEXP env)
{
    DevDesc *dd = CurrentDevice();
 
    checkArity(op, args);
    if (dd->newDevStruct) {
	return GEcreateSnapshot((GEDevDesc*) dd);
    } else {
	errorcall(call, "can't take snapshot of old-style device");
	return R_NilValue;
    }
}
 
SEXP do_playSnapshot(SEXP call, SEXP op, SEXP args, SEXP env)
{
    DevDesc *dd = CurrentDevice();
 
    checkArity(op, args);
    if (dd->newDevStruct)
	GEplaySnapshot(CAR(args), (GEDevDesc*) dd);
    else
	errorcall(call, "can't play snapshot on old-style device");
    return R_NilValue;
}
 
/* I don't think this gets called in any base R code
 */
SEXP do_replay(SEXP call, SEXP op, SEXP args, SEXP env)
{
    if (!NoDevices()) {
	GEDevDesc *dd = GEcurrentDevice();
	checkArity(op, args);
	/*     Rf_dpptr(dd)->resize(); */
	GEplayDisplayList(dd);
    }
    return R_NilValue;
}
 
SEXP do_playDL(SEXP call, SEXP op, SEXP args, SEXP env)
{
    DevDesc *dd = CurrentDevice();
    SEXP theList;
    int ask;
 
    checkArity(op, args);
    if(!isList(theList = CAR(args)))
       errorcall(call, "invalid argument");
    if (dd->newDevStruct)
	((GEDevDesc*) dd)->dev->displayList = theList;
    else
	dd->displayList = theList;
    if (theList != R_NilValue) {
	ask = Rf_gpptr(dd)->ask;
	Rf_gpptr(dd)->ask = 1;
	GReset(dd);
	while (theList != R_NilValue) {
	    SEXP theOperation = CAR(theList);
	    SEXP l_op = CAR(theOperation);
	    SEXP l_args = CDR(theOperation);
	    PRIMFUN(l_op) (R_NilValue, l_op, l_args, R_NilValue);
	    if (!Rf_gpptr(dd)->valid) break;
	    theList = CDR(theList);
	}
	Rf_gpptr(dd)->ask = ask;
    }
    return R_NilValue;
}
 
SEXP do_setGPar(SEXP call, SEXP op, SEXP args, SEXP env)
{
    DevDesc *dd = CurrentDevice();
    int lGPar = 1 + sizeof(GPar) / sizeof(int);
    SEXP GP;
 
    checkArity(op, args);
    GP = CAR(args);
    if (!isInteger(GP) || length(GP) != lGPar)
	errorcall(call, "invalid graphics parameter list");
    copyGPar((GPar *) INTEGER(GP), Rf_dpSavedptr(dd)); /* &dd->Rf_dpSaved); */
    return R_NilValue;
}
 
/* symbols(..) in ../library/base/R/symbols.R  : */
 
/* utility just computing range() */
static Rboolean SymbolRange(double *x, int n, double *xmax, double *xmin)
{
    int i;
    *xmax = -DBL_MAX;
    *xmin =  DBL_MAX;
    for(i = 0; i < n; i++)
	if (R_FINITE(x[i])) {
	    if (*xmax < x[i]) *xmax = x[i];
	    if (*xmin > x[i]) *xmin = x[i];
	}
    return(*xmax >= *xmin && *xmin >= 0);
}
 
static void CheckSymbolPar(SEXP call, SEXP p, int *nr, int *nc)
{
    SEXP dim = getAttrib(p, R_DimSymbol);
    switch(length(dim)) {
    case 0:
	*nr = LENGTH(p);
	*nc = 1;
	break;
    case 1:
	*nr = INTEGER(dim)[0];
	*nc = 1;
	break;
    case 2:
	*nr = INTEGER(dim)[0];
	*nc = INTEGER(dim)[1];
	break;
    default:
	*nr = 0;
	*nc = 0;
    }
    if (*nr == 0 || *nc == 0)
	errorcall(call, "invalid symbol parameter vector");
}
 
/* Internal  symbols(x, y, type, data, inches, bg, fg, ...) */
SEXP do_symbols(SEXP call, SEXP op, SEXP args, SEXP env)
{
    SEXP x, y, p, fg, bg;
    int i, j, nr, nc, nbg, nfg, type;
    double pmax, pmin, inches, rx, ry;
    double xx, yy, p0, p1, p2, p3, p4;
    double *pp, *xp, *yp;
    char *vmax;
 
    SEXP originalArgs = args;
    DevDesc *dd = CurrentDevice();
    GCheckState(dd);
 
    if (length(args) < 7)
	errorcall(call, "insufficient arguments");
 
    PROTECT(x = coerceVector(CAR(args), REALSXP)); args = CDR(args);
    PROTECT(y = coerceVector(CAR(args), REALSXP)); args = CDR(args);
    if (!isNumeric(x) || !isNumeric(y) || length(x) <= 0 || LENGTH(x) <= 0)
	errorcall(call, "invalid symbol coordinates");
 
    type = asInteger(CAR(args)); args = CDR(args);
 
    /* data: */
    p = PROTECT(coerceVector(CAR(args), REALSXP)); args = CDR(args);
    CheckSymbolPar(call, p, &nr, &nc);
    if (LENGTH(x) != nr || LENGTH(y) != nr)
	errorcall(call, "x/y/parameter length mismatch");
 
    inches = asReal(CAR(args)); args = CDR(args);
    if (!R_FINITE(inches) || inches < 0)
	inches = 0;
 
    PROTECT(bg = FixupCol(CAR(args), NA_INTEGER)); args = CDR(args);
    nbg = LENGTH(bg);
 
    PROTECT(fg = FixupCol(CAR(args), NA_INTEGER)); args = CDR(args);
    nfg = LENGTH(fg);
 
    GSavePars(dd);
    ProcessInlinePars(args, dd, call);
 
    GMode(1, dd);
    switch (type) {
    case 1: /* circles */
	if (nc != 1)
	    errorcall(call, "invalid circles data");
	if (!SymbolRange(REAL(p), nr, &pmax, &pmin))
	    errorcall(call, "invalid symbol parameter");
	for (i = 0; i < nr; i++) {
	    if (R_FINITE(REAL(x)[i]) && R_FINITE(REAL(y)[i]) &&
		R_FINITE(REAL(p)[i])) {
		rx = REAL(p)[i];
		if (inches > 0)
		    rx *= inches / pmax;/* <-- no GConvert ?? */
		else/* FIXME : INCHES in the following ??? */
		    rx = GConvertXUnits(rx, USER, INCHES, dd);
		GCircle(REAL(x)[i], REAL(y)[i],	USER, rx,
			INTEGER(bg)[i%nbg], INTEGER(fg)[i%nfg],	dd);
	    }
	}
	break;
    case 2: /* squares */
	if(nc != 1)
	    errorcall(call, "invalid squares data");
	if(!SymbolRange(REAL(p), nr, &pmax, &pmin))
	    errorcall(call, "invalid symbol parameter");
	for (i = 0; i < nr; i++) {
	    if (R_FINITE(REAL(x)[i]) && R_FINITE(REAL(y)[i]) &&
		R_FINITE(REAL(p)[i])) {
		p0 = REAL(p)[i];
		xx = REAL(x)[i];
		yy = REAL(y)[i];
		GConvert(&xx, &yy, USER, DEVICE, dd);
		if (inches > 0) {
		    p0 = p0 / pmax * inches;
		    rx = GConvertXUnits(0.5 * p0, INCHES, DEVICE, dd);
		    ry = GConvertYUnits(0.5 * p0, INCHES, DEVICE, dd);
		}
		else {
		    rx = GConvertXUnits(0.5 * p0, USER, DEVICE, dd);
		    ry = GConvertYUnits(0.5 * p0, USER, DEVICE, dd);
		}
		GRect(xx - rx, yy - rx, xx + rx, yy + rx, DEVICE,
		      INTEGER(bg)[i%nbg], INTEGER(fg)[i%nfg], dd);
	    }
	}
	break;
    case 3: /* rectangles */
	if (nc != 2)
	    errorcall(call, "invalid rectangles data (need 2 columns)");
	if (!SymbolRange(REAL(p), 2 * nr, &pmax, &pmin))
	    errorcall(call, "invalid symbol parameter");
	for (i = 0; i < nr; i++) {
	    if (R_FINITE(REAL(x)[i]) && R_FINITE(REAL(y)[i]) &&
		R_FINITE(REAL(p)[i]) && R_FINITE(REAL(p)[i+nr])) {
		xx = REAL(x)[i];
		yy = REAL(y)[i];
		GConvert(&xx, &yy, USER, DEVICE, dd);
		p0 = REAL(p)[i];
		p1 = REAL(p)[i+nr];
		if (inches > 0) {
		    p0 *= inches / pmax;
		    p1 *= inches / pmax;
		    rx = GConvertXUnits(0.5 * p0, INCHES, DEVICE, dd);
		    ry = GConvertYUnits(0.5 * p1, INCHES, DEVICE, dd);
		}
		else {
		    rx = GConvertXUnits(0.5 * p0, USER, DEVICE, dd);
		    ry = GConvertYUnits(0.5 * p1, USER, DEVICE, dd);
		}
		GRect(xx - rx, yy - ry, xx + rx, yy + ry, DEVICE,
		      INTEGER(bg)[i%nbg], INTEGER(fg)[i%nfg], dd);
 
	    }
	}
	break;
    case 4: /* stars */
	if (nc < 3)
	    errorcall(call, "invalid stars data");
	if (!SymbolRange(REAL(p), nc * nr, &pmax, &pmin))
	    errorcall(call, "invalid symbol parameter");
	vmax = vmaxget();
	pp = (double*)R_alloc(nc, sizeof(double));
	xp = (double*)R_alloc(nc, sizeof(double));
	yp = (double*)R_alloc(nc, sizeof(double));
	p1 = 2.0 * M_PI / nc;
	for (i = 0; i < nr; i++) {
	    xx = REAL(x)[i];
	    yy = REAL(y)[i];
	    if (R_FINITE(xx) && R_FINITE(yy)) {
		GConvert(&xx, &yy, USER, NDC, dd);
		if (inches > 0) {
		    for(j = 0; j < nc; j++) {
			p0 = REAL(p)[i + j * nr];
			if (!R_FINITE(p0)) p0 = 0;
			pp[j] = (p0 / pmax) * inches;
		    }
		}
		else {
		    for(j = 0; j < nc; j++) {
			p0 = REAL(p)[i + j * nr];
			if (!R_FINITE(p0)) p0 = 0;
			pp[j] =	 GConvertXUnits(p0, USER, INCHES, dd);
		    }
		}
		for(j = 0; j < nc; j++) {
		    xp[j] = GConvertXUnits(pp[j] * cos(j * p1),
					   INCHES, NDC, dd) + xx;
		    yp[j] = GConvertYUnits(pp[j] * sin(j * p1),
					   INCHES, NDC, dd) + yy;
		}
		GPolygon(nc, xp, yp, NDC,
			 INTEGER(bg)[i%nbg], INTEGER(fg)[i%nfg], dd);
	    }
	}
	vmaxset(vmax);
	break;
    case 5: /* thermometers */
	if (nc != 3 && nc != 4)
	    errorcall(call, "invalid thermometers data (need 3 or 4 columns)");
	SymbolRange(REAL(p)+2*nr/* <-- pointer arith*/, nr, &pmax, &pmin);
	if (pmax < pmin)
	    errorcall(call, "invalid thermometers[,%s]",(nc == 4)? "3:4" : "3");
	if (pmin < 0. || pmax > 1.) /* S-PLUS has an error here */
	    warningcall(call,"thermometers[,%s] not in [0,1] -- may look funny",
			(nc == 4)? "3:4" : "3");
	if (!SymbolRange(REAL(p), 2 * nr, &pmax, &pmin))
	    errorcall(call, "invalid thermometers[,1:2]");
	for (i = 0; i < nr; i++) {
	    xx = REAL(x)[i];
	    yy = REAL(y)[i];
	    if (R_FINITE(xx) && R_FINITE(yy)) {
		p0 = REAL(p)[i];
		p1 = REAL(p)[i + nr];
		p2 = REAL(p)[i + 2 * nr];
		p3 = (nc == 4)? REAL(p)[i + 3 * nr] : 0.;
		if (R_FINITE(p0) && R_FINITE(p1) &&
		    R_FINITE(p2) && R_FINITE(p3)) {
		    if (p2 < 0) p2 = 0; else if (p2 > 1) p2 = 1;
		    if (p3 < 0) p3 = 0; else if (p3 > 1) p3 = 1;
		    GConvert(&xx, &yy, USER, NDC, dd);
		    if (inches > 0) {
			p0 *= inches / pmax;
			p1 *= inches / pmax;
			rx = GConvertXUnits(0.5 * p0, INCHES, NDC, dd);
			ry = GConvertYUnits(0.5 * p1, INCHES, NDC, dd);
		    }
		    else {
			rx = GConvertXUnits(0.5 * p0, USER, NDC, dd);
			ry = GConvertYUnits(0.5 * p1, USER, NDC, dd);
		    }
		    GRect(xx - rx, yy - ry, xx + rx, yy + ry, NDC,
			  INTEGER(bg)[i%nbg], INTEGER(fg)[i%nfg], dd);
		    GRect(xx - rx,  yy - (1 - 2 * p2) * ry,
			  xx + rx,  yy - (1 - 2 * p3) * ry,
			  NDC,
			  INTEGER(fg)[i%nfg], INTEGER(fg)[i%nfg], dd);
		    GLine(xx - rx, yy, xx - 1.5 * rx, yy, NDC, dd);
		    GLine(xx + rx, yy, xx + 1.5 * rx, yy, NDC, dd);
 
		}
	    }
	}
	break;
    case 6: /* boxplots (wid, hei, loWhsk, upWhsk, medProp) */
	if (nc != 5)
	    errorcall(call, "invalid boxplots data (need 5 columns)");
	pmax = -DBL_MAX;
	pmin =	DBL_MAX;
	for(i = 0; i < nr; i++) {
	    p4 = REAL(p)[i + 4 * nr];	/* median proport. in [0,1] */
	    if (pmax < p4) pmax = p4;
	    if (pmin > p4) pmin = p4;
	}
	if (pmin < 0. || pmax > 1.) /* S-PLUS has an error here */
	    warningcall(call, "boxplots[,5] outside [0,1] -- may look funny");
	if (!SymbolRange(REAL(p), 4 * nr, &pmax, &pmin))
	    errorcall(call, "invalid boxplots[, 1:4]");
	for (i = 0; i < nr; i++) {
	    xx = REAL(x)[i];
	    yy = REAL(y)[i];
	    if (R_FINITE(xx) && R_FINITE(yy)) {
		p0 = REAL(p)[i];	/* width */
		p1 = REAL(p)[i + nr];	/* height */
		p2 = REAL(p)[i + 2 * nr];/* lower whisker */
		p3 = REAL(p)[i + 3 * nr];/* upper whisker */
		p4 = REAL(p)[i + 4 * nr];/* median proport. in [0,1] */
		if (R_FINITE(p0) && R_FINITE(p1) &&
		    R_FINITE(p2) && R_FINITE(p3) && R_FINITE(p4)) {
		    GConvert(&xx, &yy, USER, NDC, dd);
		    if (inches > 0) {
			p0 *= inches / pmax;
			p1 *= inches / pmax;
			p2 *= inches / pmax;
			p3 *= inches / pmax;
			p0 = GConvertXUnits(p0, INCHES, NDC, dd);
			p1 = GConvertYUnits(p1, INCHES, NDC, dd);
			p2 = GConvertYUnits(p2, INCHES, NDC, dd);
			p3 = GConvertYUnits(p3, INCHES, NDC, dd);
		    }
		    else {
			p0 = GConvertXUnits(p0, USER, NDC, dd);
			p1 = GConvertYUnits(p1, USER, NDC, dd);
			p2 = GConvertYUnits(p2, USER, NDC, dd);
			p3 = GConvertYUnits(p3, USER, NDC, dd);
		    }
		    rx = 0.5 * p0;
		    ry = 0.5 * p1;
		    p4 = (1 - p4) * (yy - ry) + p4 * (yy + ry);
		    /* Box */
		    GRect(xx - rx, yy - ry, xx + rx, yy + ry, NDC,
			  INTEGER(bg)[i%nbg], INTEGER(fg)[i%nfg], dd);
		    /* Median */
		    GLine(xx - rx, p4, xx + rx, p4, NDC, dd);
		    /* Lower Whisker */
		    GLine(xx, yy - ry, xx, yy - ry - p2, NDC, dd);
		    /* Upper Whisker */
		    GLine(xx, yy + ry, xx, yy + ry + p3, NDC, dd);
		}
	    }
	}
	break;
    default:
	errorcall(call, "invalid symbol type");
    }
    GMode(0, dd);
    GRestorePars(dd);
    if (GRecording(call))
	recordGraphicOperation(op, originalArgs, dd);
    UNPROTECT(5);
    return R_NilValue;
}

Rev 29864	Rev 30691
1	`/*`	1	`/*`
2	`* R : A Computer Language for Statistical Data Analysis`	2	`* R : A Computer Language for Statistical Data Analysis`
3	`* Copyright (C) 1995, 1996 Robert Gentleman and Ross Ihaka`	3	`* Copyright (C) 1995, 1996 Robert Gentleman and Ross Ihaka`
4	`* Copyright (C) 1997--2001 Robert Gentleman, Ross Ihaka and the`	4	`* Copyright (C) 1997--2001 Robert Gentleman, Ross Ihaka and the`
5	`* R Development Core Team`	5	`* R Development Core Team`
6	`* Copyright (C) 2002--2003 The R Foundation`	6	`* Copyright (C) 2002--2004 The R Foundation`
7	`*`	7	`*`
8	`* This program is free software; you can redistribute it and/or modify`	8	`* This program is free software; you can redistribute it and/or modify`
9	`* it under the terms of the GNU General Public License as published by`	9	`* it under the terms of the GNU General Public License as published by`
10	`* the Free Software Foundation; either version 2 of the License, or`	10	`* the Free Software Foundation; either version 2 of the License, or`
11	`* (at your option) any later version.`	11	`* (at your option) any later version.`
12	`*`	12	`*`
13	`* This program is distributed in the hope that it will be useful,`	13	`* This program is distributed in the hope that it will be useful,`
14	`* but WITHOUT ANY WARRANTY; without even the implied warranty of`	14	`* but WITHOUT ANY WARRANTY; without even the implied warranty of`
15	`* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the`	15	`* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the`
16	`* GNU General Public License for more details.`	16	`* GNU General Public License for more details.`
17	`*`	17	`*`
18	`* A copy of the GNU General Public License is available via WWW at`	18	`* A copy of the GNU General Public License is available via WWW at`
19	`* http://www.gnu.org/copyleft/gpl.html. You can also obtain it by`	19	`* http://www.gnu.org/copyleft/gpl.html. You can also obtain it by`
20	`* writing to the Free Software Foundation, Inc., 59 Temple Place,`	20	`* writing to the Free Software Foundation, Inc., 59 Temple Place,`
21	`* Suite 330, Boston, MA 02111-1307 USA.`	21	`* Suite 330, Boston, MA 02111-1307 USA.`
22	`*/`	22	`*/`
23		23
24	`#ifdef HAVE_CONFIG_H`	24	`#ifdef HAVE_CONFIG_H`
25	`#include <config.h>`	25	`#include <config.h>`
26	`#endif`	26	`#endif`
27		27
28	`#include <Defn.h>`	28	`#include <Defn.h>`
29	`#include <Rmath.h>`	29	`#include <Rmath.h>`
30	`#include <Graphics.h>`	30	`#include <Graphics.h>`
31	`#include <Rdevices.h>`	31	`#include <Rdevices.h>`
32	`#include <Print.h>`	32	`#include <Print.h>`
33		33
34	`#ifndef HAVE_HYPOT`	34	`#ifndef HAVE_HYPOT`
35	`# define hypot pythag`	35	`# define hypot pythag`
36	`#endif`	36	`#endif`
37		37
38		38
39	`void NewFrameConfirm(void)`	39	`void NewFrameConfirm(void)`
40	`{`	40	`{`
41	`unsigned char buf[16];`	41	`unsigned char buf[16];`
42	`R_ReadConsole("Hit <Return> to see next plot: ", buf, 16, 0);`	42	`R_ReadConsole("Hit <Return> to see next plot: ", buf, 16, 0);`
43	`}`	43	`}`
44		44
45	`/* Remember: +1 and/or -1 because C arrays are */`	45	`/* Remember: +1 and/or -1 because C arrays are */`
46	`/* zero-based and R-vectors are one-based. */`	46	`/* zero-based and R-vectors are one-based. */`
47		47
48	`#define checkArity_length \`	48	`#define checkArity_length \`
49	`checkArity(op, args); \`	49	`checkArity(op, args); \`
50	`if(!LENGTH(CAR(args))) \`	50	`if(!LENGTH(CAR(args))) \`
51	`errorcall(call, "argument must have positive length")`	51	`errorcall(call, "argument must have positive length")`
52		52
53		53
54	`SEXP do_devcontrol(SEXP call, SEXP op, SEXP args, SEXP env)`	54	`SEXP do_devcontrol(SEXP call, SEXP op, SEXP args, SEXP env)`
55	`{`	55	`{`
56	`int listFlag;`	56	`int listFlag;`
57		57
58	`checkArity(op, args);`	58	`checkArity(op, args);`
59	`listFlag = asLogical(CAR(args));`	59	`listFlag = asLogical(CAR(args));`
60	`if(listFlag == NA_LOGICAL) errorcall(call, "invalid argument");`	60	`if(listFlag == NA_LOGICAL) errorcall(call, "invalid argument");`
61	`if(listFlag)`	61	`if(listFlag)`
62	`enableDisplayList(CurrentDevice());`	62	`enableDisplayList(CurrentDevice());`
63	`else`	63	`else`
64	`inhibitDisplayList(CurrentDevice());`	64	`inhibitDisplayList(CurrentDevice());`
65	`return ScalarLogical(listFlag);`	65	`return ScalarLogical(listFlag);`
66	`}`	66	`}`
67		67
68	`SEXP do_devcopy(SEXP call, SEXP op, SEXP args, SEXP env)`	68	`SEXP do_devcopy(SEXP call, SEXP op, SEXP args, SEXP env)`
69	`{`	69	`{`
70	`checkArity_length;`	70	`checkArity_length;`
71	`GEcopyDisplayList(INTEGER(CAR(args))[0] - 1);`	71	`GEcopyDisplayList(INTEGER(CAR(args))[0] - 1);`
72	`return R_NilValue;`	72	`return R_NilValue;`
73	`}`	73	`}`
74		74
75	`SEXP do_devcur(SEXP call, SEXP op, SEXP args, SEXP env)`	75	`SEXP do_devcur(SEXP call, SEXP op, SEXP args, SEXP env)`
76	`{`	76	`{`
77	`SEXP cd = allocVector(INTSXP, 1);`	77	`SEXP cd = allocVector(INTSXP, 1);`
78	`checkArity(op, args);`	78	`checkArity(op, args);`
79	`INTEGER(cd)[0] = curDevice() + 1;`	79	`INTEGER(cd)[0] = curDevice() + 1;`
80	`return cd;`	80	`return cd;`
81	`}`	81	`}`
82		82
83	`SEXP do_devnext(SEXP call, SEXP op, SEXP args, SEXP env)`	83	`SEXP do_devnext(SEXP call, SEXP op, SEXP args, SEXP env)`
84	`{`	84	`{`
85	`SEXP nd = allocVector(INTSXP, 1);`	85	`SEXP nd = allocVector(INTSXP, 1);`
86	`checkArity_length;`	86	`checkArity_length;`
87	`INTEGER(nd)[0] = nextDevice(INTEGER(CAR(args))[0] - 1) + 1;`	87	`INTEGER(nd)[0] = nextDevice(INTEGER(CAR(args))[0] - 1) + 1;`
88	`return nd;`	88	`return nd;`
89	`}`	89	`}`
90		90
91	`SEXP do_devprev(SEXP call, SEXP op, SEXP args, SEXP env)`	91	`SEXP do_devprev(SEXP call, SEXP op, SEXP args, SEXP env)`
92	`{`	92	`{`
93	`SEXP pd = allocVector(INTSXP, 1);`	93	`SEXP pd = allocVector(INTSXP, 1);`
94	`checkArity_length;`	94	`checkArity_length;`
95	`INTEGER(pd)[0] = prevDevice(INTEGER(CAR(args))[0] - 1) + 1;`	95	`INTEGER(pd)[0] = prevDevice(INTEGER(CAR(args))[0] - 1) + 1;`
96	`return pd;`	96	`return pd;`
97	`}`	97	`}`
98		98
99	`SEXP do_devset(SEXP call, SEXP op, SEXP args, SEXP env)`	99	`SEXP do_devset(SEXP call, SEXP op, SEXP args, SEXP env)`
100	`{`	100	`{`
101	`int devNum = INTEGER(CAR(args))[0] - 1;`	101	`int devNum = INTEGER(CAR(args))[0] - 1;`
102	`SEXP sd = allocVector(INTSXP, 1);`	102	`SEXP sd = allocVector(INTSXP, 1);`
103	`checkArity(op, args);`	103	`checkArity(op, args);`
104	`INTEGER(sd)[0] = selectDevice(devNum) + 1;`	104	`INTEGER(sd)[0] = selectDevice(devNum) + 1;`
105	`return sd;`	105	`return sd;`
106	`}`	106	`}`
107		107
108	`SEXP do_devoff(SEXP call, SEXP op, SEXP args, SEXP env)`	108	`SEXP do_devoff(SEXP call, SEXP op, SEXP args, SEXP env)`
109	`{`	109	`{`
110	`checkArity_length;`	110	`checkArity_length;`
111	`killDevice(INTEGER(CAR(args))[0] - 1);`	111	`killDevice(INTEGER(CAR(args))[0] - 1);`
112	`return R_NilValue;`	112	`return R_NilValue;`
113	`}`	113	`}`
114		114
115		115
116	`/* P A R A M E T E R U T I L I T I E S */`	116	`/* P A R A M E T E R U T I L I T I E S */`
117		117
118		118
119	`/* ProcessInLinePars handles inline par specifications in graphics functions.`	119	`/* ProcessInLinePars handles inline par specifications in graphics functions.`
120	`* It does this by calling Specify2() from ./par.c */`	120	`* It does this by calling Specify2() from ./par.c */`
121		121
122	`void ProcessInlinePars(SEXP s, DevDesc *dd, SEXP call)`	122	`void ProcessInlinePars(SEXP s, DevDesc *dd, SEXP call)`
123	`{`	123	`{`
124	`if (isList(s)) {`	124	`if (isList(s)) {`
125	`while (s != R_NilValue) {`	125	`while (s != R_NilValue) {`
126	`if (isList(CAR(s)))`	126	`if (isList(CAR(s)))`
127	`ProcessInlinePars(CAR(s), dd, call);`	127	`ProcessInlinePars(CAR(s), dd, call);`
128	`else if (TAG(s) != R_NilValue)`	128	`else if (TAG(s) != R_NilValue)`
129	`Specify2(CHAR(PRINTNAME(TAG(s))), CAR(s), dd, call);`	129	`Specify2(CHAR(PRINTNAME(TAG(s))), CAR(s), dd, call);`
130	`s = CDR(s);`	130	`s = CDR(s);`
131	`}`	131	`}`
132	`}`	132	`}`
133	`}`	133	`}`
134		134
-		135	`/*`
-		136	`* Extract specified par from list of inline pars`
-		137	`*/`
-		138	`SEXP getInlinePar(SEXP s, char *name)`
-		139	`{`
-		140	`SEXP result = R_NilValue;`
-		141	`int found = 0;`
-		142	`if (isList(s) && !found) {`
-		143	`while (s != R_NilValue) {`
-		144	`if (isList(CAR(s))) {`
-		145	`result = getInlinePar(CAR(s), name);`
-		146	`if (result)`
-		147	`found = 1;`
-		148	`} else`
-		149	`if (TAG(s) != R_NilValue)`
-		150	`if (!strcmp(CHAR(PRINTNAME(TAG(s))), name)) {`
-		151	`result = CAR(s);`
-		152	`found = 1;`
-		153	`}`
-		154	`s = CDR(s);`
-		155	`}`
-		156	`}`
-		157	`return result;`
-		158	`}`
-		159
135	`SEXP FixupPch(SEXP pch, int dflt)`	160	`SEXP FixupPch(SEXP pch, int dflt)`
136	`{`	161	`{`
137	`int i, n;`	162	`int i, n;`
138	`SEXP ans = R_NilValue;/* -Wall*/`	163	`SEXP ans = R_NilValue;/* -Wall*/`
139		164
140	`n = length(pch);`	165	`n = length(pch);`
141	`if (n == 0) {`	166	`if (n == 0) {`
142	`ans = allocVector(INTSXP, 1);`	167	`ans = allocVector(INTSXP, 1);`
143	`INTEGER(ans)[0] = dflt;`	168	`INTEGER(ans)[0] = dflt;`
144	`}`	169	`}`
145	`else if (isList(pch)) {`	170	`else if (isList(pch)) {`
146	`ans = allocVector(INTSXP, n);`	171	`ans = allocVector(INTSXP, n);`
147	`for (i = 0; pch != R_NilValue; pch = CDR(pch))`	172	`for (i = 0; pch != R_NilValue; pch = CDR(pch))`
148	`INTEGER(ans)[i++] = asInteger(CAR(pch));`	173	`INTEGER(ans)[i++] = asInteger(CAR(pch));`
149	`}`	174	`}`
150	`else if (isInteger(pch)) {`	175	`else if (isInteger(pch)) {`
151	`ans = allocVector(INTSXP, n);`	176	`ans = allocVector(INTSXP, n);`
152	`for (i = 0; i < n; i++)`	177	`for (i = 0; i < n; i++)`
153	`INTEGER(ans)[i] = INTEGER(pch)[i];`	178	`INTEGER(ans)[i] = INTEGER(pch)[i];`
154	`}`	179	`}`
155	`else if (isReal(pch)) {`	180	`else if (isReal(pch)) {`
156	`ans = allocVector(INTSXP, n);`	181	`ans = allocVector(INTSXP, n);`
157	`for (i = 0; i < n; i++)`	182	`for (i = 0; i < n; i++)`
158	`INTEGER(ans)[i] = R_FINITE(REAL(pch)[i]) ?`	183	`INTEGER(ans)[i] = R_FINITE(REAL(pch)[i]) ?`
159	`REAL(pch)[i] : NA_INTEGER;`	184	`REAL(pch)[i] : NA_INTEGER;`
160	`}`	185	`}`
161	`else if (isString(pch)) {`	186	`else if (isString(pch)) {`
162	`ans = allocVector(INTSXP, n);`	187	`ans = allocVector(INTSXP, n);`
163	`for (i = 0; i < n; i++)`	188	`for (i = 0; i < n; i++)`
164	`INTEGER(ans)[i] = STRING_ELT(pch, i) != NA_STRING ?`	189	`INTEGER(ans)[i] = STRING_ELT(pch, i) != NA_STRING ?`
165	`CHAR(STRING_ELT(pch, i))[0] : NA_INTEGER;`	190	`CHAR(STRING_ELT(pch, i))[0] : NA_INTEGER;`
166	`}`	191	`}`
167	`else if (isLogical(pch)) {/* NA, but not TRUE/FALSE */`	192	`else if (isLogical(pch)) {/* NA, but not TRUE/FALSE */`
168	`ans = allocVector(INTSXP, n);`	193	`ans = allocVector(INTSXP, n);`
169	`for (i = 0; i < n; i++)`	194	`for (i = 0; i < n; i++)`
170	`if(LOGICAL(pch)[i] == NA_LOGICAL)`	195	`if(LOGICAL(pch)[i] == NA_LOGICAL)`
171	`INTEGER(ans)[i] = NA_INTEGER;`	196	`INTEGER(ans)[i] = NA_INTEGER;`
172	`else error("only NA allowed in logical plotting symbol");`	197	`else error("only NA allowed in logical plotting symbol");`
173	`}`	198	`}`
174	`else error("invalid plotting symbol");`	199	`else error("invalid plotting symbol");`
175	`for (i = 0; i < n; i++) {`	200	`for (i = 0; i < n; i++) {`
176	`if (INTEGER(ans)[i] < 0 && INTEGER(ans)[i] != NA_INTEGER)`	201	`if (INTEGER(ans)[i] < 0 && INTEGER(ans)[i] != NA_INTEGER)`
177	`INTEGER(ans)[i] = dflt;`	202	`INTEGER(ans)[i] = dflt;`
178	`}`	203	`}`
179	`return ans;`	204	`return ans;`
180	`}`	205	`}`
181		206
182	`SEXP FixupLty(SEXP lty, int dflt)`	207	`SEXP FixupLty(SEXP lty, int dflt)`
183	`{`	208	`{`
184	`int i, n;`	209	`int i, n;`
185	`SEXP ans;`	210	`SEXP ans;`
186	`n = length(lty);`	211	`n = length(lty);`
187	`if (n == 0) {`	212	`if (n == 0) {`
188	`ans = allocVector(INTSXP, 1);`	213	`ans = allocVector(INTSXP, 1);`
189	`INTEGER(ans)[0] = dflt;`	214	`INTEGER(ans)[0] = dflt;`
190	`}`	215	`}`
191	`else {`	216	`else {`
192	`ans = allocVector(INTSXP, n);`	217	`ans = allocVector(INTSXP, n);`
193	`for (i = 0; i < n; i++)`	218	`for (i = 0; i < n; i++)`
194	`INTEGER(ans)[i] = LTYpar(lty, i);`	219	`INTEGER(ans)[i] = LTYpar(lty, i);`
195	`}`	220	`}`
196	`return ans;`	221	`return ans;`
197	`}`	222	`}`
198		223
199	`SEXP FixupLwd(SEXP lwd, double dflt)`	224	`SEXP FixupLwd(SEXP lwd, double dflt)`
200	`{`	225	`{`
201	`int i, n;`	226	`int i, n;`
202	`double w;`	227	`double w;`
203	`SEXP ans = NULL;`	228	`SEXP ans = NULL;`
204		229
205	`n = length(lwd);`	230	`n = length(lwd);`
206	`if (n == 0) {`	231	`if (n == 0) {`
207	`ans = allocVector(REALSXP, 1);`	232	`ans = allocVector(REALSXP, 1);`
208	`REAL(ans)[0] = dflt;`	233	`REAL(ans)[0] = dflt;`
209	`}`	234	`}`
210	`else {`	235	`else {`
211	`PROTECT(lwd = coerceVector(lwd, REALSXP));`	236	`PROTECT(lwd = coerceVector(lwd, REALSXP));`
212	`n = length(lwd);`	237	`n = length(lwd);`
213	`ans = allocVector(REALSXP, n);`	238	`ans = allocVector(REALSXP, n);`
214	`for (i = 0; i < n; i++) {`	239	`for (i = 0; i < n; i++) {`
215	`w = REAL(lwd)[i];`	240	`w = REAL(lwd)[i];`
216	`if (w < 0) w = NA_REAL;`	241	`if (w < 0) w = NA_REAL;`
217	`REAL(ans)[i] = w;`	242	`REAL(ans)[i] = w;`
218		243
219	`}`	244	`}`
220	`UNPROTECT(1);`	245	`UNPROTECT(1);`
221	`}`	246	`}`
222	`return ans;`	247	`return ans;`
223	`}`	248	`}`
224		249
225	`SEXP FixupFont(SEXP font, int dflt)`	250	`SEXP FixupFont(SEXP font, int dflt)`
226	`{`	251	`{`
227	`int i, k, n;`	252	`int i, k, n;`
228	`SEXP ans = R_NilValue;/* -Wall*/`	253	`SEXP ans = R_NilValue;/* -Wall*/`
229	`n = length(font);`	254	`n = length(font);`
230	`if (n == 0) {`	255	`if (n == 0) {`
231	`ans = allocVector(INTSXP, 1);`	256	`ans = allocVector(INTSXP, 1);`
232	`INTEGER(ans)[0] = dflt;`	257	`INTEGER(ans)[0] = dflt;`
233	`}`	258	`}`
234	`else if (isInteger(font) \|\| isLogical(font)) {`	259	`else if (isInteger(font) \|\| isLogical(font)) {`
235	`ans = allocVector(INTSXP, n);`	260	`ans = allocVector(INTSXP, n);`
236	`for (i = 0; i < n; i++) {`	261	`for (i = 0; i < n; i++) {`
237	`k = INTEGER(font)[i];`	262	`k = INTEGER(font)[i];`
238	`#ifndef Win32`	263	`#ifndef Win32`
239	`if (k < 1 \|\| k > 4) k = NA_INTEGER;`	264	`if (k < 1 \|\| k > 4) k = NA_INTEGER;`
240	`#else`	265	`#else`
241	`if (k < 1 \|\| k > 32) k = NA_INTEGER;`	266	`if (k < 1 \|\| k > 32) k = NA_INTEGER;`
242	`#endif`	267	`#endif`
243	`INTEGER(ans)[i] = k;`	268	`INTEGER(ans)[i] = k;`
244	`}`	269	`}`
245	`}`	270	`}`
246	`else if (isReal(font)) {`	271	`else if (isReal(font)) {`
247	`ans = allocVector(INTSXP, n);`	272	`ans = allocVector(INTSXP, n);`
248	`for (i = 0; i < n; i++) {`	273	`for (i = 0; i < n; i++) {`
249	`k = REAL(font)[i];`	274	`k = REAL(font)[i];`
250	`#ifndef Win32`	275	`#ifndef Win32`
251	`if (k < 1 \|\| k > 4) k = NA_INTEGER;`	276	`if (k < 1 \|\| k > 4) k = NA_INTEGER;`
252	`#else`	277	`#else`
253	`if (k < 1 \|\| k > 32) k = NA_INTEGER;`	278	`if (k < 1 \|\| k > 32) k = NA_INTEGER;`
254	`#endif`	279	`#endif`
255	`INTEGER(ans)[i] = k;`	280	`INTEGER(ans)[i] = k;`
256	`}`	281	`}`
257	`}`	282	`}`
258	`else error("invalid font specification");`	283	`else error("invalid font specification");`
259	`return ans;`	284	`return ans;`
260	`}`	285	`}`
261		286
262	`SEXP FixupCol(SEXP col, unsigned int dflt)`	287	`SEXP FixupCol(SEXP col, unsigned int dflt)`
263	`{`	288	`{`
264	`int i, n;`	289	`int i, n;`
265	`SEXP ans;`	290	`SEXP ans;`
266	`n = length(col);`	291	`n = length(col);`
267	`if (n == 0) {`	292	`if (n == 0) {`
268	`ans = allocVector(INTSXP, 1);`	293	`ans = allocVector(INTSXP, 1);`
269	`INTEGER(ans)[0] = dflt;`	294	`INTEGER(ans)[0] = dflt;`
270	`}`	295	`}`
271	`else {`	296	`else {`
272	`ans = allocVector(INTSXP, n);`	297	`ans = allocVector(INTSXP, n);`
273	`if (isList(col))`	298	`if (isList(col))`
274	`for (i = 0; i < n; i++) {`	299	`for (i = 0; i < n; i++) {`
275	`INTEGER(ans)[i] = RGBpar(CAR(col), 0);`	300	`INTEGER(ans)[i] = RGBpar(CAR(col), 0);`
276	`col = CDR(col);`	301	`col = CDR(col);`
277	`}`	302	`}`
278	`else`	303	`else`
279	`for (i = 0; i < n; i++)`	304	`for (i = 0; i < n; i++)`
280	`INTEGER(ans)[i] = RGBpar(col, i);`	305	`INTEGER(ans)[i] = RGBpar(col, i);`
281	`}`	306	`}`
282	`return ans;`	307	`return ans;`
283	`}`	308	`}`
284		309
285	`SEXP FixupCex(SEXP cex, double dflt)`	310	`SEXP FixupCex(SEXP cex, double dflt)`
286	`{`	311	`{`
287	`SEXP ans;`	312	`SEXP ans;`
288	`int i, n;`	313	`int i, n;`
289	`n = length(cex);`	314	`n = length(cex);`
290	`if (n == 0) {`	315	`if (n == 0) {`
291	`ans = allocVector(REALSXP, 1);`	316	`ans = allocVector(REALSXP, 1);`
292	`if (R_FINITE(dflt) && dflt > 0)`	317	`if (R_FINITE(dflt) && dflt > 0)`
293	`REAL(ans)[0] = dflt;`	318	`REAL(ans)[0] = dflt;`
294	`else`	319	`else`
295	`REAL(ans)[0] = NA_REAL;`	320	`REAL(ans)[0] = NA_REAL;`
296	`}`	321	`}`
297	`else {`	322	`else {`
298	`double c;`	323	`double c;`
299	`ans = allocVector(REALSXP, n);`	324	`ans = allocVector(REALSXP, n);`
300	`if (isReal(cex))`	325	`if (isReal(cex))`
301	`for (i = 0; i < n; i++) {`	326	`for (i = 0; i < n; i++) {`
302	`c = REAL(cex)[i];`	327	`c = REAL(cex)[i];`
303	`if (R_FINITE(c) && c > 0)`	328	`if (R_FINITE(c) && c > 0)`
304	`REAL(ans)[i] = c;`	329	`REAL(ans)[i] = c;`
305	`else`	330	`else`
306	`REAL(ans)[i] = NA_REAL;`	331	`REAL(ans)[i] = NA_REAL;`
307	`}`	332	`}`
308	`else if (isInteger(cex) \|\| isLogical(cex))`	333	`else if (isInteger(cex) \|\| isLogical(cex))`
309		334
310	`for (i = 0; i < n; i++) {`	335	`for (i = 0; i < n; i++) {`
311	`c = INTEGER(cex)[i];`	336	`c = INTEGER(cex)[i];`
312	`if (c == NA_INTEGER \|\| c <= 0)`	337	`if (c == NA_INTEGER \|\| c <= 0)`
313	`c = NA_REAL;`	338	`c = NA_REAL;`
314	`REAL(ans)[i] = c;`	339	`REAL(ans)[i] = c;`
315	`}`	340	`}`
316	`}`	341	`}`
317	`return ans;`	342	`return ans;`
318	`}`	343	`}`
319		344
320	`SEXP FixupVFont(SEXP vfont) {`	345	`SEXP FixupVFont(SEXP vfont) {`
321	`SEXP ans = R_NilValue;`	346	`SEXP ans = R_NilValue;`
322	`if (!isNull(vfont)) {`	347	`if (!isNull(vfont)) {`
323	`SEXP vf;`	348	`SEXP vf;`
324	`int typeface, fontindex;`	349	`int typeface, fontindex;`
325	`int minindex, maxindex=0;/* -Wall*/`	350	`int minindex, maxindex=0;/* -Wall*/`
326	`int i;`	351	`int i;`
327	`PROTECT(vf = coerceVector(vfont, INTSXP));`	352	`PROTECT(vf = coerceVector(vfont, INTSXP));`
328	`if (length(vf) != 2)`	353	`if (length(vf) != 2)`
329	`error("Invalid vfont value");`	354	`error("Invalid vfont value");`
330	`typeface = INTEGER(vf)[0];`	355	`typeface = INTEGER(vf)[0];`
331	`if (typeface < 0 \|\| typeface > 7)`	356	`if (typeface < 0 \|\| typeface > 7)`
332	`error("Invalid vfont value [typeface]");`	357	`error("Invalid vfont value [typeface]");`
333	`/* For each of the typefaces {0..7}, there are several fontindices`	358	`/* For each of the typefaces {0..7}, there are several fontindices`
334	`available; how many depends on the typeface.`	359	`available; how many depends on the typeface.`
335	`The possible combinations are "given" in ./g_fontdb.c`	360	`The possible combinations are "given" in ./g_fontdb.c`
336	`and also listed in help(Hershey).`	361	`and also listed in help(Hershey).`
337	`*/`	362	`*/`
338	`minindex = 1;`	363	`minindex = 1;`
339	`switch (typeface) {`	364	`switch (typeface) {`
340	`case 0: /* serif */`	365	`case 0: /* serif */`
341	`maxindex = 7; break;`	366	`maxindex = 7; break;`
342	`case 1: /* sans serif */`	367	`case 1: /* sans serif */`
343	`case 6: /* serif symbol */`	368	`case 6: /* serif symbol */`
344	`maxindex = 4; break;`	369	`maxindex = 4; break;`
345	`case 2: /* script */`	370	`case 2: /* script */`
346	`maxindex = 3; break;`	371	`maxindex = 3; break;`
347	`case 3: /* gothic english */`	372	`case 3: /* gothic english */`
348	`case 4: /* gothic german */`	373	`case 4: /* gothic german */`
349	`case 5: /* gothic italian */`	374	`case 5: /* gothic italian */`
350	`maxindex = 1; break;`	375	`maxindex = 1; break;`
351	`case 7: /* sans serif symbol */`	376	`case 7: /* sans serif symbol */`
352	`maxindex = 2;`	377	`maxindex = 2;`
353	`}`	378	`}`
354	`fontindex = INTEGER(vf)[1];`	379	`fontindex = INTEGER(vf)[1];`
355	`if (fontindex < minindex \|\| fontindex > maxindex)`	380	`if (fontindex < minindex \|\| fontindex > maxindex)`
356	`error("Invalid vfont value [fontindex]");`	381	`error("Invalid vfont value [fontindex]");`
357	`ans = allocVector(INTSXP, 2);`	382	`ans = allocVector(INTSXP, 2);`
358	`for (i=0; i<2; i++)`	383	`for (i=0; i<2; i++)`
359	`INTEGER(ans)[i] = INTEGER(vf)[i];`	384	`INTEGER(ans)[i] = INTEGER(vf)[i];`
360	`UNPROTECT(1);`	385	`UNPROTECT(1);`
361	`}`	386	`}`
362	`return ans;`	387	`return ans;`
363	`}`	388	`}`
364		389
365	`/* GetTextArg() : extract from call and possibly set text arguments`	390	`/* GetTextArg() : extract from call and possibly set text arguments`
366	`* ("label", col=, cex=, font=, vfont=)`	391	`* ("label", col=, cex=, font=, vfont=)`
367	`*`	392	`*`
368	`* Main purpose: Treat things like title(main = list("This Title", font= 4))`	393	`* Main purpose: Treat things like title(main = list("This Title", font= 4))`
369	`*`	394	`*`
370	`* Called from do_title() [only, currently]`	395	`* Called from do_title() [only, currently]`
371	`*/`	396	`*/`
372	`static void`	397	`static void`
373	`GetTextArg(SEXP call, SEXP spec, SEXP *ptxt,`	398	`GetTextArg(SEXP call, SEXP spec, SEXP *ptxt,`
374	`int pcol, double pcex, int pfont, SEXP pvfont)`	399	`int pcol, double pcex, int pfont, SEXP pvfont)`
375	`{`	400	`{`
376	`int i, n, col, font;`	401	`int i, n, col, font, colspecd;`
377	`double cex;`	402	`double cex;`
378	`SEXP txt, vfont, nms;`	403	`SEXP txt, vfont, nms;`
379		404
380	`txt = R_NilValue;`	405	`txt = R_NilValue;`
381	`vfont = R_NilValue;`	406	`vfont = R_NilValue;`
382	`cex = NA_REAL;`	407	`cex = NA_REAL;`
383	`col = NA_INTEGER;`	408	`col = NA_INTEGER;`
-		409	`colspecd = 0;`
384	`font = NA_INTEGER;`	410	`font = NA_INTEGER;`
385	`PROTECT(txt);`	411	`PROTECT(txt);`
386		412
387	`switch (TYPEOF(spec)) {`	413	`switch (TYPEOF(spec)) {`
388	`case LANGSXP:`	414	`case LANGSXP:`
389	`case SYMSXP:`	415	`case SYMSXP:`
390	`UNPROTECT(1);`	416	`UNPROTECT(1);`
391	`PROTECT(txt = coerceVector(spec, EXPRSXP));`	417	`PROTECT(txt = coerceVector(spec, EXPRSXP));`
392	`break;`	418	`break;`
393	`case VECSXP:`	419	`case VECSXP:`
394	`if (length(spec) == 0) {`	420	`if (length(spec) == 0) {`
395	`*ptxt = R_NilValue;`	421	`*ptxt = R_NilValue;`
396	`}`	422	`}`
397	`else {`	423	`else {`
398	`nms = getAttrib(spec, R_NamesSymbol);`	424	`nms = getAttrib(spec, R_NamesSymbol);`
399	`if (nms==R_NilValue){ /* PR#1939 */`	425	`if (nms==R_NilValue){ /* PR#1939 */`
400	`txt = VECTOR_ELT(spec, 0);`	426	`txt = VECTOR_ELT(spec, 0);`
401	`if (TYPEOF(txt) == LANGSXP \|\|TYPEOF(txt) == SYMSXP ) {`	427	`if (TYPEOF(txt) == LANGSXP \|\|TYPEOF(txt) == SYMSXP ) {`
402	`UNPROTECT(1);`	428	`UNPROTECT(1);`
403	`PROTECT(txt = coerceVector(txt, EXPRSXP));`	429	`PROTECT(txt = coerceVector(txt, EXPRSXP));`
404	`}`	430	`}`
405	`else if (!isExpression(txt)) {`	431	`else if (!isExpression(txt)) {`
406	`UNPROTECT(1);`	432	`UNPROTECT(1);`
407	`PROTECT(txt = coerceVector(txt, STRSXP));`	433	`PROTECT(txt = coerceVector(txt, STRSXP));`
408	`}`	434	`}`
409	`} else {`	435	`} else {`
410	`n = length(nms);`	436	`n = length(nms);`
411	`for (i = 0; i < n; i++) {`	437	`for (i = 0; i < n; i++) {`
412	`if (!strcmp(CHAR(STRING_ELT(nms, i)), "cex")) {`	438	`if (!strcmp(CHAR(STRING_ELT(nms, i)), "cex")) {`
413	`cex = asReal(VECTOR_ELT(spec, i));`	439	`cex = asReal(VECTOR_ELT(spec, i));`
414	`}`	440	`}`
415	`else if (!strcmp(CHAR(STRING_ELT(nms, i)), "col")) {`	441	`else if (!strcmp(CHAR(STRING_ELT(nms, i)), "col")) {`
-		442	`SEXP colsxp = VECTOR_ELT(spec, i);`
-		443	`if (!isNAcol(colsxp, 0, LENGTH(colsxp))) {`
416	`col = asInteger(FixupCol(VECTOR_ELT(spec, i), NA_INTEGER));`	444	`col = asInteger(FixupCol(colsxp, NA_INTEGER));`
-		445	`colspecd = 1;`
-		446	`}`
417	`}`	447	`}`
418	`else if (!strcmp(CHAR(STRING_ELT(nms, i)), "font")) {`	448	`else if (!strcmp(CHAR(STRING_ELT(nms, i)), "font")) {`
419	`font = asInteger(FixupFont(VECTOR_ELT(spec, i), NA_INTEGER));`	449	`font = asInteger(FixupFont(VECTOR_ELT(spec, i), NA_INTEGER));`
420	`}`	450	`}`
421	`else if (!strcmp(CHAR(STRING_ELT(nms, i)), "vfont")) {`	451	`else if (!strcmp(CHAR(STRING_ELT(nms, i)), "vfont")) {`
422	`vfont = FixupVFont(VECTOR_ELT(spec, i));`	452	`vfont = FixupVFont(VECTOR_ELT(spec, i));`
423	`}`	453	`}`
424	`else if (!strcmp(CHAR(STRING_ELT(nms, i)), "")) {`	454	`else if (!strcmp(CHAR(STRING_ELT(nms, i)), "")) {`
425	`txt = VECTOR_ELT(spec, i);`	455	`txt = VECTOR_ELT(spec, i);`
426	`if (TYPEOF(txt) == LANGSXP \|\| TYPEOF(txt)==SYMSXP) {`	456	`if (TYPEOF(txt) == LANGSXP \|\| TYPEOF(txt)==SYMSXP) {`
427	`UNPROTECT(1);`	457	`UNPROTECT(1);`
428	`PROTECT(txt = coerceVector(txt, EXPRSXP));`	458	`PROTECT(txt = coerceVector(txt, EXPRSXP));`
429	`}`	459	`}`
430	`else if (!isExpression(txt)) {`	460	`else if (!isExpression(txt)) {`
431	`UNPROTECT(1);`	461	`UNPROTECT(1);`
432	`PROTECT(txt = coerceVector(txt, STRSXP));`	462	`PROTECT(txt = coerceVector(txt, STRSXP));`
433	`}`	463	`}`
434	`}`	464	`}`
435	`else errorcall(call, "invalid graphics parameter");`	465	`else errorcall(call, "invalid graphics parameter");`
436	`}`	466	`}`
437	`}`	467	`}`
438	`}`	468	`}`
439	`break;`	469	`break;`
440	`case STRSXP:`	470	`case STRSXP:`
441	`case EXPRSXP:`	471	`case EXPRSXP:`
442	`txt = spec;`	472	`txt = spec;`
443	`break;`	473	`break;`
444	`default:`	474	`default:`
445	`txt = coerceVector(spec, STRSXP);`	475	`txt = coerceVector(spec, STRSXP);`
446	`break;`	476	`break;`
447	`}`	477	`}`
448	`UNPROTECT(1);`	478	`UNPROTECT(1);`
449	`if (txt != R_NilValue) {`	479	`if (txt != R_NilValue) {`
450	`*ptxt = txt;`	480	`*ptxt = txt;`
451	`if (R_FINITE(cex)) *pcex = cex;`	481	`if (R_FINITE(cex)) *pcex = cex;`
452	`if (col != NA_INTEGER) *pcol = col;`	482	`if (colspecd) *pcol = col;`
453	`if (font != NA_INTEGER) *pfont = font;`	483	`if (font != NA_INTEGER) *pfont = font;`
454	`if (vfont != R_NilValue) *pvfont = vfont;`	484	`if (vfont != R_NilValue) *pvfont = vfont;`
455	`}`	485	`}`
456	`}/* GetTextArg */`	486	`}/* GetTextArg */`
457		487
458		488
459	`/* GRAPHICS FUNCTION ENTRY POINTS */`	489	`/* GRAPHICS FUNCTION ENTRY POINTS */`
460		490
461		491
462	`SEXP do_plot_new(SEXP call, SEXP op, SEXP args, SEXP env)`	492	`SEXP do_plot_new(SEXP call, SEXP op, SEXP args, SEXP env)`
463	`{`	493	`{`
464	`/* plot.new() - create a new plot "frame" */`	494	`/* plot.new() - create a new plot "frame" */`
465		495
466	`DevDesc *dd;`	496	`DevDesc *dd;`
467		497
468	`checkArity(op, args);`	498	`checkArity(op, args);`
469		499
470	`dd = GNewPlot(GRecording(call));`	500	`dd = GNewPlot(GRecording(call));`

The R Project SVN R

(root)/branches/djm-source/src/main/plot.c @ 34412 – Rev 29864 → 30691