1*14486Ssam #ifndef lint 2*14486Ssam static char sccsid[] = "@(#)matrix.c 4.2 08/11/83"; 3*14486Ssam #endif 411011Sshannon 511011Sshannon #include "e.h" 611011Sshannon 711011Sshannon column(type, p1) int type, p1; { 811011Sshannon int i; 911011Sshannon 1011011Sshannon lp[p1] = ct - p1 - 1; 1111011Sshannon if( dbg ){ 1211011Sshannon printf(".\t%d column of", type); 1311011Sshannon for( i=p1+1; i<ct; i++ ) 1411011Sshannon printf(" S%d", lp[i]); 1511011Sshannon printf(", rows=%d\n",lp[p1]); 1611011Sshannon } 1711011Sshannon lp[ct++] = type; 1811011Sshannon } 1911011Sshannon 2011011Sshannon matrix(p1) int p1; { 2111011Sshannon int nrow, ncol, i, j, k, hb, b, val[100]; 2211011Sshannon char *space; 2311011Sshannon 2411011Sshannon space = "\\ \\ "; 2511011Sshannon nrow = lp[p1]; /* disaster if rows inconsistent */ 2611011Sshannon ncol = 0; 2711011Sshannon for( i=p1; i<ct; i += lp[i]+2 ){ 2811011Sshannon ncol++; 2911011Sshannon if(dbg)printf(".\tcolct=%d\n",lp[i]); 3011011Sshannon } 3111011Sshannon for( k=1; k<=nrow; k++ ) { 3211011Sshannon hb = b = 0; 3311011Sshannon j = p1 + k; 3411011Sshannon for( i=0; i<ncol; i++ ) { 3511011Sshannon hb = max(hb, eht[lp[j]]-ebase[lp[j]]); 3611011Sshannon b = max(b, ebase[lp[j]]); 3711011Sshannon j += nrow + 2; 3811011Sshannon } 3911011Sshannon if(dbg)printf(".\trow %d: b=%d, hb=%d\n", k, b, hb); 4011011Sshannon j = p1 + k; 4111011Sshannon for( i=0; i<ncol; i++ ) { 4211011Sshannon ebase[lp[j]] = b; 4311011Sshannon eht[lp[j]] = b + hb; 4411011Sshannon j += nrow + 2; 4511011Sshannon } 4611011Sshannon } 4711011Sshannon j = p1; 4811011Sshannon for( i=0; i<ncol; i++ ) { 4911011Sshannon lpile(lp[j+lp[j]+1], j+1, j+lp[j]+1); 5011011Sshannon val[i] = yyval; 5111011Sshannon j += nrow + 2; 5211011Sshannon } 5311011Sshannon yyval = oalloc(); 5411011Sshannon eht[yyval] = eht[val[0]]; 5511011Sshannon ebase[yyval] = ebase[val[0]]; 5611011Sshannon lfont[yyval] = rfont[yyval] = 0; 5711011Sshannon if(dbg)printf(".\tmatrix S%d: r=%d, c=%d, h=%d, b=%d\n", 5811011Sshannon yyval,nrow,ncol,eht[yyval],ebase[yyval]); 5911011Sshannon printf(".ds %d \"", yyval); 6011011Sshannon for( i=0; i<ncol; i++ ) { 6111011Sshannon printf("\\*(%d%s", val[i], i==ncol-1 ? "" : space); 6211011Sshannon ofree(val[i]); 6311011Sshannon } 6411011Sshannon printf("\n"); 6511011Sshannon ct = p1; 6611011Sshannon } 67