🚀 go-pugleaf

In article <1990Feb10.160605.25254@ux1.cso.uiuc.edu> mcdonald@aries.scs.uiuc.edu (Doug McDonald) writes:

>	subroutine mm(A, B, C, n)
>        double precision a(n,n), b(n,n), c(n,n), t
>	do 1 i = 1, n
>	  do 2 j = 1, n
>	    t = 0.0
>	    do 3 k = 1, n
>	      t = t + B(i, k) * C(k, j)
>              A(i,j) = t
>3           continue 
>2         continue
>1       continue
>	end

>Surprised? I sure am. Why should the MIPS Fortran results be so slow?

The Fortran source is incorrect.  Correct is:

	subroutine mm(A, B, C, n)
        double precision a(n,n), b(n,n), c(n,n), t
	do 1 i = 1, n
	  do 2 j = 1, n
	    t = 0.0
	    do 3 k = 1, n
	      t = t + B(i, k) * C(k, j)
3           continue 
            A(i,j) = t
2         continue
1       continue
	end

Note the interchange of the store to A(i,j) and the continue.

Preston