Handle += and -= updates in sparse-sparse products

Submitted by Gael Guennebaud

Assigned to Nobody

As noticed there:

it should be easy to extend the current sparse-sparse products implementation to natively support operations like:

C += alpha * A * B;

C -= alpha * A * B;

without extra temporaries.