Does the definition of alphaNumerators[threadId] should be (float*)malloc(K * K * sizeof(float)) ?

[nk2014yj]

the dot product of two vectors(p-sk and tj) is a float, and in the function computeOptimalAlphaSubset, using memset(alphaNumerators[threadId], 0, K * K * sizeof(float)).