纯净版fasm 浮点矩阵乘SSE_fasm浮点计算程序-优快云博客

这是一个使用fasm编写的浮点矩阵乘法示例，通过SSE指令集进行优化。代码中包含了一个1*4向量乘以4*4矩阵的函数MatrixMul，实现了矩阵乘法的计算过程。程序调用了自定义的MatrixMul函数，输出了浮点运算结果，并使用了printf和system函数进行输出和暂停操作。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
;*--==--* fasm console float mul example. By G-Spider
;*--==--* 编译: fasm.exe MatrixMul.asm MatrixMul.exe
;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
format PE console                    ;指定为PE控制台程序
entry start                          ;指定入口点为start

;---------------------------------------------
;//代码节，命名为'.text'(可读可执行)
section '.text' code readable executable

start :
;--------函数调用---------------
     push     fc                      ;地址入栈
     push     fb
     push     fa
     call     MatrixMul               ;自实现函数调用

;-------输出浮点结果------------
     cvtps2pd xmm0 , qword [fc + 4 * 0 ]     ;转成双精度浮点
     cvtps2pd xmm1 , qword [fc + 4 * 2 ]
     sub       esp , 16 * 2

     movups    [ esp ], xmm0              ;参数入栈
     movups    [ esp + 16 ], xmm1
     push     szFmt
     call      [printf ]                ;导入函数调用
     add       esp , 16 * 2 + 4              ;栈平衡
;-------暂停--------------------
     push     szPause
     call      [system ]
     add       esp , 8
     jmp       dword [ esp - 4 ]           ;栈平衡，并结束程序

;===================================================
;MatrixMul(a14,b44,c14) 函数
; 功能:
;      a14*b44=c14 (1*4的向量乘以4*4的矩阵=1*4的向量)
; 参数:
;      a14:1*4的单精度浮点向量指针，输入参数
;      b44:4*4的单精度浮点矩阵指针，输入参数
;      c14:1*4的单精度浮点向量指针，输出参数
; 说明:
;      如果能保存数据对齐，可将movups替换成movaps
;===================================================
MatrixMul :
     ;定义三个标签
     label   .a14 dword at esp + 4      ;arg1
     label   .b44 dword at esp + 8      ;arg2
     label   .c14 dword at esp + 12     ;arg3

     mov      eax ,[.a14 ]
     mov      edx ,[.b44 ]
     mov      ecx ,[.c14 ]
     movups xmm0 ,[ eax ]               ;d c b a

     movaps xmm1 , xmm0                ;d c b a
     movaps xmm2 , xmm0                ;d c b a
     movaps xmm3 , xmm0                ;d c b a

     shufps xmm0 , xmm0 , 0              ;a a a a
     shufps xmm1 , xmm1 , 01010101b      ;b b b b
     shufps xmm2 , xmm2 , 10101010b      ;c c c c
     shufps xmm3 , xmm3 , 11111111b      ;d d d d

     movups xmm4 ,[ edx + 16 * 0 ]
     movups xmm5 ,[ edx + 16 * 1 ]
     movups xmm6 ,[ edx + 16 * 2 ]
   &