scalar_mult_add
scalar_mult_add(out, ++lhs, ++rhs);
scalar_mult_add(out, m++, ++ar);
scalar_mult_add(out, mr = mc++, ar = a);
scalar_mult_add(out, (mr += rank), ++ar);