Skip to content

nogmp: remove GMP #2

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
xurzie wants to merge 7 commits into
base: master
Choose a base branch
from
Open

nogmp: remove GMP #2

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
f2cf1a7
nogmp: remove GMP
xurzie Jan 12, 2026
d2c2449
nogmp: update fft.cpp, alt_bn128_test.cpp w/ U256 API
xurzie Jan 18, 2026
920fc8b
nogmp: add testAltBn128 to taskfile.js
xurzie Jan 18, 2026
73ce049
nogmp: update
xurzie Jan 27, 2026
9eaffd0
U256 API removed, using Limb API instead
xurzie Feb 2, 2026
68f2a46
nogmp: alt_bn128_test, fft fix up to new version of mp API
xurzie Feb 8, 2026
387a902
nogmp: mp-API fix
xurzie Feb 20, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
165 changes: 66 additions & 99 deletions c/alt_bn128_test.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,16 +1,17 @@
#include <gmp.h>
#include <iostream>
#include <cstdint>
#include <cstring>

#include "gtest/gtest.h"
#include "alt_bn128.hpp"
#include "fft.hpp"
#include "../../../build/mp.hpp"

using namespace AltBn128;

namespace {

TEST(altBn128, f2_simpleMul) {

F2Element e1;
F2.fromString(e1, "(2,2)");

Expand All @@ -23,8 +24,6 @@ TEST(altBn128, f2_simpleMul) {
F2Element e33;
F2.fromString(e33, "(0,12)");

// std::cout << F2.toString(e3) << std::endl;

ASSERT_TRUE(F2.eq(e3, e33));
}

Expand Down Expand Up @@ -94,17 +93,13 @@ TEST(altBn128, f12_inv) {

TEST(altBn128, g1_PlusZero) {
G1Point p1;

G1.add(p1, G1.one(), G1.zero());

ASSERT_TRUE(G1.eq(p1, G1.one()));
}

TEST(altBn128, g1_minus_g1) {
G1Point p1;

G1.sub(p1, G1.one(), G1.one());

ASSERT_TRUE(G1.isZero(p1));
}

Expand All @@ -121,7 +116,6 @@ TEST(altBn128, g1_times_4) {
ASSERT_TRUE(G1.eq(p1,p2));
}


TEST(altBn128, g1_times_3) {
G1Point p1;
G1.add(p1, G1.one(), G1.one());
Expand All @@ -135,21 +129,18 @@ TEST(altBn128, g1_times_3) {
ASSERT_TRUE(G1.eq(p1,p2));
}

TEST(altBn128, g1_times_3_exp) {
TEST(altBn128, g1_times_3_exp) {
G1Point p1;
G1.add(p1, G1.one(), G1.one());
G1.add(p1, p1, G1.one());

mpz_t e;
mpz_init_set_str(e, "3", 10);

uint8_t scalar[32];
for (int i=0;i<32;i++) scalar[i] = 0;
mpz_export((void *)scalar, NULL, -1, 8, -1, 0, e);
mpz_clear(e);
mp_uint_t scalar;
mp_uint_t x;
mp_set(x, 3);
mp_copy(scalar, x);

G1Point p2;
G1.mulByScalar(p2, G1.one(), scalar, 32);
G1.mulByScalar(p2, G1.one(), (uint8_t*)scalar, MP_N);

ASSERT_TRUE(G1.eq(p1,p2));
}
Expand All @@ -174,8 +165,7 @@ TEST(altBn128, g1_times_5) {
ASSERT_TRUE(G1.eq(p1,p6));
}

TEST(altBn128, g1_times_65_exp) {

TEST(altBn128, g1_times_65_exp) {
G1Point p1;
G1.dbl(p1, G1.one());
G1.dbl(p1, p1);
Expand All @@ -185,89 +175,79 @@ TEST(altBn128, g1_times_65_exp) {
G1.dbl(p1, p1);
G1.add(p1, p1, G1.one());

mpz_t e;
mpz_init_set_str(e, "65", 10);

uint8_t scalar[32];
for (int i=0;i<32;i++) scalar[i] = 0;
mpz_export((void *)scalar, NULL, -1, 8, -1, 0, e);
mpz_clear(e);
mp_uint_t scalar;
mp_uint_t x;
mp_set(x, 65);
mp_copy(scalar, x);

G1Point p2;
G1.mulByScalar(p2, G1.one(), scalar, 32);
G1.mulByScalar(p2, G1.one(), (uint8_t*)scalar, MP_N);

ASSERT_TRUE(G1.eq(p1,p2));
}

TEST(altBn128, g1_expToOrder) {
mpz_t e;
mpz_init_set_str(e, "21888242871839275222246405745257275088548364400416034343698204186575808495617", 10);

uint8_t scalar[32];

for (int i=0;i<32;i++) scalar[i] = 0;
mpz_export((void *)scalar, NULL, -1, 8, -1, 0, e);
mpz_clear(e);
TEST(altBn128, g1_expToOrder) {
mp_uint_t scalar;
mp_uint_t x;
ASSERT_EQ(mp_set(x,
"21888242871839275222246405745257275088548364400416034343698204186575808495617",
10
), 0);
mp_copy(scalar, x);

G1Point p1;

G1.mulByScalar(p1, G1.one(), scalar, 32);
G1.mulByScalar(p1, G1.one(), (uint8_t *)scalar, MP_N);

ASSERT_TRUE(G1.isZero(p1));
}

TEST(altBn128, g2_expToOrder) {
mpz_t e;
mpz_init_set_str(e, "21888242871839275222246405745257275088548364400416034343698204186575808495617", 10);

uint8_t scalar[32];

for (int i=0;i<32;i++) scalar[i] = 0;
mpz_export((void *)scalar, NULL, -1, 8, -1, 0, e);
mpz_clear(e);
TEST(altBn128, g2_expToOrder) {
mp_uint_t scalar;
mp_uint_t x;
ASSERT_EQ(mp_set(x,
"21888242871839275222246405745257275088548364400416034343698204186575808495617",
10
), 0);
mp_copy(scalar, x);

Curve<F2Field<RawFq>>::Point p1;

G2.mulByScalar(p1, G2.one(), scalar, 32);
G2.mulByScalar(p1, G2.one(), (uint8_t *)scalar, MP_N);

ASSERT_TRUE(G2.isZero(p1));
}

TEST(altBn128, multiExp) {

int NMExp = 40000;

typedef uint8_t Scalar[32];
typedef mp_uint_t Scalar;

Scalar *scalars = new Scalar[NMExp];
G1PointAffine *bases = new G1PointAffine[NMExp];

uint64_t acc=0;
for (int i=0; i<NMExp; i++) {
if (i==0) {
uint64_t acc = 0;
for (int i = 0; i < NMExp; i++) {
if (i == 0) {
G1.copy(bases[0], G1.one());
} else {
G1.add(bases[i], bases[i-1], G1.one());
}
for (int j=0; j<32; j++) scalars[i][j] = 0;
*(int *)&scalars[i][0] = i+1;
acc += (i+1)*(i+1);

mp_uint_t x;
mp_set(x, (uint64_t)(i + 1));
mp_copy(scalars[i], x);
acc += (uint64_t)(i + 1) * (uint64_t)(i + 1);
}

G1Point p1;
G1.multiMulByScalar(p1, bases, (uint8_t *)scalars, 32, NMExp);

mpz_t e;
mpz_init_set_ui(e, acc);

Scalar sAcc;
G1.multiMulByScalar(p1, bases, (uint8_t *)scalars, MP_N, NMExp);

for (int i=0;i<32;i++) sAcc[i] = 0;
mpz_export((void *)sAcc, NULL, -1, 8, -1, 0, e);
mpz_clear(e);
mp_uint_t sAcc;
mp_uint_t x;
mp_set(x, acc);
mp_copy(sAcc, x);

G1Point p2;
G1.mulByScalar(p2, G1.one(), sAcc, 32);
G1.mulByScalar(p2, G1.one(), (uint8_t *)sAcc, MP_N);

ASSERT_TRUE(G1.eq(p1, p2));

Expand All @@ -276,40 +256,37 @@ TEST(altBn128, multiExp) {
}

TEST(altBn128, multiExpMSM) {

int NMExp = 40000;

typedef uint8_t Scalar[32];
typedef mp_uint_t Scalar;

Scalar *scalars = new Scalar[NMExp];
G1PointAffine *bases = new G1PointAffine[NMExp];

uint64_t acc=0;
for (int i=0; i<NMExp; i++) {
if (i==0) {
uint64_t acc = 0;
for (int i = 0; i < NMExp; i++) {
if (i == 0) {
G1.copy(bases[0], G1.one());
} else {
G1.add(bases[i], bases[i-1], G1.one());
}
for (int j=0; j<32; j++) scalars[i][j] = 0;
*(int *)&scalars[i][0] = i+1;
acc += (i+1)*(i+1);

mp_uint_t x;
mp_set(x, (uint64_t)(i + 1));
mp_copy(scalars[i], x);
acc += (uint64_t)(i + 1) * (uint64_t)(i + 1);
}

G1Point p1;
G1.multiMulByScalarMSM(p1, bases, (uint8_t *)scalars, 32, NMExp);

mpz_t e;
mpz_init_set_ui(e, acc);
G1.multiMulByScalarMSM(p1, bases, (uint8_t *)scalars, MP_N, NMExp);

Scalar sAcc;

for (int i=0;i<32;i++) sAcc[i] = 0;
mpz_export((void *)sAcc, NULL, -1, 8, -1, 0, e);
mpz_clear(e);
mp_uint_t sAcc;
mp_uint_t x;
mp_set(x, acc);
mp_copy(sAcc, x);

G1Point p2;
G1.mulByScalar(p2, G1.one(), sAcc, 32);
G1.mulByScalar(p2, G1.one(), (uint8_t *)sAcc, MP_N);

ASSERT_TRUE(G1.eq(p1, p2));

Expand All @@ -318,7 +295,6 @@ TEST(altBn128, multiExpMSM) {
}

TEST(altBn128, multiExp2) {

int NMExp = 2;

AltBn128::FrElement *scalars = new AltBn128::FrElement[NMExp];
Expand All @@ -341,19 +317,16 @@ TEST(altBn128, multiExp2) {
F1.fromString(ref.x, "9163953212624378696742080269971059027061360176019470242548968584908855004282");
F1.fromString(ref.y, "20922060990592511838374895951081914567856345629513259026540392951012456141360");

G1.multiMulByScalar(r, bases, (uint8_t *)scalars, 32, 2);
G1.multiMulByScalar(r, bases, (uint8_t *)scalars, MP_N, 2);
G1.copy(ra, r);

// std::cout << G1.toString(r, 10);

ASSERT_TRUE(G1.eq(ra, ref));

delete[] bases;
delete[] scalars;
}

TEST(altBn128, multiExp2MSM) {

int NMExp = 2;

AltBn128::FrElement *scalars = new AltBn128::FrElement[NMExp];
Expand All @@ -376,19 +349,16 @@ TEST(altBn128, multiExp2MSM) {
F1.fromString(ref.x, "9163953212624378696742080269971059027061360176019470242548968584908855004282");
F1.fromString(ref.y, "20922060990592511838374895951081914567856345629513259026540392951012456141360");

G1.multiMulByScalarMSM(r, bases, (uint8_t *)scalars, 32, 2);
G1.multiMulByScalarMSM(r, bases, (uint8_t *)scalars, MP_N, 2);
G1.copy(ra, r);

// std::cout << G1.toString(r, 10);

ASSERT_TRUE(G1.eq(ra, ref));

delete[] bases;
delete[] scalars;
}

TEST(altBn128, multiExp8MSM) {

int NMExp = 8;

AltBn128::FrElement *scalars = new AltBn128::FrElement[NMExp];
Expand Down Expand Up @@ -435,11 +405,9 @@ TEST(altBn128, multiExp8MSM) {
F1.fromString(ref.x, "17747920359253913546551417160303297937542312574889904290131615776238588901697");
F1.fromString(ref.y, "8815119438581789680513912776342567599606944899217792926373871775002956510503");

G1.multiMulByScalarMSM(r, bases, (uint8_t *)scalars, 32, NMExp);
G1.multiMulByScalarMSM(r, bases, (uint8_t *)scalars, MP_N, NMExp);
G1.copy(ra, r);

// std::cout << G1.toString(r, 10);

ASSERT_TRUE(G1.eq(ra, ref));

delete[] bases;
Expand Down Expand Up @@ -469,7 +437,6 @@ TEST(altBn128, fft) {
delete[] a;
}


} // namespace

int main(int argc, char **argv) {
Expand Down
Loading