快速sin 和 cos 函数计算

CEO    发表于2007年12月26日 15:02 阅读(27) 评论( 1) 分类： 技术资料

举报

void _SSE2_SinCos(float x, float* s, float* c)  // any x { __asm {   movss xmm0, x    movaps xmm7, xmm0    movss xmm1, _ps_am_inv_sign_mask    movss xmm2, _ps_am_sign_mask    movss xmm3, _ps_am_2_o_pi    andps xmm0, xmm1    andps xmm7, xmm2    mulss xmm0, xmm3       pxor xmm3, xmm3    movd xmm5, _epi32_1    movss xmm4, _ps_am_1       cvttps2dq xmm2, xmm0    pand xmm5, xmm2    movd xmm1, _epi32_2    pcmpeqd xmm5, xmm3    movd xmm3, _epi32_1    cvtdq2ps xmm6, xmm2    paddd xmm3, xmm2    pand xmm2, xmm1    pand xmm3, xmm1    subss xmm0, xmm6    pslld xmm2, (31 - 1)    minss xmm0, xmm4       mov  eax, s     // mov eax, [esp + 4 + 16]    mov  edx, c    // mov edx, [esp + 4 + 16 + 4]       subss xmm4, xmm0    pslld xmm3, (31 - 1)       movaps xmm6, xmm4    xorps xmm2, xmm7    movaps xmm7, xmm5    andps xmm6, xmm7    andnps xmm7, xmm0    andps xmm0, xmm5    andnps xmm5, xmm4    movss xmm4, _ps_sincos_p3    orps xmm6, xmm7    orps xmm0, xmm5    movss xmm5, _ps_sincos_p2       movaps xmm1, xmm0    movaps xmm7, xmm6    mulss xmm0, xmm0    mulss xmm6, xmm6    orps xmm1, xmm2    orps xmm7, xmm3    movaps xmm2, xmm0    movaps xmm3, xmm6    mulss xmm0, xmm4    mulss xmm6, xmm4    movss xmm4, _ps_sincos_p1    addss xmm0, xmm5    addss xmm6, xmm5    movss xmm5, _ps_sincos_p0    mulss xmm0, xmm2    mulss xmm6, xmm3    addss xmm0, xmm4    addss xmm6, xmm4    mulss xmm0, xmm2    mulss xmm6, xmm3    addss xmm0, xmm5    addss xmm6, xmm5    mulss xmm0, xmm1    mulss xmm6, xmm7       // use full stores since caller might reload with full loads    movss [eax], xmm0    movss [edx], xmm6    } } 比标准库的sin , cos 速度提高三倍