1
0
Fork 0
eschac/src/magics.rs
2025-09-03 20:57:14 +02:00

326 lines
10 KiB
Rust

use crate::bitboard::*;
use crate::board::*;
use crate::rays::Rays;
const BISHOP_SHR: u8 = 55;
const ROOK_SHR: u8 = 52;
pub(crate) struct Magics {
bishop: BySquare<Magic>,
rook: BySquare<Magic>,
table: Box<[Bitboard]>,
}
#[derive(Clone, Copy)]
struct Magic {
premask: Bitboard,
factor: u64,
offset: isize,
}
impl Magics {
pub(crate) fn compute(rays: &Rays) -> Self {
let mut data = Vec::new();
let mut aux =
|shr,
factors: fn(Square) -> u64,
make_table: fn(&Rays, Square) -> (Bitboard, Vec<(Bitboard, Bitboard)>)| {
BySquare::new(|square| {
let (premask, table) = make_table(rays, square);
let factor = factors(square);
let offset = fill_table(&mut data, shr, factor, table);
Magic {
premask,
factor,
offset,
}
})
};
let bishop = aux(BISHOP_SHR, bishop_factors, make_bishop_table);
let rook = aux(ROOK_SHR, rook_factors, make_rook_table);
let mut table = Box::new_uninit_slice(data.len());
for (i, entry) in data.into_iter().enumerate() {
table[i].write(entry);
}
Self {
bishop,
rook,
table: unsafe { table.assume_init() },
}
}
#[inline]
pub(crate) fn bishop(&self, square: Square, blockers: Bitboard) -> Bitboard {
unsafe { self.get_unchecked(BISHOP_SHR, *self.bishop.get(square), blockers) }
}
#[inline]
pub(crate) fn rook(&self, square: Square, blockers: Bitboard) -> Bitboard {
unsafe { self.get_unchecked(ROOK_SHR, *self.rook.get(square), blockers) }
}
#[inline]
unsafe fn get_unchecked(&self, shr: u8, magic: Magic, blockers: Bitboard) -> Bitboard {
let Magic {
premask,
factor,
offset,
} = magic;
*self.table.get_unchecked(
((hash(shr, factor, blockers | premask) as isize).unchecked_add(offset)) as usize,
)
}
}
fn fill_table(
data: &mut Vec<Bitboard>,
shr: u8,
factor: u64,
table: Vec<(Bitboard, Bitboard)>,
) -> isize {
let offset = data.len() as isize
- table
.iter()
.map(|(x, _)| hash(shr, factor, *x) as isize)
.min()
.unwrap();
for (x, y) in &table {
let i = (hash(shr, factor, *x) as isize + offset) as usize;
while data.len() <= i {
data.push(Bitboard::new());
}
if data[i] != Bitboard::new() && data[i] != *y {
panic!();
}
data[i] = *y;
}
offset
}
fn make_bishop_table(rays: &Rays, square: Square) -> (Bitboard, Vec<(Bitboard, Bitboard)>) {
let mut premask = Bitboard::new();
for direction in [
Direction::NorthWest,
Direction::SouthWest,
Direction::SouthEast,
Direction::NorthEast,
] {
premask |= rays.ray(square, direction);
}
premask &= !Rank::First.bitboard();
premask &= !Rank::Eighth.bitboard();
premask &= !File::A.bitboard();
premask &= !File::H.bitboard();
let mut table = make_table(premask, |blockers| {
let mut res = Bitboard::new();
for direction in [
Direction::NorthWest,
Direction::SouthWest,
Direction::SouthEast,
Direction::NorthEast,
] {
res |= rays.blocked(square, direction, blockers);
}
res
});
premask = !premask;
for (x, _) in &mut table {
*x |= premask;
}
(premask, table)
}
fn make_rook_table(rays: &Rays, square: Square) -> (Bitboard, Vec<(Bitboard, Bitboard)>) {
let mut premask = Bitboard::new();
premask |= rays.ray(square, Direction::North) & !Rank::Eighth.bitboard();
premask |= rays.ray(square, Direction::West) & !File::A.bitboard();
premask |= rays.ray(square, Direction::South) & !Rank::First.bitboard();
premask |= rays.ray(square, Direction::East) & !File::H.bitboard();
let mut table = make_table(premask, |blockers| {
let mut res = Bitboard::new();
for direction in [
Direction::North,
Direction::West,
Direction::South,
Direction::East,
] {
res |= rays.blocked(square, direction, blockers);
}
res
});
premask = !premask;
for (x, _) in &mut table {
*x |= premask;
}
(premask, table)
}
fn make_table<T, F>(premask: Bitboard, f: F) -> Vec<(Bitboard, T)>
where
F: Fn(Bitboard) -> T,
{
let mut res = Vec::new();
let mut subset: u64 = 0;
loop {
subset = subset.wrapping_sub(premask.0) & premask.0;
let x = Bitboard(subset);
let y = f(x);
res.push((x, y));
if subset == 0 {
break;
}
}
res
}
fn hash(shr: u8, factor: u64, x: Bitboard) -> usize {
(x.0.wrapping_mul(factor) >> shr) as usize
}
fn bishop_factors(square: Square) -> u64 {
match square {
Square::A1 => 0x0000404040404040,
Square::B1 => 0x0040C100081000E8,
Square::C1 => 0x0000401020200000,
Square::D1 => 0x0040802004000000,
Square::E1 => 0x10403C0180000000,
Square::F1 => 0x0040210100800000,
Square::G1 => 0x0068104002008000,
Square::H1 => 0x0048082080040080,
Square::A2 => 0x0000004040404040,
Square::B2 => 0x0000002020202020,
Square::C2 => 0x00040080184001E4,
Square::D2 => 0x0040008020040000,
Square::E2 => 0x1040003C01800000,
Square::F2 => 0x0078002001008000,
Square::G2 => 0x0068001040020080,
Square::H2 => 0x0068000820010040,
Square::A3 => 0x0000400080808080,
Square::B3 => 0x0000200040404040,
Square::C3 => 0x0000400080808080,
Square::D3 => 0x0000200200801000,
Square::E3 => 0x0060200100080000,
Square::F3 => 0x0000100021C60021,
Square::G3 => 0x0000040010410040,
Square::H3 => 0x0000020008208020,
Square::A4 => 0x0000804000810100,
Square::B4 => 0x0000402000408080,
Square::C4 => 0x0000040800802080,
Square::D4 => 0x000020100C010020,
Square::E4 => 0x0000840000802000,
Square::F4 => 0x0001801800240010,
Square::G4 => 0x0000080800104100,
Square::H4 => 0x0000040400082080,
Square::A5 => 0x0000010278010040,
Square::B5 => 0x000000813C004040,
Square::C5 => 0x000001027A010040,
Square::D5 => 0x0000018180280200,
Square::E5 => 0x0000204018003080,
Square::F5 => 0x0000202040008040,
Square::G5 => 0x0000101010002080,
Square::H5 => 0x0000080808001040,
Square::A6 => 0x0000004100F90080,
Square::B6 => 0x0000002080BC0040,
Square::C6 => 0x0000004103440080,
Square::D6 => 0x0000000080FD0080,
Square::E6 => 0x0000020040100100,
Square::F6 => 0x0000404040400080,
Square::G6 => 0x000000206027D010,
Square::H6 => 0x00000008400DE806,
Square::A7 => 0x0000002101007200,
Square::B7 => 0x0000001041003900,
Square::C7 => 0x080000000F8080A0,
Square::D7 => 0x0000000008003FC0,
Square::E7 => 0x0000000100202000,
Square::F7 => 0x0000004040802000,
Square::G7 => 0x00000060401043D0,
Square::H7 => 0x00000020200413F0,
Square::A8 => 0x1400000F00410088,
Square::B8 => 0x0000000010410039,
Square::C8 => 0x000080000800807E,
Square::D8 => 0x000C69003008003F,
Square::E8 => 0x0000000001002020,
Square::F8 => 0x0000000040408020,
Square::G8 => 0x001980004010801F,
Square::H8 => 0x0000404040404040,
}
}
fn rook_factors(square: Square) -> u64 {
match square {
Square::A1 => 0x002000A28110000C,
Square::B1 => 0x0018000C01060001,
Square::C1 => 0x0040080010004004,
Square::D1 => 0x0028004084200028,
Square::E1 => 0x0030018000900300,
Square::F1 => 0x0020008020010202,
Square::G1 => 0x001800410080001F,
Square::H1 => 0x0068006801040004,
Square::A2 => 0x000028010114000A,
Square::B2 => 0x00000C0083000600,
Square::C2 => 0x0000080401020008,
Square::D2 => 0x0000200200040020,
Square::E2 => 0x0000200100020020,
Square::F2 => 0x00001800C0006018,
Square::G2 => 0x0000180070400018,
Square::H2 => 0x0000180030640018,
Square::A3 => 0x00300018010C0004,
Square::B3 => 0x0004001000080010,
Square::C3 => 0x0001000804020008,
Square::D3 => 0x0002002004002002,
Square::E3 => 0x0001002002002001,
Square::F3 => 0x0001001000801040,
Square::G3 => 0x0000004040008001,
Square::H3 => 0x0000802000200040,
Square::A4 => 0x0040200010080008,
Square::B4 => 0x0000080010040010,
Square::C4 => 0x0001020008040008,
Square::D4 => 0x0000020020040020,
Square::E4 => 0x0000010020020020,
Square::F4 => 0x0000008020010020,
Square::G4 => 0x0000008020200040,
Square::H4 => 0x0000200020004081,
Square::A5 => 0x0000081C00380020,
Square::B5 => 0x0000080400100010,
Square::C5 => 0x0000400880410010,
Square::D5 => 0x0000200200200400,
Square::E5 => 0x0000200100200200,
Square::F5 => 0x0000200080200100,
Square::G5 => 0x0000008000404001,
Square::H5 => 0x0000802000200040,
Square::A6 => 0x0000010B14002800,
Square::B6 => 0x0000030086000C00,
Square::C6 => 0x0000084040804200,
Square::D6 => 0x0000020004002020,
Square::E6 => 0x0000009001803003,
Square::F6 => 0x0000004001004002,
Square::G6 => 0x000000100800A804,
Square::H6 => 0x000000082800D002,
Square::A7 => 0x00000109040200A8,
Square::B7 => 0x000000808A050014,
Square::C7 => 0x0000004048038018,
Square::D7 => 0x0000020020040020,
Square::E7 => 0x0000002030018030,
Square::F7 => 0x0000001800E08018,
Square::G7 => 0x0000000810580050,
Square::H7 => 0x0000000C04600050,
Square::A8 => 0x0000001020891046,
Square::B8 => 0x00000080090015C1,
Square::C8 => 0x0000004005489101,
Square::D8 => 0x0000040810204002,
Square::E8 => 0x000C040810002022,
Square::F8 => 0x0008000404883002,
Square::G8 => 0x0008000400548802,
Square::H8 => 0x0000000224104486,
}
}