GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-04-20 16:49:38 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix noise in sparse_basic_3 (numerical cancellation)

Browse Source
This commit is contained in:
Gael Guennebaud 2018-12-08 00:13:37 +01:00
parent efaf03bf96
commit cd25b538ab
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
test/sparse_basic.cpp
View File

@ -658,7 +658,8 @@ void big_sparse_triplet(Index rows, Index cols, double density) {
{ {
Index r = internal::random<Index>(0,rows-1); Index r = internal::random<Index>(0,rows-1);
Index c = internal::random<Index>(0,cols-1); Index c = internal::random<Index>(0,cols-1);
Scalar v = internal::random<Scalar>(); // use positive values to prevent numerical cancellation errors in sum
Scalar v = numext::abs(internal::random<Scalar>());
triplets.push_back(TripletType(r,c,v)); triplets.push_back(TripletType(r,c,v));
sum += v; sum += v;
} }