nb_kernel214_ppc_altivec.c

来自「最著名最快的分子模拟软件」· C语言 代码 · 共 1,653 行 · 第 1/5 页

C
1,653 行
字号
12345
				jnrc            = jjnr[k+2];				jnrd            = jjnr[k+3];				j3a             = 3*jnra;				j3b             = 3*jnrb;				j3c             = 3*jnrc;				j3d             = 3*jnrd;				load_4_4atoms(pos+j3a,pos+j3b,pos+j3c,pos+j3d,							  &jx1,&jy1,&jz1,&jx2,&jy2,&jz2,&jx3,&jy3,&jz3,							  &jx4,&jy4,&jz4);				dx11            = vec_sub(ix1,jx1);				dy11            = vec_sub(iy1,jy1);				dz11            = vec_sub(iz1,jz1);				dx22            = vec_sub(ix2,jx2);				dy22            = vec_sub(iy2,jy2);				dz22            = vec_sub(iz2,jz2);				dx23            = vec_sub(ix2,jx3);				dy23            = vec_sub(iy2,jy3);				dz23            = vec_sub(iz2,jz3);				dx24            = vec_sub(ix2,jx4);				dy24            = vec_sub(iy2,jy4);				dz24            = vec_sub(iz2,jz4);				dx32            = vec_sub(ix3,jx2);				dy32            = vec_sub(iy3,jy2);				dz32            = vec_sub(iz3,jz2);				dx33            = vec_sub(ix3,jx3);				dy33            = vec_sub(iy3,jy3);				dz33            = vec_sub(iz3,jz3);				dx34            = vec_sub(ix3,jx4);				dy34            = vec_sub(iy3,jy4);				dz34            = vec_sub(iz3,jz4);				dx42            = vec_sub(ix4,jx2);				dy42            = vec_sub(iy4,jy2);				dz42            = vec_sub(iz4,jz2);				dx43            = vec_sub(ix4,jx3);				dy43            = vec_sub(iy4,jy3);				dz43            = vec_sub(iz4,jz3);				dx44            = vec_sub(ix4,jx4);				dy44            = vec_sub(iy4,jy4);				dz44            = vec_sub(iz4,jz4);				rsq11           = vec_madd(dx11,dx11,nul);				rsq22           = vec_madd(dx22,dx22,nul);				rsq23           = vec_madd(dx23,dx23,nul);				rsq24           = vec_madd(dx24,dx24,nul);				rsq32           = vec_madd(dx32,dx32,nul);				rsq33           = vec_madd(dx33,dx33,nul);				rsq34           = vec_madd(dx34,dx34,nul);				rsq42           = vec_madd(dx42,dx42,nul);				rsq43           = vec_madd(dx43,dx43,nul);				rsq44           = vec_madd(dx44,dx44,nul);				rsq11           = vec_madd(dy11,dy11,rsq11);				rsq22           = vec_madd(dy22,dy22,rsq22);				rsq23           = vec_madd(dy23,dy23,rsq23);				rsq24           = vec_madd(dy24,dy24,rsq24);				rsq32           = vec_madd(dy32,dy32,rsq32);				rsq33           = vec_madd(dy33,dy33,rsq33);				rsq34           = vec_madd(dy34,dy34,rsq34);				rsq42           = vec_madd(dy42,dy42,rsq42);				rsq43           = vec_madd(dy43,dy43,rsq43);				rsq44           = vec_madd(dy44,dy44,rsq44);     				rsq11           = vec_madd(dz11,dz11,rsq11);				rsq22           = vec_madd(dz22,dz22,rsq22);				rsq23           = vec_madd(dz23,dz23,rsq23);				rsq24           = vec_madd(dz24,dz24,rsq24);				rsq32           = vec_madd(dz32,dz32,rsq32);				rsq33           = vec_madd(dz33,dz33,rsq33);				rsq34           = vec_madd(dz34,dz34,rsq34);				rsq42           = vec_madd(dz42,dz42,rsq42);				rsq43           = vec_madd(dz43,dz43,rsq43);				rsq44           = vec_madd(dz44,dz44,rsq44);				rinvsq11        = do_recip(rsq11);				do_9_invsqrt(rsq22,rsq23,rsq24,							 rsq32,rsq33,rsq34,							 rsq42,rsq43,rsq44,							 &rinv22,&rinv23,&rinv24,							 &rinv32,&rinv33,&rinv34,							 &rinv42,&rinv43,&rinv44);      				rinvsix         = vec_madd(rinvsq11,rinvsq11,nul);				rinvsix         = vec_madd(rinvsix,rinvsq11,nul);				krsq22          = vec_madd(vkrf,rsq22,nul);				krsq23          = vec_madd(vkrf,rsq23,nul);				krsq24          = vec_madd(vkrf,rsq24,nul);				krsq32          = vec_madd(vkrf,rsq32,nul);				krsq33          = vec_madd(vkrf,rsq33,nul);				krsq34          = vec_madd(vkrf,rsq34,nul);				krsq42          = vec_madd(vkrf,rsq42,nul);				krsq43          = vec_madd(vkrf,rsq43,nul);				krsq44          = vec_madd(vkrf,rsq44,nul);				rinv22          = vec_add(rinv22,krsq22);				rinv23          = vec_add(rinv23,krsq23);				rinv24          = vec_add(rinv24,krsq24);				rinv32          = vec_add(rinv32,krsq32);				rinv33          = vec_add(rinv33,krsq33);				rinv34          = vec_add(rinv34,krsq34);				rinv42          = vec_add(rinv42,krsq42);				rinv43          = vec_add(rinv43,krsq43);				rinv44          = vec_add(rinv44,krsq44);				Vvdwtot          = vec_nmsub(c6,rinvsix,Vvdwtot);				Vvdwtot          = vec_madd(c12,vec_madd(rinvsix,rinvsix,nul),										   Vvdwtot);				rinv22          = vec_sub(rinv22,vcrf);				rinv23          = vec_sub(rinv23,vcrf);				rinv24          = vec_sub(rinv24,vcrf);				rinv32          = vec_sub(rinv32,vcrf);				rinv33          = vec_sub(rinv33,vcrf);				rinv34          = vec_sub(rinv34,vcrf);				rinv42          = vec_sub(rinv42,vcrf);				rinv43          = vec_sub(rinv43,vcrf);				rinv44          = vec_sub(rinv44,vcrf);      				vctot           = vec_madd(qqHH,rinv22,vctot);				vctot           = vec_madd(qqHH,rinv23,vctot);				vctot           = vec_madd(qqMH,rinv24,vctot);				vctot           = vec_madd(qqHH,rinv32,vctot);				vctot           = vec_madd(qqHH,rinv33,vctot);				vctot           = vec_madd(qqMH,rinv34,vctot);				vctot           = vec_madd(qqMH,rinv42,vctot);				vctot           = vec_madd(qqMH,rinv43,vctot);				vctot           = vec_madd(qqMM,rinv44,vctot);			} 			if(k<(nj1-2)) {				jnra            = jjnr[k];				jnrb            = jjnr[k+1];				jnrc            = jjnr[k+2];				j3a             = 3*jnra;				j3b             = 3*jnrb;				j3c             = 3*jnrc;				load_3_4atoms(pos+j3a,pos+j3b,pos+j3c,							  &jx1,&jy1,&jz1,&jx2,&jy2,&jz2,&jx3,&jy3,&jz3,							  &jx4,&jy4,&jz4);				qqMMt           = vec_sld(qqMM,nul,4);				qqMHt           = vec_sld(qqMH,nul,4);				qqHHt           = vec_sld(qqHH,nul,4);				c6t             = vec_sld(c6,nul,4);				c12t            = vec_sld(c12,nul,4);				dx11            = vec_sub(ix1,jx1);				dy11            = vec_sub(iy1,jy1);				dz11            = vec_sub(iz1,jz1);				dx22            = vec_sub(ix2,jx2);				dy22            = vec_sub(iy2,jy2);				dz22            = vec_sub(iz2,jz2);				dx23            = vec_sub(ix2,jx3);				dy23            = vec_sub(iy2,jy3);				dz23            = vec_sub(iz2,jz3);				dx24            = vec_sub(ix2,jx4);				dy24            = vec_sub(iy2,jy4);				dz24            = vec_sub(iz2,jz4);				dx32            = vec_sub(ix3,jx2);				dy32            = vec_sub(iy3,jy2);				dz32            = vec_sub(iz3,jz2);				dx33            = vec_sub(ix3,jx3);				dy33            = vec_sub(iy3,jy3);				dz33            = vec_sub(iz3,jz3);				dx34            = vec_sub(ix3,jx4);				dy34            = vec_sub(iy3,jy4);				dz34            = vec_sub(iz3,jz4);				dx42            = vec_sub(ix4,jx2);				dy42            = vec_sub(iy4,jy2);				dz42            = vec_sub(iz4,jz2);				dx43            = vec_sub(ix4,jx3);				dy43            = vec_sub(iy4,jy3);				dz43            = vec_sub(iz4,jz3);				dx44            = vec_sub(ix4,jx4);				dy44            = vec_sub(iy4,jy4);				dz44            = vec_sub(iz4,jz4);				rsq11           = vec_madd(dx11,dx11,nul);				rsq22           = vec_madd(dx22,dx22,nul);				rsq23           = vec_madd(dx23,dx23,nul);				rsq24           = vec_madd(dx24,dx24,nul);				rsq32           = vec_madd(dx32,dx32,nul);				rsq33           = vec_madd(dx33,dx33,nul);				rsq34           = vec_madd(dx34,dx34,nul);				rsq42           = vec_madd(dx42,dx42,nul);				rsq43           = vec_madd(dx43,dx43,nul);				rsq44           = vec_madd(dx44,dx44,nul);				rsq11           = vec_madd(dy11,dy11,rsq11);				rsq22           = vec_madd(dy22,dy22,rsq22);				rsq23           = vec_madd(dy23,dy23,rsq23);				rsq24           = vec_madd(dy24,dy24,rsq24);				rsq32           = vec_madd(dy32,dy32,rsq32);				rsq33           = vec_madd(dy33,dy33,rsq33);				rsq34           = vec_madd(dy34,dy34,rsq34);				rsq42           = vec_madd(dy42,dy42,rsq42);				rsq43           = vec_madd(dy43,dy43,rsq43);				rsq44           = vec_madd(dy44,dy44,rsq44);     				rsq11           = vec_madd(dz11,dz11,rsq11);				rsq22           = vec_madd(dz22,dz22,rsq22);				rsq23           = vec_madd(dz23,dz23,rsq23);				rsq24           = vec_madd(dz24,dz24,rsq24);				rsq32           = vec_madd(dz32,dz32,rsq32);				rsq33           = vec_madd(dz33,dz33,rsq33);				rsq34           = vec_madd(dz34,dz34,rsq34);				rsq42           = vec_madd(dz42,dz42,rsq42);				rsq43           = vec_madd(dz43,dz43,rsq43);				rsq44           = vec_madd(dz44,dz44,rsq44);				rinvsq11        = do_recip(rsq11);				do_9_invsqrt(rsq22,rsq23,rsq24,							 rsq32,rsq33,rsq34,							 rsq42,rsq43,rsq44,							 &rinv22,&rinv23,&rinv24,							 &rinv32,&rinv33,&rinv34,							 &rinv42,&rinv43,&rinv44);				zero_highest_element_in_vector(&rinvsq11);				zero_highest_element_in_9_vectors(&rinv22,&rinv23,&rinv24,												  &rinv32,&rinv33,&rinv34,												  &rinv42,&rinv43,&rinv44);            				rinvsix         = vec_madd(rinvsq11,rinvsq11,nul);				rinvsix         = vec_madd(rinvsix,rinvsq11,nul);				krsq22          = vec_madd(vkrf,rsq22,nul);				krsq23          = vec_madd(vkrf,rsq23,nul);				krsq24          = vec_madd(vkrf,rsq24,nul);				krsq32          = vec_madd(vkrf,rsq32,nul);				krsq33          = vec_madd(vkrf,rsq33,nul);				krsq34          = vec_madd(vkrf,rsq34,nul);				krsq42          = vec_madd(vkrf,rsq42,nul);				krsq43          = vec_madd(vkrf,rsq43,nul);				krsq44          = vec_madd(vkrf,rsq44,nul);				rinv22          = vec_add(rinv22,krsq22);				rinv23          = vec_add(rinv23,krsq23);				rinv24          = vec_add(rinv24,krsq24);				rinv32          = vec_add(rinv32,krsq32);				rinv33          = vec_add(rinv33,krsq33);				rinv34          = vec_add(rinv34,krsq34);				rinv42          = vec_add(rinv42,krsq42);				rinv43          = vec_add(rinv43,krsq43);				rinv44          = vec_add(rinv44,krsq44);				Vvdwtot          = vec_nmsub(c6t,rinvsix,Vvdwtot);				Vvdwtot          = vec_madd(c12t,vec_madd(rinvsix,rinvsix,nul),										   Vvdwtot);				rinv22          = vec_sub(rinv22,vcrf);				rinv23          = vec_sub(rinv23,vcrf);				rinv24          = vec_sub(rinv24,vcrf);				rinv32          = vec_sub(rinv32,vcrf);				rinv33          = vec_sub(rinv33,vcrf);				rinv34          = vec_sub(rinv34,vcrf);				rinv42          = vec_sub(rinv42,vcrf);				rinv43          = vec_sub(rinv43,vcrf);				rinv44          = vec_sub(rinv44,vcrf);      				vctot           = vec_madd(qqHHt,rinv22,vctot);				vctot           = vec_madd(qqHHt,rinv23,vctot);				vctot           = vec_madd(qqMHt,rinv24,vctot);				vctot           = vec_madd(qqHHt,rinv32,vctot);				vctot           = vec_madd(qqHHt,rinv33,vctot);				vctot           = vec_madd(qqMHt,rinv34,vctot);				vctot           = vec_madd(qqMHt,rinv42,vctot);				vctot           = vec_madd(qqMHt,rinv43,vctot);				vctot           = vec_madd(qqMMt,rinv44,vctot);			} else if(k<(nj1-1)) {				jnra            = jjnr[k];				jnrb            = jjnr[k+1];				j3a             = 3*jnra;				j3b             = 3*jnrb;				load_2_4atoms(pos+j3a,pos+j3b,							  &jx1,&jy1,&jz1,&jx2,&jy2,&jz2,&jx3,&jy3,&jz3,							  &jx4,&jy4,&jz4);				qqMMt           = vec_sld(qqMM,nul,8);				qqMHt           = vec_sld(qqMH,nul,8);				qqHHt           = vec_sld(qqHH,nul,8);				c6t             = vec_sld(c6,nul,8);				c12t            = vec_sld(c12,nul,8);				dx11            = vec_sub(ix1,jx1);				dy11            = vec_sub(iy1,jy1);				dz11            = vec_sub(iz1,jz1);				dx22            = vec_sub(ix2,jx2);				dy22            = vec_sub(iy2,jy2);				dz22            = vec_sub(iz2,jz2);				dx23            = vec_sub(ix2,jx3);				dy23            = vec_sub(iy2,jy3);				dz23            = vec_sub(iz2,jz3);				dx24            = vec_sub(ix2,jx4);				dy24            = vec_sub(iy2,jy4);				dz24            = vec_sub(iz2,jz4);				dx32            = vec_sub(ix3,jx2);				dy32            = vec_sub(iy3,jy2);				dz32            = vec_sub(iz3,jz2);				dx33            = vec_sub(ix3,jx3);				dy33            = vec_sub(iy3,jy3);				dz33            = vec_sub(iz3,jz3);				dx34            = vec_sub(ix3,jx4);				dy34            = vec_sub(iy3,jy4);				dz34            = vec_sub(iz3,jz4);				dx42            = vec_sub(ix4,jx2);				dy42            = vec_sub(iy4,jy2);				dz42            = vec_sub(iz4,jz2);				dx43            = vec_sub(ix4,jx3);				dy43            = vec_sub(iy4,jy3);				dz43            = vec_sub(iz4,jz3);				dx44            = vec_sub(ix4,jx4);				dy44            = vec_sub(iy4,jy4);				dz44            = vec_sub(iz4,jz4);				rsq11           = vec_madd(dx11,dx11,nul);				rsq22           = vec_madd(dx22,dx22,nul);				rsq23           = vec_madd(dx23,dx23,nul);				rsq24           = vec_madd(dx24,dx24,nul);				rsq32           = vec_madd(dx32,dx32,nul);				rsq33           = vec_madd(dx33,dx33,nul);				rsq34           = vec_madd(dx34,dx34,nul);				rsq42           = vec_madd(dx42,dx42,nul);				rsq43           = vec_madd(dx43,dx43,nul);				rsq44           = vec_madd(dx44,dx44,nul);				rsq11           = vec_madd(dy11,
12345

nb_kernel214_ppc_altivec.c - 源码说明

本页面展示了「最著名最快的分子模拟软件」中的 nb_kernel214_ppc_altivec.c 源码文件，采用 C语言 编程语言编写，共 1,653 行代码。您可以在线阅读完整代码内容，也可以返回资源详情页下载完整源码包进行本地学习和开发。

虫虫开发者社区收录了大量与分子模拟相关的技术资源，包括源代码、技术文档、电路图等，是电子工程师和嵌入式开发者的专业学习平台。

⌨️ 快捷键说明

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