咨询自编VMIPS代码与指定C代码在双精度向量场景下的等价性

整体逻辑匹配性 #

从功能意图来看，这段VMIPS代码的核心目标和C代码一致：处理35个双精度元素，计算X、Y、Z三个向量的对应值，但存在两处关键错误，会导致实际运行结果和C代码不一致：

错误1：Y向量的计算逻辑错误 #

C代码中Y的计算式为：

Y[i] = A[i]*D[i] + B[i]*C[i];

但你的VMIPS代码里对运算数的引用出现了混淆：

MULVV.D   V7, V1, V3     #A*D  <-- 错误：V3是C向量，实际计算的是A*C
MULVV.D   V8, V2, V4     #B*C  <-- 错误：V4是D向量，实际计算的是B*D
ADDVV.D   V10, V7, V8    #Y = A*D + B*C

正确的写法应修正为：

MULVV.D   V7, V1, V4     #A*D（V4对应D向量）
MULVV.D   V8, V2, V3     #B*C（V3对应C向量）
ADDVV.D   V10, V7, V8    #Y = A*D + B*C

错误2：SV存储指令的操作数顺序错误 #

VMIPS中SV（Store Vector）指令的格式为**SV 向量寄存器, 内存地址**，作用是将向量寄存器内容存入指定内存地址。你的代码里操作数顺序颠倒了：

SV        Rx, V9         #Store X
SV        Ry, V10        #Store Y
SV        Rz, V13        #Store Z

正确的写法应为：

SV        V9, Rx         #将X向量存入Rx指向的内存
SV        V10, Ry        #将Y向量存入Ry指向的内存
SV        V13, Rz        #将Z向量存入Rz指向的内存

其他正确的实现细节 #

• 正确设置向量长度VLR为35，匹配C代码中循环35次的逻辑（目标向量机长度为64，设置VLR=35会仅处理前35个元素）
• X和Z向量的计算逻辑完全符合C代码的运算步骤

修正上述两处错误后，这段VMIPS代码就能和C代码完全等价。

内容的提问来源于stack exchange，提问作者maffffff