Le calcul numérique de haute performance - Université de Laval
Le calcul numérique de haute performance - Université de Laval
Le calcul numérique de haute performance - Université de Laval
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
subroutine shortrange(r,pairs,src,npairs,nmax,np)<br />
parameter (npt=10000)<br />
integer pairs(nmax,2)<br />
dimension src(nmax), r(np,3)<br />
common /short/ chi(0:npt), dchi(0:npt-1), dr<br />
do m=1,npairs<br />
i=pairs(m,1)<br />
j=pairs(m,2)<br />
dx=r(i,1)-r(j,1)<br />
dy=r(i,2)-r(j,2)<br />
dz=r(i,3)-r(j,3)<br />
r2=dx*dx+dy*dy+dz*dz<br />
r=sqrt(r2)<br />
d=r/dr<br />
k=int(d)<br />
c=d-k<br />
src(m)=chi(k)+c*dchi(k)<br />
enddo<br />
return<br />
end<br />
+ + �