topic ippsWinKaiser is not thread safe in Intel® Integrated Performance Primitives

ippsWinKaiser is not thread safe

Oliver1971 — Mon, 15 Apr 2024 14:37:16 GMT

Hello everyone,

We found that the ippsWinKaiser_32f_I() method is not thread safe, an example is attached to this message.

In our example we create a sinus wave as input vector and run the ippsWinKaiser_32f_I() method and store the result as "expected" result. We then launch multiple threads running the ippsWinKaiser_32f_I() method again with the same input vector and compare the "expected" result against the new computation. We abort the test routine if any result does not match.

The test program fails every time until we protect the call to ippsWinKaiser_32f_I() with a mutex. Without mutex the example always fails, with mutex all tests were good.

Any suggestions how to fix this or is this a bug?

Best regards,

Oliver

P.S.: ippsWinHann_32f_I() and ippsWinHamming_32f_I() are working as expected

Re:ippsWinKaiser is not thread safe

Gennady_F_Intel — Tue, 16 Apr 2024 12:40:38 GMT

This is an unexpected behavior as IPP is not internally threaded library.

I didn't run this example, but I see that you compare like as follows: computed != expected

which is not correct in the case of floating point data.

Re: ippsWinKaiser is not thread safe

Oliver1971 — Tue, 16 Apr 2024 13:19:56 GMT

Ok, I asked Microsofts KI for better comparision methods of two std vectors of floats:

bool compareVectors(const std::vector<float>&vec1, const std::vector<float>&vec2) {
    if (vec1.size() != vec2.size()) {
        return false;
    }

    float epsilon = 0.1f;
    for (size_t i = 0; i < vec1.size(); ++i) {
       if (std::abs(vec1[i] - vec2[i]) > epsilon) {
          return false;
       }
    }

return true;
}

Now I observe the following difference (which disappears once I mutex lock the calls):

Left expected (correct), right when error occurs. I even once saw an INF

Re: ippsWinKaiser is not thread safe

Oliver1971 — Tue, 16 Apr 2024 13:52:56 GMT

Ok, I asked the KI how to properly compare two std::vectors of type float:

bool compareVectors(const std::vector<float>&vec1, const std::vector<float>&vec2) {
    if (vec1.size() != vec2.size()) {
       return false;
    }

    float epsilon = 0.1f;
    for (size_t i = 0; i < vec1.size(); ++i) {
       if (std::abs(vec1[i] - vec2[i]) > epsilon) {
          return false;
       }
    }

return true;
}

I still get the errors (left expected, right new computation):

Re: ippsWinKaiser is not thread safe

Gennady_F_Intel — Wed, 17 Apr 2024 07:47:02 GMT

Oliver,

yes, we have to confirmed that we see the same behavior as you reported.

The problem will be escalated internally and we will keep you informed with the status.

--Gennady

Re:ippsWinKaiser is not thread safe

Ruqiu_C_Intel — Wed, 09 Oct 2024 01:45:16 GMT

The fix will be available in an upcoming release.

We will no longer follow this post, if you still experience the same issue in a new release, please post a new one and link this one.