The R Project SVN – R – Diff – /branches/R-exp-R5/src/main/complex.c

    against MSVCRT's pow, and gets (0+0i)^Y as 0+0i for all Y */
 static double complex mycpow (double complex X, double complex Y)
+{
     double complex Res; int k;
     if (X == 0.0) {
+	if(cimag(Y) == 0.0){
-	__real__ Res = R_pow(0.0, __real__ Y);
+	    __real__ Res = R_pow(0.0, __real__ Y);
-	__imag__ Res = 0.0;
+	    __imag__ Res = 0.0;
+	} else {
+	    __real__ Res = R_NaN;
+	    __imag__ Res = R_NaN;
+	}
     } else if (__imag__ Y == 0.0 && __real__ Y == (k = (int)__real__ Y)
 	       && abs(k) <= 65536) {
 	return R_cpow_n(X, k);
     } else {
 	double rho, r,i, theta;
 */
 static double complex mycpow (double complex X, double complex Y)
+{
     double iY = cimag(Y);
-    if (iY == 0) {/* X ^ <real> */
+    if (iY == 0.0) {/* X ^ <real> */
 	double complex Z; int k;
-	if (X == 0.) {
+	if (X == 0.0) {
 	    Z = R_pow(0., creal(Y));
 	} else if(creal(Y) == (k = (int)creal(Y)) && abs(k) <= 65536) {
 	    Z = R_cpow_n(X, k);
 	} else {
 	    Z = cpow(X, Y);

The R Project SVN R