Consider a semiring concept

[BenBrock]

This could simplify the definition of multiply, ewise_add, and ewise_multiply. Users could also define semirings in-place using a monoid and a binary operator.

multiply(a, b, c, {grb::plus(), grb::times()});