The R Project SVN R

Rev

Rev 7002 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 7002 Rev 7192
Line 768... Line 768...
768 
	    if(!R_FINITE(COMPLEX(z)[i].r) || !R_FINITE(COMPLEX(z)[i].i))
 768 
	    if(!R_FINITE(COMPLEX(z)[i].r) || !R_FINITE(COMPLEX(z)[i].i))
769 
		errorcall(call, "invalid polynomial coefficient");
 769 
		errorcall(call, "invalid polynomial coefficient");
770 
	    REAL(zr)[degree-i] = COMPLEX(z)[i].r;
 770 
	    REAL(zr)[degree-i] = COMPLEX(z)[i].r;
771 
	    REAL(zi)[degree-i] = COMPLEX(z)[i].i;
 771 
	    REAL(zi)[degree-i] = COMPLEX(z)[i].i;
772 
	}
 772 
	}
773 
	F77_SYMBOL(cpoly)(REAL(zr), REAL(zi), &degree,
 - 
 
774 
			  REAL(rr), REAL(ri), &fail);
 773 
	cpoly(REAL(zr), REAL(zi), &degree, REAL(rr), REAL(ri), &fail);
775 
	if(fail) errorcall(call, "root finding code failed");
 774 
	if(fail) errorcall(call, "root finding code failed");
776 
	UNPROTECT(2);
 775 
	UNPROTECT(2);
777 
	r = allocVector(CPLXSXP, degree);
 776 
	r = allocVector(CPLXSXP, degree);
778 
	for(i=0 ; i<n ; i++) {
 777 
	for(i=0 ; i<n ; i++) {
779 
	    COMPLEX(r)[i].r = REAL(rr)[i];
 778 
	    COMPLEX(r)[i].r = REAL(rr)[i];