1 subroutine tce_sortacc_2(unsorted,sorted,a,b,i,j,factor) 2 implicit none 3 integer a,b 4 integer i,j 5 integer id(2),jd(2),ia,ib,j1,j2 6 double precision sorted(a*b) 7 double precision unsorted(a*b) 8 double precision factor 9 ia = 0 10 jd(1) = a 11 jd(2) = b 12 do j1 = 1,a 13 id(1) = j1 14 do j2 = 1,b 15 id(2) = j2 16 ia = ia + 1 17 ib = id(j)+jd(j)*(id(i)-1) 18 sorted(ib) = sorted(ib) + unsorted(ia) * factor 19 enddo 20 enddo 21 return 22 end 23