xref: /plan9/sys/src/cmd/eqn/shift.c (revision 3e12c5d1bb89fc02707907988834ef147769ddaf)
1*3e12c5d1SDavid du Colombier #include "e.h"
2*3e12c5d1SDavid du Colombier #include "y.tab.h"
3*3e12c5d1SDavid du Colombier 
subsup(int p1,int p2,int p3)4*3e12c5d1SDavid du Colombier void subsup(int p1, int p2, int p3)
5*3e12c5d1SDavid du Colombier {
6*3e12c5d1SDavid du Colombier 	if (p2 != 0 && p3 != 0)
7*3e12c5d1SDavid du Colombier 		shift2(p1, p2, p3);
8*3e12c5d1SDavid du Colombier 	else if (p2 != 0)
9*3e12c5d1SDavid du Colombier 		bshiftb(p1, SUB, p2);
10*3e12c5d1SDavid du Colombier 	else if (p3 != 0)
11*3e12c5d1SDavid du Colombier 		bshiftb(p1, SUP, p3);
12*3e12c5d1SDavid du Colombier }
13*3e12c5d1SDavid du Colombier 
14*3e12c5d1SDavid du Colombier extern double Subbase, Supshift;
15*3e12c5d1SDavid du Colombier extern char *Sub1space, *Sup1space, *Sub2space;
16*3e12c5d1SDavid du Colombier extern char *SS1space, *SS2space;
17*3e12c5d1SDavid du Colombier 
bshiftb(int p1,int dir,int p2)18*3e12c5d1SDavid du Colombier void bshiftb(int p1, int dir, int p2)
19*3e12c5d1SDavid du Colombier {
20*3e12c5d1SDavid du Colombier 	int subps, n;
21*3e12c5d1SDavid du Colombier 	double shval, d1, h1, b1, h2, b2;
22*3e12c5d1SDavid du Colombier 	char *sh1, *sh2;
23*3e12c5d1SDavid du Colombier 
24*3e12c5d1SDavid du Colombier 	yyval = p1;
25*3e12c5d1SDavid du Colombier 	h1 = eht[p1];
26*3e12c5d1SDavid du Colombier 	b1 = ebase[p1];
27*3e12c5d1SDavid du Colombier 	h2 = eht[p2];
28*3e12c5d1SDavid du Colombier 	b2 = ebase[p2];
29*3e12c5d1SDavid du Colombier 	subps = ps;
30*3e12c5d1SDavid du Colombier 	ps += deltaps;
31*3e12c5d1SDavid du Colombier 	if (dir == SUB) {
32*3e12c5d1SDavid du Colombier 		/* base .2m below bottom of main box */
33*3e12c5d1SDavid du Colombier 		shval = b1 + EM(Subbase, ps);
34*3e12c5d1SDavid du Colombier 		ebase[yyval] = shval + b2;
35*3e12c5d1SDavid du Colombier 		eht[yyval] = max(h1-b1+shval+b2, h2);
36*3e12c5d1SDavid du Colombier 		if (rfont[p1] == ITAL && lfont[p2] == ROM)
37*3e12c5d1SDavid du Colombier 			n = 2;		/* Sub1space */
38*3e12c5d1SDavid du Colombier 		else
39*3e12c5d1SDavid du Colombier 			n = max(2, class[rclass[p1]][lclass[p2]]);
40*3e12c5d1SDavid du Colombier 		sh1 = pad(n);
41*3e12c5d1SDavid du Colombier 		rclass[p1] = OTHER;	/* OTHER leaves too much after sup */
42*3e12c5d1SDavid du Colombier 	} else {	/* superscript */
43*3e12c5d1SDavid du Colombier 		/* 4/10 up main box */
44*3e12c5d1SDavid du Colombier 		d1 = EM(Subbase, subps);
45*3e12c5d1SDavid du Colombier 		ebase[yyval] = b1;
46*3e12c5d1SDavid du Colombier 		shval = -(Supshift * (h1-b1)) - b2;
47*3e12c5d1SDavid du Colombier 		if (Supshift*(h1-b1) + h2 < h1-b1)	/* raise little super */
48*3e12c5d1SDavid du Colombier 			shval = -(h1-b1) + h2-b2 - d1;
49*3e12c5d1SDavid du Colombier 		eht[yyval] = h1 + max(0, h2 - (1-Supshift)*(h1-b1));
50*3e12c5d1SDavid du Colombier 		if (rclass[p1] == ILETF)
51*3e12c5d1SDavid du Colombier 			n = 4;
52*3e12c5d1SDavid du Colombier 		else if (rfont[p1] == ITAL)
53*3e12c5d1SDavid du Colombier 			n = 2;		/* Sup1space */
54*3e12c5d1SDavid du Colombier 		else
55*3e12c5d1SDavid du Colombier 			n = max(1, class[rclass[p1]][lclass[p2]]);
56*3e12c5d1SDavid du Colombier 		sh1 = pad(n);
57*3e12c5d1SDavid du Colombier 		rclass[p1] = rclass[p2];	/* OTHER leaves too much after sup */
58*3e12c5d1SDavid du Colombier 	}
59*3e12c5d1SDavid du Colombier 	dprintf(".\tS%d <- %d shift %g %d; b=%g, h=%g, ps=%d, subps=%d\n",
60*3e12c5d1SDavid du Colombier 		yyval, p1, shval, p2, ebase[yyval], eht[yyval], ps, subps);
61*3e12c5d1SDavid du Colombier 	sh2 = Sub2space;	/* was Sub2space; */
62*3e12c5d1SDavid du Colombier 	printf(".as %d \\v'%gm'%s%s\\*(%d%s%s\\v'%gm'\n",
63*3e12c5d1SDavid du Colombier 		yyval, REL(shval,ps), DPS(ps,subps), sh1, p2,
64*3e12c5d1SDavid du Colombier 		DPS(subps,ps), sh2, REL(-shval,ps));
65*3e12c5d1SDavid du Colombier 	rfont[p1] = 0;
66*3e12c5d1SDavid du Colombier 	sfree(p2);
67*3e12c5d1SDavid du Colombier }
68*3e12c5d1SDavid du Colombier 
shift2(int p1,int p2,int p3)69*3e12c5d1SDavid du Colombier void shift2(int p1, int p2, int p3)
70*3e12c5d1SDavid du Colombier {
71*3e12c5d1SDavid du Colombier 	int subps;
72*3e12c5d1SDavid du Colombier 	double h1, h2, h3, b1, b2, b3, subsh, d2, supsh;
73*3e12c5d1SDavid du Colombier 	int treg;
74*3e12c5d1SDavid du Colombier 	char *sh2;
75*3e12c5d1SDavid du Colombier 
76*3e12c5d1SDavid du Colombier 	treg = salloc();
77*3e12c5d1SDavid du Colombier 	yyval = p1;
78*3e12c5d1SDavid du Colombier 	subps = ps;	/* sub and sup at this size */
79*3e12c5d1SDavid du Colombier 	ps += deltaps;	/* outer size */
80*3e12c5d1SDavid du Colombier 	h1 = eht[p1]; b1 = ebase[p1];
81*3e12c5d1SDavid du Colombier 	h2 = eht[p2]; b2 = ebase[p2];
82*3e12c5d1SDavid du Colombier 	h3 = eht[p3]; b3 = ebase[p3];
83*3e12c5d1SDavid du Colombier 	subsh = EM(Subbase, ps);
84*3e12c5d1SDavid du Colombier 	if (b1 > b2 + subsh) /* move little sub down */
85*3e12c5d1SDavid du Colombier 		subsh += b1;
86*3e12c5d1SDavid du Colombier 	eht[yyval] = max(subsh+b2-b1+h1, h2);
87*3e12c5d1SDavid du Colombier 	supsh = -Supshift*(h1-b1) - b3;
88*3e12c5d1SDavid du Colombier 	d2 = EM(Subbase, subps);
89*3e12c5d1SDavid du Colombier 	if (h3 < (1-Supshift)*(h1-b1))
90*3e12c5d1SDavid du Colombier 		supsh = -(h1-b1) + (h3-b3) - d2;
91*3e12c5d1SDavid du Colombier 	ebase[yyval] = subsh + b2 - b1;
92*3e12c5d1SDavid du Colombier 	eht[yyval] = h1 + subsh+b2-b1 + max(0, h3-(1-Supshift)*(h1-b1));
93*3e12c5d1SDavid du Colombier 	dprintf(".\tS%d <- %d sub %d sup %d, ps=%d, subps=%d, h=%g, b=%g\n",
94*3e12c5d1SDavid du Colombier 		yyval, p1, p2, p3, ps, subps, eht[yyval], ebase[yyval]);
95*3e12c5d1SDavid du Colombier 	if (rclass[p1] == ILETF)
96*3e12c5d1SDavid du Colombier 		sh2 = "\\|\\|";
97*3e12c5d1SDavid du Colombier 	else
98*3e12c5d1SDavid du Colombier 		sh2 = SS2space;
99*3e12c5d1SDavid du Colombier 	/*n = max(class[rclass[p1]][lclass[p2]], class[rclass[p1]][lclass[p3]]);
100*3e12c5d1SDavid du Colombier 	/*sh2 = pad(max(2, n));
101*3e12c5d1SDavid du Colombier 	*/
102*3e12c5d1SDavid du Colombier 	printf(".ds %d %s\\*(%d\n", p2, SS1space, p2);
103*3e12c5d1SDavid du Colombier 	nrwid(p2, subps, p2);
104*3e12c5d1SDavid du Colombier 	printf(".ds %d %s\\*(%d\n", p3, sh2, p3);
105*3e12c5d1SDavid du Colombier 	nrwid(p3, subps, p3);
106*3e12c5d1SDavid du Colombier 	printf(".nr %d \\n(%d\n", treg, p3);
107*3e12c5d1SDavid du Colombier 	printf(".if \\n(%d>\\n(%d .nr %d \\n(%d\n", p2, treg, treg, p2);
108*3e12c5d1SDavid du Colombier 	printf(".as %d %s\\v'%gm'\\*(%d\\v'%gm'\\h'-\\n(%du'\\\n",
109*3e12c5d1SDavid du Colombier 		p1, DPS(ps,subps), REL(subsh,subps), p2, REL(-subsh,subps), p2);
110*3e12c5d1SDavid du Colombier 	printf("\\v'%gm'\\*(%d\\v'%gm'\\h'-\\n(%du+\\n(%du'%s%s\n",
111*3e12c5d1SDavid du Colombier 		REL(supsh,subps), p3, REL(-supsh,subps), p3, treg, DPS(subps,ps), Sub2space);
112*3e12c5d1SDavid du Colombier 	if (rfont[p2] == ITAL)
113*3e12c5d1SDavid du Colombier 		rfont[yyval] = 0;	/* lie */
114*3e12c5d1SDavid du Colombier 	rclass[yyval] = rclass[p3];	/* was OTHER */
115*3e12c5d1SDavid du Colombier 	sfree(p2); sfree(p3); sfree(treg);
116*3e12c5d1SDavid du Colombier }
117