- matrix_t *data=(matrix_t*)payload_data;
- matrix_t result;
- int i,j,k;
-
- /* 2. Make some room to return the result */
- result.rows = data[0].rows;
- result.cols = data[1].cols;
- result.ctn = xbt_malloc0(sizeof(double) * result.rows * result.cols);
-
- /* 3. Do the computation */
- for (i=0; i<result.rows; i++)
- for (j=0; j<result.cols; j++)
- for (k=0; k<data[1].rows; k++)
- result.ctn[i*result.cols + j] += data[0].ctn[i*result.cols +k] *data[1].ctn[k*result.cols +j];
+ xbt_matrix_t *request = (xbt_matrix_t*)payload_data;
+ xbt_matrix_t result;
+
+ /* 2. Do the computation */
+ result = xbt_matrix_double_new_mult(request[0], request[1]);