Subject: [ublas] Function 'equals' in matrix_assign.hpp
From: Shangtong Zhang (zhangshangtong.cpp_at_[hidden])
Date: 2016-03-04 00:46:00
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 <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?
School of Computer Science,
Fudan University, PRC.