+/** \brief Copy the values from the matrix src into the matrix dst
+ * \param dest: 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
+ */
+void xbt_matrix_copy_values(xbt_matrix_t dst, xbt_matrix_t src,
+ int lsize, int rsize,
+ int lpos_dst,int rpos_dst,
+ int lpos_src,int rpos_src,
+ pvoid_f_pvoid_t *const cpy_f) {
+ int i,j;
+
+ DEBUG10("Copy a %dx%d submatrix from %dx%d(of %dx%d) to %dx%d (of %dx%d)",
+ 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 (particulary 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 (i=0;i<rsize;i++) {
+ if (cpy_f) {
+ for (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);
+ }
+ }
+
+}