Fixes the Faure sequence when PrimeSieve is used

Closes openturns#2653
mbaudin47 · May 14, 2024 · fae5686 · fae5686
1 parent 501d3dc
commit fae5686
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/ChangeLog b/ChangeLog
@@ -27,6 +27,8 @@
 
 === Miscellaneous ===
 
+=== Bug fixes (total 45) ===
+ * #2653 (The Faure sequence is wrong)
 
 
 == 1.23rc1 release (2024-05-03) == #release-1.23rc1

diff --git a/lib/src/Base/Stat/LowDiscrepancySequenceImplementation.cxx b/lib/src/Base/Stat/LowDiscrepancySequenceImplementation.cxx
@@ -381,9 +381,9 @@ Unsigned64BitsInteger LowDiscrepancySequenceImplementation::GetNextPrimeNumber(c
 #ifdef OPENTURNS_HAVE_PRIMESIEVE
   primesieve::iterator it;
 #if PRIMESIEVE_VERSION_MAJOR >= 11
-  it.jump_to(n - 1);
+  it.jump_to(n);
 #else
-  it.skipto(n - 1);
+  it.skipto(n);
 #endif
   return it.next_prime();
 #else