nb_kernel314_ia32_sse2.intel_syntax.s

来自「最著名最快的分子模拟软件」· S 代码 · 共 2,193 行 · 第 1/5 页

S
2,193
字号
	mulpd  xmm2, [esp + nb314_c12]	movapd xmm3, xmm2	subpd  xmm3, xmm1	;# Vvdw=Vvdw12-Vvdw6 			addpd  xmm3, [esp + nb314_Vvdwtot]	mulpd  xmm1, [esp + nb314_six]	mulpd  xmm2, [esp + nb314_twelve]	subpd  xmm2, xmm1	mulpd  xmm2, xmm0	movapd [esp + nb314_Vvdwtot], xmm3	movapd xmm0, xmm2	movapd xmm1, xmm2 ;# fscal	xorpd xmm3, xmm3	movapd xmm4, xmm3	movapd xmm5, xmm3	mulpd xmm0, [esp + nb314_dxOO]	mulpd xmm1, [esp + nb314_dyOO]	mulpd xmm2, [esp + nb314_dzOO]	subpd xmm3, xmm0	subpd xmm4, xmm1	subpd xmm5, xmm2	addpd xmm0, [esp + nb314_fixO]	addpd xmm1, [esp + nb314_fiyO]	addpd xmm2, [esp + nb314_fizO]	movapd [esp + nb314_fjxO], xmm3	movapd [esp + nb314_fjyO], xmm4	movapd [esp + nb314_fjzO], xmm5	movapd [esp + nb314_fixO], xmm0	movapd [esp + nb314_fiyO], xmm1	movapd [esp + nb314_fizO], xmm2	;# H1-H1 interaction 	movapd xmm0, [esp + nb314_rinvH1H1]	movapd xmm1, xmm0	mulpd  xmm1, [esp + nb314_rsqH1H1] ;# xmm1=r 	mulpd  xmm1, [esp + nb314_tsc]	movd mm0, eax	movd mm1, ebx	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 + nb314_VFtab]	movd eax, mm6	psrlq mm6, 32	movd ebx, mm6		;# indices in eax/ebx 	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 + nb314_two]	;# two*Heps2 	movapd xmm3, [esp + nb314_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 + nb314_vctot]    	movapd [esp + nb314_vctot], xmm5	xorpd  xmm1, xmm1	mulpd  xmm3,  [esp + nb314_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 + nb314_dxH1H1]	mulpd xmm1, [esp + nb314_dyH1H1]	mulpd xmm2, [esp + nb314_dzH1H1]	subpd xmm3, xmm0	subpd xmm4, xmm1	subpd xmm5, xmm2	addpd xmm0, [esp + nb314_fixH1]	addpd xmm1, [esp + nb314_fiyH1]	addpd xmm2, [esp + nb314_fizH1]	movapd [esp + nb314_fjxH1], xmm3	movapd [esp + nb314_fjyH1], xmm4	movapd [esp + nb314_fjzH1], xmm5	movapd [esp + nb314_fixH1], xmm0	movapd [esp + nb314_fiyH1], xmm1	movapd [esp + nb314_fizH1], xmm2	;# H1-H2 interaction  	movapd xmm0, [esp + nb314_rinvH1H2]	movapd xmm1, xmm0	mulpd  xmm1, [esp + nb314_rsqH1H2] ;# xmm1=r 	mulpd  xmm1, [esp + nb314_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 + nb314_VFtab]	movd eax, mm6	psrlq mm6, 32	movd ebx, mm6		;# indices in eax/ebx 	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 + nb314_two]	;# two*Heps2 	movapd xmm3, [esp + nb314_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 + nb314_vctot]    	movapd [esp + nb314_vctot], xmm5	xorpd  xmm1, xmm1	mulpd  xmm3,  [esp + nb314_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 + nb314_dxH1H2]	mulpd xmm1, [esp + nb314_dyH1H2]	mulpd xmm2, [esp + nb314_dzH1H2]	subpd xmm3, xmm0	subpd xmm4, xmm1	subpd xmm5, xmm2	addpd xmm0, [esp + nb314_fixH1]	addpd xmm1, [esp + nb314_fiyH1]	addpd xmm2, [esp + nb314_fizH1]	movapd [esp + nb314_fjxH2], xmm3	movapd [esp + nb314_fjyH2], xmm4	movapd [esp + nb314_fjzH2], xmm5	movapd [esp + nb314_fixH1], xmm0	movapd [esp + nb314_fiyH1], xmm1	movapd [esp + nb314_fizH1], xmm2	;# H1-M interaction 	movapd xmm0, [esp + nb314_rinvH1M]	movapd xmm1, xmm0	mulpd  xmm1, [esp + nb314_rsqH1M] ;# xmm1=r 	mulpd  xmm1, [esp + nb314_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 + nb314_VFtab]	movd eax, mm6	psrlq mm6, 32	movd ebx, mm6		;# indices in eax/ebx 	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 + nb314_two]	;# two*Heps2 	movapd xmm3, [esp + nb314_qqMH]	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 + nb314_vctot]    	movapd [esp + nb314_vctot], xmm5	xorpd  xmm1, xmm1	mulpd  xmm3,  [esp + nb314_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 + nb314_dxH1M]	mulpd xmm1, [esp + nb314_dyH1M]	mulpd xmm2, [esp + nb314_dzH1M]	subpd xmm3, xmm0	subpd xmm4, xmm1	subpd xmm5, xmm2	addpd xmm0, [esp + nb314_fixH1]	addpd xmm1, [esp + nb314_fiyH1]	addpd xmm2, [esp + nb314_fizH1]	movapd [esp + nb314_fjxM], xmm3	movapd [esp + nb314_fjyM], xmm4	movapd [esp + nb314_fjzM], xmm5	movapd [esp + nb314_fixH1], xmm0	movapd [esp + nb314_fiyH1], xmm1	movapd [esp + nb314_fizH1], xmm2	;# H2-H1 interaction 	movapd xmm0, [esp + nb314_rinvH2H1]	movapd xmm1, xmm0	mulpd  xmm1, [esp + nb314_rsqH2H1] ;# xmm1=r 	mulpd  xmm1, [esp + nb314_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 + nb314_VFtab]	movd eax, mm6	psrlq mm6, 32	movd ebx, mm6		;# indices in eax/ebx 	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 + nb314_two]	;# two*Heps2 	movapd xmm3, [esp + nb314_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 + nb314_vctot]    	movapd [esp + nb314_vctot], xmm5	xorpd  xmm1, xmm1	mulpd  xmm3,  [esp + nb314_tsc]	mulpd  xmm3, xmm0	subpd  xmm1, xmm3	movapd xmm0, xmm1	movapd xmm2, xmm1		movapd xmm3, [esp + nb314_fjxH1]	movapd xmm4, [esp + nb314_fjyH1]	movapd xmm5, [esp + nb314_fjzH1]	mulpd xmm0, [esp + nb314_dxH2H1]	mulpd xmm1, [esp + nb314_dyH2H1]	mulpd xmm2, [esp + nb314_dzH2H1]	subpd xmm3, xmm0	subpd xmm4, xmm1	subpd xmm5, xmm2	addpd xmm0, [esp + nb314_fixH2]	addpd xmm1, [esp + nb314_fiyH2]	addpd xmm2, [esp + nb314_fizH2]	movapd [esp + nb314_fjxH1], xmm3	movapd [esp + nb314_fjyH1], xmm4	movapd [esp + nb314_fjzH1], xmm5	movapd [esp + nb314_fixH2], xmm0	movapd [esp + nb314_fiyH2], xmm1	movapd [esp + nb314_fizH2], xmm2	;# H2-H2 interaction 	movapd xmm0, [esp + nb314_rinvH2H2]	movapd xmm1, xmm0	mulpd  xmm1, [esp + nb314_rsqH2H2] ;# xmm1=r 	mulpd  xmm1, [esp + nb314_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 + nb314_VFtab]	movd eax, mm6	psrlq mm6, 32	movd ebx, mm6		;# indices in eax/ebx 	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 + nb314_two]	;# two*Heps2 	movapd xmm3, [esp + nb314_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 + nb314_vctot]    	movapd [esp + nb314_vctot], xmm5	xorpd  xmm1, xmm1	mulpd  xmm3,  [esp + nb314_tsc]	mulpd  xmm3, xmm0	subpd  xmm1, xmm3	movapd xmm0, xmm1	movapd xmm2, xmm1		movapd xmm3, [esp + nb314_fjxH2]	movapd xmm4, [esp + nb314_fjyH2]	movapd xmm5, [esp + nb314_fjzH2]	mulpd xmm0, [esp + nb314_dxH2H2]	mulpd xmm1, [esp + nb314_dyH2H2]	mulpd xmm2, [esp + nb314_dzH2H2]	subpd xmm3, xmm0	subpd xmm4, xmm1	subpd xmm5, xmm2	addpd xmm0, [esp + nb314_fixH2]	addpd xmm1, [esp + nb314_fiyH2]	addpd xmm2, [esp + nb314_fizH2]	movapd [esp + nb314_fjxH2], xmm3	movapd [esp + nb314_fjyH2], xmm4	movapd [esp + nb314_fjzH2], xmm5	movapd [esp + nb314_fixH2], xmm0	movapd [esp + nb314_fiyH2], xmm1	movapd [esp + nb314_fizH2], xmm2	;# H2-M interaction 	movapd xmm0, [esp + nb314_rinvH2M]	movapd xmm1, xmm0	mulpd  xmm1, [esp + nb314_rsqH2M] ;# xmm1=r 	mulpd  xmm1, [esp + nb314_tsc]	cvttpd2pi mm6, xmm1	;# mm6 = lu idx 	cvtpi2pd xmm6, mm6	subpd xmm1, xmm6	;# xmm1=eps 	movapd xmm2, xmm1		mulpd  xmm2, xmm2	;# xmm2=eps2 

⌨️ 快捷键说明

复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?