MSWS seems incorrect

[TheIronBorn]

This lib outputs the full 64-bits while the reference implementation only outputs 32-bits. (It also adds the Weyl sequence after squaring instead of before squaring)

inline static uint32_t msws() {

   x *= x; x += (w += s); return x = (x>>32) | (x<<32);

}

[TheIronBorn]

It also sets x to the rotated x rather than merely outputting it