bsnes-libretro/nall/galois-field.hpp
Tim Allen d87a0f633d Update to bsnes v107r4 beta release.
byuu says:

  - bsnes: added video filters from bsnes v082
  - bsnes: added ZSNES snow effect option when games paused or unloaded
    (no, I'm not joking)
  - bsnes: added 7-zip support (LZMA 19.00 SDK)

[Recent higan WIPs have also mentioned bsnes changes, although the higan code
no longer includes the bsnes code. These changes include:

  - higan, bsnes: added EXLOROM, EXLOROM-RAM, EXHIROM mappings
  - higan, bsnes: focus the viewport after leaving fullscreen exclusive
    mode
  - bsnes: re-added mightymo's cheat code database
  - bsnes: improved make install rules for the game and cheat code
    databases
  - bsnes: delayed construction of hiro::Window objects to properly show
    bsnes window icons

- Ed.]
2019-07-07 19:44:09 +10:00

71 lines
2.0 KiB
C++

#pragma once
//table-driven galois field modulo 2
//do not use with GF(2^17) or larger
namespace nall {
template<typename field, uint Elements, uint Polynomial>
struct GaloisField {
using type = GaloisField;
GaloisField(uint x = 0) : x(x) {}
operator field() const { return x; }
auto operator^(field y) const -> type { return x ^ y; }
auto operator+(field y) const -> type { return x ^ y; }
auto operator-(field y) const -> type { return x ^ y; }
auto operator*(field y) const -> type { return x && y ? exp(log(x) + log(y)) : 0; }
auto operator/(field y) const -> type { return x && y ? exp(log(x) + Elements - log(y)) : 0; }
auto& operator =(field y) { return x = y, *this; }
auto& operator^=(field y) { return x = operator^(y), *this; }
auto& operator+=(field y) { return x = operator^(y), *this; }
auto& operator-=(field y) { return x = operator^(y), *this; }
auto& operator*=(field y) { return x = operator*(y), *this; }
auto& operator/=(field y) { return x = operator/(y), *this; }
auto pow(field y) const -> type { return exp(log(x) * y); }
auto inv() const -> type { return exp(Elements - log(x)); } // 1/x
static auto log(uint x) -> uint {
enum : uint { Size = bit::round(Elements), Mask = Size - 1 };
static array<field[Size]> log = [] {
uint shift = 0, polynomial = Polynomial;
while(polynomial >>= 1) shift++;
shift--;
array<field[Size]> log;
field x = 1;
for(uint n : range(Elements)) {
log[x] = n;
x = x << 1 ^ (x >> shift ? Polynomial : 0);
}
log[0] = 0; //-inf (undefined)
return log;
}();
return log[x & Mask];
}
static auto exp(uint x) -> uint {
static array<field[Elements]> exp = [] {
uint shift = 0, polynomial = Polynomial;
while(polynomial >>= 1) shift++;
shift--;
array<field[Elements]> exp;
field x = 1;
for(uint n : range(Elements)) {
exp[n] = x;
x = x << 1 ^ (x >> shift ? Polynomial : 0);
}
return exp;
}();
return exp[x % Elements];
}
field x;
};
}