home *** CD-ROM | disk | FTP | other *** search
- ;
- ; This program originally available on the Motorola DSP bulletin board.
- ; It is provided under a DISCLAMER OF WARRANTY available from
- ; Motorola DSP Operation, 6501 Wm. Cannon Drive W., Austin, Tx., 78735.
- ;
- ; [1x3][3x3]=[1x3] Matrix Multiplication
- ;
- ; Last Update 04 Feb 87 Version 1.0
- ;
- opt cex
- page 132,66,0,0
- ;
- ; matrix multiply
- ;
- ;
- org x:0
- ;
- mat_a
- dc .4 ;a(1,1)
- dc .3 ;a(1,2)
- dc -.6 ;a(1,3)
-
- org y:0
- mat_b
- dc -.4 ;b(1,1)
- dc .5 ;b(2,1)
- dc .6 ;b(3,1)
- dc .35 ;b(1,2)
- dc -.3 ;b(2,2)
- dc -.35 ;b(3,2)
- dc .15 ;b(1,3)
- dc .4 ;b(2,3)
- dc .45 ;b(3,3)
-
- mat_x ds 3
- ;
- ;
- org p:$100
- matmul move #mat_a,r0 ;point to A matrix
- move #mat_b,r4 ;point to B matrix
- move #2,m0 ;mod 3
- move #mat_x,r1 ;output X matrix
- ;
- move x:(r0)+,x0 y:(r4)+,y0
- mpy x0,y0,a x:(r0)+,x0 y:(r4)+,y0
- mac x0,y0,a x:(r0)+,x0 y:(r4)+,y0
- macr x0,y0,a x:(r0)+,x0 y:(r4)+,y0
- move a,y:(r1)+
- mpy x0,y0,a x:(r0)+,x0 y:(r4)+,y0
- mac x0,y0,a x:(r0)+,x0 y:(r4)+,y0
- macr x0,y0,a x:(r0)+,x0 y:(r4)+,y0
- move a,y:(r1)+
- mpy x0,y0,a x:(r0)+,x0 y:(r4)+,y0
- mac x0,y0,a x:(r0)+,x0 y:(r4)+,y0
- macr x0,y0,a x:(r0)+,x0 y:(r4)+,y0
- move a,y:(r1)+
- end
-