1
Download Global Arrays with hands-on time
Transcript
Matrix Multiply Example – C (cont.)
/* Multiply a and b to get c */
for(i=0; i < hi[0] - lo[0] + 1; i++) {
for(j=0; j < hi[1] - lo[1] + 1; j++) {
c[i][j] = 0.0;
for(k=0; k<dims[0]; k++)
c[i][j] = c[i][j] + a[i][k]*btrns[j][k];
}
}
/* Copy c back to g_c */
NGA_Put(g_c, lo, hi, c, ld);
verify(g_a, g_b, g_c, lo1, hi1, ld);
/* Deallocate arrays */
GA_Destroy(g_a);
GA_Destroy(g_b);
GA_Destroy(g_c);
