我正在编写一个SSE方法,用于音频处理,但遇到了问题。我已经根据英特尔的论文实现了基于SSE的随机函数:
如果
此外,我还有一个使用SSE进行从Float到S16转换的方法,转换方法如下:
unsigned int Float_S16LE(float *data, const unsigned int samples, uint8_t *dest)
{
int16_t *dst = (int16_t*)dest;
const __m128 mul = _mm_set_ps1((float)INT16_MAX);
__m128 rand;
const uint32_t even = count & ~0x3;
for(uint32_t i = 0; i < even; i += 4, data += 4, dst += 4)
{
/* random round to dither */
FloatRand4(-0.5f, 0.5f, NULL, &rand);
__m128 rmul = _mm_add_ps(mul, rand);
__m128 in = _mm_mul_ps(_mm_load_ps(data),rmul);
__m64 con = _mm_cvtps_pi16(in);
memcpy(dst, &con, sizeof(int16_t) * 4);
}
}
FloatRand4的定义如下:
static inline void FloatRand4(const float min, const float max, float result[4], __m128 *sseresult = NULL)
{
const float delta = (max - min) / 2.0f;
const float factor = delta / (float)INT32_MAX;
...
}
如果
sseresult != NULL
,则返回__m128
结果,而result
未使用。
这在第一个循环中表现得完美,但在下一个循环中,delta
变成了-1.#INF
,而不是1.0
。如果我注释掉__m64 con = _mm_cvtps_pi16(in);
这一行,问题就解决了。
我认为FPU进入了一个未知状态或其他什么状态。