PPT Slide
mydgemm:
do kb = 1, kk, blk
ke = min(kb+blk-1,kk)
do ib = 1, ii, blk
ie = min(ib+blk-1, ii)
do i = ib,ie
do k = kb, ke
T(k-kb+1, i-ib+1) = B(i,k)
enddo
enddo
do jb = 1, jj, blk
je = min(jb+blk-1, jj)
do j = jb, je, 2
do i = ib, ie, 2
T1 = 0.0d0
T2 = 0.0d0
T3 = 0.0d0
T4 = 0.0d0
do k = kb, ke
T1 = T1 + T (k-kb+1,i-ib+1)*C(k,j)
T2 = T2 + T(k-kb+1, i-ib+2)*C(k,j)
T3 = T3 + T(k-kb+1, i-ib+1)*C(k,j+1)
T4 = T4 +T(k-kb+1, i-ib+2)*C(k,j+1)
enddo
A(i,j) = A(i,j)+T1
A(i+1,j) = A(i+1, j)+T2
A(i,j+1) = A(i, j+1)+T3
A(i+1, j+1) = A(i+1, j+1) +T4
enddo
enddo
enddo
enddo
enddo