-
-/** \brief Display the content of a matrix (debugging purpose)
- * \param coords: boolean indicating whether we should add the coords of each cell to the output*/
-void xbt_matrix_dump(xbt_matrix_t matrix, const char*name, int coords,
- void_f_pvoid_t display_fun) {
- int i,j;
-
- fprintf(stderr,">>> Matrix %s dump (%d x %d)\n",name,matrix->lines,matrix->rows);
- for (i=0; i<matrix->lines; i++) {
- fprintf(stderr," ");
- for (j=0; j<matrix->rows; j++) {
- if (coords)
- fprintf(stderr," (%d,%d)=",i,j);
- else
- fprintf(stderr," ");
- display_fun(xbt_matrix_get_ptr(matrix,i,j));
+/** @brief Copy the values from the matrix src into the matrix dst
+ * @param dst: destination
+ * @param src: source
+ * @param lsize: number of lines to copy
+ * @param rsize: number of rows to copy
+ * @param lpos_dst: line offset on destination matrix
+ * @param rpos_dst: row offset on destination matrix
+ * @param lpos_src: line offset on destination matrix
+ * @param rpos_src: row offset on destination matrix
+ * @param cpy_f: the function to use to copy the elements over
+ */
+void xbt_matrix_copy_values(xbt_matrix_t dst, xbt_matrix_t src, unsigned int lsize, unsigned int rsize,
+ unsigned int lpos_dst, unsigned int rpos_dst,unsigned int lpos_src, unsigned int rpos_src,
+ pvoid_f_pvoid_t const cpy_f)
+{
+ XBT_DEBUG ("Copy a %ux%u submatrix from %ux%u(of %ux%u) to %ux%u (of %ux%u)",
+ lsize, rsize, lpos_src, rpos_src, src->lines, src->rows, lpos_dst, rpos_dst, dst->lines, dst->rows);
+
+ /* everybody knows that issue is between the chair and the screen (particularly in my office) */
+ xbt_assert(src->elmsize == dst->elmsize);
+ /* don't check free_f since the user may play weird games with this */
+
+ xbt_assert(lpos_src + lsize <= src->lines);
+ xbt_assert(rpos_src + rsize <= src->rows);
+
+ xbt_assert(lpos_dst + lsize <= dst->lines);
+ xbt_assert(rpos_dst + rsize <= dst->rows);
+
+ /* Lets get serious here */
+ for (unsigned i = 0; i < rsize; i++) {
+ if (cpy_f) {
+ for (unsigned j = 0; j < lsize; j++)
+ xbt_matrix_get_as(dst, j + lpos_dst, i + rpos_dst, void *) =
+ cpy_f(xbt_matrix_get_ptr(src, j + rpos_src, i + lpos_src));
+ } else {
+ memcpy(xbt_matrix_get_ptr(dst, lpos_dst, i + rpos_dst), xbt_matrix_get_ptr(src, lpos_src, i + rpos_src),
+ dst->elmsize * lsize);