portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

RNG fix

Browse Source
This commit is contained in:
Peter Boyle 2022-06-27 12:09:52 -04:00
parent 8208a6214f
commit e1e7b1e224
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Grid/lattice/Lattice_rng.h
View File

@ -428,18 +428,20 @@ public:
thread_for( lidx, _grid->lSites(), { thread_for( lidx, _grid->lSites(), {
int gidx; int gidx;
int o_idx;
int i_idx;
int rank;
Coordinate pcoor; Coordinate pcoor;
Coordinate lcoor; Coordinate lcoor;
Coordinate gcoor; Coordinate gcoor;
_grid->LocalIndexToLocalCoor(lidx,lcoor); _grid->LocalIndexToLocalCoor(lidx,lcoor);
pcoor=_grid->ThisProcessorCoor(); pcoor=_grid->ThisProcessorCoor();
_grid->ProcessorCoorLocalCoorToGlobalCoor(pcoor,lcoor,gcoor); _grid->ProcessorCoorLocalCoorToGlobalCoor(pcoor,lcoor,gcoor);
_grid->GlobalCoorToGlobalIndex(gcoor,gidx);
int o_idx;
int i_idx;
int rank;
_grid->GlobalCoorToRankIndex(rank,o_idx,i_idx,gcoor); _grid->GlobalCoorToRankIndex(rank,o_idx,i_idx,gcoor);
assert(rank == _grid->ThisRank() ); assert(rank == _grid->ThisRank() );
int l_idx=generator_idx(o_idx,i_idx); int l_idx=generator_idx(o_idx,i_idx);
_generators[l_idx] = master_engine; _generators[l_idx] = master_engine;
Skip(_generators[l_idx],gidx); // Skip to next RNG sequence Skip(_generators[l_idx],gidx); // Skip to next RNG sequence