I’m trying to contribute to toeplitz matrix, here is my prototype

However I have a problem.

When I execute code 

where m is a toeplitz matrix, it will raise exception.

When matrix_assign is executed, it calls ‘equals' in matrix_assign.hpp

https://github.com/ShangtongZhang/ublas/blob/toeplitz/include/boost/numeric/ublas/detail/matrix_assign.hpp line 32

However the parameter epsilon is 0, norm_inf (e1 - e2 ) < 0 will never succeed.

So it will return false and raise exception.

Same things happen even m is a banded matrix.

Why is it ‘<' rather than '<= '? Is it by design?

Thanks,