O.K i see. Is ist possible that domainSize doesn't work with a reduce kernel? Because the following program always calculates the whole sum of the array not only the first column.
reduce kernel void SumRed(float a<>, reduce float b<>
{
b+=a;
}
const int MAX = 4;
float *zahlen = new float[MAX*MAX];
for(int i=0;i<MAX;i++)
{
for(int j=0;j<MAX;j++)
{
zahlen[i*MAX+j] = i*MAX+j;
cout << zahlen[i*MAX+j] << " ";
}
cout << "\n";
}
const int rank = 2;
unsigned int dims[] = { MAX, MAX };
::brook::Stream<float>s1(rank,dims);
s1.read(zahlen);
unsigned int dimsR[] = { 1, 1 };
::brook::Stream<float>s2(rank,dimsR);
SumRed.domainSize(uint4(1,MAX,0,0));
SumRed(s1,s2);
float summe = 0.0f;
s2.write(&summe);
cout << "Sum: " << summe << "\n";
delete[] zahlen;
Hi,
i've another problem, what i try to do is to run a kernel only for the first column, after that i need all the data for a another kernel call. If i run the following program the values of the first column are set to 1.0, this is what i want. But the other values are changed after the kernel call. I thought if i use domainOffset and domainSize the whole Data is copied to the GPU and the kernel is called for the domain and the other data is untouched.
This is the data how it looks like in the beginning:
0 0 0
0 0 0
0 0 0
this is how it should look like after the kernel call
1 0 0
1 0 0
1 0 0
this is how it looks,
1 0.5e 0.6e
1 45 42e
1 4e 4e
kernel void SumRed(float a<>, out float b<>
{
b=1.0f;
}
const int MAX = 3;
float *zahlen = new float[MAX*MAX];
for(int i=0;i<MAX;i++)
{
for(int j=0;j<MAX;j++)
{
zahlen[i*MAX+j] =0.0f;
cout << zahlen[i*MAX+j] << " ";
}
cout << "\n";
}
cout << "\n";
const int rank = 2;
unsigned int dims[] = { MAX, MAX };
::brook::Stream<float>s1(rank,dims);
::brook::Stream<float>s2(rank,dims);
SumRed.domainOffset(uint4(0,0,0,0));
SumRed.domainSize(uint4(1,MAX,0,0));
s1.read(zahlen);
SumRed(s1,s2);
s2.write(zahlen);
for(int i=0;i<MAX;i++)
{
for(int j=0;j<MAX;j++)
{
cout << zahlen[i*MAX+j] << " ";
}
cout << "\n";
}
