void mult(const DMatZZpFFPACK& A,
const DMatZZpFFPACK& B,
DMatZZpFFPACK& C)
{
// We assume that C is set to the correct size, and is the zero matrix here.
addMultipleTo(C,A,B);
}
void subtractMultipleTo(DMatZZpFFPACK& C,
const DMatZZpFFPACK& A,
const DMatZZpFFPACK& B)
{
DMatZZpFFPACK::ElementType minus_one;
A.ring().set_from_long(minus_one, -1);
addMultipleTo(C, minus_one, A, B);
}
void addMultipleTo(DMatZZpFFPACK& C,
const DMatZZpFFPACK& A,
const DMatZZpFFPACK& B)
{
DMatZZpFFPACK::ElementType one;
A.ring().set_from_long(one, 1);
addMultipleTo(C, one, A, B);
}
