revm_precompile/bls12_381/
g1_msm.rs

1//! BLS12-381 G1 msm precompile. More details in [`g1_msm`]
2use crate::bls12_381::utils::{pad_g1_point, remove_g1_padding};
3use crate::bls12_381::G1Point;
4use crate::bls12_381_const::{
5    DISCOUNT_TABLE_G1_MSM, G1_MSM_ADDRESS, G1_MSM_BASE_GAS_FEE, G1_MSM_INPUT_LENGTH,
6    PADDED_G1_LENGTH, SCALAR_LENGTH,
7};
8use crate::bls12_381_utils::msm_required_gas;
9use crate::{
10    crypto, Precompile, PrecompileError, PrecompileId, PrecompileOutput, PrecompileResult,
11};
12
13/// [EIP-2537](https://eips.ethereum.org/EIPS/eip-2537#specification) BLS12_G1MSM precompile.
14pub const PRECOMPILE: Precompile =
15    Precompile::new(PrecompileId::Bls12G1Msm, G1_MSM_ADDRESS, g1_msm);
16
17/// Implements EIP-2537 G1MSM precompile.
18/// G1 multi-scalar-multiplication call expects `160*k` bytes as an input that is interpreted
19/// as byte concatenation of `k` slices each of them being a byte concatenation
20/// of encoding of G1 point (`128` bytes) and encoding of a scalar value (`32`
21/// bytes).
22/// Output is an encoding of multi-scalar-multiplication operation result - single G1
23/// point (`128` bytes).
24/// See also: <https://eips.ethereum.org/EIPS/eip-2537#abi-for-g1-multiexponentiation>
25pub fn g1_msm(input: &[u8], gas_limit: u64) -> PrecompileResult {
26    let input_len = input.len();
27    if input_len == 0 || !input_len.is_multiple_of(G1_MSM_INPUT_LENGTH) {
28        return Err(PrecompileError::Bls12381G1MsmInputLength);
29    }
30
31    let k = input_len / G1_MSM_INPUT_LENGTH;
32    let required_gas = msm_required_gas(k, &DISCOUNT_TABLE_G1_MSM, G1_MSM_BASE_GAS_FEE);
33    if required_gas > gas_limit {
34        return Err(PrecompileError::OutOfGas);
35    }
36
37    let mut valid_pairs_iter = (0..k).map(|i| {
38        let start = i * G1_MSM_INPUT_LENGTH;
39        let padded_g1 = &input[start..start + PADDED_G1_LENGTH];
40        let scalar_bytes = &input[start + PADDED_G1_LENGTH..start + G1_MSM_INPUT_LENGTH];
41
42        // Remove padding from G1 point - this validates padding format
43        let [x, y] = remove_g1_padding(padded_g1)?;
44        let scalar_array: [u8; SCALAR_LENGTH] = scalar_bytes.try_into().unwrap();
45
46        let point: G1Point = (*x, *y);
47        Ok((point, scalar_array))
48    });
49
50    let unpadded_result = crypto().bls12_381_g1_msm(&mut valid_pairs_iter)?;
51
52    // Pad the result for EVM compatibility
53    let padded_result = pad_g1_point(&unpadded_result);
54
55    Ok(PrecompileOutput::new(required_gas, padded_result.into()))
56}
57
58#[cfg(test)]
59mod test {
60    use super::*;
61    use primitives::{hex, Bytes};
62
63    #[test]
64    fn bls_g1multiexp_g1_not_on_curve_but_in_subgroup() {
65        let input = Bytes::from(hex!("000000000000000000000000000000000a2833e497b38ee3ca5c62828bf4887a9f940c9e426c7890a759c20f248c23a7210d2432f4c98a514e524b5184a0ddac00000000000000000000000000000000150772d56bf9509469f9ebcd6e47570429fd31b0e262b66d512e245c38ec37255529f2271fd70066473e393a8bead0c30000000000000000000000000000000000000000000000000000000000000000"));
66        let fail = g1_msm(&input, G1_MSM_BASE_GAS_FEE);
67        assert_eq!(fail, Err(PrecompileError::Bls12381G1NotOnCurve));
68    }
69}