binius_field/arch/portable/
packed_128.rs

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
// Copyright 2024-2025 Irreducible Inc.

use super::{
	packed::{impl_broadcast, impl_ops_for_zero_height, PackedPrimitiveType},
	packed_arithmetic::{alphas, impl_tower_constants},
};
use crate::{
	arch::{PackedStrategy, PairwiseRecursiveStrategy, PairwiseStrategy, PairwiseTableStrategy},
	arithmetic_traits::{
		impl_invert_with, impl_mul_alpha_with, impl_mul_with, impl_square_with,
		impl_transformation_with_strategy,
	},
	BinaryField128b, BinaryField16b, BinaryField1b, BinaryField2b, BinaryField32b, BinaryField4b,
	BinaryField64b, BinaryField8b,
};

// Define 128 bit packed field types
pub type PackedBinaryField128x1b = PackedPrimitiveType<u128, BinaryField1b>;
pub type PackedBinaryField64x2b = PackedPrimitiveType<u128, BinaryField2b>;
pub type PackedBinaryField32x4b = PackedPrimitiveType<u128, BinaryField4b>;
pub type PackedBinaryField16x8b = PackedPrimitiveType<u128, BinaryField8b>;
pub type PackedBinaryField8x16b = PackedPrimitiveType<u128, BinaryField16b>;
pub type PackedBinaryField4x32b = PackedPrimitiveType<u128, BinaryField32b>;
pub type PackedBinaryField2x64b = PackedPrimitiveType<u128, BinaryField64b>;
pub type PackedBinaryField1x128b = PackedPrimitiveType<u128, BinaryField128b>;

// Define broadcast
impl_broadcast!(u128, BinaryField1b);
impl_broadcast!(u128, BinaryField2b);
impl_broadcast!(u128, BinaryField4b);
impl_broadcast!(u128, BinaryField8b);
impl_broadcast!(u128, BinaryField16b);
impl_broadcast!(u128, BinaryField32b);
impl_broadcast!(u128, BinaryField64b);
impl_broadcast!(u128, BinaryField128b);

// Define operations for height 0
impl_ops_for_zero_height!(PackedBinaryField128x1b);

// Define constants
impl_tower_constants!(BinaryField1b, u128, { alphas!(u128, 0) });
impl_tower_constants!(BinaryField2b, u128, { alphas!(u128, 1) });
impl_tower_constants!(BinaryField4b, u128, { alphas!(u128, 2) });
impl_tower_constants!(BinaryField8b, u128, { alphas!(u128, 3) });
impl_tower_constants!(BinaryField16b, u128, { alphas!(u128, 4) });
impl_tower_constants!(BinaryField32b, u128, { alphas!(u128, 5) });
impl_tower_constants!(BinaryField64b, u128, { alphas!(u128, 6) });

// Define multiplication
impl_mul_with!(PackedBinaryField64x2b @ PackedStrategy);
impl_mul_with!(PackedBinaryField32x4b @ PackedStrategy);
impl_mul_with!(PackedBinaryField16x8b @ PackedStrategy);
impl_mul_with!(PackedBinaryField8x16b @ PairwiseStrategy);
impl_mul_with!(PackedBinaryField4x32b @ PairwiseStrategy);
impl_mul_with!(PackedBinaryField2x64b @ PairwiseRecursiveStrategy);
impl_mul_with!(PackedBinaryField1x128b @ PairwiseRecursiveStrategy);

// Define square
impl_square_with!(PackedBinaryField64x2b @ PackedStrategy);
impl_square_with!(PackedBinaryField32x4b @ PackedStrategy);
impl_square_with!(PackedBinaryField16x8b @ PackedStrategy);
impl_square_with!(PackedBinaryField8x16b @ PairwiseRecursiveStrategy);
impl_square_with!(PackedBinaryField4x32b @ PairwiseStrategy);
impl_square_with!(PackedBinaryField2x64b @ PairwiseStrategy);
impl_square_with!(PackedBinaryField1x128b @ PairwiseRecursiveStrategy);

// Define invert
impl_invert_with!(PackedBinaryField64x2b @ PackedStrategy);
impl_invert_with!(PackedBinaryField32x4b @ PackedStrategy);
impl_invert_with!(PackedBinaryField16x8b @ PairwiseTableStrategy);
impl_invert_with!(PackedBinaryField8x16b @ PairwiseRecursiveStrategy);
impl_invert_with!(PackedBinaryField4x32b @ PairwiseStrategy);
impl_invert_with!(PackedBinaryField2x64b @ PairwiseRecursiveStrategy);
impl_invert_with!(PackedBinaryField1x128b @ PairwiseRecursiveStrategy);

// Define multiply by alpha
impl_mul_alpha_with!(PackedBinaryField64x2b @ PackedStrategy);
impl_mul_alpha_with!(PackedBinaryField32x4b @ PackedStrategy);
impl_mul_alpha_with!(PackedBinaryField16x8b @ PackedStrategy);
impl_mul_alpha_with!(PackedBinaryField8x16b @ PackedStrategy);
impl_mul_alpha_with!(PackedBinaryField4x32b @ PackedStrategy);
impl_mul_alpha_with!(PackedBinaryField2x64b @ PackedStrategy);
impl_mul_alpha_with!(PackedBinaryField1x128b @ PairwiseRecursiveStrategy);

// Define linear transformations
impl_transformation_with_strategy!(PackedBinaryField128x1b, PackedStrategy);
impl_transformation_with_strategy!(PackedBinaryField64x2b, PackedStrategy);
impl_transformation_with_strategy!(PackedBinaryField32x4b, PackedStrategy);
impl_transformation_with_strategy!(PackedBinaryField16x8b, PackedStrategy);
impl_transformation_with_strategy!(PackedBinaryField8x16b, PackedStrategy);
impl_transformation_with_strategy!(PackedBinaryField4x32b, PackedStrategy);
impl_transformation_with_strategy!(PackedBinaryField2x64b, PairwiseStrategy);
impl_transformation_with_strategy!(PackedBinaryField1x128b, PairwiseStrategy);