#include #include #define MAT_SIZE 1024 int main(int argc, char **argv) { double **a, **b, **c; int i,j,k; a = (double **)malloc(sizeof(double *)*MAT_SIZE); b = (double **)malloc(sizeof(double *)*MAT_SIZE); c = (double **)malloc(sizeof(double *)*MAT_SIZE); for (i = 0; i < MAT_SIZE; i++) { a[i] = (double *)malloc(sizeof(double)*MAT_SIZE); b[i] = (double *)malloc(sizeof(double)*MAT_SIZE); c[i] = (double *)malloc(sizeof(double)*MAT_SIZE); for (j = 0; j < MAT_SIZE; j++) { a[i][j] = i+j; b[i][j] = i-j; c[i][j] = 0.0; } } for (i = 0; i < MAT_SIZE; i++) { for (j = 0; j < MAT_SIZE; j++) { for (k = 0; k < MAT_SIZE; k++) { c[i][j] = c[i][j] + a[i][k] * b[k][j]; } } } return 0; }