anchor.diff:  -spec exs1024_calc(uint64(), uint64()) -> {uint64(), uint64()}.
   
   exs1024_calc(S0, S1) ->
  -    S11 = S1 bxor ((S1 bsl 31) band ?UINT64MASK),
  +    S11 = S1 bxor ((S1 band ?UINT33MASK) bsl 31),
       S12 = S11 bxor (S11 bsr 11),
       S01 = S0 bxor (S0 bsr 30),
       NS1 = S01 bxor S12,
  
positive.code: 
  exs64_next(R) ->
      R1 = R bxor (R bsr 12),
      R2 = R1 bxor ((R1 bsl 25) band ?UINT64MASK),
      R3 = R2 bxor (R2 bsr 27),
      {(R3 * 2685821657736338717) band ?UINT64MASK, R3}.
  
  
negative.code: 
  -type exs1024_state() :: {list(uint64()), list(uint64())}.
  
  -define(UINT33MASK, 16#1ffffffff).
  
  %% Calculation of xorshift1024star.
  %% exs1024_calc(S0, S1) -> {X, NS1}.
  %% X: random number output
  