#include <string.h>
#define OUTFILENAMEFMT "OUT.%04d"

FILE* makeOutfile(int myID) {
    char name[strlen(OUTFILENAMEFMT)+1];
    snprintf(name, sizeof(name), OUTFILENAMEFMT, myID);
    return fopen(name, "w");
}

/* print A
   A is dimN by dimM (FIX THIS)
   strideA is as described for matmul 
 */
void printMatrix(FILE *stream, char *name, double *A, int dimNb, int dimMb, 
        int strideA, int myID_i, int myID_j) {
    int i, j;
    for (i = 0; i < dimNb; ++i) 
        for (j = 0; j < dimMb; ++j) {
            /* print A[i][j] */
            fprintf(stream, "%s[%04d][%04d] = %g\n", name, 
                    global(myID_i,dimNb,i), global(myID_j,dimMb,j), 
                    A[i*strideA + j]);
        }
}