[Starpu-devel] Question about matrix splicing

[christian eder <Christian.Eder@inria.fr>]

Dear all,

I am new to starpu and try to implement an easy matrix multiplication
to learn the interface. I have two matrices A and B and I assume that
B is already transposed. Trying to make 2-dimensional splicings of A
and B works, but when I try a 1-dimensional splicing. i.e.

task->handles[0] = starpu_data_get_sub_data(A_hdl, 1, j);

where j goes over all rows of A resp. the same for B I have problems
achieving the correct offsets in the tasks. Assume the pretty easy
(and of course not optimized) splicing in the following
fa.filter_func    = starpu_matrix_filter_block;
  fa.nchildren      = l;
  fb.filter_func    = starpu_matrix_filter_block;
  fb.nchildren      = n;
  fc.filter_func    = starpu_matrix_filter_vertical_block;
  fc.nchildren      = n;

  starpu_data_partition(a_hdl, &fa);
  starpu_data_partition(b_hdl, &fb);

  //starpu_data_map_filters(a_hdl, 2, &fa, &fc);
  //starpu_data_map_filters(b_hdl, 2, &fb, &fc);
  printf("- - - \n");
  starpu_data_map_filters(c_hdl, 2, &fa, &fc);
  printf("! ! !\n");


--
gpg: http://pgp.mit.edu:11371/pks/lookup?op=get&search=0x6E0D034B959FDE6A