#48 sparse_multipolynom::normalize doesn't work

[Sergey S. Lyalin]

Where: trunk #417, Version 2.2.9.412.

When sparse_multipolynom::normalize is called for 2*x^(0, 1, 1, 0, 0)+x^(0, 0, 0, 1, 2)+x^(1, 0, 0, 0, 0), it doesn't sort anything inside the polynom, but it should give x^(1, 0, 0, 0, 0)+2*x^(0, 1, 1, 0, 0)+x^(0, 0, 0, 1, 2). Lexicographical monomial order is used. I think the problem is in function sparse_multipolynom_config_default::normalize at src/arageli/sparse_multipolynom.hpp : 1145.

Current workaround: do not call sparse_multipolynom::normalize. Instead do the following: form a new sparse_multipolynom value v_new from v_old by v_new += m, for all m in a set of v_old monomials, then use v_new as normalized v_old.

[Sergey S. Lyalin]

Applied workaround at #452. Not tested much, so still open.