dm11=tmp(1,1) dm12=tmp(1,2) dm13=tmp(1,3) dm21=tmp(2,1) dm22=tmp(2,2) dm23 = tmp(2,3) dm31 =tmp(3,1) dm32 = tmp(3,2) dm33 = tmp(3,3) do j = 1,mn_9 l = ls(j) s1=dm11*p(j,1)+dm12*p(j,2)+dm13*p(j,3) s2=dm12*p(j,1)+dm22*p(j,2)+dm23*p(j,3) s3=dm13*p(j,1)+dm23*p(j,2)+dm33*p(j,3) do i= 1,j k=ls(i) sm(k,l)=sm(k,l)+s1*p(i,1)+s2*p(i,2)+s3*p(i,3) sm(l,k)=sm(k,l) end do end do
This is from some 80's code that ran on a Sun, does anyone see any good reason for the dm variables that are lost on me?