The R Project SVN R

Rev

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

Rev 72021 Rev 72022
Line 161... Line 161...
161 
	n = XLENGTH(t);
 161 
	n = XLENGTH(t);
162 
	PROTECT(s = allocVector(CPLXSXP, n));
 162 
	PROTECT(s = allocVector(CPLXSXP, n));
163 
	setAttrib(s, R_NamesSymbol, getAttrib(t, R_NamesSymbol));
 163 
	setAttrib(s, R_NamesSymbol, getAttrib(t, R_NamesSymbol));
164 
	UNPROTECT(1);
 164 
	UNPROTECT(1);
165 
	if(n == 0) return s;
 165 
	if(n == 0) return s;
166 
	for (i = 0 ; i < n ; i++) {
 166 
	/* no need to initialize s, ccum* set all elements */
167 
	    COMPLEX(s)[i].r = NA_REAL;
 - 
 
168 
	    COMPLEX(s)[i].i = NA_REAL;
 - 
 
169 
	}
 - 
 
170 
	switch (PRIMVAL(op) ) {
 167 
	switch (PRIMVAL(op) ) {
171 
	case 1:	/* cumsum */
 168 
	case 1:	/* cumsum */
172 
	    return ccumsum(t, s);
 169 
	    return ccumsum(t, s);
173 
	    break;
 170 
	    break;
174 
	case 2: /* cumprod */
 171 
	case 2: /* cumprod */
Line 215... Line 212...
215 
	n = XLENGTH(t);
 212 
	n = XLENGTH(t);
216 
	PROTECT(s = allocVector(REALSXP, n));
 213 
	PROTECT(s = allocVector(REALSXP, n));
217 
	setAttrib(s, R_NamesSymbol, getAttrib(t, R_NamesSymbol));
 214 
	setAttrib(s, R_NamesSymbol, getAttrib(t, R_NamesSymbol));
218 
	UNPROTECT(2);
 215 
	UNPROTECT(2);
219 
	if(n == 0) return s;
 216 
	if(n == 0) return s;
220 
	for(i = 0 ; i < n ; i++) REAL(s)[i] = NA_REAL;
 217 
	/* no need to initialize s, cum* set all elements */
221 
	switch (PRIMVAL(op) ) {
 218 
	switch (PRIMVAL(op) ) {
222 
	case 1:	/* cumsum */
 219 
	case 1:	/* cumsum */
223 
	    return cumsum(t,s);
 220 
	    return cumsum(t,s);
224 
	    break;
 221 
	    break;
225 
	case 2: /* cumprod */
 222 
	case 2: /* cumprod */