nb_kernel332_ia32_sse2.intel_syntax.s
来自「最著名最快的分子模拟软件」· S 代码 · 共 2,183 行 · 第 1/5 页
S
2,183 行
movapd xmm7, xmm6 unpcklpd xmm6, xmm3 ;# G1 G2 unpckhpd xmm7, xmm3 ;# H1 H2 ;# Dispersion table ready, in xmm4-xmm7 mulpd xmm6, xmm1 ;# xmm6=Geps mulpd xmm7, xmm2 ;# xmm7=Heps2 addpd xmm5, xmm6 addpd xmm5, xmm7 ;# xmm5=Fp mulpd xmm7, [esp + nb332_two] ;# two*Heps2 addpd xmm7, xmm6 addpd xmm7, xmm5 ;# xmm7=FF mulpd xmm5, xmm1 ;# xmm5=eps*Fp addpd xmm5, xmm4 ;# xmm5=VV movapd xmm4, [esp + nb332_c12] mulpd xmm7, xmm4 ;# fijR mulpd xmm5, xmm4 ;# Vvdw12 addpd xmm7, [esp + nb332_fscal] addpd xmm5, [esp + nb332_Vvdwtot] movapd [esp + nb332_Vvdwtot], xmm5 xorpd xmm4, xmm4 mulpd xmm7, [esp + nb332_tsc] mulpd xmm7, [esp + nb332_rinvOO] subpd xmm4, xmm7 movapd xmm0, xmm4 movapd xmm1, xmm4 movapd xmm2, xmm4 xorpd xmm3, xmm3 movapd xmm4, xmm3 movapd xmm5, xmm3 mulpd xmm0, [esp + nb332_dxOO] mulpd xmm1, [esp + nb332_dyOO] mulpd xmm2, [esp + nb332_dzOO] subpd xmm3, xmm0 subpd xmm4, xmm1 subpd xmm5, xmm2 addpd xmm0, [esp + nb332_fixO] addpd xmm1, [esp + nb332_fiyO] addpd xmm2, [esp + nb332_fizO] movapd [esp + nb332_fjxO], xmm3 movapd [esp + nb332_fjyO], xmm4 movapd [esp + nb332_fjzO], xmm5 movapd [esp + nb332_fixO], xmm0 movapd [esp + nb332_fiyO], xmm1 movapd [esp + nb332_fizO], xmm2 ;# O-H1 interaction movapd xmm0, [esp + nb332_rinvOH1] movapd xmm1, xmm0 mulpd xmm1, [esp + nb332_rsqOH1] ;# xmm1=r mulpd xmm1, [esp + nb332_tsc] cvttpd2pi mm6, xmm1 ;# mm6 = lu idx cvtpi2pd xmm6, mm6 subpd xmm1, xmm6 ;# xmm1=eps movapd xmm2, xmm1 mulpd xmm2, xmm2 ;# xmm2=eps2 pslld mm6, 2 ;# idx *= 4 mov esi, [ebp + nb332_VFtab] movd eax, mm6 psrlq mm6, 32 movd ebx, mm6 ;# indices in eax/ebx lea eax, [eax + eax*2] ;# idx*=3 (12 total now) lea ebx, [ebx + ebx*2] ;# idx*=3 (12 total now) movlpd xmm4, [esi + eax*8] ;# Y1 movlpd xmm3, [esi + ebx*8] ;# Y2 movhpd xmm4, [esi + eax*8 + 8] ;# Y1 F1 movhpd xmm3, [esi + ebx*8 + 8] ;# Y2 F2 movapd xmm5, xmm4 unpcklpd xmm4, xmm3 ;# Y1 Y2 unpckhpd xmm5, xmm3 ;# F1 F2 movlpd xmm6, [esi + eax*8 + 16] ;# G1 movlpd xmm3, [esi + ebx*8 + 16] ;# G2 movhpd xmm6, [esi + eax*8 + 24] ;# G1 H1 movhpd xmm3, [esi + ebx*8 + 24] ;# G2 H2 movapd xmm7, xmm6 unpcklpd xmm6, xmm3 ;# G1 G2 unpckhpd xmm7, xmm3 ;# H1 H2 ;# coulomb table ready, in xmm4-xmm7 mulpd xmm6, xmm1 ;# xmm6=Geps mulpd xmm7, xmm2 ;# xmm7=Heps2 addpd xmm5, xmm6 addpd xmm5, xmm7 ;# xmm5=Fp mulpd xmm7, [esp + nb332_two] ;# two*Heps2 movapd xmm3, [esp + nb332_qqOH] addpd xmm7, xmm6 addpd xmm7, xmm5 ;# xmm7=FF mulpd xmm5, xmm1 ;# xmm5=eps*Fp addpd xmm5, xmm4 ;# xmm5=VV mulpd xmm5, xmm3 ;# vcoul=qq*VV mulpd xmm3, xmm7 ;# fijC=FF*qq ;# at this point mm5 contains vcoul and xmm3 fijC addpd xmm5, [esp + nb332_vctot] movapd [esp + nb332_vctot], xmm5 xorpd xmm1, xmm1 mulpd xmm3, [esp + nb332_tsc] mulpd xmm3, xmm0 subpd xmm1, xmm3 movapd xmm0, xmm1 movapd xmm2, xmm1 xorpd xmm3, xmm3 movapd xmm4, xmm3 movapd xmm5, xmm3 mulpd xmm0, [esp + nb332_dxOH1] mulpd xmm1, [esp + nb332_dyOH1] mulpd xmm2, [esp + nb332_dzOH1] subpd xmm3, xmm0 subpd xmm4, xmm1 subpd xmm5, xmm2 addpd xmm0, [esp + nb332_fixO] addpd xmm1, [esp + nb332_fiyO] addpd xmm2, [esp + nb332_fizO] movapd [esp + nb332_fjxH1], xmm3 movapd [esp + nb332_fjyH1], xmm4 movapd [esp + nb332_fjzH1], xmm5 movapd [esp + nb332_fixO], xmm0 movapd [esp + nb332_fiyO], xmm1 movapd [esp + nb332_fizO], xmm2 ;# O-H2 interaction movapd xmm0, [esp + nb332_rinvOH2] movapd xmm1, xmm0 mulpd xmm1, [esp + nb332_rsqOH2] ;# xmm1=r mulpd xmm1, [esp + nb332_tsc] cvttpd2pi mm6, xmm1 ;# mm6 = lu idx cvtpi2pd xmm6, mm6 subpd xmm1, xmm6 ;# xmm1=eps movapd xmm2, xmm1 mulpd xmm2, xmm2 ;# xmm2=eps2 pslld mm6, 2 ;# idx *= 4 mov esi, [ebp + nb332_VFtab] movd eax, mm6 psrlq mm6, 32 movd ebx, mm6 ;# indices in eax/ebx lea eax, [eax + eax*2] ;# idx*=3 (12 total now) lea ebx, [ebx + ebx*2] ;# idx*=3 (12 total now) movlpd xmm4, [esi + eax*8] ;# Y1 movlpd xmm3, [esi + ebx*8] ;# Y2 movhpd xmm4, [esi + eax*8 + 8] ;# Y1 F1 movhpd xmm3, [esi + ebx*8 + 8] ;# Y2 F2 movapd xmm5, xmm4 unpcklpd xmm4, xmm3 ;# Y1 Y2 unpckhpd xmm5, xmm3 ;# F1 F2 movlpd xmm6, [esi + eax*8 + 16] ;# G1 movlpd xmm3, [esi + ebx*8 + 16] ;# G2 movhpd xmm6, [esi + eax*8 + 24] ;# G1 H1 movhpd xmm3, [esi + ebx*8 + 24] ;# G2 H2 movapd xmm7, xmm6 unpcklpd xmm6, xmm3 ;# G1 G2 unpckhpd xmm7, xmm3 ;# H1 H2 ;# coulomb table ready, in xmm4-xmm7 mulpd xmm6, xmm1 ;# xmm6=Geps mulpd xmm7, xmm2 ;# xmm7=Heps2 addpd xmm5, xmm6 addpd xmm5, xmm7 ;# xmm5=Fp mulpd xmm7, [esp + nb332_two] ;# two*Heps2 movapd xmm3, [esp + nb332_qqOH] addpd xmm7, xmm6 addpd xmm7, xmm5 ;# xmm7=FF mulpd xmm5, xmm1 ;# xmm5=eps*Fp addpd xmm5, xmm4 ;# xmm5=VV mulpd xmm5, xmm3 ;# vcoul=qq*VV mulpd xmm3, xmm7 ;# fijC=FF*qq ;# at this point mm5 contains vcoul and xmm3 fijC addpd xmm5, [esp + nb332_vctot] movapd [esp + nb332_vctot], xmm5 xorpd xmm1, xmm1 mulpd xmm3, [esp + nb332_tsc] mulpd xmm3, xmm0 subpd xmm1, xmm3 movapd xmm0, xmm1 movapd xmm2, xmm1 xorpd xmm3, xmm3 movapd xmm4, xmm3 movapd xmm5, xmm3 mulpd xmm0, [esp + nb332_dxOH2] mulpd xmm1, [esp + nb332_dyOH2] mulpd xmm2, [esp + nb332_dzOH2] subpd xmm3, xmm0 subpd xmm4, xmm1 subpd xmm5, xmm2 addpd xmm0, [esp + nb332_fixO] addpd xmm1, [esp + nb332_fiyO] addpd xmm2, [esp + nb332_fizO] movapd [esp + nb332_fjxH2], xmm3 movapd [esp + nb332_fjyH2], xmm4 movapd [esp + nb332_fjzH2], xmm5 movapd [esp + nb332_fixO], xmm0 movapd [esp + nb332_fiyO], xmm1 movapd [esp + nb332_fizO], xmm2 ;# H1-O interaction movapd xmm0, [esp + nb332_rinvH1O] movapd xmm1, xmm0 mulpd xmm1, [esp + nb332_rsqH1O] ;# xmm1=r mulpd xmm1, [esp + nb332_tsc] cvttpd2pi mm6, xmm1 ;# mm6 = lu idx cvtpi2pd xmm6, mm6 subpd xmm1, xmm6 ;# xmm1=eps movapd xmm2, xmm1 mulpd xmm2, xmm2 ;# xmm2=eps2 pslld mm6, 2 ;# idx *= 4 mov esi, [ebp + nb332_VFtab] movd eax, mm6 psrlq mm6, 32 movd ebx, mm6 ;# indices in eax/ebx lea eax, [eax + eax*2] ;# idx*=3 (12 total now) lea ebx, [ebx + ebx*2] ;# idx*=3 (12 total now) movlpd xmm4, [esi + eax*8] ;# Y1 movlpd xmm3, [esi + ebx*8] ;# Y2 movhpd xmm4, [esi + eax*8 + 8] ;# Y1 F1 movhpd xmm3, [esi + ebx*8 + 8] ;# Y2 F2 movapd xmm5, xmm4 unpcklpd xmm4, xmm3 ;# Y1 Y2 unpckhpd xmm5, xmm3 ;# F1 F2 movlpd xmm6, [esi + eax*8 + 16] ;# G1 movlpd xmm3, [esi + ebx*8 + 16] ;# G2 movhpd xmm6, [esi + eax*8 + 24] ;# G1 H1 movhpd xmm3, [esi + ebx*8 + 24] ;# G2 H2 movapd xmm7, xmm6 unpcklpd xmm6, xmm3 ;# G1 G2 unpckhpd xmm7, xmm3 ;# H1 H2 ;# coulomb table ready, in xmm4-xmm7 mulpd xmm6, xmm1 ;# xmm6=Geps mulpd xmm7, xmm2 ;# xmm7=Heps2 addpd xmm5, xmm6 addpd xmm5, xmm7 ;# xmm5=Fp mulpd xmm7, [esp + nb332_two] ;# two*Heps2 movapd xmm3, [esp + nb332_qqOH] addpd xmm7, xmm6 addpd xmm7, xmm5 ;# xmm7=FF mulpd xmm5, xmm1 ;# xmm5=eps*Fp addpd xmm5, xmm4 ;# xmm5=VV mulpd xmm5, xmm3 ;# vcoul=qq*VV mulpd xmm3, xmm7 ;# fijC=FF*qq ;# at this point mm5 contains vcoul and xmm3 fijC addpd xmm5, [esp + nb332_vctot] movapd [esp + nb332_vctot], xmm5 xorpd xmm1, xmm1 mulpd xmm3, [esp + nb332_tsc] mulpd xmm3, xmm0 subpd xmm1, xmm3 movapd xmm0, xmm1 movapd xmm2, xmm1 movapd xmm3, [esp + nb332_fjxO] movapd xmm4, [esp + nb332_fjyO] movapd xmm5, [esp + nb332_fjzO] mulpd xmm0, [esp + nb332_dxH1O] mulpd xmm1, [esp + nb332_dyH1O] mulpd xmm2, [esp + nb332_dzH1O] subpd xmm3, xmm0 subpd xmm4, xmm1 subpd xmm5, xmm2 addpd xmm0, [esp + nb332_fixH1] addpd xmm1, [esp + nb332_fiyH1] addpd xmm2, [esp + nb332_fizH1] movapd [esp + nb332_fjxO], xmm3 movapd [esp + nb332_fjyO], xmm4 movapd [esp + nb332_fjzO], xmm5 movapd [esp + nb332_fixH1], xmm0 movapd [esp + nb332_fiyH1], xmm1 movapd [esp + nb332_fizH1], xmm2 ;# H1-H1 interaction movapd xmm0, [esp + nb332_rinvH1H1] movapd xmm1, xmm0 mulpd xmm1, [esp + nb332_rsqH1H1] ;# xmm1=r mulpd xmm1, [esp + nb332_tsc] cvttpd2pi mm6, xmm1 ;# mm6 = lu idx cvtpi2pd xmm6, mm6 subpd xmm1, xmm6 ;# xmm1=eps movapd xmm2, xmm1 mulpd xmm2, xmm2 ;# xmm2=eps2 pslld mm6, 2 ;# idx *= 4 mov esi, [ebp + nb332_VFtab] movd eax, mm6 psrlq mm6, 32 movd ebx, mm6 ;# indices in eax/ebx lea eax, [eax + eax*2] ;# idx*=3 (12 total now) lea ebx, [ebx + ebx*2] ;# idx*=3 (12 total now) movlpd xmm4, [esi + eax*8] ;# Y1 movlpd xmm3, [esi + ebx*8] ;# Y2 movhpd xmm4, [esi + eax*8 + 8] ;# Y1 F1 movhpd xmm3, [esi + ebx*8 + 8] ;# Y2 F2 movapd xmm5, xmm4 unpcklpd xmm4, xmm3 ;# Y1 Y2 unpckhpd xmm5, xmm3 ;# F1 F2 movlpd xmm6, [esi + eax*8 + 16] ;# G1 movlpd xmm3, [esi + ebx*8 + 16] ;# G2 movhpd xmm6, [esi + eax*8 + 24] ;# G1 H1 movhpd xmm3, [esi + ebx*8 + 24] ;# G2 H2 movapd xmm7, xmm6 unpcklpd xmm6, xmm3 ;# G1 G2 unpckhpd xmm7, xmm3 ;# H1 H2 ;# coulomb table ready, in xmm4-xmm7 mulpd xmm6, xmm1 ;# xmm6=Geps mulpd xmm7, xmm2 ;# xmm7=Heps2 addpd xmm5, xmm6 addpd xmm5, xmm7 ;# xmm5=Fp mulpd xmm7, [esp + nb332_two] ;# two*Heps2 movapd xmm3, [esp + nb332_qqHH] addpd xmm7, xmm6 addpd xmm7, xmm5 ;# xmm7=FF mulpd xmm5, xmm1 ;# xmm5=eps*Fp addpd xmm5, xmm4 ;# xmm5=VV mulpd xmm5, xmm3 ;# vcoul=qq*VV mulpd xmm3, xmm7 ;# fijC=FF*qq ;# at this point mm5 contains vcoul and xmm3 fijC addpd xmm5, [esp + nb332_vctot] movapd [esp + nb332_vctot], xmm5 xorpd xmm1, xmm1 mulpd xmm3, [esp + nb332_tsc] mulpd xmm3, xmm0 subpd xmm1, xmm3 movapd xmm0, xmm1 movapd xmm2, xmm1 movapd xmm3, [esp + nb332_fjxH1] movapd xmm4, [esp + nb332_fjyH1] movapd xmm5, [esp + nb332_fjzH1] mulpd xmm0, [esp + nb332_dxH1H1] mulpd xmm1, [esp + nb332_dyH1H1] mulpd xmm2, [esp + nb332_dzH1H1] subpd xmm3, xmm0 subpd xmm4, xmm1 subpd xmm5, xmm2 addpd xmm0, [esp + nb332_fixH1] addpd xmm1, [esp + nb332_fiyH1] addpd xmm2, [esp + nb332_fizH1] movapd [esp + nb332_fjxH1], xmm3 movapd [esp + nb332_fjyH1], xmm4 movapd [esp + nb332_fjzH1], xmm5 movapd [esp + nb332_fixH1], xmm0 movapd [esp + nb332_fiyH1], xmm1 movapd [esp + nb332_fizH1], xmm2 ;# H1-H2 interaction movapd xmm0, [esp + nb332_rinvH1H2] movapd xmm1, xmm0 mulpd xmm1, [esp + nb332_rsqH1H2] ;# xmm1=r mulpd xmm1, [esp + nb332_tsc] cvttpd2pi mm6, xmm1 ;# mm6 = lu idx cvtpi2pd xmm6, mm6 subpd xmm1, xmm6 ;# xmm1=eps movapd xmm2, xmm1 mulpd xmm2, xmm2 ;# xmm2=eps2 pslld mm6, 2 ;# idx *= 4 mov esi, [ebp + nb332_VFtab] movd eax, mm6 psrlq mm6, 32 movd ebx, mm6 ;# indices in eax/ebx lea eax, [eax + eax*2] ;# idx*=3 (12 total now) lea ebx, [ebx + ebx*2] ;# idx*=3 (12 total now) movlpd xmm4, [esi + eax*8] ;# Y1 movlpd xmm3, [esi + ebx*8] ;# Y2 movhpd xmm4, [esi + eax*8 + 8] ;# Y1 F1 movhpd xmm3, [esi + ebx*8 + 8] ;# Y2 F2 movapd xmm5, xmm4 unpcklpd xmm4, xmm3 ;# Y1 Y2 unpckhpd xmm5, xmm3 ;# F1 F2 movlpd xmm6, [esi + eax*8 + 16] ;# G1 movlpd xmm3, [esi + ebx*8 + 16] ;# G2 movhpd xmm6, [esi + eax*8 + 24] ;# G1 H1 movhpd xmm3, [esi + ebx*8 + 24] ;# G2 H2 movapd xmm7, xmm6 unpcklpd xmm6, xmm3 ;# G1 G2 unpckhpd xmm7, xmm3 ;# H1 H2 ;# coulomb table ready, in xmm4-xmm7 mulpd xmm6, xmm1 ;# xmm6=Geps mulpd xmm7, xmm2 ;# xmm7=Heps2 addpd xmm5, xmm6 addpd xmm5, xmm7 ;# xmm5=Fp mulpd xmm7, [esp + nb332_two] ;# two*Heps2 movapd xmm3, [esp + nb332_qqHH] addpd xmm7, xmm6 addpd xmm7, xmm5 ;# xmm7=FF mulpd xmm5, xmm1 ;# xmm5=eps*Fp addpd xmm5, xmm4 ;# xmm5=VV mulpd xmm5, xmm3 ;# vcoul=qq*VV mulpd xmm3, xmm7 ;# fijC=FF*qq ;# at this point mm5 contains vcoul and xmm3 fijC addpd xmm5, [esp + nb332_vctot] movapd [esp + nb332_vctot], xmm5 xorpd xmm1, xmm1 mulpd xmm3, [esp + nb332_tsc] mulpd xmm3, xmm0 subpd xmm1, xmm3 movapd xmm0, xmm1 movapd xmm2, xmm1 movapd xmm3, [esp + nb332_fjxH2] movapd xmm4, [esp + nb332_fjyH2] movapd xmm5, [esp + nb332_fjzH2]
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?