diff options
| author | MITSUNARI Shigeo <herumi@nifty.com> | 2018-10-21 15:47:03 +0800 |
|---|---|---|
| committer | MITSUNARI Shigeo <herumi@nifty.com> | 2018-10-21 15:47:03 +0800 |
| commit | 0cc0dfee05b74c99d4ad5f85c58b2dbb118bb668 (patch) | |
| tree | 9cff87669ebf99fb3bc0c8802e7c735f9e21bbc7 | |
| parent | b6e26d80d8c20d899e0d0aff4b0bb5c4e4a175ec (diff) | |
| download | dexon-mcl-0cc0dfee05b74c99d4ad5f85c58b2dbb118bb668.tar.gz dexon-mcl-0cc0dfee05b74c99d4ad5f85c58b2dbb118bb668.tar.zst dexon-mcl-0cc0dfee05b74c99d4ad5f85c58b2dbb118bb668.zip | |
add test for bls12-381 prime
| -rw-r--r-- | test/mont_fp_test.cpp | 69 |
1 files changed, 69 insertions, 0 deletions
diff --git a/test/mont_fp_test.cpp b/test/mont_fp_test.cpp index 582bae2..b96cc14 100644 --- a/test/mont_fp_test.cpp +++ b/test/mont_fp_test.cpp @@ -1,12 +1,19 @@ #define PUT(x) std::cout << #x "=" << (x) << std::endl #include <cybozu/test.hpp> #include <cybozu/benchmark.hpp> +#include <cybozu/xorshift.hpp> #include <time.h> +#if 0 +#include <mcl/bls12_381.hpp> +using namespace mcl::bls12; +typedef Fr Zn; +#else #include <mcl/fp.hpp> struct ZnTag; typedef mcl::FpT<ZnTag> Zn; typedef mcl::FpT<> Fp; +#endif struct Montgomery { typedef mcl::fp::Unit Unit; @@ -67,6 +74,25 @@ struct Montgomery { } #endif } + void mod(mpz_class& z, const mpz_class& xy) const + { + z = xy; + for (size_t i = 0; i < pn_; i++) { +//printf("i=%zd\n", i); +//std::cout << "z=" << std::hex << z << std::endl; + Unit q = mcl::gmp::getUnit(z, 0) * rp_; +//std::cout << "q=" << q << std::endl; + mpz_class t; + mcl::gmp::set(t, q); + z += p_ * t; + z >>= sizeof(Unit) * 8; +//std::cout << "z=" << std::hex << z << std::endl; + } + if (z >= p_) { + z -= p_; + } +//std::cout << "z=" << std::hex << z << std::endl; + } }; template<class T> @@ -217,6 +243,7 @@ CYBOZU_TEST_AUTO(test) { Test test; const char *tbl[] = { +#if 1 // N = 2 "0x0000000000000001000000000000000d", "0x7fffffffffffffffffffffffffffffff", @@ -241,9 +268,11 @@ CYBOZU_TEST_AUTO(test) "0x7523648240000001ba344d80000000086121000000000013a700000000000017", "0x800000000000000000000000000000000000000000000000000000000000005f", "0xffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff43", // max prime +#endif #if MCL_MAX_BIT_SIZE >= 384 // N = 6 + "0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffaaab", "0xfffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffeffffffff0000000000000000ffffffff", #endif @@ -254,6 +283,46 @@ CYBOZU_TEST_AUTO(test) }; for (size_t i = 0; i < CYBOZU_NUM_OF_ARRAY(tbl); i++) { printf("prime=%s\n", tbl[i]); +#if 0 + mpz_class p(tbl[i]); + initPairing(mcl::BLS12_381); +printf("FpDbl::mod=%p\n", FpDbl::mod); +#if 1 + cybozu::XorShift rg; + for (int i = 0; i < 1000; i++) { + Fp x, y, z; + FpDbl xy; + x.setByCSPRNG(rg); + y.setByCSPRNG(rg); + FpDbl::mulPre(xy, x, y); + FpDbl::mod(z, xy); + if (z != x * y) { + puts("ERR"); + std::cout << std::hex; + PUT(x); + PUT(y); + PUT(z); + PUT(x * y); + exit(1); + } + } +#else + Montgomery mont(p); + mpz_class xy1("00783239517141aafb2ffc854bebe27a47655610fa048c1bb6482ba78a79f232c9043ac4a898274ced5b2f65cf281204f4266605b5602ca09cacd1ea9e650c7936a436836ecfd3861d601ef424dbd6e8dbc8c71e9bc7da343f640cf5dffb903f", 16); +// mpz_class xy1 = (mpz_class(1) << (64 * 12 - 4)) - 1; + mpz_class z1; + mont.mod(z1, xy1); +std::cout << "1=" << z1 << std::endl; + FpDbl xy2, z2; + Fp z3; + z3.clear(); + xy2.setMpz(xy1); +xy2.dump(); + FpDbl::mod(z3, xy2); + z3.dump(); +#endif + exit(1); +#endif test.run(tbl[i]); } } |