nb_kernel334_ppc_altivec.c
来自「最著名最快的分子模拟软件」· C语言 代码 · 共 1,610 行 · 第 1/5 页
C
1,610 行
fiz3 = vec_madd(fs32,dz32,fiz3); fix4 = vec_madd(fs42,dx42,fix4); fiy4 = vec_madd(fs42,dy42,fiy4); fiz4 = vec_madd(fs42,dz42,fiz4); fix2 = vec_madd(fs23,dx23,fix2); fiy2 = vec_madd(fs23,dy23,fiy2); fiz2 = vec_madd(fs23,dz23,fiz2); fix3 = vec_madd(fs33,dx33,fix3); fiy3 = vec_madd(fs33,dy33,fiy3); fiz3 = vec_madd(fs33,dz33,fiz3); fix4 = vec_madd(fs43,dx43,fix4); fiy4 = vec_madd(fs43,dy43,fiy4); fiz4 = vec_madd(fs43,dz43,fiz4); fix2 = vec_madd(fs24,dx24,fix2); fiy2 = vec_madd(fs24,dy24,fiy2); fiz2 = vec_madd(fs24,dz24,fiz2); fix3 = vec_madd(fs34,dx34,fix3); fiy3 = vec_madd(fs34,dy34,fiy3); fiz3 = vec_madd(fs34,dz34,fiz3); fix4 = vec_madd(fs44,dx44,fix4); fiy4 = vec_madd(fs44,dy44,fiy4); fiz4 = vec_madd(fs44,dz44,fiz4); fjx1 = vec_nmsub(fs11,dx11,nul); fjy1 = vec_nmsub(fs11,dy11,nul); fjz1 = vec_nmsub(fs11,dz11,nul); fjx2 = vec_nmsub(fs22,dx22,nul); fjy2 = vec_nmsub(fs22,dy22,nul); fjz2 = vec_nmsub(fs22,dz22,nul); fjx3 = vec_nmsub(fs23,dx23,nul); fjy3 = vec_nmsub(fs23,dy23,nul); fjz3 = vec_nmsub(fs23,dz23,nul); fjx4 = vec_nmsub(fs24,dx24,nul); fjy4 = vec_nmsub(fs24,dy24,nul); fjz4 = vec_nmsub(fs24,dz24,nul); fjx2 = vec_nmsub(fs32,dx32,fjx2); fjy2 = vec_nmsub(fs32,dy32,fjy2); fjz2 = vec_nmsub(fs32,dz32,fjz2); fjx3 = vec_nmsub(fs33,dx33,fjx3); fjy3 = vec_nmsub(fs33,dy33,fjy3); fjz3 = vec_nmsub(fs33,dz33,fjz3); fjx4 = vec_nmsub(fs34,dx34,fjx4); fjy4 = vec_nmsub(fs34,dy34,fjy4); fjz4 = vec_nmsub(fs34,dz34,fjz4); fjx2 = vec_nmsub(fs42,dx42,fjx2); fjy2 = vec_nmsub(fs42,dy42,fjy2); fjz2 = vec_nmsub(fs42,dz42,fjz2); fjx3 = vec_nmsub(fs43,dx43,fjx3); fjy3 = vec_nmsub(fs43,dy43,fjy3); fjz3 = vec_nmsub(fs43,dz43,fjz3); fjx4 = vec_nmsub(fs44,dx44,fjx4); fjy4 = vec_nmsub(fs44,dy44,fjy4); fjz4 = vec_nmsub(fs44,dz44,fjz4); add_force_to_4_4atoms(faction+j3a,faction+j3b, faction+j3c,faction+j3d, fjx1,fjy1,fjz1,fjx2,fjy2,fjz2, fjx3,fjy3,fjz3,fjx4,fjy4,fjz4); } 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); zero_highest_element_in_vector(&rsq11); zero_highest_element_in_9_vectors(&rsq22,&rsq23,&rsq24, &rsq32,&rsq33,&rsq34, &rsq42,&rsq43,&rsq44); rinv11 = do_invsqrt(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(&rinv11); zero_highest_element_in_9_vectors(&rinv22,&rinv23,&rinv24, &rinv32,&rinv33,&rinv34, &rinv42,&rinv43,&rinv44); r11 = vec_madd(rsq11,rinv11,nul); r22 = vec_madd(rsq22,rinv22,nul); r23 = vec_madd(rsq23,rinv23,nul); r24 = vec_madd(rsq24,rinv24,nul); r32 = vec_madd(rsq32,rinv32,nul); r33 = vec_madd(rsq33,rinv33,nul); r34 = vec_madd(rsq34,rinv34,nul); r42 = vec_madd(rsq42,rinv42,nul); r43 = vec_madd(rsq43,rinv43,nul); r44 = vec_madd(rsq44,rinv44,nul); do_3_ljctable_lj(VFtab,vec_madd(r11,tsc,nul), &VVd,&FFd,&VVr,&FFr); do_3_ljctable_coul(VFtab,vec_madd(r22,tsc,nul),&VV22c,&FF22c); do_3_ljctable_coul(VFtab,vec_madd(r23,tsc,nul),&VV23c,&FF23c); do_3_ljctable_coul(VFtab,vec_madd(r24,tsc,nul),&VV24c,&FF24c); do_3_ljctable_coul(VFtab,vec_madd(r32,tsc,nul),&VV32c,&FF32c); do_3_ljctable_coul(VFtab,vec_madd(r33,tsc,nul),&VV33c,&FF33c); do_3_ljctable_coul(VFtab,vec_madd(r34,tsc,nul),&VV34c,&FF34c); do_3_ljctable_coul(VFtab,vec_madd(r42,tsc,nul),&VV42c,&FF42c); do_3_ljctable_coul(VFtab,vec_madd(r43,tsc,nul),&VV43c,&FF43c); do_3_ljctable_coul(VFtab,vec_madd(r44,tsc,nul),&VV44c,&FF44c); fs11 = vec_nmsub(c6t,FFd,nul); fs11 = vec_nmsub(c12t,FFr,fs11); fs22 = vec_nmsub(qqHHt,FF22c,nul); fs23 = vec_nmsub(qqHHt,FF23c,nul); fs24 = vec_nmsub(qqMHt,FF24c,nul); fs32 = vec_nmsub(qqHHt,FF32c,nul); fs33 = vec_nmsub(qqHHt,FF33c,nul); fs34 = vec_nmsub(qqMHt,FF34c,nul); fs42 = vec_nmsub(qqMHt,FF42c,nul); fs43 = vec_nmsub(qqMHt,FF43c,nul); fs44 = vec_nmsub(qqMMt,FF44c,nul); Vvdwtot = vec_madd(c6t,VVd,Vvdwtot); Vvdwtot = vec_madd(c12t,VVr,Vvdwtot); vctot = vec_madd(qqHHt,VV22c,vctot); vctot = vec_madd(qqHHt,VV23c,vctot); vctot = vec_madd(qqMHt,VV24c,vctot); vctot = vec_madd(qqHHt,VV32c,vctot); vctot = vec_madd(qqHHt,VV33c,vctot); vctot = vec_madd(qqMHt,VV34c,vctot); vctot = vec_madd(qqMHt,VV42c,vctot); vctot = vec_madd(qqMHt,VV43c,vctot); vctot = vec_madd(qqMMt,VV44c,vctot); fs11 = vec_madd(fs11,tsc,nul); fs22 = vec_madd(fs22,tsc,nul); fs23 = vec_madd(fs23,tsc,nul); fs24 = vec_madd(fs24,tsc,nul); fs32 = vec_madd(fs32,tsc,nul); fs33 = vec_madd(fs33,tsc,nul); fs34 = vec_madd(fs34,tsc,nul); fs42 = vec_madd(fs42,tsc,nul); fs43 = vec_madd(fs43,tsc,nul); fs44 = vec_madd(fs44,tsc,nul); fs11 = vec_madd(fs11,rinv11,nul); fs22 = vec_madd(fs22,rinv22,nul); fs23 = vec_madd(fs23,rinv23,nul); fs24 = vec_madd(fs24,rinv24,nul); fs32 = vec_madd(fs32,rinv32,nul); fs33 = vec_madd(fs33,rinv33,nul); fs34 = vec_madd(fs34,rinv34,nul); fs42 = vec_madd(fs42,rinv42,nul); fs43 = vec_madd(fs43,rinv43,nul); fs44 = vec_madd(fs44,rinv44,nul); fix1 = vec_madd(fs11,dx11,fix1); fiy1 = vec_madd(fs11,dy11,fiy1); fiz1 = vec_madd(fs11,dz11,fiz1); fix2 = vec_madd(fs22,dx22,fix2); fiy2 = vec_madd(fs22,dy22,fiy2); fiz2 = vec_madd(fs22,dz22,fiz2); fix3 = vec_madd(fs32,dx32,fix3); fiy3 = vec_madd(fs32,dy32,fiy3); fiz3 = vec_madd(fs32,dz32,fiz3); fix4 = vec_madd(fs42,dx42,fix4); fiy4 = vec_madd(fs42,dy42,fiy4); fiz4 = vec_madd(fs42,dz42,fiz4); fix2 = vec_madd(fs23,dx23,fix2); fiy2 = vec_madd(fs23,dy23,fiy2); fiz2 = vec_madd(fs23,dz23,fiz2); fix3 = vec_madd(fs33,dx33,fix3); fiy3 = vec_madd(fs33,dy33,fiy3); fiz3 = vec_madd(fs33,dz33,fiz3); fix4 = vec_madd(fs43,dx43,fix4); fiy4 = vec_madd(fs43,dy43,fiy4); fiz4 = vec_madd(fs43,dz43,fiz4); fix2 = vec_madd(fs24,dx24,fix2); fiy2 = vec_madd(fs24,dy24,fiy2); fiz2 = vec_madd(fs24,dz24,fiz2); fix3 = vec_madd(fs34,dx34,fix3); fiy3 = vec_madd(fs34,dy34,fiy3); fiz3 = vec_madd(fs34,dz34,fiz3); fix4 = vec_madd(fs44,dx44,fix4); fiy4 = vec_madd(fs44,dy44,fiy4); fiz4 = vec_madd(fs44,dz44,fiz4); fjx1 = vec_nmsub(fs11,dx11,nul); fjy1 = vec_nmsub(fs11,dy11,nul); fjz1 = vec_nmsub(fs11,dz11,nul); fjx2 = vec_nmsub(fs22,dx22,nul); fjy2 = vec_nmsub(fs22,dy22,nul); fjz2 = vec_nmsub(fs22,dz22,nul); fjx3 = vec_nmsub(fs23,dx23,nul); fjy3 = vec_nmsub(fs23,dy23,nul); fjz3 = vec_nmsub(fs23,dz23,nul); fjx4 = vec_nmsub(fs24,dx24,nul); fjy4 = vec_nmsub(fs24,dy24,nul); fjz4 = vec_nmsub(fs24,dz24,nul); fjx2 = vec_nmsub(fs32,dx32,fjx2); fjy2 = vec_nmsub(fs32,dy32,fjy2); fjz2 = vec_nmsub(fs32,dz32,fjz2); fjx3 = vec_nmsub(fs33,dx33,fjx3); fjy3 = vec_nmsub(fs33,dy33,fjy3); fjz3 = vec_nmsub(fs33,dz33,fjz3); fjx4 = vec_nmsub(fs34,dx34,fjx4); fjy4 = vec_nmsub(fs34,dy34,fjy4); fjz4 = vec_nmsub(fs34,dz34,fjz4); fjx2 = vec_nmsub(fs42,dx42,fjx2); fjy2 = vec_nmsub(fs42,dy42,fjy2); fjz2 = vec_nmsub(fs42,dz42,fjz2); fjx3 = vec_nmsub(fs43,dx43,fjx3); fjy3 = vec_nmsub(fs43,dy43,fjy3); fjz3 = vec_nmsub(fs43,dz43,fjz3); fjx4 = vec_nmsub(fs44,dx44,fjx4); fjy4 = vec_nmsub(fs44,dy44,fjy4); fjz4 = vec_nmsub(fs44,dz44,fjz4); add_force_to_3_4atoms(faction+j3a,faction+j3b,faction+j3c, fjx1,fjy1,fjz1,fjx2,fjy2,fjz2, fjx3,fjy3,fjz3,fjx4,fjy4,fjz4); } 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);⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?