windlx实验-矩阵相乘

 ;********multiply an array to an array*****
;*******    Written by heiden     *********
;三个矩阵按行存储
;第一个矩阵的大小为：5*4

;第二个矩阵的大小为：4*3

;存放结果的矩阵大小为：5*3
;l 代表了第一个矩阵的当前被扫描的那行
;r4<=r  相对于result的偏移量
;r5<-l 
;r10<-i  第一个矩阵某一行的某个元素的索引
;r7<-j  第二个矩阵某一列的某个元素的索引

  .data
 

;**buffer of first array
  .global dat1
dat1: .space 20
;**buffer of second array
  .global dat2
dat2: .space 12
;**buffer of result
  .global result
result: .space 15

  .text
  .global main
main:
  addi r11,r0,0 ;temp
  addi r4,r0,0   ;r
  addi r5,r0,0   ;l
  addi r1,r0,5   ;

  movi2fp f1,r1
  addi r1,r0,4
  movi2fp f2,r1   ;n
  addi r1,r0,3
  movi2fp f3,r1   ;k
 
 
line: movi2fp f4,r5  ;判断是否扫描完
  ltf f4,f1  ;l<m?
  bfpf finish
  addi r6,r0,0  ;col
column: movi2fp f4,r6
  ltf f4,f3  ;col<k
  bfpf leveladd
 
  movfp2i r1,f2 
  multu r10,r5,r1   ;i<-l*n
  addi r7,r6,0  ;mov col to j
  addi r11,r0,0    ;temp=0
  addi r9,r5,1  ;r9<-l+1
  movfp2i r1,f2
  multu r9,r1,r9 ;p<-n*(l+1)
calculate: movi2fp f4,r10
  movi2fp f9,r9
  ltf f4,f9  ;compare i to p
  bfpf asign
  addi r1,r10,dat1
  lbu r2,0(r1)  ;??
  addi r1,r7,dat2
  lbu r3,0(r1)   ;??
  multu r1,r2,r3 ;temp1<-dat1[i]+dat2[j]
  add r11,r11,r1 ;temp<-temp+temp1
  addi r10,r10,1  ;i++
  movfp2i r1,f3  
  add r7,r7,r1   ;j<-j+k
  j calculate
asign:addi r1,r4,result
  sb 0(r1),r11  ;store result ??
  addi r4,r4,1
  addi r6,r6,1
  j column
leveladd:addi r5,r5,1
  j line
finish: trap 0  ;**end
 
 

此程序为实现5*4的矩阵与4*3的矩阵相乘.尚未使用循环展开技术