| 1 | --- lagged_fibonacci.hpp.orig 2008-11-27 02:10:55.000000000 +0800
|
|---|
| 2 | +++ lagged_fibonacci.hpp 2009-02-01 23:41:14.000000000 +0800
|
|---|
| 3 | @@ -317,7 +317,7 @@
|
|---|
| 4 | unsigned long mask = ~((~0u) << (w%32)); // now lowest w bits set
|
|---|
| 5 | RealType two32 = pow(RealType(2), 32);
|
|---|
| 6 | unsigned int j;
|
|---|
| 7 | - for(j = 0; j < long_lag && first != last; ++j, ++first) {
|
|---|
| 8 | + for(j = 0; j < long_lag && first != last; ++j) {
|
|---|
| 9 | x[j] = RealType(0);
|
|---|
| 10 | for(int k = 0; k < w/32 && first != last; ++k, ++first)
|
|---|
| 11 | x[j] += *first / pow(two32,k+1);
|
|---|