const char
*EncodeComplex(Rcomplex x, int wr, int dr, int er, int wi, int di, int ei,
const char *dec)
{
static char buff[NB];
/* IEEE allows signed zeros; strip these here */
if (x.r == 0.0) x.r = 0.0;
if (x.i == 0.0) x.i = 0.0;
if (ISNA(x.r) || ISNA(x.i)) {
snprintf(buff, NB,
"%*s", /* was "%*s%*s", R_print.gap, "", */
min(wr+wi+2, (NB-1)), CHAR(R_print.na_string));
} else {
char Re[NB];
const char *Im, *tmp;
int flagNegIm = 0;
Rcomplex y;
/* formatComplex rounded, but this does not, and we need to
keep it that way so we don't get strange trailing zeros.
But we do want to avoid printing small exponentials that
are probably garbage.
*/
z_prec_r(&y, &x, R_print.digits);
/* EncodeReal has static buffer, so copy */
tmp = EncodeReal0(y.r == 0. ? y.r : x.r, wr, dr, er, dec);
strcpy(Re, tmp);
if ( (flagNegIm = (x.i < 0)) ) x.i = -x.i;
Im = EncodeReal0(y.i == 0. ? y.i : x.i, wi, di, ei, dec);
snprintf(buff, NB, "%s%s%si", Re, flagNegIm ? "-" : "+", Im);
}
buff[NB-1] = '\0';
return buff;
}
static void z_prec(Rcomplex *r, Rcomplex *x, Rcomplex *p)
{
z_prec_r(r, x, p->r);
}
