Modified C++ code from
http://www.devarticles.com/c/a/Cplus...g-in-C-plus/4/
I haven't tested this, but it looks like it should work to me.
Calculates l * m, stores result in dest. m_numRows = number of rows in m and number of columns in l. numRows is number of rows in dest, numCols is number of cols in dest. Dimensions are assumed to be correct.
Code:
void matrixmul(double ** dest, double ** l, double ** m, int m_numRows, int numRows, int numCols) {
int r,c,i;
for (r=0;r<numRows;r++) {
for (c=0;c<m_numCols;c++) {
for (i=0;i<numCols;i++) {
dest[r][c]+=l[r][i]*m[i][c];
}
}
}
}
Parallel matrix multiplication
