Demo entry 4090640

md

   

Submitted by anonymous on Mar 18, 2016 at 15:01
Language: Fortran. Code size: 698 Bytes.

      !$omp parallel do private(i,j,k,xx,r2,r,fc), schedule(runtime)
      do 100 i=myrank,n-1,nprocs
         a(:,i)=0.0d0
         do 90 j=0,n-1
            if(i.ne.j) then
               r2=0.0d0
               do k=1,3
                  xx(k)=x(k,i)-x(k,j)
                  if(xx(k).gt.+halfl) xx(k)=xx(k)-xl
                  if(xx(k).lt.-halfl) xx(k)=xx(k)+xl
                  r2=r2+xx(k)*xx(k)
               enddo
               r=sqrt(r2)
               fc = zii(j)*exp(-xmuc*r)*(1./r+xmuc)/r2
               a(:,i) = a(:,i)+fc*xx(:)
            endif
   90    continue
         a(:,i) = (frp*zii(i)*vc*a(:,i))/(aii(i)*xmass)
  100 continue
      !$omp end parallel do

This snippet took 0.00 seconds to highlight.

Back to the Entry List or Home.

Delete this entry (admin only).