The R Project SVN R

Rev

Rev 4928 | Rev 5187 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 4928 Rev 5055
Line 271... Line 271...
271
    PROTECT(ser = allocVector(INTSXP, c));
271
    PROTECT(ser = allocVector(INTSXP, c));
272
    PROTECT(swr = allocVector(INTSXP, c));
272
    PROTECT(swr = allocVector(INTSXP, c));
273
    PROTECT(sdi = allocVector(INTSXP, c));
273
    PROTECT(sdi = allocVector(INTSXP, c));
274
    PROTECT(sei = allocVector(INTSXP, c));
274
    PROTECT(sei = allocVector(INTSXP, c));
275
    PROTECT(swi = allocVector(INTSXP, c));
275
    PROTECT(swi = allocVector(INTSXP, c));
276
    PROTECT(sw = allocVector(INTSXP, c));
276
    PROTECT(sw  = allocVector(INTSXP, c));
277
    UNPROTECT(7);
277
    UNPROTECT(7);
278
    x = COMPLEX(sx) + offset;
278
    x = COMPLEX(sx) + offset;
279
    dr = INTEGER(sdr);
279
    dr = INTEGER(sdr);
280
    er = INTEGER(ser);
280
    er = INTEGER(ser);
281
    wr = INTEGER(swr);
281
    wr = INTEGER(swr);
Line 285... Line 285...
285
    w = INTEGER(sw);
285
    w = INTEGER(sw);
286
 
286
 
287
    /* Determine the column widths */
287
    /* Determine the column widths */
288
 
288
 
289
    for (j = 0; j < c; j++) {
289
    for (j = 0; j < c; j++) {
290
	formatComplex(&x[j * r], r,
290
	formatComplex(&x[j * r], r, 
291
		      &wr[j], &dr[j], &er[j],
291
		      &wr[j], &dr[j], &er[j],
292
		      &wi[j], &di[j], &ei[j]);
292
		      &wi[j], &di[j], &ei[j]);
293
	if (!isNull(cl))
293
	if (!isNull(cl))
294
	    clabw = strlen(CHAR(STRING(cl)[j]));
294
	    clabw = strlen(CHAR(STRING(cl)[j]));
295
	else
295
	else
Line 324... Line 324...
324
	    MatrixRowLabel(rl, i, rlabw, lbloff);
324
	    MatrixRowLabel(rl, i, rlabw, lbloff);
325
	    for (j = jmin; j < jmax; j++) {
325
	    for (j = jmin; j < jmax; j++) {
326
		if (ISNA(x[i + j * r].r) || ISNA(x[i + j * r].i))
326
		if (ISNA(x[i + j * r].r) || ISNA(x[i + j * r].i))
327
		    Rprintf("%s", EncodeReal(NA_REAL, w[j], 0, 0));
327
		    Rprintf("%s", EncodeReal(NA_REAL, w[j], 0, 0));
328
		else
328
		else
-
 
329
		    Rprintf("%s", 
329
		    Rprintf("%s", EncodeComplex(x[i + j * r],
330
			    EncodeComplex(x[i + j * r],
330
						wr[j] + R_print.gap, dr[j], er[j],
331
					  wr[j] + R_print.gap, dr[j], er[j],
331
						wi[j], di[j], ei[j]));
332
					  wi[j], di[j], ei[j]));
332
	    }
333
	    }
333
	}
334
	}
334
	Rprintf("\n");
335
	Rprintf("\n");
335
	jmin = jmax;
336
	jmin = jmax;
336
    }
337
    }