Function AtomicMultiplyAdd

src/cryptlib/integer.cpp:275–320

}

static word AtomicMultiplyAdd(word *C, word A0, word A1, word B0, word B1)
{
    word s;
    dword_union d;

    if (A1 >= A0)
        if (B0 >= B1)
        {
            s = 0;
            d.dw = (dword)(A1-A0)*(B0-B1);
        }
        else
        {
            s = (A1-A0);
            d.dw = (dword)s*(word)(B0-B1);
        }
    else
        if (B0 > B1)
        {
            s = (B0-B1);
            d.dw = (word)(A1-A0)*(dword)s;
        }
        else
        {
            s = 0;
            d.dw = (dword)(A0-A1)*(B1-B0);
        }

    dword_union A0B0;
    A0B0.dw = (dword)A0*B0;
    dword_union t;
    t.dw = A0B0.dw + C[0];
    C[0] = t.low;

    dword_union A1B1;
        A1B1.dw = (dword)A1*B1;
    t.dw = (dword) t.high + A0B0.low + d.low + A1B1.low + C[1];
    C[1] = t.low;

    t.dw = (dword) t.high + A1B1.low + A0B0.high + d.high + A1B1.high - s + C[2];
    C[2] = t.low;

    t.dw = (dword) t.high + A1B1.high + C[3];
    C[3] = t.low;
    return t.high;
}

static inline void AtomicSquare(word *C, word A, word B)
{