The R Project SVN – R – Diff – /branches/R-gcontext/src/main/array.c

     int i, j;
     if (nr > 0 && nc > 0) {
         F77_CALL(dsyrk)(uplo, trans, &nc, &nr, &one, x, &nr, &zero, z, &nc);
 	for (i = 1; i < nc; i++)
 	    for (j = 0; j < i; j++) z[i + nc *j] = z[j + nc * i];
+    } else { /* zero-extent operations should return zeroes */
+	for(i = 0; i < nc*nc; i++) z[i] = 0;
+    }
+}
 static void crossprod(double *x, int nrx, int ncx,
 		      double *y, int nry, int ncy, double *z)
+{
     double one = 1.0, zero = 0.0;
     if (nrx > 0 && ncx > 0 && nry > 0 && ncy > 0) {
         F77_CALL(dgemm)(transa, transb, &ncx, &ncy, &nrx, &one,
 			x, &nrx, y, &nry, &zero, z, &ncx);
+    }
+    else { /* zero-extent operations should return zeroes */
+	int i;
+	for(i = 0; i < ncx*ncy; i++) z[i] = 0;
+    }
 #else
     int i, j, k;
     double xji, yjk, sum;
     for (i = 0; i < ncx; i++)
     one.r = 1.0; one.i = zero.r = zero.i = 0.0;
     if (nrx > 0 && ncx > 0 && nry > 0 && ncy > 0) {
         F77_CALL(zgemm)(transa, transb, &ncx, &ncy, &nrx, &one,
 			x, &nrx, y, &nry, &zero, z, &ncx);
+    }
+    else { /* zero-extent operations should return zeroes */
+	int i;
+	for(i = 0; i < ncx*ncy; i++) z[i].r = z[i].i = 0;
+    }
 #else
     int i, j, k;
     double xji_r, xji_i, yjk_r, yjk_i, sum_r, sum_i;
     for (i = 0; i < ncx; i++)

The R Project SVN R