I think it is a problem with matrices in R and C being stored in different row/column order
> void gowsim ( double *mat, int *OBJ, int *MATDESC) { > double x [*MATDESC][*OBJ]; > int i, j, nrow, ncol; > nrow = *OBJ; > ncol = *MATDESC; > > /* Rebuild Matrix */ > for (j=0; j < ncol; j++) { > for (i=0; i < nrow; i++) { > x[i][j] = *mat; Swapping x[i][j] to x[j][i] should fix the problem. Sorry, I did not test it. > Rprintf("row %d col %d value %f\n", i, j, x[i][j]); > mat++; > } > } > for (i=0; i< nrow; i++) { > Rprintf("%f %f %f %f\n", x[i][0], x[i][1], x[i][2], x[i][3]); > } >} Jarek ====================================================\======= Jarek Tuszynski, PhD. o / \ Science Applications International Corporation <\__,| (703) 676-4192 "> \ [EMAIL PROTECTED] ` \ ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html