xref: /plan9/sys/src/cmd/eqn/matrix.c (revision 3e12c5d1bb89fc02707907988834ef147769ddaf)
1*3e12c5d1SDavid du Colombier #include "e.h"
2*3e12c5d1SDavid du Colombier 
startcol(int type)3*3e12c5d1SDavid du Colombier startcol(int type)	/* mark start of column in lp[] array */
4*3e12c5d1SDavid du Colombier {
5*3e12c5d1SDavid du Colombier 	int oct = ct;
6*3e12c5d1SDavid du Colombier 
7*3e12c5d1SDavid du Colombier 	lp[ct++] = type;
8*3e12c5d1SDavid du Colombier 	lp[ct++] = 0;	/* count, to come */
9*3e12c5d1SDavid du Colombier 	lp[ct++] = 0;	/* separation, to come */
10*3e12c5d1SDavid du Colombier 	return oct;
11*3e12c5d1SDavid du Colombier }
12*3e12c5d1SDavid du Colombier 
column(int oct,int sep)13*3e12c5d1SDavid du Colombier void column(int oct, int sep)	/* remember end of column that started at lp[oct] */
14*3e12c5d1SDavid du Colombier {
15*3e12c5d1SDavid du Colombier 	int i, type;
16*3e12c5d1SDavid du Colombier 
17*3e12c5d1SDavid du Colombier 	lp[oct+1] = ct - oct - 3;
18*3e12c5d1SDavid du Colombier 	lp[oct+2] = sep;
19*3e12c5d1SDavid du Colombier 	type = lp[oct];
20*3e12c5d1SDavid du Colombier 	if (dbg) {
21*3e12c5d1SDavid du Colombier 		printf(".\t%d column of", type);
22*3e12c5d1SDavid du Colombier 		for (i = oct+3; i < ct; i++ )
23*3e12c5d1SDavid du Colombier 			printf(" S%d", lp[i]);
24*3e12c5d1SDavid du Colombier 		printf(", rows=%d, sep=%d\n", lp[oct+1], lp[oct+2]);
25*3e12c5d1SDavid du Colombier 	}
26*3e12c5d1SDavid du Colombier }
27*3e12c5d1SDavid du Colombier 
matrix(int oct)28*3e12c5d1SDavid du Colombier void matrix(int oct)	/* matrix is list of columns */
29*3e12c5d1SDavid du Colombier {
30*3e12c5d1SDavid du Colombier 	int nrow, ncol, i, j, k, val[100];
31*3e12c5d1SDavid du Colombier 	double b, hb;
32*3e12c5d1SDavid du Colombier 	char *space;
33*3e12c5d1SDavid du Colombier 	extern char *Matspace;
34*3e12c5d1SDavid du Colombier 
35*3e12c5d1SDavid du Colombier 	space = Matspace;	/* between columns of matrix */
36*3e12c5d1SDavid du Colombier 	nrow = lp[oct+1];	/* disaster if rows inconsistent */
37*3e12c5d1SDavid du Colombier 				/* also assumes just columns */
38*3e12c5d1SDavid du Colombier 				/* fix when add other things */
39*3e12c5d1SDavid du Colombier 	ncol = 0;
40*3e12c5d1SDavid du Colombier 	for (i = oct+1; i < ct; i += lp[i]+3 ) {
41*3e12c5d1SDavid du Colombier 		ncol++;
42*3e12c5d1SDavid du Colombier 		dprintf(".\tcolct=%d\n", lp[i]);
43*3e12c5d1SDavid du Colombier 	}
44*3e12c5d1SDavid du Colombier 	for (k=1; k <= nrow; k++) {
45*3e12c5d1SDavid du Colombier 		hb = b = 0;
46*3e12c5d1SDavid du Colombier 		j = oct + k + 2;
47*3e12c5d1SDavid du Colombier 		for (i=0; i < ncol; i++) {
48*3e12c5d1SDavid du Colombier 			hb = max(hb, eht[lp[j]]-ebase[lp[j]]);
49*3e12c5d1SDavid du Colombier 			b = max(b, ebase[lp[j]]);
50*3e12c5d1SDavid du Colombier 			j += nrow + 3;
51*3e12c5d1SDavid du Colombier 		}
52*3e12c5d1SDavid du Colombier 		dprintf(".\trow %d: b=%g, hb=%g\n", k, b, hb);
53*3e12c5d1SDavid du Colombier 		j = oct + k + 2;
54*3e12c5d1SDavid du Colombier 		for (i=0; i<ncol; i++) {
55*3e12c5d1SDavid du Colombier 			ebase[lp[j]] = b;
56*3e12c5d1SDavid du Colombier 			eht[lp[j]] = b + hb;
57*3e12c5d1SDavid du Colombier 			j += nrow + 3;
58*3e12c5d1SDavid du Colombier 		}
59*3e12c5d1SDavid du Colombier 	}
60*3e12c5d1SDavid du Colombier 	j = oct;
61*3e12c5d1SDavid du Colombier 	for (i=0; i<ncol; i++) {
62*3e12c5d1SDavid du Colombier 		pile(j);
63*3e12c5d1SDavid du Colombier 		val[i] = yyval;
64*3e12c5d1SDavid du Colombier 		j += nrow + 3;
65*3e12c5d1SDavid du Colombier 	}
66*3e12c5d1SDavid du Colombier 	yyval = salloc();
67*3e12c5d1SDavid du Colombier 	eht[yyval] = eht[val[0]];
68*3e12c5d1SDavid du Colombier 	ebase[yyval] = ebase[val[0]];
69*3e12c5d1SDavid du Colombier 	lfont[yyval] = rfont[yyval] = 0;
70*3e12c5d1SDavid du Colombier 	dprintf(".\tmatrix S%d: r=%d, c=%d, h=%g, b=%g\n",
71*3e12c5d1SDavid du Colombier 		yyval,nrow,ncol,eht[yyval],ebase[yyval]);
72*3e12c5d1SDavid du Colombier 	printf(".ds %d \"", yyval);
73*3e12c5d1SDavid du Colombier 	for( i=0; i<ncol; i++ )  {
74*3e12c5d1SDavid du Colombier 		printf("\\*(%d%s", val[i], i==ncol-1 ? "" : space);
75*3e12c5d1SDavid du Colombier 		sfree(val[i]);
76*3e12c5d1SDavid du Colombier 	}
77*3e12c5d1SDavid du Colombier 	printf("\n");
78*3e12c5d1SDavid du Colombier }
79