exvector tv;
tv.reserve(vars.size() + 1);
for (std::size_t j = 0; j < vars.size(); ++j) {
- if (ec[i].first[j] != 0)
- tv.push_back(power(vars[j], ec[i].first[j]));
+ const exp_vector_t& exp_vector(ec[i].first);
+
+ bug_on(exp_vector.size() != vars.size(),
+ "expected " << vars.size() << " variables, "
+ "expression has " << exp_vector.size() << " instead");
+
+ if (exp_vector[j] != 0)
+ tv.push_back(power(vars[j], exp_vector[j]));
}
tv.push_back(ec[i].second);
ex tmp = (new mul(tv))->setflag(status_flags::dynallocated);
// p_1(x_n) p_2(x_0, \ldots, x_{n-1})
c = ec.rbegin()->second;
ec.rbegin()->second = ex1;
- pp = ex_collect_to_ex(ec, vars).expand().smod(numeric(p));
+ pp = ex_collect_to_ex(ec, rest_vars).expand().smod(numeric(p));
return;
}