1*11011Sshannon /* matrix.c 4.1 83/02/11 */ 2*11011Sshannon 3*11011Sshannon #include "e.h" 4*11011Sshannon 5*11011Sshannon column(type, p1) int type, p1; { 6*11011Sshannon int i; 7*11011Sshannon 8*11011Sshannon lp[p1] = ct - p1 - 1; 9*11011Sshannon if( dbg ){ 10*11011Sshannon printf(".\t%d column of", type); 11*11011Sshannon for( i=p1+1; i<ct; i++ ) 12*11011Sshannon printf(" S%d", lp[i]); 13*11011Sshannon printf(", rows=%d\n",lp[p1]); 14*11011Sshannon } 15*11011Sshannon lp[ct++] = type; 16*11011Sshannon } 17*11011Sshannon 18*11011Sshannon matrix(p1) int p1; { 19*11011Sshannon int nrow, ncol, i, j, k, hb, b, val[100]; 20*11011Sshannon char *space; 21*11011Sshannon 22*11011Sshannon space = "\\ \\ "; 23*11011Sshannon nrow = lp[p1]; /* disaster if rows inconsistent */ 24*11011Sshannon ncol = 0; 25*11011Sshannon for( i=p1; i<ct; i += lp[i]+2 ){ 26*11011Sshannon ncol++; 27*11011Sshannon if(dbg)printf(".\tcolct=%d\n",lp[i]); 28*11011Sshannon } 29*11011Sshannon for( k=1; k<=nrow; k++ ) { 30*11011Sshannon hb = b = 0; 31*11011Sshannon j = p1 + k; 32*11011Sshannon for( i=0; i<ncol; i++ ) { 33*11011Sshannon hb = max(hb, eht[lp[j]]-ebase[lp[j]]); 34*11011Sshannon b = max(b, ebase[lp[j]]); 35*11011Sshannon j += nrow + 2; 36*11011Sshannon } 37*11011Sshannon if(dbg)printf(".\trow %d: b=%d, hb=%d\n", k, b, hb); 38*11011Sshannon j = p1 + k; 39*11011Sshannon for( i=0; i<ncol; i++ ) { 40*11011Sshannon ebase[lp[j]] = b; 41*11011Sshannon eht[lp[j]] = b + hb; 42*11011Sshannon j += nrow + 2; 43*11011Sshannon } 44*11011Sshannon } 45*11011Sshannon j = p1; 46*11011Sshannon for( i=0; i<ncol; i++ ) { 47*11011Sshannon lpile(lp[j+lp[j]+1], j+1, j+lp[j]+1); 48*11011Sshannon val[i] = yyval; 49*11011Sshannon j += nrow + 2; 50*11011Sshannon } 51*11011Sshannon yyval = oalloc(); 52*11011Sshannon eht[yyval] = eht[val[0]]; 53*11011Sshannon ebase[yyval] = ebase[val[0]]; 54*11011Sshannon lfont[yyval] = rfont[yyval] = 0; 55*11011Sshannon if(dbg)printf(".\tmatrix S%d: r=%d, c=%d, h=%d, b=%d\n", 56*11011Sshannon yyval,nrow,ncol,eht[yyval],ebase[yyval]); 57*11011Sshannon printf(".ds %d \"", yyval); 58*11011Sshannon for( i=0; i<ncol; i++ ) { 59*11011Sshannon printf("\\*(%d%s", val[i], i==ncol-1 ? "" : space); 60*11011Sshannon ofree(val[i]); 61*11011Sshannon } 62*11011Sshannon printf("\n"); 63*11011Sshannon ct = p1; 64*11011Sshannon } 65