The R Project SVN – R – Diff – /branches/R-exp-R5/src/main/plotmath.c

 #ifdef max
 #undef max
 #endif
 #ifdef Unix
+static DevDesc *mathDevice;
 		/* return maximum of two doubles */
 static double max(double x, double y)
+{
 	if (x > y) return x;
 	/* initialisation code for mathematical notation */
 static double ratioScale = 0.8;
 static double scriptScale = 0.65;
 static int ratioDepth = 0;
-static int metricUnit = 3;
+static int metricUnit = INCHES;
 static SEXP plusSymbol;
 static SEXP minusSymbol;
 static SEXP timesSymbol;
 static SEXP divideSymbol;
 static float OperatorSpace[] = { 0.1, 0.15, 0.2, 0.6, 0.1 };
 static double fontHeight()
+{
 	double height, depth, width;
-	GMetricInfo(0, &height, &depth, &width, metricUnit);
+	GMetricInfo(0, &height, &depth, &width, metricUnit, mathDevice);
 	return height + depth;
+}
 static double xHeight()
+{
 	int x = 'x';
 	double xheight, depth, width;
-	GMetricInfo(x, &xheight, &depth, &width, metricUnit);
+	GMetricInfo(x, &xheight, &depth, &width, metricUnit, mathDevice);
 	return xheight;
+}
 static double axisHeight()
+{
 	int plus = '+';
 	double plusHeight, depth, width;
-	GMetricInfo(plus, &plusHeight, &depth, &width, metricUnit);
+	GMetricInfo(plus, &plusHeight, &depth, &width, metricUnit, mathDevice);
 	return 0.5 * plusHeight;
+}
 static double superscriptDrop()
+{
 static double radicalExWidth()
+{
 	int radicalEx = 96;
 	double height, depth, REWidth;
-	GMetricInfo(radicalEx, &height, &depth, &REWidth, metricUnit);
+	GMetricInfo(radicalEx, &height, &depth, &REWidth, metricUnit, mathDevice);
 	return REWidth;
+}
 static double superscriptShift(SEXP body, SEXP sup)
+{
 	BBOX bodyBBox = elementBBox(body);
 	BBOX superscriptBBox;
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
 	double temp1 = bboxHeight(bodyBBox) - superscriptDrop();
 	double temp2 = superscript();
 	double temp3;
-	GP->cex = GP->cex * scriptScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
 	superscriptBBox = elementBBox(sup);
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 	temp3 = bboxDepth(superscriptBBox) + 0.25 * xHeight();
 	return max(temp1, max(temp2, temp3));
+}
 static double subscriptShift(SEXP body, SEXP sub, int subOnly)
+{
 	BBOX bodyBBox = elementBBox(body);
 	BBOX subscriptBBox;
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
 	double temp1 = bboxDepth(bodyBBox) + subscriptDrop();
 	double temp2 = subscript();
 	double temp3;
-	GP->cex = GP->cex * scriptScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
 	subscriptBBox = elementBBox(sub);
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 	temp3 = bboxHeight(subscriptBBox) - (4 * xHeight() / 5);
 	if (subOnly)
 		return max(temp1, max(temp2, temp3));
 	else
 static void supsubShift(SEXP body, SEXP sup, SEXP sub,
 			double *supShift, double *subShift)
+{
 	BBOX superscriptBBox, subscriptBBox;
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
 	double temp1, temp2;
-	GP->cex = GP->cex * scriptScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
 	superscriptBBox = elementBBox(sup);
 	subscriptBBox = elementBBox(sub);
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 	*supShift = superscriptShift(body, sup);
 	*subShift = subscriptShift(body, sub, 0);
 	temp1 = (*supShift - bboxDepth(superscriptBBox)) -
 	    (bboxHeight(subscriptBBox) - *subShift);
+}
 static double numeratorVShift(SEXP num)
+{
 	BBOX numBBox;
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
 	double theShift, theClearance, minClearance;
-	GP->cex = GP->cex * scriptScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
 	numBBox = elementBBox(num);
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 	theShift = numeratorShift();
 	theClearance = numeratorShift() - bboxDepth(numBBox) -
 	    (axisHeight() + 0.5 * lineWidth());
 	minClearance = 3 * lineWidth();
 /* RATIO */
 static double denominatorVShift(SEXP denom)
+{
 	BBOX denomBBox;
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
 	double theShift, theClearance, minClearance;
 #ifdef OLD
-	GP->cex = GP->cex * scriptScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
 #else
-	GP->cex = GP->cex * ratioScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * ratioScale;
 #endif
 	denomBBox = elementBBox(denom);
 #ifdef OLD
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 #else
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 #endif
 	theShift = denominatorShift();
 	theClearance = axisHeight() - 0.5 * lineWidth() -
 	    (bboxHeight(denomBBox) - denominatorShift());
 /* RATIO */
 static double fractionWidth(SEXP num, SEXP denom)
+{
 	BBOX numBBox;
 	BBOX denomBBox;
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
 	double temp1, temp2;
 #ifdef OLD
-	GP->cex = GP->cex * scriptScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
 #else
-	GP->cex = GP->cex * ratioScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * ratioScale;
 #endif
 	numBBox = elementBBox(num);
 	denomBBox = elementBBox(denom);
 #ifdef OLD
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 #else
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 #endif
 	temp1 = bboxWidth(numBBox);
 	temp2 = bboxWidth(denomBBox);
 	return max(temp1, temp2);
 static void numdenomHShift(SEXP num, SEXP denom,
 			   double *numShift, double *denomShift)
+{
 	BBOX numBBox;
 	BBOX denomBBox;
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
 	double temp1, temp2;
 #ifdef OLD
-	GP->cex = GP->cex * scriptScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
 #else
-	GP->cex = GP->cex * ratioScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * ratioScale;
 #endif
 	numBBox = elementBBox(num);
 	denomBBox = elementBBox(denom);
 #ifdef OLD
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 #else
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 #endif
 	temp1 = bboxWidth(numBBox);
 	temp2 = bboxWidth(denomBBox);
 	if (temp1 > temp2) {
+}
 static BBOX makeBBoxFromChar(int chr)
+{
 	double height, depth, width;
-	GMetricInfo(chr, &height, &depth, &width, metricUnit);
+	GMetricInfo(chr, &height, &depth, &width, metricUnit, mathDevice);
 	return makeBBox(height, depth, width);
+}
 static BBOX shiftBBox(BBOX bbox, double shiftV)
+{
 static double currentAngle;
 static double cosAngle;
 static double sinAngle;
 		/*
-		   // convert currentX and currentY from inches and 0 angle
+		   // convert currentX and currentY from 0 angle
-		   // to figure units and currentAngle
+		   // to and currentAngle
 		 */
 static double convertedX()
+{
 	double rotatedX = referenceX +
 	(currentX - referenceX) * cosAngle -
 	(currentY - referenceY) * sinAngle;
-	return xInchtoFig(rotatedX);
+	return rotatedX;
+}
 static double convertedY()
+{
 	double rotatedY = referenceY +
 	(currentY - referenceY) * cosAngle +
 	(currentX - referenceX) * sinAngle;
-	return yInchtoFig(rotatedY);
+	return rotatedY;
+}
 static void moveAcross(double xamount)
+{
 	currentX += xamount;
 	/* been converted to ascii are in the symbol font */
 static BBOX asciiBBox(int ascii)
+{
 	if ((ascii == hatAscii()) || (ascii == tildeAscii()))
-		GP->font = 1;
+		mathDevice->gp.font = 1;
 	else
-		GP->font = 5;
+		mathDevice->gp.font = 5;
 	return makeBBoxFromChar(ascii);
+}
 static void drawAscii(int ascii)
+{
 	char asciiStr[2];
 	if ((ascii == hatAscii()) || (ascii == tildeAscii()))
-		GP->font = 1;
+		mathDevice->gp.font = 1;
 	else
-		GP->font = 5;
+		mathDevice->gp.font = 5;
 	asciiStr[0] = ascii;
 	asciiStr[1] = '\0';
-	GText(convertedX(), convertedY(), asciiStr, 0.0, 0.0, currentAngle);
+	GText(convertedX(), convertedY(), INCHES, asciiStr,
+.0, 0.0, currentAngle, mathDevice);
-	moveAcross(GStrWidth(asciiStr, metricUnit));
+	moveAcross(GStrWidth(asciiStr, metricUnit, mathDevice));
+}
 /* code for character atoms */
 static BBOX charBBox(char *str, SEXP expr)
+{
 	BBOX resultBBox = nullBBox();
 	int i;
-	GP->font = atomFontFace(expr);
+	mathDevice->gp.font = atomFontFace(expr);
 	for (i = 0; i < strlen(str); i++)
 		resultBBox = combineBBoxes(resultBBox, makeBBoxFromChar(str[i]));
 	return resultBBox;
+}
 static void drawChar(char *str, SEXP expr)
+{
-	GP->font = atomFontFace(expr);
+	mathDevice->gp.font = atomFontFace(expr);
-	GText(convertedX(), convertedY(), str, 0.0, 0.0, currentAngle);
+	GText(convertedX(), convertedY(), INCHES, str,
+.0, 0.0, currentAngle, mathDevice);
-	moveAcross(GStrWidth(str, metricUnit));
+	moveAcross(GStrWidth(str, metricUnit, mathDevice));
+}
 /* code for symbol atoms */
 static BBOX symbolBBox(SEXP expr)
 static int cexGap = 1;
 static void setGapCEX()
+{
-  cexGap = GP->cex;
+  cexGap = mathDevice->gp.cex;
+}
 static BBOX gapBBox(double gap)
+{
-  double cexSaved = GP->cex;
+  double cexSaved = mathDevice->gp.cex;
   BBOX theBBox;
-  GP->cex = cexGap;
+  mathDevice->gp.cex = cexGap;
   theBBox = makeBBox(0, 0, gap * fontHeight());
-  GP->cex = cexSaved;
+  mathDevice->gp.cex = cexSaved;
   return theBBox;
+}
 static BBOX smallgapBBox(double gap)
+{
-  double cexSaved = GP->cex;
+  double cexSaved = mathDevice->gp.cex;
   BBOX theBBox;
-  GP->cex = cexGap;
+  mathDevice->gp.cex = cexGap;
   theBBox = makeBBox(0, 0, 0.5 * gap * fontHeight());
-  GP->cex = cexSaved;
+  mathDevice->gp.cex = cexSaved;
   return theBBox;
+}
 static void drawGap(double gap)
+{
-  double cexSaved = GP->cex;
+  double cexSaved = mathDevice->gp.cex;
-  GP->cex = cexGap;
+  mathDevice->gp.cex = cexGap;
   moveAcross(gap * fontHeight());
-  GP->cex = cexSaved;
+  mathDevice->gp.cex = cexSaved;
+}
 static void drawSmallGap(double gap)
+{
-  double cexSaved = GP->cex;
+  double cexSaved = mathDevice->gp.cex;
-  GP->cex = cexGap;
+  mathDevice->gp.cex = cexGap;
   moveAcross(0.5 * gap * fontHeight());
-  GP->cex = cexSaved;
+  mathDevice->gp.cex = cexSaved;
+}
 /* code for binary operator (+, -, *, /) expressions */
 /* NOTE that gaps are specified as proportions of the current font height */
 static BBOX supsubBBox(SEXP body, SEXP superscript, SEXP subscript);
 static BBOX superscriptBBox(SEXP superscript)
+{
 	BBOX result;
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
-	GP->cex = GP->cex * scriptScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
 	result = elementBBox(superscript);
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 	return result;
+}
 static BBOX supBBox(SEXP expr)
+}
 static BBOX subscriptBBox(SEXP subscript)
+{
 	BBOX result;
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
-	GP->cex = GP->cex * scriptScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
 	result = elementBBox(subscript);
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 	return result;
+}
 static BBOX subBBox(SEXP expr)
 		   shiftBBox(subscriptBBox(subscript), subShift)));
+}
 static void drawScriptElement(SEXP expr)
+{
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
-	GP->cex = GP->cex * scriptScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
 	drawElement(expr);
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
+}
 static void drawSupSub(SEXP body, SEXP superscript, SEXP subscript);
 static void drawSuper(SEXP expr)
+{
 	BBOX bodyBBox = elementBBox(body);
 	double width = bboxWidth(bodyBBox);
 	double savedX = currentX;
 	double savedY = currentY;
+	double x[3], y[3];
 	moveUp(bboxHeight(bodyBBox) + customAccentGap());
-	GStartPath();
-	GMoveTo(convertedX(), convertedY());
+	x[0] = convertedX(); y[0] = convertedY();
 	moveUp(customHatHeight());
 	moveAcross(width / 2);
-	GLineTo(convertedX(), convertedY());
+	x[1] = convertedX(); y[1] = convertedY();
 	moveUp(-customHatHeight());
 	moveAcross(width / 2);
-	GLineTo(convertedX(), convertedY());
+	x[2] = convertedX(); y[2] = convertedY();
-	GEndPath();
+	GPolyline(3, x, y, INCHES, mathDevice);
 	moveTo(savedX, savedY);
 	drawElement(body);
+}
 static void drawBar(SEXP body)
+{
 	BBOX bodyBBox = elementBBox(body);
 	double savedX = currentX;
 	double savedY = currentY;
+	double x[2], y[2];
 	moveUp(bboxHeight(bodyBBox) + customAccentGap());
-	GStartPath();
-	GMoveTo(convertedX(), convertedY());
+	x[0] = convertedX(); y[0] = convertedY();
 	moveAcross(bboxWidth(bodyBBox));
-	GLineTo(convertedX(), convertedY());
+	x[1] = convertedX(); y[1] = convertedY();
-	GEndPath();
+	GPolyline(2, x, y, INCHES, mathDevice);
 	moveTo(savedX, savedY);
 	drawElement(body);
+}
 static void drawAccent(SEXP expr)
+{
 	SEXP numerator = CADR(expr);
 	SEXP denominator = CADDR(expr);
 	BBOX numBBox, denomBBox;
 	double numHShift, denomHShift;
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
 #ifdef OLD
-	GP->cex = GP->cex * scriptScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
 #else
-	GP->cex = GP->cex * ratioScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * ratioScale;
 #endif
 	numBBox = elementBBox(numerator);
 	denomBBox = elementBBox(denominator);
 #ifdef OLD
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 #else
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
 #endif
 	numdenomHShift(numerator, denominator, &numHShift, &denomHShift);
 	return combineAlignedBBoxes(
 	     shiftBBox(combineBBoxes(makeBBox(numHShift, 0, 0), numBBox),
 		   -denominatorVShift(denominator)));
+}
 static void drawRatioElement(SEXP expr)
+{
-	float cexSaved = GP->cex;
+	float cexSaved = mathDevice->gp.cex;
-	GP->cex = GP->cex * ratioScale;
+	mathDevice->gp.cex = mathDevice->gp.cex * ratioScale;
 	drawElement(expr);
-	GP->cex = cexSaved;
+	mathDevice->gp.cex = cexSaved;
+}
 static void drawFraction(SEXP expr)
+{
 	SEXP numerator = CADR(expr);
 	SEXP denominator = CADDR(expr);
 	double savedX = currentX;
 	double savedY = currentY;
 	double fWidth = fractionWidth(numerator, denominator);
 	double numHShift, denomHShift;
+	double x[2], y[2];
 	numdenomHShift(numerator, denominator, &numHShift, &denomHShift);
 	moveAcross(numHShift);
 	moveUp(numeratorVShift(numerator));
 #ifdef OLD
 #else
 	drawRatioElement(numerator);
 #endif
 	moveTo(savedX, savedY);
 	moveUp(axisHeight());
-	GStartPath();
-	GMoveTo(convertedX(), convertedY());
+	x[0] = convertedX(); y[0] = convertedY();
 	moveAcross(fWidth);
-	GLineTo(convertedX(), convertedY());
+	x[1] = convertedX(); y[1] = convertedY();
-	GEndPath();
+	GPolyline(2, x, y, INCHES, mathDevice);
 	moveTo(savedX, savedY);
 	moveAcross(denomHShift);
 	moveUp(-denominatorVShift(denominator));
 #ifdef OLD
 	drawScriptElement(denominator);
 static useRelGap = 1;
 static BBOX operatorLimitBBox(SEXP limit)
+{
-  float cexSaved = GP->cex;
+  float cexSaved = mathDevice->gp.cex;
   BBOX limitBBox;
-  GP->cex = GP->cex * scriptScale;
+  mathDevice->gp.cex = mathDevice->gp.cex * scriptScale;
   useRelGap = 0;
   limitBBox = elementBBox(limit);
   useRelGap = 1;
-  GP->cex = cexSaved;
+  mathDevice->gp.cex = cexSaved;
   return limitBBox;
+}
 static BBOX operatorBBox(SEXP expr)
 	double depth = bboxDepth(bodyBBox);
 	double width = bboxWidth(bodyBBox);
 	double twiddleHeight = (height - depth) / 2;
 	double savedX = currentX;
 	double savedY = currentY;
+	double x[5], y[5];
 	moveUp(0.8 * twiddleHeight);
-	GStartPath();
-	GMoveTo(convertedX(), convertedY());
+	x[0] = convertedX(); y[0] = convertedY();
 	moveUp(0.2 * twiddleHeight);
 	moveAcross(0.3 * customRadicalWidth());
-	GLineTo(convertedX(), convertedY());
+	x[1] = convertedX(); y[1] = convertedY();
 	moveUp(-(twiddleHeight + depth));
 	moveAcross(0.3 * customRadicalWidth());
-	GLineTo(convertedX(), convertedY());
+	x[2] = convertedX(); y[2] = convertedY();
 	moveUp(depth + height + customRadicalGap());
 	moveAcross(0.4 * customRadicalWidth());
-	GLineTo(convertedX(), convertedY());
+	x[3] = convertedX(); y[3] = convertedY();
 	moveAcross(customRadicalSpace() + width);
-	GLineTo(convertedX(), convertedY());
+	x[4] = convertedX(); y[4] = convertedY();
-	GEndPath();
+	GPolyline(5, x, y, INCHES, mathDevice);
 	moveTo(savedX, savedY);
 	moveAcross(customRadicalWidth() + customRadicalSpace());
 	drawElement(body);
+}
+{
 	SEXP body = CADR(expr);
 	BBOX bodyBBox = elementBBox(expr);
 	double height = bboxHeight(bodyBBox);
 	double depth = bboxDepth(bodyBBox);
+	double width = bboxWidth(bodyBBox);
+	double savedX = currentX;
+	double savedY = currentY;
+	double x[2], y[2];
 	moveUp(-depth);
-	GStartPath();
-	GMoveTo(convertedX(), convertedY());
+	x[0] = convertedX(); y[0] = convertedY();
 	moveUp(depth + height);
-	GLineTo(convertedX(), convertedY());
+	x[1] = convertedX(); y[1] = convertedY();
+	GPolyline(2, x, y, INCHES, mathDevice);
 	moveUp(-height);
 	moveAcross(absSpace());
 	drawElement(body);
 	moveAcross(absSpace());
 	moveUp(-depth);
-	GMoveTo(convertedX(), convertedY());
+	x[0] = convertedX(); y[0] = convertedY();
 	moveUp(depth + height);
-	GLineTo(convertedX(), convertedY());
+	x[1] = convertedX(); y[1] = convertedY();
-	GEndPath();
+	GPolyline(2, x, y, INCHES, mathDevice);
 	moveUp(-height);
+}
 /* code for general expressions with no special meaning in */
 /* mathematical notation syntax (e.g., f(x)) */
 		drawAtom(expr);
+}
         /* calculate width of expression */
         /* BBOXes are in INCHES (see metricUnit) */
-double GExpressionWidth(SEXP expr, int units)
+double GExpressionWidth(SEXP expr, int units, DevDesc *dd)
+{
         BBOX exprBBox = elementBBox(expr);
         double w  = exprBBox.width;
-        switch(units) {
+	if (units == INCHES)
-                case 1: /* user == world */
-                        w = ((exprBBox.width / GP->ipr[0]) / GP->fig2dev.bx) / GP->win2fig.bx;
-                        break;
-                case 2: /* figure */
-                        w = (exprBBox.width / GP->ipr[0]) / GP->fig2dev.bx;
-                        break;
-                case 3: /* inches */
-                        w = exprBBox.width;
-                        break;
-                case 4: /* rasters */
-                        w = exprBBox.width / GP->ipr[0];
-                        break;
+		return w;
-        }
+	else
-        return w;
+		return GConvertXUnits(w, INCHES, units, dd);
+}
 #define ABS(a)  ((a)>=0 ? (a) : -(a))
-double GExpressionHeight(SEXP expr, int units)
+double GExpressionHeight(SEXP expr, int units, DevDesc *dd)
+{
         BBOX exprBBox = elementBBox(expr);
         double h = exprBBox.height + exprBBox.depth;
-        switch(units) {
+	if (units == INCHES)
-                case 1: /* user == world */
-                        h = ((h / GP->ipr[1]) / ABS(GP->fig2dev.by)) / GP->win2fig.by;
-                        break;
+		return h;
-                case 2: /* figure */
-                        h = (h / GP->ipr[1]) / ABS(GP->fig2dev.by);
-                        break;
-                case 3: /* inches */
-                        break;
-                case 4: /* rasters */
-                        h = h / GP->ipr[1];
-                        break;
-        }
+	else
-        return h;
+		return GConvertYUnits(h, INCHES, units, dd);
+}
 		/* functions forming the API */
-void GMathText(double x, double y, SEXP expr, double xc, double yc, double rot)
+void GMathText(double x, double y, int coords, SEXP expr,
+	       double xc, double yc, double rot, DevDesc *dd)
+{
 	BBOX expressionBBox;
+	mathDevice = dd;
 	initFormulaSymbols();
 	expressionBBox = elementBBox(expr);
-	/* NOTE that x and y are in Figure coordinates */
-	referenceX = xFigtoInch(x);
+	referenceX = x;
-	referenceY = yFigtoInch(y);
+	referenceY = y;
+	GConvert(&referenceX, &referenceY, coords, INCHES, dd);
 	currentX = referenceX - xc * bboxWidth(expressionBBox);
 	currentY = referenceY - yc * bboxHeight(expressionBBox);
 	currentAngle = rot;
 	cosAngle = cos(rot / 90 * half_pi);
 	sinAngle = sin(rot / 90 * half_pi);
 	drawElement(expr);
+}
-#define XINVFMAP(x) ((x - GP->fig2dev.ax)/GP->fig2dev.bx)
-#define YINVFMAP(y) ((y - GP->fig2dev.ay)/GP->fig2dev.by)
-void GMMathText(SEXP str, int side, double line, int outer, double at, int las)
+void GMMathText(SEXP str, int side, double line, int outer, double at, int las,
+		DevDesc *dd)
+{
+	int coords;
-	double a, x, y, xadj, yadj;
+	double a, xadj, yadj;
+	mathDevice = dd;
 	if (outer) {
 		switch (side) {
 		case 1:
-			x = at;
+			line = line + 1;
-			y = yChartoNDC(GP->cexbase * GP->mex * (GP->oma[0] - line + 1));
-			x = DP->inner2dev.ax + DP->inner2dev.bx * x;
-			y = DP->ndc2dev.ay + DP->ndc2dev.by * y;
+			coords = MAR1;
 			a = 0.0;
-			xadj = GP->adj;
+			xadj = mathDevice->gp.adj;
 			yadj = 0.0;
 			break;
 		case 2:
-			x = xChartoNDC(GP->cexbase * GP->mex * (GP->oma[1] - line));
-			y = at;
+			coords = MAR2;
-			x = DP->ndc2dev.ax + DP->ndc2dev.bx * x;
-			y = DP->inner2dev.ay + DP->inner2dev.by * y;
 			a = 90.0;
-			xadj = GP->adj;
+			xadj = mathDevice->gp.adj;
 			yadj = 0.0;
 			break;
 		case 3:
-			x = at;
+			coords = MAR3;
-			y = 1.0 - yChartoNDC(GP->cexbase * GP->mex * (GP->oma[2] - line));
-			x = DP->inner2dev.ax + DP->inner2dev.bx * x;
-			y = DP->ndc2dev.ay + DP->ndc2dev.by * y;
 			a = 0.0;
-			xadj = GP->adj;
+			xadj = mathDevice->gp.adj;
 			yadj = 0.0;
 			break;
 		case 4:
-			x = 1.0 - xChartoNDC(GP->cexbase * GP->mex * (GP->oma[3] - line));
-			y = at;
+			line = line + 1;
-			x = DP->ndc2dev.ax + DP->ndc2dev.bx * x;
+			coords = MAR4;
-			y = DP->inner2dev.ay + DP->inner2dev.by * y;
 			a = 90.0;
-			xadj = GP->adj;
+			xadj = mathDevice->gp.adj;
 			yadj = 0.0;
 			break;
+		}
-		x = XINVFMAP(x);
-		y = YINVFMAP(y);
-		GMathText(x, y, str, xadj, yadj, a);
+		GMathText(at, line, coords, str, xadj, yadj, a, dd);
+	}
 	else {
 		switch (side) {
 		case 1:
 			if (las == 2) {
-				y = GP->plt[2] - yInchtoFig(yChartoInch(GP->cexbase * GP->mex * (line + GP->yLineBias)));
+				at = at - GConvertXUnits(dd->gp.yLineBias,
+							 LINES, USER, dd);
-				x = XMAP(at) - xInchtoFig(xChartoInch(GP->cexbase * GP->mex * GP->yLineBias));
+				line = line + dd->gp.yLineBias;
 				a = 90.0;
 				xadj = 1.0;
 				yadj = 0.5;
+			}
 			else {
-				y = GP->plt[2] - yInchtoFig(yChartoInch(GP->cexbase * GP->mex * (line + 1 - GP->yLineBias)));
+				line = line + 1 - dd->gp.yLineBias;
-				x = XMAP(at);
 				a = 0.0;
-				xadj = GP->adj;
+				xadj = mathDevice->gp.adj;
 				yadj = 0.0;
+			}
+			coords = MAR1;
 			break;
 		case 2:
 			if (las == 1 || las == 2) {
-				x = GP->plt[0] - xInchtoFig(xChartoInch(GP->cexbase * GP->mex * (line + GP->yLineBias)));
+				at = at + GConvertYUnits(dd->gp.yLineBias,
+							 LINES, USER, dd);
-				y = YMAP(at) + yInchtoFig(yChartoInch(GP->cexbase * GP->mex * GP->yLineBias));
+				line = line + dd->gp.yLineBias;
 				a = 0.0;
 				xadj = 1.0;
 				yadj = 0.5;
+			}
 			else {
-				x = GP->plt[0] - xInchtoFig(xChartoInch(GP->cexbase * GP->mex * (line + GP->yLineBias)));
-				y = YMAP(at);
+				line = line + dd->gp.yLineBias;
 				a = 90.0;
-				xadj = GP->adj;
+				xadj = mathDevice->gp.adj;
 				yadj = 0.0;
+			}
+			coords = MAR2;
 			break;
 		case 3:
 			if (las == 2) {
-				y = GP->plt[3] + yInchtoFig(yChartoInch(GP->cexbase * GP->mex * (line + GP->yLineBias)));
+				at = at - GConvertXUnits(dd->gp.yLineBias,
+							 LINES, USER, dd);
-				x = XMAP(at) - xInchtoFig(xChartoInch(GP->cexbase * GP->mex * GP->yLineBias));
+				line = line + dd->gp.yLineBias;
 				a = 90.0;
 				xadj = 0.0;
 				yadj = 0.5;
+			}
 			else {
-				y = GP->plt[3] + yInchtoFig(yChartoInch(GP->cexbase * GP->mex * (line + GP->yLineBias)));
-				x = XMAP(at);
+				line = line + dd->gp.yLineBias;
 				a = 0.0;
-				xadj = GP->adj;
+				xadj = mathDevice->gp.adj;
 				yadj = 0.0;
+			}
+			coords = MAR3;
 			break;
 		case 4:
 			if (las == 1 || las == 2) {
-				x = GP->plt[1] + xInchtoFig(xChartoInch(GP->cexbase * GP->mex * (line + GP->yLineBias)));
+				at = at + GConvertYUnits(dd->gp.yLineBias,
+							 LINES, USER, dd);
-				y = YMAP(at) + yInchtoFig(yChartoInch(GP->cexbase * GP->mex * GP->yLineBias));
+				line = line + dd->gp.yLineBias;
 				a = 0.0;
 				xadj = 0.0;
 				yadj = 0.5;
+			}
 			else {
-				x = GP->plt[1] + xInchtoFig(xChartoInch(GP->cexbase * GP->mex * (line + 1 - GP->yLineBias)));
+				line = line + 1 - dd->gp.yLineBias;
-				y = YMAP(at);
 				a = 90.0;
-				xadj = GP->adj;
+				xadj = mathDevice->gp.adj;
 				yadj = 0.0;
+			}
+			coords = MAR4;
 			break;
+		}
-		GMathText(x, y, str, xadj, yadj, a);
+		GMathText(at, line, coords, str, xadj, yadj, a, dd);
+	}
+}
 #else

The R Project SVN R

(root)/branches/R-exp-R5/src/main/plotmath.c @ 63243 – Rev 526 → 581