SSE組み込み関数は、通常のfloatオペレーションが-1を返すようにします。＃INV

Question

私は、音声処理を実行する書き込み中のSSEメソッドに問題があります。私はここでインテルの論文に基づいてSSEランダムな機能を実装しました：

次のように私もまた、SSEを使用して、S16にフロートから変換を実行されているメソッドを持って

は、変換は非常に単純に行われ

http://software.intel.com/en-us/articles/fast-random-number-generator-on-the-intel-pentiumr-4-processor/

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; 
    ... 
} 

場合結果が返され、resultは未使用です。 これは最初のループで完全に実行されますが、次のループでdeltaは1.0ではなく-1.#INFになります。私が__m64 con = _mm_cvtps_pi16(in);という行をコメントアウトすると、問題は解消されます。

私はFPUが未知の状態になっていると思います。

Answer 1

。両方が同じレジスタ上で動作しているので、奇妙な結果を生むことがあります。

_mm_empty() 

FPUは正しい状態にリセットされます。マイクロソフトはGuidelines for When to Use EMMS

Answer 2

• _mm_load_psは、整列された負荷を行うことは保証されていません。 float *データは16バイトではなく4バイトにアライメントすることができます_ => _mm_loadu_ps
• memcpyはおそらくSSEで得られる利点を無効にしますが、__m64のストアコマンドを使用する必要があります。 アラインされていないストリームまたは__m64のストアを実行できない場合は、_m128iの内部に保持し、_mm_maskmoveu_si128でマスクされた書き込みを実行するか、手動で8バイトを格納します。 SSE整数演算及び（レギュラー）浮動小数点演算を混合

http://msdn.microsoft.com/en-us/library/bytwczae.aspx