The R Project SVN R

Rev

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

Rev 45771 Rev 45774
Line 977... Line 977...
977 
    for(i = 0; i < n; i++) {
 977 
    for(i = 0; i < n; i++) {
978 
	if(COMPLEX(z)[i].r!= 0.0 || COMPLEX(z)[i].i != 0.0) degree = i;
 978 
	if(COMPLEX(z)[i].r!= 0.0 || COMPLEX(z)[i].i != 0.0) degree = i;
979 
    }
 979 
    }
980 
    n = degree + 1; /* omit trailing zeroes */
 980 
    n = degree + 1; /* omit trailing zeroes */
981 
    if(degree >= 1) {
 981 
    if(degree >= 1) {
982 
	if(n > 50) error(_("polynomial degree too high (49 max)"));
 - 
 
983 
	/* <==>	 #define NMAX 50  in  ../appl/cpoly.c */
 - 
 
984 
 
 - 
 
985 
	/* if(COMPLEX(z)[n-1].r == 0.0 && COMPLEX(z)[n-1].i == 0.0)
 - 
 
986 
	   error("highest power has coefficient 0");*/
 - 
 
987 
 
 - 
 
988 
	PROTECT(rr = allocVector(REALSXP, n));
 982 
	PROTECT(rr = allocVector(REALSXP, n));
989 
	PROTECT(ri = allocVector(REALSXP, n));
 983 
	PROTECT(ri = allocVector(REALSXP, n));
990 
	PROTECT(zr = allocVector(REALSXP, n));
 984 
	PROTECT(zr = allocVector(REALSXP, n));
991 
	PROTECT(zi = allocVector(REALSXP, n));
 985 
	PROTECT(zi = allocVector(REALSXP, n));
992 
 
 986