xref: /plan9/sys/src/cmd/eqn/fromto.c (revision 3e12c5d1bb89fc02707907988834ef147769ddaf)
1*3e12c5d1SDavid du Colombier # include "e.h"
2*3e12c5d1SDavid du Colombier 
fromto(int p1,int p2,int p3)3*3e12c5d1SDavid du Colombier void fromto(int p1, int p2, int p3)
4*3e12c5d1SDavid du Colombier {
5*3e12c5d1SDavid du Colombier 	double b, h1, b1, t;
6*3e12c5d1SDavid du Colombier 	int subps;
7*3e12c5d1SDavid du Colombier 
8*3e12c5d1SDavid du Colombier 	yyval = salloc();
9*3e12c5d1SDavid du Colombier 	lfont[yyval] = rfont[yyval] = 0;
10*3e12c5d1SDavid du Colombier 	h1 = eht[yyval] = eht[p1];
11*3e12c5d1SDavid du Colombier 	b1 = ebase[p1];
12*3e12c5d1SDavid du Colombier 	b = 0;
13*3e12c5d1SDavid du Colombier 	subps = ps;
14*3e12c5d1SDavid du Colombier 	ps += deltaps;
15*3e12c5d1SDavid du Colombier 	nrwid(p1, ps, p1);
16*3e12c5d1SDavid du Colombier 	printf(".nr %d \\n(%d\n", yyval, p1);
17*3e12c5d1SDavid du Colombier 	if (p2 > 0) {
18*3e12c5d1SDavid du Colombier 		nrwid(p2, subps, p2);
19*3e12c5d1SDavid du Colombier 		printf(".if \\n(%d>\\n(%d .nr %d \\n(%d\n", p2, yyval, yyval, p2);
20*3e12c5d1SDavid du Colombier 		eht[yyval] += eht[p2];
21*3e12c5d1SDavid du Colombier 		b = eht[p2];
22*3e12c5d1SDavid du Colombier 	}
23*3e12c5d1SDavid du Colombier 	if (p3 > 0) {
24*3e12c5d1SDavid du Colombier 		nrwid(p3, subps, p3);
25*3e12c5d1SDavid du Colombier 		printf(".if \\n(%d>\\n(%d .nr %d \\n(%d\n", p3, yyval, yyval, p3);
26*3e12c5d1SDavid du Colombier 		eht[yyval] += eht[p3];
27*3e12c5d1SDavid du Colombier 	}
28*3e12c5d1SDavid du Colombier 	printf(".ds %d ", yyval);	/* bottom of middle box */
29*3e12c5d1SDavid du Colombier 	if (p2 > 0) {
30*3e12c5d1SDavid du Colombier 		t = eht[p2]-ebase[p2]+b1;
31*3e12c5d1SDavid du Colombier 		printf("\\v'%gm'\\h'\\n(%du-\\n(%du/2u'%s\\*(%d%s",
32*3e12c5d1SDavid du Colombier 			REL(t,ps), yyval, p2, DPS(ps,subps), p2, DPS(subps,ps));
33*3e12c5d1SDavid du Colombier 		printf("\\h'-\\n(%du-\\n(%du/2u'\\v'%gm'\\\n",
34*3e12c5d1SDavid du Colombier 			yyval, p2, REL(-t,ps));
35*3e12c5d1SDavid du Colombier 	}
36*3e12c5d1SDavid du Colombier 	printf("\\h'\\n(%du-\\n(%du/2u'\\*(%d\\h'\\n(%du-\\n(%du/2u'\\\n",
37*3e12c5d1SDavid du Colombier 		yyval, p1, p1, yyval, p1);
38*3e12c5d1SDavid du Colombier 	if (p3  >0) {
39*3e12c5d1SDavid du Colombier 		t = h1-b1+ebase[p3];
40*3e12c5d1SDavid du Colombier 		printf("\\v'%gm'\\h'-\\n(%du-\\n(%du/2u'%s\\*(%d%s\\h'\\n(%du-\\n(%du/2u'\\v'%gm'\\\n",
41*3e12c5d1SDavid du Colombier 			REL(-t,ps), yyval, p3, DPS(ps,subps), p3, DPS(subps,ps), yyval, p3, REL(t,ps));
42*3e12c5d1SDavid du Colombier 	}
43*3e12c5d1SDavid du Colombier 	printf("\n");
44*3e12c5d1SDavid du Colombier 	ebase[yyval] = b + b1;
45*3e12c5d1SDavid du Colombier 	dprintf(".\tS%d <- %d from %d to %d; h=%g b=%g\n",
46*3e12c5d1SDavid du Colombier 		yyval, p1, p2, p3, eht[yyval], ebase[yyval]);
47*3e12c5d1SDavid du Colombier 	sfree(p1);
48*3e12c5d1SDavid du Colombier 	if (p2 > 0)
49*3e12c5d1SDavid du Colombier 		sfree(p2);
50*3e12c5d1SDavid du Colombier 	if (p3 > 0)
51*3e12c5d1SDavid du Colombier 		sfree(p3);
52*3e12c5d1SDavid du Colombier }
53