c c ------------------- c subroutine atxb2c(a,b,c) c multiply a-TRANSPOSE by b to get c real a(3,3),b(3,3),c(3,3) c do 100, i=1,3 do 100, j=1,3 c(i,j)=0.0 do 100,ijk=1,3 c(i,j)=c(i,j)+a(ijk,i)*b(ijk,j) 100 continue return end c c ------------------- c