--- lagged_fibonacci.hpp.orig 2008-11-27 02:10:55.000000000 +0800 +++ lagged_fibonacci.hpp 2009-01-31 22:05:44.000000000 +0800 @@ -319,7 +319,7 @@ unsigned int j; for(j = 0; j < long_lag && first != last; ++j, ++first) { x[j] = RealType(0); - for(int k = 0; k < w/32 && first != last; ++k, ++first) + for(int k = 0; k < w/32 && first != last; ++k) x[j] += *first / pow(two32,k+1); if(first != last && mask != 0) x[j] += fmod((*first & mask) / _modulus, RealType(1));