- matrix_t *data=(matrix_t*)payload_data;
- matrix_t result;
- int i,j,k;
-
- /* 2. Make some room to return the result */
- result.lines = data[0].lines;
- result.rows = data[1].rows;
- result.ctn = xbt_malloc0(sizeof(double) * result.lines * result.rows);
-
- /* 3. Do the computation */
- for (i=0; i<result.lines; i++)
- for (j=0; j<result.rows; j++)
- for (k=0; k<data[1].lines; k++)
- result.ctn[i*result.rows + j] += data[0].ctn[i*result.rows +k] *data[1].ctn[k*result.rows +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]);