diff --git a/MMVII/src/Matrix/cL1BarrodaleSover.cpp b/MMVII/src/Matrix/cL1BarrodaleSover.cpp index 9c87fcd269..b73bb97b65 100755 --- a/MMVII/src/Matrix/cL1BarrodaleSover.cpp +++ b/MMVII/src/Matrix/cL1BarrodaleSover.cpp @@ -225,12 +225,15 @@ template void cCraig_Barrodale_Roberts_l1::Bench() { int aNbEq = 3 + aDim * 3; cCraig_Barrodale_Roberts_l1 aSys(aDim); - for (int aKEq = 0 ; aKEq::RanGenerate(aDim),RandInInterval(-10.0,10.0)); - } - cDenseVect aVec = aSys.Solve(); - aSys.Bench1Sol(aVec); + for (int aKEq = 0 ; aKEq::RanGenerate(aDim); + auto v3 = RandInInterval(-10.0,10.0); + aSys.PublicAddObservation(v1,v2,v3); // use variable to force evaluation order + } + cDenseVect aVec = aSys.Solve(); + aSys.Bench1Sol(aVec); } } }