使用FMA(融合乘数）指令进行复杂的乘法

我想利用可用的融合乘以添加/减去CPU指令，以协助通过尺寸较大的数组进行复杂的乘法。本质上，基本数学看起来像：

void ComplexMultiplyAddToArray(float* pDstR, float* pDstI, const float* pSrc1R, const float* pSrc1I, const float* pSrc2R, const float* pSrc2I, int len)
{
    for (int i = 0; i < len; ++i)
    {
        const float fSrc1R = pSrc1R[i];
        const float fSrc1I = pSrc1I[i];
        const float fSrc2R = pSrc2R[i];
        const float fSrc2I = pSrc2I[i];
        //  Perform complex multiplication on the input and accumulate with the output
        pDstR[i] += fSrc1R*fSrc2R - fSrc1I*fSrc2I;
        pDstI[i] += fSrc1R*fSrc2I + fSrc2R*fSrc1I;
    }
}

您可能会看到，数据是在我们具有独立数量和虚构数字的单独数组的情况下结构的。现在，假设我具有以下功能作为内在的单个指令，分别执行A b c和a b-c：

float fmadd(float a, float b, float c);
float fmsub(float a, float b, float c);

天真地，我可以看到我可以替换2次，一个添加，一个添加和一个fmadd和一个fmsub，例如：

//  Perform complex multiplication on the input and accumulate with the output
pDstR[i] += fmsub(fSrc1R, fSrc2R, fSrc1I*fSrc2I);
pDstI[i] += fmadd(fSrc1R, fSrc2I, fSrc2R*fSrc1I);

我认为这会导致非常适中的性能改进，但我认为我确实缺少可以通过代数修改数学的东西，以便我可以替换几个多/add/add/mult/sub组合。在每一行中，都有一个额外的添加，我觉得我可以将其转换为单个FMA，但是令人沮丧的是，我不知道如何在不更改操作顺序并获得错误的结果的情况下弄清楚如何做。有想法的数学专家吗？

出于问题的目的，目标平台可能并不重要，因为我知道各种平台上都存在这些说明。