#include #include void InitBlock(int *a, int *b, int *c, int dim) { int len, ind; int i,j; srand(time(NULL)); len = dim*dim; for (ind = 0; ind < len; ind++) { a[ind] = (int)(rand()%1000)/100.0; c[ind] = 0.0; } for (i = 0; i < dim; i++) { for (j = 0; j < dim; j++) b[i*dim+j] = (i==j)? 1.0 : 0.0; } } void BlockMult(int* c, int* a, int* b, int dim) { int i,j,k; for (i = 0; i < dim; i++) for (j = 0; j < dim; j ++) for (k = 0; k < dim; k++) c[i*dim+j] += (a[i*dim+k] * b[k*dim+j]); }