topic ippiFilter wrong output 64bit version (IPP 8.2) in Intel® Integrated Performance Primitives

ippiFilter wrong output 64bit version (IPP 8.2)

Alexandr_S_1 — Thu, 11 Feb 2016 14:28:04 GMT

Hi,

My problem is wrong output from function ippiFilter_32f_C1R for 64-bit version, for 32-bit version it works correctly

Here is my code:

Ipp32f* convolution(float **a, float **filter, int width, int height, int coreSize, int newWidth, int newHeight)
{
   Ipp32f *src1 = new Ipp32f[width * height];
   Ipp32f *src2 = new Ipp32f[coreSize * coreSize];
   Ipp32f *dst = new Ipp32f[newWidth * newHeight];
   IppiSize src1Size = { width, height };
   IppiSize src2Size = { coreSize, coreSize };
   IppiSize dstSize = { newWidth, newHeight };
   IppiPoint anchor = { coreSize - 1 , coreSize - 1 };

   int k = 0;
   for (int i = 0; i < height; i++)
   {
       for (int j = 0; j < width; j++)
       {
           src1 = a;
           k++;
       }
   }

   k = 0;
   for (int i = coreSize - 1; i >= 0; i--)
   {
       for (int j = coreSize - 1; j >= 0; j--)
       {
           src2 = filter;
           k++;
       }
   }

   ippiFilter_32f_C1R(src1, width * sizeof(Ipp32f), dst, newWidth * sizeof(Ipp32f), dstSize, src2, src2Size, anchor);
   delete []src1;
   delete []src2;
   return dst;
}

The difference between the outputs of this function appears, when width = 19, height = 24, coreSize = 3, newWidth = 17, newHeight = 22

Thanks,

Alexandr

Alexander, which version of

Gennady_F_Intel — Fri, 12 Feb 2016 14:08:50 GMT

Alexander, which version of IPP has been used? is that sequential or threaded mode? are there any specific CPU type where the problem happens?

ipp version 8.2 from Composer

Alexandr_S_1 — Fri, 12 Feb 2016 14:48:53 GMT

ipp version 8.2 from Composer XE 2015, sequential mode, microsoft visual c++ compliler, processor Intel Core i5-3450 CPU 3.10 GHz x64, 64bit operation system Windows 10.

Can you test code, which I sent on Win32 and x64 platform with function parameters width = 19, height = 24, coreSize = 3, newWidth = 17, newHeight = 22?

Hi Alexander.

Andrey_B_Intel — Mon, 15 Feb 2016 12:38:28 GMT

Hi Alexander.

Thanks for using IPP library and your test case.

Could you provide please input data values to understand difference? The small calculation error for float data and different CPU letters is normal situation for IPP. If difference is very big - it can be a bug.

Hi Andrey.

Alexandr_S_1 — Mon, 15 Feb 2016 14:06:13 GMT

Hi Andrey.

The most simple test:

Look at the function, I sent (first post):

Ipp32f* convolution(float **a, float **filter, int width, int height, int coreSize, int newWidth, int newHeight)

Let width = 19, height = 24, coreSize = 3, newWidth = 17, newHeight = 22, each element of "a" equals 1.0f, each element of "filter" equals 1.0f,

then result of ippiFilter_32f_C1R must be array of 17x22(374) elements, each element equals 9.

I attached result of this function, when I built my test application on x86 configuration and x64.

As you can see, when application was built on x64, some elements equals 9 and some elements uninitialized.

Alexandr, we see the problem

Gennady_F_Intel — Mon, 15 Feb 2016 16:06:35 GMT

Alexandr, we see the problem with IPP v.8.2.3 between 32 and 64 bits code. This problem has gone with the latest 9.0. update 1. You need slightly modify code because of some ippiFilter_32f_C1R from 8.2. has been replaced by ippiFilterBorderInit_32f+ippiFilterBorder_32f_C1R functions. See the code you need to modify:

#if 0 status = ippiFilter_32f_C1R(src1, width * sizeof(Ipp32f), dst, newWidth * sizeof(Ipp32f), dstSize, src2, src2Size, anchor);

#else

{ Ipp32f borderValue[1] = { 0.0f }; IppiFilterBorderSpec* pSpec; Ipp8u* pBuffer;

int specSize; int bufferSize; Ipp32f* pSrc; Ipp32f* src2Dir;

status = ippiFilterBorderGetSize(src2Size, dstSize, ipp32f, ipp32f, 1, &specSize, &bufferSize);

pSpec = (IppiFilterBorderSpec*)ippsMalloc_8u(specSize);

pBuffer = ippsMalloc_8u(bufferSize); src2Dir = ippsMalloc_32f(src2Size.width*src2Size.height);

for (k = 0; k < src2Size.width*src2Size.height; k++) { src2Dir = src2[src2Size.width*src2Size.height - 1 - k]; }

status = ippiFilterBorderInit_32f(src2Dir, src2Size, ipp32f, 1, ippRndNear, pSpec);

pSrc = src1 + ((src2Size.height - 1) / 2)*width + (src2Size.width - 1) / 2;

status = ippiFilterBorder_32f_C1R(pSrc, width * sizeof(Ipp32f), dst, newWidth * sizeof(Ipp32f), dstSize, ippBorderInMem, borderValue, pSpec,pBuffer);

ippsFree(pSpec); ippsFree(pBuffer); ippsFree(src2Dir);

}

#endif

Ok, Gennady, thanks.

Alexandr_S_1 — Tue, 16 Feb 2016 10:48:10 GMT

Ok, Gennady, thanks.