topic integral problem in Intel® Integrated Performance Primitives

integral problem

Youcef_K_ — Mon, 03 Jun 2013 15:19:37 GMT

Good morning dear all,

I'm currently using ippiSqrIntegral_8u32f64f_C1R to compute integral and square integral of an image. The square integral double array resulted seem to be right but it's not the case for the float integral array. So i decided to experiment the function in a simple 512x512 image entirely filled with 255 value. With this case, I expect that every value of integral array is gonna finish with "5" or "0" as it is a sum of 255 values... but i notice that some values can finish with "4" or "7", etc... errors occur between (126,508) and (513,508) .Here is my problem, because even if i'm wrong with the function paramters, it cannot gives float finishing with value other than "0" or "5". Thanks for any help

Here is my code :

#include <vector>

int main(){

std::vector<unsigned char> white(512*512,255);

Ipp8u* src = &white[0];

int step = 0;

Ipp32f* integral = ippiMalloc_32f_C1(512+1,521+1,&step);
Ipp64f* integral2 = (double*)malloc((512+1)*(512+1)*sizeof(double));

IppiSize roi = {w+2*R,h+2*R};
st = ippiSqrIntegral_8u32f64f_C1R(src,512*sizeof(Ipp8u),integral,(512+1)*sizeof(Ipp32f),integral2,(512+1)*sizeof(double),roi,0,0);

if(st!=ippStsNoErr)
{
std::cout << ippGetStatusString(st) << std::endl;
return 1;
}

for(int i=0 ; i<512+1 ; i++)

for(int j=0 ; j<512+1 ; j++)

{

if(*(integral + i*(512+1) + j)%5 != 0)

{

std::cout << i << "," << j << std::endl;

system("PAUSE");

}

return 0;

}

Sorry it is :

Youcef_K_ — Mon, 03 Jun 2013 15:21:25 GMT

Sorry it is :

IppiSize roi = {512,512};

What version ( update ) of

SergeyKostrov — Tue, 04 Jun 2013 01:03:37 GMT

What version ( update ) of IPP are you using? >>...some values can finish with "4" or "7", etc... errors occur between (126,508) and (513,508)... Please attach a txt-file with a complete output of the test application. Thanks.

I run w_ipp_7.1.1.119.exe to

Youcef_K_ — Tue, 04 Jun 2013 08:06:12 GMT

I run w_ipp_7.1.1.119.exe to install ipp.

Here attached, you'll find the .cpp and corresponding command window output in output.txt. I notice that the problem appears for square matrices bigger than 256x256.

thankls

thanks

Thanks for the output.txt

SergeyKostrov — Tue, 04 Jun 2013 12:56:07 GMT

Thanks for the output.txt file. However, I can not see your problem: >>...some values can finish with "4" or "7", etc... errors occur between (126,508) and (513,508)... and there are no any substrings 126,508 or 513,508 in the file.

I just copy-pasted the ouput

Youcef_K_ — Tue, 04 Jun 2013 13:24:18 GMT

I just copy-pasted the ouput of command window but the buffer of this window is not big enough to contain all value (it begins only from i=493,j=310). So, i now write directly values in .txt. Please find attached the corresponding output.txt and the .cpp file (where you can change the size of matrix and regenerate output.txt). I don't know why but now , first wrong value occurs at i=129,j=511 (not 126,508) and last one occurs at i=512,j=321 (not 513,508).

Thanks for you help.

Hi Sergey, any news about

Youcef_K_ — Fri, 07 Jun 2013 09:14:21 GMT

Hi Sergey, any news about integral function... thanks a lot for your help

Hi Youcef,

Sergey_K_Intel — Mon, 10 Jun 2013 06:39:47 GMT

Hi Youcef,

Due to limited width of Ipp32f fraction we cannot expect, that adding 255 to Ipp32f value we'll be getting multiple of 5 indefinitely. Look at the code:

[cpp]
Ipp32f sum = 0.;
int i;
for(i = 0; i < 512*512; i++) {
sum += (Ipp32f) 255;
if((sum/5 - (int)(sum/5)) != 0) {
printf("i=%d, sum=%f\n", i, sum);
break;
}
}
[/cpp]
We'll get:
i=65794, sum=16777724.000000
65794/512 = 128.5. I.e. that starting on 129 row you will be not getting 5 or 0 at the end of integral value. That's limitation. If you initialize your src with 5 (instead of 255), you'll get more or less good result.

Regards,
Sergey

Thanks Sergey, it sounds

Youcef_K_ — Mon, 10 Jun 2013 12:29:20 GMT

Thanks Sergey, it sounds clear for me now.

>>...Due to limited width of

SergeyKostrov — Fri, 28 Jun 2013 01:20:57 GMT

>>...Due to limited width of Ipp32f fraction we cannot expect, that adding 255 to Ipp32f value we'll be getting >>multiple of 5 indefinitely... I did a quick test and I changed 'float a' to 'double a' in a test case 'integral-ipp-test.cpp'. Take a look at results if interested and use Windiff to compare output files as follows: Windiff.exe

This is again truncation

Sergey_K_Intel — Fri, 28 Jun 2013 06:16:01 GMT

This is again truncation issue. The expression "a/5.0 - (int)(a/5.0)" is performed with different precision for 'float' and 'double' variables. So, the differences can take place.

Regards,
Sergey

>>...This is again truncation

SergeyKostrov — Fri, 28 Jun 2013 14:17:07 GMT

>>...This is again truncation issue... Sergey, This is Not a truncation issue and this is a precision issue ( 24-bit precision for 'float' vs. 53-bit precision for 'double' ). I simply wanted to stress on the right usage of that term and I understand that there are limitations with accuracy of floating-point calculations. Thanks.