CFD/eigen
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fix noise in sparse_basic_3 (numerical cancellation)

Browse Source 
(grafted from cd25b538ab
)
This commit is contained in:
Gael Guennebaud 2018-12-08 00:13:37 +01:00
parent c285ed1033
commit 204d1f1456
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

@ -630,7 +630,8 @@ void big_sparse_triplet(Index rows, Index cols, double density) {
{
Index r = internal::random<Index>(0,rows-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));
sum += v;
}