[uBlas] compressed_matrix::erase_element() segfaults

[Niko Vuokko <niko.vuokko@…>]

The following example segfaults, system is GCC 3.4.4 on Cygwin, Boost 1.35.0.

compressed_matrix<int> m(2,2);
m(0,0) = 1;
m.erase_element(1,0);

The problem is that inside erase_element std::lower_bound oddly returns it_begin, even though simple

while(*itt < elem2)
   ++itt;

returns it_begin+2.

[(none)]

Milestone Boost 1.35.1 deleted

[Gunter]

bug still exists with gcc version 4.3.2 (Debian 4.3.2-1.1)

[Gunter]

example to provoke segfault

[Gunter]

(In [54230]) comp_mat_erase.cpp - unit test to verify compressed matrix layout, insertion and erasure of elements, see #2195

CMakeLists.txt - added new test: comp_mat_erase

Jamfile.v2 - added new test: comp_mat_erase

[Gunter]

(In [54231]) matrix_sparse.hpp - see #2195, fixed compressed_matrix::erase_element(); must use >= instead of > because the last element of index1_data is the "end of list" marker. Segfault was caused by *(itv + 1) which pointed to an element _after_ index1_data.end()

[Gunter]

(In [56605]) fix #2195:

• merged changes into release branch

[Gunter]

merged into 1.41.0