2018-01-07 18:55:52 -05:00
|
|
|
// Copyright (c) 2017-2018, The Monero Project
|
2017-11-30 14:59:10 -05:00
|
|
|
//
|
|
|
|
// All rights reserved.
|
|
|
|
//
|
|
|
|
// Redistribution and use in source and binary forms, with or without modification, are
|
|
|
|
// permitted provided that the following conditions are met:
|
|
|
|
//
|
|
|
|
// 1. Redistributions of source code must retain the above copyright notice, this list of
|
|
|
|
// conditions and the following disclaimer.
|
|
|
|
//
|
|
|
|
// 2. Redistributions in binary form must reproduce the above copyright notice, this list
|
|
|
|
// of conditions and the following disclaimer in the documentation and/or other
|
|
|
|
// materials provided with the distribution.
|
|
|
|
//
|
|
|
|
// 3. Neither the name of the copyright holder nor the names of its contributors may be
|
|
|
|
// used to endorse or promote products derived from this software without specific
|
|
|
|
// prior written permission.
|
|
|
|
//
|
|
|
|
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY
|
|
|
|
// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF
|
|
|
|
// MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL
|
|
|
|
// THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
|
|
|
|
// SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
|
|
|
|
// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
|
|
|
|
// INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
|
|
|
|
// STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF
|
|
|
|
// THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
|
|
|
//
|
|
|
|
// Adapted from Java code by Sarang Noether
|
|
|
|
|
|
|
|
#include <stdlib.h>
|
|
|
|
#include <openssl/ssl.h>
|
|
|
|
#include <boost/thread/mutex.hpp>
|
|
|
|
#include "misc_log_ex.h"
|
|
|
|
#include "common/perf_timer.h"
|
|
|
|
extern "C"
|
|
|
|
{
|
|
|
|
#include "crypto/crypto-ops.h"
|
|
|
|
}
|
|
|
|
#include "rctOps.h"
|
|
|
|
#include "bulletproofs.h"
|
|
|
|
|
|
|
|
#undef MONERO_DEFAULT_LOG_CATEGORY
|
|
|
|
#define MONERO_DEFAULT_LOG_CATEGORY "bulletproofs"
|
|
|
|
|
|
|
|
//#define DEBUG_BP
|
|
|
|
|
|
|
|
#define PERF_TIMER_START_BP(x) PERF_TIMER_START_UNIT(x, 1000000)
|
|
|
|
|
|
|
|
namespace rct
|
|
|
|
{
|
|
|
|
|
|
|
|
static rct::key vector_exponent(const rct::keyV &a, const rct::keyV &b);
|
2018-01-03 16:37:18 -05:00
|
|
|
static rct::keyV vector_powers(const rct::key &x, size_t n);
|
|
|
|
static rct::keyV vector_dup(const rct::key &x, size_t n);
|
2017-11-30 14:59:10 -05:00
|
|
|
static rct::key inner_product(const rct::keyV &a, const rct::keyV &b);
|
|
|
|
|
|
|
|
static constexpr size_t maxN = 64;
|
2018-01-03 16:37:18 -05:00
|
|
|
static constexpr size_t maxM = 16;
|
|
|
|
static rct::key Hi[maxN*maxM], Gi[maxN*maxM];
|
|
|
|
static ge_dsmp Gprecomp[maxN*maxM], Hprecomp[maxN*maxM];
|
2017-11-30 14:59:10 -05:00
|
|
|
static const rct::key TWO = { {0x02, 0x00, 0x00,0x00 , 0x00, 0x00, 0x00,0x00 , 0x00, 0x00, 0x00,0x00 , 0x00, 0x00, 0x00,0x00 , 0x00, 0x00, 0x00,0x00 , 0x00, 0x00, 0x00,0x00 , 0x00, 0x00, 0x00,0x00 , 0x00, 0x00, 0x00,0x00 } };
|
2018-01-03 16:37:18 -05:00
|
|
|
static const rct::keyV oneN = vector_dup(rct::identity(), maxN);
|
2017-11-30 14:59:10 -05:00
|
|
|
static const rct::keyV twoN = vector_powers(TWO, maxN);
|
|
|
|
static const rct::key ip12 = inner_product(oneN, twoN);
|
|
|
|
static boost::mutex init_mutex;
|
|
|
|
|
2018-01-05 17:39:59 -05:00
|
|
|
//addKeys3acc_p3
|
|
|
|
//aAbB += a*A + b*B where a, b are scalars, A, B are curve points
|
|
|
|
//A and B must be input after applying "precomp"
|
|
|
|
static void addKeys3acc_p3(ge_p3 *aAbB, const key &a, const ge_dsmp A, const key &b, const ge_dsmp B)
|
|
|
|
{
|
|
|
|
ge_p3 rv;
|
|
|
|
ge_p1p1 p1;
|
|
|
|
ge_p2 p2;
|
|
|
|
ge_double_scalarmult_precomp_vartime2_p3(&rv, a.bytes, A, b.bytes, B);
|
|
|
|
ge_cached cached;
|
|
|
|
ge_p3_to_cached(&cached, aAbB);
|
|
|
|
ge_add(&p1, &rv, &cached);
|
|
|
|
ge_p1p1_to_p3(aAbB, &p1);
|
|
|
|
}
|
|
|
|
|
|
|
|
static void addKeys_acc_p3(ge_p3 *acc_p3, const rct::key &a, const rct::key &point)
|
|
|
|
{
|
|
|
|
ge_p3 p3;
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(ge_frombytes_vartime(&p3, point.bytes) == 0, "ge_frombytes_vartime failed");
|
|
|
|
ge_scalarmult_p3(&p3, a.bytes, &p3);
|
|
|
|
ge_cached cached;
|
|
|
|
ge_p3_to_cached(&cached, acc_p3);
|
|
|
|
ge_p1p1 p1;
|
|
|
|
ge_add(&p1, &p3, &cached);
|
|
|
|
ge_p1p1_to_p3(acc_p3, &p1);
|
|
|
|
}
|
|
|
|
|
2017-11-30 14:59:10 -05:00
|
|
|
static rct::key get_exponent(const rct::key &base, size_t idx)
|
|
|
|
{
|
|
|
|
static const std::string salt("bulletproof");
|
|
|
|
std::string hashed = std::string((const char*)base.bytes, sizeof(base)) + salt + tools::get_varint_data(idx);
|
|
|
|
return rct::hashToPoint(rct::hash2rct(crypto::cn_fast_hash(hashed.data(), hashed.size())));
|
|
|
|
}
|
|
|
|
|
|
|
|
static void init_exponents()
|
|
|
|
{
|
|
|
|
boost::lock_guard<boost::mutex> lock(init_mutex);
|
|
|
|
|
|
|
|
static bool init_done = false;
|
|
|
|
if (init_done)
|
|
|
|
return;
|
2018-01-03 16:37:18 -05:00
|
|
|
for (size_t i = 0; i < maxN*maxM; ++i)
|
2017-11-30 14:59:10 -05:00
|
|
|
{
|
|
|
|
Hi[i] = get_exponent(rct::H, i * 2);
|
|
|
|
rct::precomp(Hprecomp[i], Hi[i]);
|
|
|
|
Gi[i] = get_exponent(rct::H, i * 2 + 1);
|
|
|
|
rct::precomp(Gprecomp[i], Gi[i]);
|
|
|
|
}
|
|
|
|
init_done = true;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Given two scalar arrays, construct a vector commitment */
|
|
|
|
static rct::key vector_exponent(const rct::keyV &a, const rct::keyV &b)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(a.size() == b.size(), "Incompatible sizes of a and b");
|
2018-01-03 16:37:18 -05:00
|
|
|
CHECK_AND_ASSERT_THROW_MES(a.size() <= maxN*maxM, "Incompatible sizes of a and maxN");
|
2018-01-05 17:39:59 -05:00
|
|
|
ge_p3 res_p3 = ge_p3_identity;
|
2017-11-30 14:59:10 -05:00
|
|
|
for (size_t i = 0; i < a.size(); ++i)
|
|
|
|
{
|
2018-01-05 17:39:59 -05:00
|
|
|
rct::addKeys3acc_p3(&res_p3, a[i], Gprecomp[i], b[i], Hprecomp[i]);
|
2017-11-30 14:59:10 -05:00
|
|
|
}
|
2018-01-05 17:39:59 -05:00
|
|
|
rct::key res;
|
|
|
|
ge_p3_tobytes(res.bytes, &res_p3);
|
2017-11-30 14:59:10 -05:00
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Compute a custom vector-scalar commitment */
|
|
|
|
static rct::key vector_exponent_custom(const rct::keyV &A, const rct::keyV &B, const rct::keyV &a, const rct::keyV &b)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(A.size() == B.size(), "Incompatible sizes of A and B");
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(a.size() == b.size(), "Incompatible sizes of a and b");
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(a.size() == A.size(), "Incompatible sizes of a and A");
|
2018-01-03 16:37:18 -05:00
|
|
|
CHECK_AND_ASSERT_THROW_MES(a.size() <= maxN*maxM, "Incompatible sizes of a and maxN");
|
2018-01-05 17:39:59 -05:00
|
|
|
ge_p3 res_p3 = ge_p3_identity;
|
2017-11-30 14:59:10 -05:00
|
|
|
for (size_t i = 0; i < a.size(); ++i)
|
|
|
|
{
|
|
|
|
#if 0
|
2018-01-05 17:39:59 -05:00
|
|
|
rct::key term;
|
2017-11-30 14:59:10 -05:00
|
|
|
// we happen to know where A and B might fall, so don't bother checking the rest
|
|
|
|
ge_dsmp *Acache = NULL, *Bcache = NULL;
|
|
|
|
ge_dsmp Acache_custom[1], Bcache_custom[1];
|
|
|
|
if (Gi[i] == A[i])
|
|
|
|
Acache = Gprecomp + i;
|
|
|
|
else if (i<32 && Gi[i+32] == A[i])
|
|
|
|
Acache = Gprecomp + i + 32;
|
|
|
|
else
|
|
|
|
{
|
|
|
|
rct::precomp(Acache_custom[0], A[i]);
|
|
|
|
Acache = Acache_custom;
|
|
|
|
}
|
|
|
|
if (i == 0 && B[i] == Hi[0])
|
|
|
|
Bcache = Hprecomp;
|
|
|
|
else
|
|
|
|
{
|
|
|
|
rct::precomp(Bcache_custom[0], B[i]);
|
|
|
|
Bcache = Bcache_custom;
|
|
|
|
}
|
|
|
|
rct::addKeys3(term, a[i], *Acache, b[i], *Bcache);
|
2018-01-05 17:39:59 -05:00
|
|
|
rct::addKeys(res, res, term);
|
2017-11-30 14:59:10 -05:00
|
|
|
#else
|
|
|
|
ge_dsmp Acache, Bcache;
|
|
|
|
rct::precomp(Bcache, B[i]);
|
2018-01-05 17:39:59 -05:00
|
|
|
rct::precomp(Acache, A[i]);
|
|
|
|
addKeys3acc_p3(&res_p3, a[i], Acache, b[i], Bcache);
|
2017-11-30 14:59:10 -05:00
|
|
|
#endif
|
|
|
|
}
|
2018-01-05 17:39:59 -05:00
|
|
|
rct::key res;
|
|
|
|
ge_p3_tobytes(res.bytes, &res_p3);
|
2017-11-30 14:59:10 -05:00
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Given a scalar, construct a vector of powers */
|
2018-01-03 16:37:18 -05:00
|
|
|
static rct::keyV vector_powers(const rct::key &x, size_t n)
|
2017-11-30 14:59:10 -05:00
|
|
|
{
|
|
|
|
rct::keyV res(n);
|
|
|
|
if (n == 0)
|
|
|
|
return res;
|
|
|
|
res[0] = rct::identity();
|
|
|
|
if (n == 1)
|
|
|
|
return res;
|
|
|
|
res[1] = x;
|
|
|
|
for (size_t i = 2; i < n; ++i)
|
|
|
|
{
|
|
|
|
sc_mul(res[i].bytes, res[i-1].bytes, x.bytes);
|
|
|
|
}
|
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
2018-01-05 17:39:59 -05:00
|
|
|
/* Given a scalar, return the sum of its powers from 0 to n-1 */
|
|
|
|
static rct::key vector_power_sum(const rct::key &x, size_t n)
|
|
|
|
{
|
|
|
|
if (n == 0)
|
|
|
|
return rct::zero();
|
|
|
|
rct::key res = rct::identity();
|
|
|
|
if (n == 1)
|
|
|
|
return res;
|
|
|
|
rct::key prev = x;
|
|
|
|
for (size_t i = 1; i < n; ++i)
|
|
|
|
{
|
|
|
|
if (i > 1)
|
|
|
|
sc_mul(prev.bytes, prev.bytes, x.bytes);
|
|
|
|
sc_add(res.bytes, res.bytes, prev.bytes);
|
|
|
|
}
|
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
2017-11-30 14:59:10 -05:00
|
|
|
/* Given two scalar arrays, construct the inner product */
|
|
|
|
static rct::key inner_product(const rct::keyV &a, const rct::keyV &b)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(a.size() == b.size(), "Incompatible sizes of a and b");
|
|
|
|
rct::key res = rct::zero();
|
|
|
|
for (size_t i = 0; i < a.size(); ++i)
|
|
|
|
{
|
|
|
|
sc_muladd(res.bytes, a[i].bytes, b[i].bytes, res.bytes);
|
|
|
|
}
|
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Given two scalar arrays, construct the Hadamard product */
|
|
|
|
static rct::keyV hadamard(const rct::keyV &a, const rct::keyV &b)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(a.size() == b.size(), "Incompatible sizes of a and b");
|
|
|
|
rct::keyV res(a.size());
|
|
|
|
for (size_t i = 0; i < a.size(); ++i)
|
|
|
|
{
|
|
|
|
sc_mul(res[i].bytes, a[i].bytes, b[i].bytes);
|
|
|
|
}
|
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Given two curvepoint arrays, construct the Hadamard product */
|
|
|
|
static rct::keyV hadamard2(const rct::keyV &a, const rct::keyV &b)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(a.size() == b.size(), "Incompatible sizes of a and b");
|
|
|
|
rct::keyV res(a.size());
|
|
|
|
for (size_t i = 0; i < a.size(); ++i)
|
|
|
|
{
|
|
|
|
rct::addKeys(res[i], a[i], b[i]);
|
|
|
|
}
|
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Add two vectors */
|
|
|
|
static rct::keyV vector_add(const rct::keyV &a, const rct::keyV &b)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(a.size() == b.size(), "Incompatible sizes of a and b");
|
|
|
|
rct::keyV res(a.size());
|
|
|
|
for (size_t i = 0; i < a.size(); ++i)
|
|
|
|
{
|
|
|
|
sc_add(res[i].bytes, a[i].bytes, b[i].bytes);
|
|
|
|
}
|
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Subtract two vectors */
|
|
|
|
static rct::keyV vector_subtract(const rct::keyV &a, const rct::keyV &b)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(a.size() == b.size(), "Incompatible sizes of a and b");
|
|
|
|
rct::keyV res(a.size());
|
|
|
|
for (size_t i = 0; i < a.size(); ++i)
|
|
|
|
{
|
|
|
|
sc_sub(res[i].bytes, a[i].bytes, b[i].bytes);
|
|
|
|
}
|
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Multiply a scalar and a vector */
|
|
|
|
static rct::keyV vector_scalar(const rct::keyV &a, const rct::key &x)
|
|
|
|
{
|
|
|
|
rct::keyV res(a.size());
|
|
|
|
for (size_t i = 0; i < a.size(); ++i)
|
|
|
|
{
|
|
|
|
sc_mul(res[i].bytes, a[i].bytes, x.bytes);
|
|
|
|
}
|
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
2018-01-03 16:37:18 -05:00
|
|
|
/* Create a vector from copies of a single value */
|
|
|
|
static rct::keyV vector_dup(const rct::key &x, size_t N)
|
|
|
|
{
|
|
|
|
return rct::keyV(N, x);
|
|
|
|
}
|
|
|
|
|
2017-11-30 14:59:10 -05:00
|
|
|
/* Exponentiate a curve vector by a scalar */
|
|
|
|
static rct::keyV vector_scalar2(const rct::keyV &a, const rct::key &x)
|
|
|
|
{
|
|
|
|
rct::keyV res(a.size());
|
|
|
|
for (size_t i = 0; i < a.size(); ++i)
|
|
|
|
{
|
|
|
|
rct::scalarmultKey(res[i], a[i], x);
|
|
|
|
}
|
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
2018-01-03 16:37:18 -05:00
|
|
|
/* Get the sum of a vector's elements */
|
|
|
|
static rct::key vector_sum(const rct::keyV &a)
|
|
|
|
{
|
|
|
|
rct::key res = rct::zero();
|
|
|
|
for (size_t i = 0; i < a.size(); ++i)
|
|
|
|
{
|
|
|
|
sc_add(res.bytes, res.bytes, a[i].bytes);
|
|
|
|
}
|
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
2017-11-30 14:59:10 -05:00
|
|
|
static rct::key switch_endianness(rct::key k)
|
|
|
|
{
|
|
|
|
std::reverse(k.bytes, k.bytes + sizeof(k));
|
|
|
|
return k;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Compute the inverse of a scalar, the stupid way */
|
|
|
|
static rct::key invert(const rct::key &x)
|
|
|
|
{
|
|
|
|
rct::key inv;
|
|
|
|
|
|
|
|
BN_CTX *ctx = BN_CTX_new();
|
|
|
|
BIGNUM *X = BN_new();
|
|
|
|
BIGNUM *L = BN_new();
|
|
|
|
BIGNUM *I = BN_new();
|
|
|
|
|
|
|
|
BN_bin2bn(switch_endianness(x).bytes, sizeof(rct::key), X);
|
|
|
|
BN_bin2bn(switch_endianness(rct::curveOrder()).bytes, sizeof(rct::key), L);
|
|
|
|
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(BN_mod_inverse(I, X, L, ctx), "Failed to invert");
|
|
|
|
|
|
|
|
const int len = BN_num_bytes(I);
|
|
|
|
CHECK_AND_ASSERT_THROW_MES((size_t)len <= sizeof(rct::key), "Invalid number length");
|
|
|
|
inv = rct::zero();
|
|
|
|
BN_bn2bin(I, inv.bytes);
|
|
|
|
std::reverse(inv.bytes, inv.bytes + len);
|
|
|
|
|
|
|
|
BN_free(I);
|
|
|
|
BN_free(L);
|
|
|
|
BN_free(X);
|
|
|
|
BN_CTX_free(ctx);
|
|
|
|
|
|
|
|
#ifdef DEBUG_BP
|
|
|
|
rct::key tmp;
|
|
|
|
sc_mul(tmp.bytes, inv.bytes, x.bytes);
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(tmp == rct::identity(), "invert failed");
|
|
|
|
#endif
|
|
|
|
return inv;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Compute the slice of a vector */
|
|
|
|
static rct::keyV slice(const rct::keyV &a, size_t start, size_t stop)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(start < a.size(), "Invalid start index");
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(stop <= a.size(), "Invalid stop index");
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(start < stop, "Invalid start/stop indices");
|
|
|
|
rct::keyV res(stop - start);
|
|
|
|
for (size_t i = start; i < stop; ++i)
|
|
|
|
{
|
|
|
|
res[i - start] = a[i];
|
|
|
|
}
|
|
|
|
return res;
|
|
|
|
}
|
|
|
|
|
2017-12-18 16:58:45 -05:00
|
|
|
static rct::key hash_cache_mash(rct::key &hash_cache, const rct::key &mash0, const rct::key &mash1)
|
|
|
|
{
|
|
|
|
rct::keyV data;
|
|
|
|
data.reserve(3);
|
|
|
|
data.push_back(hash_cache);
|
|
|
|
data.push_back(mash0);
|
|
|
|
data.push_back(mash1);
|
|
|
|
return hash_cache = rct::hash_to_scalar(data);
|
|
|
|
}
|
|
|
|
|
|
|
|
static rct::key hash_cache_mash(rct::key &hash_cache, const rct::key &mash0, const rct::key &mash1, const rct::key &mash2)
|
|
|
|
{
|
|
|
|
rct::keyV data;
|
|
|
|
data.reserve(4);
|
|
|
|
data.push_back(hash_cache);
|
|
|
|
data.push_back(mash0);
|
|
|
|
data.push_back(mash1);
|
|
|
|
data.push_back(mash2);
|
|
|
|
return hash_cache = rct::hash_to_scalar(data);
|
|
|
|
}
|
|
|
|
|
|
|
|
static rct::key hash_cache_mash(rct::key &hash_cache, const rct::key &mash0, const rct::key &mash1, const rct::key &mash2, const rct::key &mash3)
|
|
|
|
{
|
|
|
|
rct::keyV data;
|
|
|
|
data.reserve(5);
|
|
|
|
data.push_back(hash_cache);
|
|
|
|
data.push_back(mash0);
|
|
|
|
data.push_back(mash1);
|
|
|
|
data.push_back(mash2);
|
|
|
|
data.push_back(mash3);
|
|
|
|
return hash_cache = rct::hash_to_scalar(data);
|
|
|
|
}
|
|
|
|
|
2017-11-30 14:59:10 -05:00
|
|
|
/* Given a value v (0..2^N-1) and a mask gamma, construct a range proof */
|
|
|
|
Bulletproof bulletproof_PROVE(const rct::key &sv, const rct::key &gamma)
|
|
|
|
{
|
|
|
|
init_exponents();
|
|
|
|
|
|
|
|
PERF_TIMER_UNIT(PROVE, 1000000);
|
|
|
|
|
|
|
|
constexpr size_t logN = 6; // log2(64)
|
|
|
|
constexpr size_t N = 1<<logN;
|
|
|
|
|
|
|
|
rct::key V;
|
|
|
|
rct::keyV aL(N), aR(N);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_v);
|
2017-12-02 16:17:06 -05:00
|
|
|
rct::addKeys2(V, gamma, sv, rct::H);
|
2017-11-30 14:59:10 -05:00
|
|
|
PERF_TIMER_STOP(PROVE_v);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_aLaR);
|
|
|
|
for (size_t i = N; i-- > 0; )
|
|
|
|
{
|
|
|
|
if (sv[i/8] & (((uint64_t)1)<<(i%8)))
|
|
|
|
{
|
|
|
|
aL[i] = rct::identity();
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
aL[i] = rct::zero();
|
|
|
|
}
|
|
|
|
sc_sub(aR[i].bytes, aL[i].bytes, rct::identity().bytes);
|
|
|
|
}
|
|
|
|
PERF_TIMER_STOP(PROVE_aLaR);
|
|
|
|
|
2017-12-18 16:58:45 -05:00
|
|
|
rct::key hash_cache = rct::hash_to_scalar(V);
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
// DEBUG: Test to ensure this recovers the value
|
|
|
|
#ifdef DEBUG_BP
|
|
|
|
uint64_t test_aL = 0, test_aR = 0;
|
|
|
|
for (size_t i = 0; i < N; ++i)
|
|
|
|
{
|
|
|
|
if (aL[i] == rct::identity())
|
|
|
|
test_aL += ((uint64_t)1)<<i;
|
|
|
|
if (aR[i] == rct::zero())
|
|
|
|
test_aR += ((uint64_t)1)<<i;
|
|
|
|
}
|
|
|
|
uint64_t v_test = 0;
|
|
|
|
for (int n = 0; n < 8; ++n) v_test |= (((uint64_t)sv[n]) << (8*n));
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(test_aL == v_test, "test_aL failed");
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(test_aR == v_test, "test_aR failed");
|
|
|
|
#endif
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_step1);
|
|
|
|
// PAPER LINES 38-39
|
|
|
|
rct::key alpha = rct::skGen();
|
|
|
|
rct::key ve = vector_exponent(aL, aR);
|
|
|
|
rct::key A;
|
2017-12-02 16:17:06 -05:00
|
|
|
rct::addKeys(A, ve, rct::scalarmultBase(alpha));
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
// PAPER LINES 40-42
|
|
|
|
rct::keyV sL = rct::skvGen(N), sR = rct::skvGen(N);
|
|
|
|
rct::key rho = rct::skGen();
|
|
|
|
ve = vector_exponent(sL, sR);
|
|
|
|
rct::key S;
|
2017-12-02 16:17:06 -05:00
|
|
|
rct::addKeys(S, ve, rct::scalarmultBase(rho));
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
// PAPER LINES 43-45
|
2017-12-18 16:58:45 -05:00
|
|
|
rct::key y = hash_cache_mash(hash_cache, A, S);
|
|
|
|
rct::key z = hash_cache = rct::hash_to_scalar(y);
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
// Polynomial construction before PAPER LINE 46
|
|
|
|
rct::key t0 = rct::zero();
|
|
|
|
rct::key t1 = rct::zero();
|
|
|
|
rct::key t2 = rct::zero();
|
|
|
|
|
|
|
|
const auto yN = vector_powers(y, N);
|
|
|
|
|
2018-01-03 16:37:18 -05:00
|
|
|
rct::key ip1y = vector_sum(yN);
|
2017-11-30 14:59:10 -05:00
|
|
|
rct::key tmp;
|
|
|
|
sc_muladd(t0.bytes, z.bytes, ip1y.bytes, t0.bytes);
|
|
|
|
|
|
|
|
rct::key zsq;
|
|
|
|
sc_mul(zsq.bytes, z.bytes, z.bytes);
|
|
|
|
sc_muladd(t0.bytes, zsq.bytes, sv.bytes, t0.bytes);
|
|
|
|
|
|
|
|
rct::key k = rct::zero();
|
|
|
|
sc_mulsub(k.bytes, zsq.bytes, ip1y.bytes, k.bytes);
|
|
|
|
|
|
|
|
rct::key zcu;
|
|
|
|
sc_mul(zcu.bytes, zsq.bytes, z.bytes);
|
|
|
|
sc_mulsub(k.bytes, zcu.bytes, ip12.bytes, k.bytes);
|
|
|
|
sc_add(t0.bytes, t0.bytes, k.bytes);
|
|
|
|
|
|
|
|
// DEBUG: Test the value of t0 has the correct form
|
|
|
|
#ifdef DEBUG_BP
|
|
|
|
rct::key test_t0 = rct::zero();
|
|
|
|
rct::key iph = inner_product(aL, hadamard(aR, yN));
|
|
|
|
sc_add(test_t0.bytes, test_t0.bytes, iph.bytes);
|
|
|
|
rct::key ips = inner_product(vector_subtract(aL, aR), yN);
|
|
|
|
sc_muladd(test_t0.bytes, z.bytes, ips.bytes, test_t0.bytes);
|
|
|
|
rct::key ipt = inner_product(twoN, aL);
|
|
|
|
sc_muladd(test_t0.bytes, zsq.bytes, ipt.bytes, test_t0.bytes);
|
|
|
|
sc_add(test_t0.bytes, test_t0.bytes, k.bytes);
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(t0 == test_t0, "t0 check failed");
|
|
|
|
#endif
|
|
|
|
PERF_TIMER_STOP(PROVE_step1);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_step2);
|
|
|
|
const auto HyNsR = hadamard(yN, sR);
|
2018-01-03 16:37:18 -05:00
|
|
|
const auto vpIz = vector_dup(z, N);
|
2017-11-30 14:59:10 -05:00
|
|
|
const auto vp2zsq = vector_scalar(twoN, zsq);
|
|
|
|
const auto aL_vpIz = vector_subtract(aL, vpIz);
|
|
|
|
const auto aR_vpIz = vector_add(aR, vpIz);
|
|
|
|
|
|
|
|
rct::key ip1 = inner_product(aL_vpIz, HyNsR);
|
|
|
|
sc_add(t1.bytes, t1.bytes, ip1.bytes);
|
|
|
|
|
|
|
|
rct::key ip2 = inner_product(sL, vector_add(hadamard(yN, aR_vpIz), vp2zsq));
|
|
|
|
sc_add(t1.bytes, t1.bytes, ip2.bytes);
|
|
|
|
|
|
|
|
rct::key ip3 = inner_product(sL, HyNsR);
|
|
|
|
sc_add(t2.bytes, t2.bytes, ip3.bytes);
|
|
|
|
|
|
|
|
// PAPER LINES 47-48
|
|
|
|
rct::key tau1 = rct::skGen(), tau2 = rct::skGen();
|
|
|
|
|
2017-12-02 16:17:06 -05:00
|
|
|
rct::key T1 = rct::addKeys(rct::scalarmultKey(rct::H, t1), rct::scalarmultBase(tau1));
|
|
|
|
rct::key T2 = rct::addKeys(rct::scalarmultKey(rct::H, t2), rct::scalarmultBase(tau2));
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
// PAPER LINES 49-51
|
2017-12-18 16:58:45 -05:00
|
|
|
rct::key x = hash_cache_mash(hash_cache, z, T1, T2);
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
// PAPER LINES 52-53
|
|
|
|
rct::key taux = rct::zero();
|
|
|
|
sc_mul(taux.bytes, tau1.bytes, x.bytes);
|
|
|
|
rct::key xsq;
|
|
|
|
sc_mul(xsq.bytes, x.bytes, x.bytes);
|
|
|
|
sc_muladd(taux.bytes, tau2.bytes, xsq.bytes, taux.bytes);
|
|
|
|
sc_muladd(taux.bytes, gamma.bytes, zsq.bytes, taux.bytes);
|
|
|
|
rct::key mu;
|
|
|
|
sc_muladd(mu.bytes, x.bytes, rho.bytes, alpha.bytes);
|
|
|
|
|
|
|
|
// PAPER LINES 54-57
|
|
|
|
rct::keyV l = vector_add(aL_vpIz, vector_scalar(sL, x));
|
|
|
|
rct::keyV r = vector_add(hadamard(yN, vector_add(aR_vpIz, vector_scalar(sR, x))), vp2zsq);
|
|
|
|
PERF_TIMER_STOP(PROVE_step2);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_step3);
|
|
|
|
rct::key t = inner_product(l, r);
|
|
|
|
|
|
|
|
// DEBUG: Test if the l and r vectors match the polynomial forms
|
|
|
|
#ifdef DEBUG_BP
|
|
|
|
rct::key test_t;
|
|
|
|
sc_muladd(test_t.bytes, t1.bytes, x.bytes, t0.bytes);
|
|
|
|
sc_muladd(test_t.bytes, t2.bytes, xsq.bytes, test_t.bytes);
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(test_t == t, "test_t check failed");
|
|
|
|
#endif
|
|
|
|
|
|
|
|
// PAPER LINES 32-33
|
2017-12-18 16:58:45 -05:00
|
|
|
rct::key x_ip = hash_cache_mash(hash_cache, x, taux, mu, t);
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
// These are used in the inner product rounds
|
|
|
|
size_t nprime = N;
|
|
|
|
rct::keyV Gprime(N);
|
|
|
|
rct::keyV Hprime(N);
|
|
|
|
rct::keyV aprime(N);
|
|
|
|
rct::keyV bprime(N);
|
|
|
|
const rct::key yinv = invert(y);
|
|
|
|
rct::key yinvpow = rct::identity();
|
|
|
|
for (size_t i = 0; i < N; ++i)
|
|
|
|
{
|
|
|
|
Gprime[i] = Gi[i];
|
|
|
|
Hprime[i] = scalarmultKey(Hi[i], yinvpow);
|
|
|
|
sc_mul(yinvpow.bytes, yinvpow.bytes, yinv.bytes);
|
|
|
|
aprime[i] = l[i];
|
|
|
|
bprime[i] = r[i];
|
|
|
|
}
|
|
|
|
rct::keyV L(logN);
|
|
|
|
rct::keyV R(logN);
|
|
|
|
int round = 0;
|
|
|
|
rct::keyV w(logN); // this is the challenge x in the inner product protocol
|
|
|
|
PERF_TIMER_STOP(PROVE_step3);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_step4);
|
|
|
|
// PAPER LINE 13
|
|
|
|
while (nprime > 1)
|
|
|
|
{
|
|
|
|
// PAPER LINE 15
|
|
|
|
nprime /= 2;
|
|
|
|
|
|
|
|
// PAPER LINES 16-17
|
|
|
|
rct::key cL = inner_product(slice(aprime, 0, nprime), slice(bprime, nprime, bprime.size()));
|
|
|
|
rct::key cR = inner_product(slice(aprime, nprime, aprime.size()), slice(bprime, 0, nprime));
|
|
|
|
|
|
|
|
// PAPER LINES 18-19
|
|
|
|
L[round] = vector_exponent_custom(slice(Gprime, nprime, Gprime.size()), slice(Hprime, 0, nprime), slice(aprime, 0, nprime), slice(bprime, nprime, bprime.size()));
|
|
|
|
sc_mul(tmp.bytes, cL.bytes, x_ip.bytes);
|
2017-12-02 16:17:06 -05:00
|
|
|
rct::addKeys(L[round], L[round], rct::scalarmultKey(rct::H, tmp));
|
2017-11-30 14:59:10 -05:00
|
|
|
R[round] = vector_exponent_custom(slice(Gprime, 0, nprime), slice(Hprime, nprime, Hprime.size()), slice(aprime, nprime, aprime.size()), slice(bprime, 0, nprime));
|
|
|
|
sc_mul(tmp.bytes, cR.bytes, x_ip.bytes);
|
2017-12-02 16:17:06 -05:00
|
|
|
rct::addKeys(R[round], R[round], rct::scalarmultKey(rct::H, tmp));
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
// PAPER LINES 21-22
|
2017-12-18 16:58:45 -05:00
|
|
|
w[round] = hash_cache_mash(hash_cache, L[round], R[round]);
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
// PAPER LINES 24-25
|
|
|
|
const rct::key winv = invert(w[round]);
|
|
|
|
Gprime = hadamard2(vector_scalar2(slice(Gprime, 0, nprime), winv), vector_scalar2(slice(Gprime, nprime, Gprime.size()), w[round]));
|
|
|
|
Hprime = hadamard2(vector_scalar2(slice(Hprime, 0, nprime), w[round]), vector_scalar2(slice(Hprime, nprime, Hprime.size()), winv));
|
|
|
|
|
|
|
|
// PAPER LINES 28-29
|
|
|
|
aprime = vector_add(vector_scalar(slice(aprime, 0, nprime), w[round]), vector_scalar(slice(aprime, nprime, aprime.size()), winv));
|
|
|
|
bprime = vector_add(vector_scalar(slice(bprime, 0, nprime), winv), vector_scalar(slice(bprime, nprime, bprime.size()), w[round]));
|
|
|
|
|
|
|
|
++round;
|
|
|
|
}
|
|
|
|
PERF_TIMER_STOP(PROVE_step4);
|
|
|
|
|
|
|
|
// PAPER LINE 58 (with inclusions from PAPER LINE 8 and PAPER LINE 20)
|
|
|
|
return Bulletproof(V, A, S, T1, T2, taux, mu, L, R, aprime[0], bprime[0], t);
|
|
|
|
}
|
|
|
|
|
|
|
|
Bulletproof bulletproof_PROVE(uint64_t v, const rct::key &gamma)
|
|
|
|
{
|
|
|
|
// vG + gammaH
|
|
|
|
PERF_TIMER_START_BP(PROVE_v);
|
|
|
|
rct::key sv = rct::zero();
|
|
|
|
sv.bytes[0] = v & 255;
|
|
|
|
sv.bytes[1] = (v >> 8) & 255;
|
|
|
|
sv.bytes[2] = (v >> 16) & 255;
|
|
|
|
sv.bytes[3] = (v >> 24) & 255;
|
|
|
|
sv.bytes[4] = (v >> 32) & 255;
|
|
|
|
sv.bytes[5] = (v >> 40) & 255;
|
|
|
|
sv.bytes[6] = (v >> 48) & 255;
|
|
|
|
sv.bytes[7] = (v >> 56) & 255;
|
|
|
|
PERF_TIMER_STOP(PROVE_v);
|
|
|
|
return bulletproof_PROVE(sv, gamma);
|
|
|
|
}
|
|
|
|
|
2018-01-03 16:37:18 -05:00
|
|
|
/* Given a set of values v (0..2^N-1) and masks gamma, construct a range proof */
|
|
|
|
Bulletproof bulletproof_PROVE(const rct::keyV &sv, const rct::keyV &gamma)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(sv.size() == gamma.size(), "Incompatible sizes of sv and gamma");
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(!sv.empty(), "sv is empty");
|
|
|
|
|
|
|
|
init_exponents();
|
|
|
|
|
|
|
|
PERF_TIMER_UNIT(PROVE, 1000000);
|
|
|
|
|
|
|
|
constexpr size_t logN = 6; // log2(64)
|
|
|
|
constexpr size_t N = 1<<logN;
|
|
|
|
size_t M, logM;
|
|
|
|
for (logM = 0; (M = 1<<logM) <= maxM && M < sv.size(); ++logM);
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(M <= maxM, "sv/gamma are too large");
|
|
|
|
const size_t logMN = logM + logN;
|
|
|
|
const size_t MN = M * N;
|
|
|
|
|
|
|
|
rct::keyV V(sv.size());
|
|
|
|
rct::keyV aL(MN), aR(MN);
|
|
|
|
rct::key tmp;
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_v);
|
|
|
|
for (size_t i = 0; i < sv.size(); ++i)
|
|
|
|
rct::addKeys2(V[i], gamma[i], sv[i], rct::H);
|
|
|
|
PERF_TIMER_STOP(PROVE_v);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_aLaR);
|
|
|
|
for (size_t j = 0; j < M; ++j)
|
|
|
|
{
|
|
|
|
for (size_t i = N; i-- > 0; )
|
|
|
|
{
|
|
|
|
if (j >= sv.size())
|
|
|
|
{
|
|
|
|
aL[j*N+i] = rct::zero();
|
|
|
|
}
|
|
|
|
else if (sv[j][i/8] & (((uint64_t)1)<<(i%8)))
|
|
|
|
{
|
|
|
|
aL[j*N+i] = rct::identity();
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
aL[j*N+i] = rct::zero();
|
|
|
|
}
|
|
|
|
sc_sub(aR[j*N+i].bytes, aL[j*N+i].bytes, rct::identity().bytes);
|
|
|
|
}
|
|
|
|
}
|
|
|
|
PERF_TIMER_STOP(PROVE_aLaR);
|
|
|
|
|
|
|
|
rct::key hash_cache = rct::hash_to_scalar(V);
|
|
|
|
|
|
|
|
// DEBUG: Test to ensure this recovers the value
|
|
|
|
#ifdef DEBUG_BP
|
|
|
|
for (size_t j = 0; j < M; ++j)
|
|
|
|
{
|
|
|
|
uint64_t test_aL = 0, test_aR = 0;
|
|
|
|
for (size_t i = 0; i < N; ++i)
|
|
|
|
{
|
|
|
|
if (aL[j*N+i] == rct::identity())
|
|
|
|
test_aL += ((uint64_t)1)<<i;
|
|
|
|
if (aR[j*N+i] == rct::zero())
|
|
|
|
test_aR += ((uint64_t)1)<<i;
|
|
|
|
}
|
|
|
|
uint64_t v_test = 0;
|
|
|
|
if (j < sv.size())
|
|
|
|
for (int n = 0; n < 8; ++n) v_test |= (((uint64_t)sv[j][n]) << (8*n));
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(test_aL == v_test, "test_aL failed");
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(test_aR == v_test, "test_aR failed");
|
|
|
|
}
|
|
|
|
#endif
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_step1);
|
|
|
|
// PAPER LINES 38-39
|
|
|
|
rct::key alpha = rct::skGen();
|
|
|
|
rct::key ve = vector_exponent(aL, aR);
|
|
|
|
rct::key A;
|
|
|
|
rct::addKeys(A, ve, rct::scalarmultBase(alpha));
|
|
|
|
|
|
|
|
// PAPER LINES 40-42
|
|
|
|
rct::keyV sL = rct::skvGen(MN), sR = rct::skvGen(MN);
|
|
|
|
rct::key rho = rct::skGen();
|
|
|
|
ve = vector_exponent(sL, sR);
|
|
|
|
rct::key S;
|
|
|
|
rct::addKeys(S, ve, rct::scalarmultBase(rho));
|
|
|
|
|
|
|
|
// PAPER LINES 43-45
|
|
|
|
rct::key y = hash_cache_mash(hash_cache, A, S);
|
|
|
|
rct::key z = hash_cache = rct::hash_to_scalar(y);
|
|
|
|
|
|
|
|
// Polynomial construction by coefficients
|
|
|
|
const auto zMN = vector_dup(z, MN);
|
|
|
|
rct::keyV l0 = vector_subtract(aL, zMN);
|
|
|
|
const rct::keyV &l1 = sL;
|
|
|
|
|
|
|
|
// This computes the ugly sum/concatenation from PAPER LINE 65
|
|
|
|
rct::keyV zero_twos(MN);
|
|
|
|
const rct::keyV zpow = vector_powers(z, M+2);
|
|
|
|
for (size_t i = 0; i < MN; ++i)
|
|
|
|
{
|
|
|
|
zero_twos[i] = rct::zero();
|
|
|
|
for (size_t j = 1; j <= M; ++j)
|
|
|
|
{
|
|
|
|
if (i >= (j-1)*N && i < j*N)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(1+j < zpow.size(), "invalid zpow index");
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(i-(j-1)*N < twoN.size(), "invalid twoN index");
|
|
|
|
sc_muladd(zero_twos[i].bytes, zpow[1+j].bytes, twoN[i-(j-1)*N].bytes, zero_twos[i].bytes);
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
rct::keyV r0 = vector_add(aR, zMN);
|
|
|
|
const auto yMN = vector_powers(y, MN);
|
|
|
|
r0 = hadamard(r0, yMN);
|
|
|
|
r0 = vector_add(r0, zero_twos);
|
|
|
|
rct::keyV r1 = hadamard(yMN, sR);
|
|
|
|
|
|
|
|
// Polynomial construction before PAPER LINE 46
|
|
|
|
rct::key t1_1 = inner_product(l0, r1);
|
|
|
|
rct::key t1_2 = inner_product(l1, r0);
|
|
|
|
rct::key t1;
|
|
|
|
sc_add(t1.bytes, t1_1.bytes, t1_2.bytes);
|
|
|
|
rct::key t2 = inner_product(l1, r1);
|
|
|
|
|
|
|
|
PERF_TIMER_STOP(PROVE_step1);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_step2);
|
|
|
|
// PAPER LINES 47-48
|
|
|
|
rct::key tau1 = rct::skGen(), tau2 = rct::skGen();
|
|
|
|
|
|
|
|
rct::key T1 = rct::addKeys(rct::scalarmultKey(rct::H, t1), rct::scalarmultBase(tau1));
|
|
|
|
rct::key T2 = rct::addKeys(rct::scalarmultKey(rct::H, t2), rct::scalarmultBase(tau2));
|
|
|
|
|
|
|
|
// PAPER LINES 49-51
|
|
|
|
rct::key x = hash_cache_mash(hash_cache, z, T1, T2);
|
|
|
|
|
|
|
|
// PAPER LINES 52-53
|
|
|
|
rct::key taux;
|
|
|
|
sc_mul(taux.bytes, tau1.bytes, x.bytes);
|
|
|
|
rct::key xsq;
|
|
|
|
sc_mul(xsq.bytes, x.bytes, x.bytes);
|
|
|
|
sc_muladd(taux.bytes, tau2.bytes, xsq.bytes, taux.bytes);
|
|
|
|
for (size_t j = 1; j <= sv.size(); ++j)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(j+1 < zpow.size(), "invalid zpow index");
|
|
|
|
sc_muladd(taux.bytes, zpow[j+1].bytes, gamma[j-1].bytes, taux.bytes);
|
|
|
|
}
|
|
|
|
rct::key mu;
|
|
|
|
sc_muladd(mu.bytes, x.bytes, rho.bytes, alpha.bytes);
|
|
|
|
|
|
|
|
// PAPER LINES 54-57
|
|
|
|
rct::keyV l = l0;
|
|
|
|
l = vector_add(l, vector_scalar(l1, x));
|
|
|
|
rct::keyV r = r0;
|
|
|
|
r = vector_add(r, vector_scalar(r1, x));
|
|
|
|
PERF_TIMER_STOP(PROVE_step2);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_step3);
|
|
|
|
rct::key t = inner_product(l, r);
|
|
|
|
|
|
|
|
// DEBUG: Test if the l and r vectors match the polynomial forms
|
|
|
|
#ifdef DEBUG_BP
|
|
|
|
rct::key test_t;
|
|
|
|
const rct::key t0 = inner_product(l0, r0);
|
|
|
|
sc_muladd(test_t.bytes, t1.bytes, x.bytes, t0.bytes);
|
|
|
|
sc_muladd(test_t.bytes, t2.bytes, xsq.bytes, test_t.bytes);
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(test_t == t, "test_t check failed");
|
|
|
|
#endif
|
|
|
|
|
|
|
|
// PAPER LINES 32-33
|
|
|
|
rct::key x_ip = hash_cache_mash(hash_cache, x, taux, mu, t);
|
|
|
|
|
|
|
|
// These are used in the inner product rounds
|
|
|
|
size_t nprime = MN;
|
|
|
|
rct::keyV Gprime(MN);
|
|
|
|
rct::keyV Hprime(MN);
|
|
|
|
rct::keyV aprime(MN);
|
|
|
|
rct::keyV bprime(MN);
|
|
|
|
const rct::key yinv = invert(y);
|
|
|
|
rct::key yinvpow = rct::identity();
|
|
|
|
for (size_t i = 0; i < MN; ++i)
|
|
|
|
{
|
|
|
|
Gprime[i] = Gi[i];
|
|
|
|
Hprime[i] = scalarmultKey(Hi[i], yinvpow);
|
|
|
|
sc_mul(yinvpow.bytes, yinvpow.bytes, yinv.bytes);
|
|
|
|
aprime[i] = l[i];
|
|
|
|
bprime[i] = r[i];
|
|
|
|
}
|
|
|
|
rct::keyV L(logMN);
|
|
|
|
rct::keyV R(logMN);
|
|
|
|
int round = 0;
|
|
|
|
rct::keyV w(logMN); // this is the challenge x in the inner product protocol
|
|
|
|
PERF_TIMER_STOP(PROVE_step3);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(PROVE_step4);
|
|
|
|
// PAPER LINE 13
|
|
|
|
while (nprime > 1)
|
|
|
|
{
|
|
|
|
// PAPER LINE 15
|
|
|
|
nprime /= 2;
|
|
|
|
|
|
|
|
// PAPER LINES 16-17
|
|
|
|
rct::key cL = inner_product(slice(aprime, 0, nprime), slice(bprime, nprime, bprime.size()));
|
|
|
|
rct::key cR = inner_product(slice(aprime, nprime, aprime.size()), slice(bprime, 0, nprime));
|
|
|
|
|
|
|
|
// PAPER LINES 18-19
|
|
|
|
L[round] = vector_exponent_custom(slice(Gprime, nprime, Gprime.size()), slice(Hprime, 0, nprime), slice(aprime, 0, nprime), slice(bprime, nprime, bprime.size()));
|
|
|
|
sc_mul(tmp.bytes, cL.bytes, x_ip.bytes);
|
|
|
|
rct::addKeys(L[round], L[round], rct::scalarmultKey(rct::H, tmp));
|
|
|
|
R[round] = vector_exponent_custom(slice(Gprime, 0, nprime), slice(Hprime, nprime, Hprime.size()), slice(aprime, nprime, aprime.size()), slice(bprime, 0, nprime));
|
|
|
|
sc_mul(tmp.bytes, cR.bytes, x_ip.bytes);
|
|
|
|
rct::addKeys(R[round], R[round], rct::scalarmultKey(rct::H, tmp));
|
|
|
|
|
|
|
|
// PAPER LINES 21-22
|
|
|
|
w[round] = hash_cache_mash(hash_cache, L[round], R[round]);
|
|
|
|
|
|
|
|
// PAPER LINES 24-25
|
|
|
|
const rct::key winv = invert(w[round]);
|
|
|
|
Gprime = hadamard2(vector_scalar2(slice(Gprime, 0, nprime), winv), vector_scalar2(slice(Gprime, nprime, Gprime.size()), w[round]));
|
|
|
|
Hprime = hadamard2(vector_scalar2(slice(Hprime, 0, nprime), w[round]), vector_scalar2(slice(Hprime, nprime, Hprime.size()), winv));
|
|
|
|
|
|
|
|
// PAPER LINES 28-29
|
|
|
|
aprime = vector_add(vector_scalar(slice(aprime, 0, nprime), w[round]), vector_scalar(slice(aprime, nprime, aprime.size()), winv));
|
|
|
|
bprime = vector_add(vector_scalar(slice(bprime, 0, nprime), winv), vector_scalar(slice(bprime, nprime, bprime.size()), w[round]));
|
|
|
|
|
|
|
|
++round;
|
|
|
|
}
|
|
|
|
PERF_TIMER_STOP(PROVE_step4);
|
|
|
|
|
|
|
|
// PAPER LINE 58 (with inclusions from PAPER LINE 8 and PAPER LINE 20)
|
|
|
|
return Bulletproof(V, A, S, T1, T2, taux, mu, L, R, aprime[0], bprime[0], t);
|
|
|
|
}
|
|
|
|
|
|
|
|
Bulletproof bulletproof_PROVE(const std::vector<uint64_t> &v, const rct::keyV &gamma)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(v.size() == gamma.size(), "Incompatible sizes of v and gamma");
|
|
|
|
|
|
|
|
// vG + gammaH
|
|
|
|
PERF_TIMER_START_BP(PROVE_v);
|
|
|
|
rct::keyV sv(v.size());
|
|
|
|
for (size_t i = 0; i < v.size(); ++i)
|
|
|
|
{
|
|
|
|
sv[i] = rct::zero();
|
|
|
|
sv[i].bytes[0] = v[i] & 255;
|
|
|
|
sv[i].bytes[1] = (v[i] >> 8) & 255;
|
|
|
|
sv[i].bytes[2] = (v[i] >> 16) & 255;
|
|
|
|
sv[i].bytes[3] = (v[i] >> 24) & 255;
|
|
|
|
sv[i].bytes[4] = (v[i] >> 32) & 255;
|
|
|
|
sv[i].bytes[5] = (v[i] >> 40) & 255;
|
|
|
|
sv[i].bytes[6] = (v[i] >> 48) & 255;
|
|
|
|
sv[i].bytes[7] = (v[i] >> 56) & 255;
|
|
|
|
}
|
|
|
|
PERF_TIMER_STOP(PROVE_v);
|
|
|
|
return bulletproof_PROVE(sv, gamma);
|
|
|
|
}
|
|
|
|
|
2017-11-30 14:59:10 -05:00
|
|
|
/* Given a range proof, determine if it is valid */
|
|
|
|
bool bulletproof_VERIFY(const Bulletproof &proof)
|
|
|
|
{
|
|
|
|
init_exponents();
|
|
|
|
|
2018-01-03 16:37:18 -05:00
|
|
|
CHECK_AND_ASSERT_MES(proof.V.size() >= 1, false, "V does not have at least one element");
|
2017-11-30 14:59:10 -05:00
|
|
|
CHECK_AND_ASSERT_MES(proof.L.size() == proof.R.size(), false, "Mismatched L and R sizes");
|
|
|
|
CHECK_AND_ASSERT_MES(proof.L.size() > 0, false, "Empty proof");
|
|
|
|
|
2018-01-03 16:37:18 -05:00
|
|
|
const size_t logN = 6;
|
2017-11-30 14:59:10 -05:00
|
|
|
const size_t N = 1 << logN;
|
2018-01-03 16:37:18 -05:00
|
|
|
rct::key tmp, tmp2;
|
|
|
|
|
|
|
|
size_t M, logM;
|
|
|
|
for (logM = 0; (M = 1<<logM) <= maxM && M < proof.V.size(); ++logM);
|
|
|
|
CHECK_AND_ASSERT_MES(proof.L.size() == 6+logM, false, "Proof is not the expected size");
|
|
|
|
const size_t MN = M*N;
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
// Reconstruct the challenges
|
|
|
|
PERF_TIMER_START_BP(VERIFY);
|
|
|
|
PERF_TIMER_START_BP(VERIFY_start);
|
2018-01-03 16:37:18 -05:00
|
|
|
rct::key hash_cache = rct::hash_to_scalar(proof.V);
|
2017-12-18 16:58:45 -05:00
|
|
|
rct::key y = hash_cache_mash(hash_cache, proof.A, proof.S);
|
|
|
|
rct::key z = hash_cache = rct::hash_to_scalar(y);
|
|
|
|
rct::key x = hash_cache_mash(hash_cache, z, proof.T1, proof.T2);
|
2017-11-30 14:59:10 -05:00
|
|
|
PERF_TIMER_STOP(VERIFY_start);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(VERIFY_line_60);
|
|
|
|
// Reconstruct the challenges
|
2017-12-18 16:58:45 -05:00
|
|
|
rct::key x_ip = hash_cache_mash(hash_cache, x, proof.taux, proof.mu, proof.t);
|
2017-11-30 14:59:10 -05:00
|
|
|
PERF_TIMER_STOP(VERIFY_line_60);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(VERIFY_line_61);
|
|
|
|
// PAPER LINE 61
|
2018-01-05 17:39:59 -05:00
|
|
|
rct::key L61Left;
|
|
|
|
rct::addKeys2(L61Left, proof.taux, proof.t, rct::H);
|
2017-11-30 14:59:10 -05:00
|
|
|
|
2018-01-03 16:37:18 -05:00
|
|
|
const rct::keyV zpow = vector_powers(z, M+3);
|
|
|
|
|
|
|
|
rct::key k;
|
2018-01-05 17:39:59 -05:00
|
|
|
const rct::key ip1y = vector_power_sum(y, MN);
|
2018-01-03 16:37:18 -05:00
|
|
|
sc_mulsub(k.bytes, zpow[2].bytes, ip1y.bytes, rct::zero().bytes);
|
|
|
|
for (size_t j = 1; j <= M; ++j)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_MES(j+2 < zpow.size(), false, "invalid zpow index");
|
|
|
|
sc_mulsub(k.bytes, zpow[j+2].bytes, ip12.bytes, k.bytes);
|
|
|
|
}
|
2017-11-30 14:59:10 -05:00
|
|
|
PERF_TIMER_STOP(VERIFY_line_61);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(VERIFY_line_61rl);
|
|
|
|
sc_muladd(tmp.bytes, z.bytes, ip1y.bytes, k.bytes);
|
2017-12-02 16:17:06 -05:00
|
|
|
rct::key L61Right = rct::scalarmultKey(rct::H, tmp);
|
2018-01-05 17:39:59 -05:00
|
|
|
ge_p3 L61Right_p3;
|
|
|
|
CHECK_AND_ASSERT_THROW_MES(ge_frombytes_vartime(&L61Right_p3, L61Right.bytes) == 0, "ge_frombytes_vartime failed");
|
|
|
|
for (size_t j = 0; j+1 < proof.V.size(); j += 2)
|
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_MES(j+2+1 < zpow.size(), false, "invalid zpow index");
|
|
|
|
ge_dsmp precomp0, precomp1;
|
|
|
|
rct::precomp(precomp0, j < proof.V.size() ? proof.V[j] : rct::identity());
|
|
|
|
rct::precomp(precomp1, j+1 < proof.V.size() ? proof.V[j+1] : rct::identity());
|
|
|
|
rct::addKeys3acc_p3(&L61Right_p3, zpow[j+2], precomp0, zpow[j+2+1], precomp1);
|
|
|
|
}
|
|
|
|
for (size_t j = proof.V.size() & 0xfffffffe; j < M; j++)
|
2018-01-03 16:37:18 -05:00
|
|
|
{
|
|
|
|
CHECK_AND_ASSERT_MES(j+2 < zpow.size(), false, "invalid zpow index");
|
2018-01-05 17:39:59 -05:00
|
|
|
// faster equivalent to:
|
|
|
|
// tmp = rct::scalarmultKey(j < proof.V.size() ? proof.V[j] : rct::identity(), zpow[j+2]);
|
|
|
|
// rct::addKeys(L61Right, L61Right, tmp);
|
|
|
|
if (j < proof.V.size())
|
|
|
|
addKeys_acc_p3(&L61Right_p3, zpow[j+2], proof.V[j]);
|
2018-01-03 16:37:18 -05:00
|
|
|
}
|
2017-11-30 14:59:10 -05:00
|
|
|
|
2018-01-05 17:39:59 -05:00
|
|
|
addKeys_acc_p3(&L61Right_p3, x, proof.T1);
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
rct::key xsq;
|
|
|
|
sc_mul(xsq.bytes, x.bytes, x.bytes);
|
2018-01-05 17:39:59 -05:00
|
|
|
addKeys_acc_p3(&L61Right_p3, xsq, proof.T2);
|
|
|
|
ge_p3_tobytes(L61Right.bytes, &L61Right_p3);
|
2017-11-30 14:59:10 -05:00
|
|
|
PERF_TIMER_STOP(VERIFY_line_61rl);
|
|
|
|
|
|
|
|
if (!(L61Right == L61Left))
|
|
|
|
{
|
|
|
|
MERROR("Verification failure at step 1");
|
|
|
|
return false;
|
|
|
|
}
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(VERIFY_line_62);
|
|
|
|
// PAPER LINE 62
|
|
|
|
rct::key P = rct::addKeys(proof.A, rct::scalarmultKey(proof.S, x));
|
|
|
|
PERF_TIMER_STOP(VERIFY_line_62);
|
|
|
|
|
|
|
|
// Compute the number of rounds for the inner product
|
2018-01-03 16:37:18 -05:00
|
|
|
const size_t rounds = logM+logN;
|
2017-11-30 14:59:10 -05:00
|
|
|
CHECK_AND_ASSERT_MES(rounds > 0, false, "Zero rounds");
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(VERIFY_line_21_22);
|
|
|
|
// PAPER LINES 21-22
|
|
|
|
// The inner product challenges are computed per round
|
|
|
|
rct::keyV w(rounds);
|
2017-12-18 16:58:45 -05:00
|
|
|
for (size_t i = 0; i < rounds; ++i)
|
2017-11-30 14:59:10 -05:00
|
|
|
{
|
2017-12-18 16:58:45 -05:00
|
|
|
w[i] = hash_cache_mash(hash_cache, proof.L[i], proof.R[i]);
|
2017-11-30 14:59:10 -05:00
|
|
|
}
|
|
|
|
PERF_TIMER_STOP(VERIFY_line_21_22);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(VERIFY_line_24_25);
|
|
|
|
// Basically PAPER LINES 24-25
|
|
|
|
// Compute the curvepoints from G[i] and H[i]
|
2018-01-05 17:39:59 -05:00
|
|
|
ge_p3 inner_prod_p3 = ge_p3_identity;
|
2017-11-30 14:59:10 -05:00
|
|
|
rct::key yinvpow = rct::identity();
|
|
|
|
rct::key ypow = rct::identity();
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(VERIFY_line_24_25_invert);
|
|
|
|
const rct::key yinv = invert(y);
|
|
|
|
rct::keyV winv(rounds);
|
|
|
|
for (size_t i = 0; i < rounds; ++i)
|
|
|
|
winv[i] = invert(w[i]);
|
|
|
|
PERF_TIMER_STOP(VERIFY_line_24_25_invert);
|
|
|
|
|
2018-01-03 16:37:18 -05:00
|
|
|
for (size_t i = 0; i < MN; ++i)
|
2017-11-30 14:59:10 -05:00
|
|
|
{
|
|
|
|
// Convert the index to binary IN REVERSE and construct the scalar exponent
|
|
|
|
rct::key g_scalar = proof.a;
|
|
|
|
rct::key h_scalar;
|
|
|
|
sc_mul(h_scalar.bytes, proof.b.bytes, yinvpow.bytes);
|
|
|
|
|
|
|
|
for (size_t j = rounds; j-- > 0; )
|
|
|
|
{
|
|
|
|
size_t J = w.size() - j - 1;
|
|
|
|
|
|
|
|
if ((i & (((size_t)1)<<j)) == 0)
|
|
|
|
{
|
|
|
|
sc_mul(g_scalar.bytes, g_scalar.bytes, winv[J].bytes);
|
|
|
|
sc_mul(h_scalar.bytes, h_scalar.bytes, w[J].bytes);
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
sc_mul(g_scalar.bytes, g_scalar.bytes, w[J].bytes);
|
|
|
|
sc_mul(h_scalar.bytes, h_scalar.bytes, winv[J].bytes);
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
// Adjust the scalars using the exponents from PAPER LINE 62
|
|
|
|
sc_add(g_scalar.bytes, g_scalar.bytes, z.bytes);
|
2018-01-03 16:37:18 -05:00
|
|
|
CHECK_AND_ASSERT_MES(2+i/N < zpow.size(), false, "invalid zpow index");
|
|
|
|
CHECK_AND_ASSERT_MES(i%N < twoN.size(), false, "invalid twoN index");
|
|
|
|
sc_mul(tmp.bytes, zpow[2+i/N].bytes, twoN[i%N].bytes);
|
2017-11-30 14:59:10 -05:00
|
|
|
sc_muladd(tmp.bytes, z.bytes, ypow.bytes, tmp.bytes);
|
|
|
|
sc_mulsub(h_scalar.bytes, tmp.bytes, yinvpow.bytes, h_scalar.bytes);
|
|
|
|
|
|
|
|
// Now compute the basepoint's scalar multiplication
|
|
|
|
// Each of these could be written as a multiexp operation instead
|
2018-01-05 17:39:59 -05:00
|
|
|
addKeys3acc_p3(&inner_prod_p3, g_scalar, Gprecomp[i], h_scalar, Hprecomp[i]);
|
2017-11-30 14:59:10 -05:00
|
|
|
|
2018-01-03 16:37:18 -05:00
|
|
|
if (i != MN-1)
|
2017-11-30 14:59:10 -05:00
|
|
|
{
|
|
|
|
sc_mul(yinvpow.bytes, yinvpow.bytes, yinv.bytes);
|
|
|
|
sc_mul(ypow.bytes, ypow.bytes, y.bytes);
|
|
|
|
}
|
|
|
|
}
|
2018-01-05 17:39:59 -05:00
|
|
|
rct::key inner_prod;
|
|
|
|
ge_p3_tobytes(inner_prod.bytes, &inner_prod_p3);
|
2017-11-30 14:59:10 -05:00
|
|
|
PERF_TIMER_STOP(VERIFY_line_24_25);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(VERIFY_line_26);
|
|
|
|
// PAPER LINE 26
|
|
|
|
rct::key pprime;
|
|
|
|
sc_sub(tmp.bytes, rct::zero().bytes, proof.mu.bytes);
|
2017-12-02 16:17:06 -05:00
|
|
|
rct::addKeys(pprime, P, rct::scalarmultBase(tmp));
|
2018-01-05 17:39:59 -05:00
|
|
|
ge_p3 pprime_p3;
|
|
|
|
CHECK_AND_ASSERT_MES(ge_frombytes_vartime(&pprime_p3, pprime.bytes) == 0, false, "ge_frombytes_vartime failed");
|
2017-11-30 14:59:10 -05:00
|
|
|
|
|
|
|
for (size_t i = 0; i < rounds; ++i)
|
|
|
|
{
|
|
|
|
sc_mul(tmp.bytes, w[i].bytes, w[i].bytes);
|
|
|
|
sc_mul(tmp2.bytes, winv[i].bytes, winv[i].bytes);
|
|
|
|
#if 1
|
|
|
|
ge_dsmp cacheL, cacheR;
|
|
|
|
rct::precomp(cacheL, proof.L[i]);
|
|
|
|
rct::precomp(cacheR, proof.R[i]);
|
2018-01-05 17:39:59 -05:00
|
|
|
addKeys3acc_p3(&pprime_p3, tmp, cacheL, tmp2, cacheR);
|
2017-11-30 14:59:10 -05:00
|
|
|
#else
|
|
|
|
rct::addKeys(pprime, pprime, rct::scalarmultKey(proof.L[i], tmp));
|
|
|
|
rct::addKeys(pprime, pprime, rct::scalarmultKey(proof.R[i], tmp2));
|
|
|
|
#endif
|
|
|
|
}
|
|
|
|
sc_mul(tmp.bytes, proof.t.bytes, x_ip.bytes);
|
2018-01-05 17:39:59 -05:00
|
|
|
addKeys_acc_p3(&pprime_p3, tmp, rct::H);
|
|
|
|
ge_p3_tobytes(pprime.bytes, &pprime_p3);
|
2017-11-30 14:59:10 -05:00
|
|
|
PERF_TIMER_STOP(VERIFY_line_26);
|
|
|
|
|
|
|
|
PERF_TIMER_START_BP(VERIFY_step2_check);
|
|
|
|
sc_mul(tmp.bytes, proof.a.bytes, proof.b.bytes);
|
|
|
|
sc_mul(tmp.bytes, tmp.bytes, x_ip.bytes);
|
2017-12-02 16:17:06 -05:00
|
|
|
tmp = rct::scalarmultKey(rct::H, tmp);
|
2017-11-30 14:59:10 -05:00
|
|
|
rct::addKeys(tmp, tmp, inner_prod);
|
|
|
|
PERF_TIMER_STOP(VERIFY_step2_check);
|
|
|
|
if (!(pprime == tmp))
|
|
|
|
{
|
|
|
|
MERROR("Verification failure at step 2");
|
|
|
|
return false;
|
|
|
|
}
|
|
|
|
|
|
|
|
PERF_TIMER_STOP(VERIFY);
|
|
|
|
return true;
|
|
|
|
}
|
|
|
|
|
|
|
|
}
|