path: root/src/bls_c_impl.hpp

#define MCLBN_DONT_EXPORT
#define BLS_DLL_EXPORT

#include <bls/bls.h>

#if 1
#include "mcl/impl/bn_c_impl.hpp"
#else
#if MCLBN_FP_UNIT_SIZE == 4 && MCLBN_FR_UNIT_SIZE == 4
#include <mcl/bn256.hpp>
#elif MCLBN_FP_UNIT_SIZE == 6 && MCLBN_FR_UNIT_SIZE == 6
#include <mcl/bn384.hpp>
#elif MCLBN_FP_UNIT_SIZE == 6 && MCLBN_FR_UNIT_SIZE == 4
#include <mcl/bls12_381.hpp>
#elif MCLBN_FP_UNIT_SIZE == 8 && MCLBN_FR_UNIT_SIZE == 8
#include <mcl/bn512.hpp>
#else
    #error "not supported size"
#endif
#include <mcl/lagrange.hpp>
using namespace mcl::bn;
inline Fr *cast(mclBnFr *p) { return reinterpret_cast<Fr*>(p); }
inline const Fr *cast(const mclBnFr *p) { return reinterpret_cast<const Fr*>(p); }

inline G1 *cast(mclBnG1 *p) { return reinterpret_cast<G1*>(p); }
inline const G1 *cast(const mclBnG1 *p) { return reinterpret_cast<const G1*>(p); }

inline G2 *cast(mclBnG2 *p) { return reinterpret_cast<G2*>(p); }
inline const G2 *cast(const mclBnG2 *p) { return reinterpret_cast<const G2*>(p); }

inline Fp12 *cast(mclBnGT *p) { return reinterpret_cast<Fp12*>(p); }
inline const Fp12 *cast(const mclBnGT *p) { return reinterpret_cast<const Fp12*>(p); }

inline Fp6 *cast(uint64_t *p) { return reinterpret_cast<Fp6*>(p); }
inline const Fp6 *cast(const uint64_t *p) { return reinterpret_cast<const Fp6*>(p); }
#endif

void Gmul(G1& z, const G1& x, const Fr& y) { G1::mul(z, x, y); }
void Gmul(G2& z, const G2& x, const Fr& y) { G2::mul(z, x, y); }
void GmulCT(G1& z, const G1& x, const Fr& y) { G1::mulCT(z, x, y); }
void GmulCT(G2& z, const G2& x, const Fr& y) { G2::mulCT(z, x, y); }

/*
    BLS signature
    e : G1 x G2 -> GT
    Q in G2 ; fixed global parameter
    H : {str} -> G1
    s : secret key
    sQ ; public key
    s H(m) ; signature of m
    verify ; e(sQ, H(m)) = e(Q, s H(m))

    swap G1 and G2 if BLS_SWAP_G is defined
    @note the current implementation does not support precomputed miller loop
*/

#ifdef BLS_SWAP_G
static G1 g_P;
inline const G1& getBasePoint() { return g_P; }
#else
static G2 g_Q;
const size_t maxQcoeffN = 128;
static mcl::FixedArray<Fp6, maxQcoeffN> g_Qcoeff; // precomputed Q
inline const G2& getBasePoint() { return g_Q; }
inline const mcl::FixedArray<Fp6, maxQcoeffN>& getQcoeff() { return g_Qcoeff; }
#endif

int blsInitNotThreadSafe(int curve, int compiledTimeVar)
{
    if (compiledTimeVar != MCLBN_COMPILED_TIME_VAR) {
        return -(compiledTimeVar | (MCLBN_COMPILED_TIME_VAR * 100));
    }
    const mcl::CurveParam& cp = mcl::getCurveParam(curve);
    bool b;
    initPairing(&b, cp);
    if (!b) return -1;

#ifdef BLS_SWAP_G
    mapToG1(&b, g_P, 1);
#else

    if (curve == MCL_BN254) {
        const char *Qx_BN254 = "11ccb44e77ac2c5dc32a6009594dbe331ec85a61290d6bbac8cc7ebb2dceb128 f204a14bbdac4a05be9a25176de827f2e60085668becdd4fc5fa914c9ee0d9a";
        const char *Qy_BN254 = "7c13d8487903ee3c1c5ea327a3a52b6cc74796b1760d5ba20ed802624ed19c8 8f9642bbaacb73d8c89492528f58932f2de9ac3e80c7b0e41f1a84f1c40182";
        g_Q.x.setStr(&b, Qx_BN254, 16);
        g_Q.y.setStr(&b, Qy_BN254, 16);
        g_Q.z = 1;
    } else {
        mapToG2(&b, g_Q, 1);
    }
    if (!b) return -100;
    if (curve == MCL_BN254) {
        #include "./qcoeff-bn254.hpp"
        g_Qcoeff.resize(BN::param.precomputedQcoeffSize);
        assert(g_Qcoeff.size() == CYBOZU_NUM_OF_ARRAY(QcoeffTblBN254));
        for (size_t i = 0; i < g_Qcoeff.size(); i++) {
            Fp6& x6 = g_Qcoeff[i];
            for (size_t j = 0; j < 6; j++) {
                Fp& x = x6.getFp0()[j];
                mcl::fp::Unit *p = const_cast<mcl::fp::Unit*>(x.getUnit());
                for (size_t k = 0; k < 4; k++) {
                    p[k] = QcoeffTblBN254[i][j][k];
                }
            }
        }
    } else {
        precomputeG2(&b, g_Qcoeff, getBasePoint());
    }
#endif
    if (!b) return -101;
    return 0;
}

#ifdef __EMSCRIPTEN__
extern "C" BLS_DLL_API void *blsMalloc(size_t n)
{
    return malloc(n);
}
extern "C" BLS_DLL_API void blsFree(void *p)
{
    free(p);
}
#endif

#if !defined(__EMSCRIPTEN__) && !defined(__wasm__)
    #if defined(CYBOZU_CPP_VERSION) && CYBOZU_CPP_VERSION >= CYBOZU_CPP_VERSION_CPP11
    #include <mutex>
        #define USE_STD_MUTEX
    #else
    #include <cybozu/mutex.hpp>
        #define USE_CYBOZU_MUTEX
    #endif
#endif

int blsInit(int curve, int compiledTimeVar)
{
    int ret = 0;
#ifdef USE_STD_MUTEX
    static std::mutex m;
    std::lock_guard<std::mutex> lock(m);
#elif defined(USE_CYBOZU_MUTEX)
    static cybozu::Mutex m;
    cybozu::AutoLock lock(m);
#endif
    static int g_curve = -1;
    if (g_curve != curve) {
        ret = blsInitNotThreadSafe(curve, compiledTimeVar);
        g_curve = curve;
    }
    return ret;
}

static inline const mclBnG1 *cast(const G1* x) { return (const mclBnG1*)x; }
static inline const mclBnG2 *cast(const G2* x) { return (const mclBnG2*)x; }

void blsIdSetInt(blsId *id, int x)
{
    *cast(&id->v) = x;
}

int blsSecretKeySetLittleEndian(blsSecretKey *sec, const void *buf, mclSize bufSize)
{
    cast(&sec->v)->setArrayMask((const char *)buf, bufSize);
    return 0;
}
int blsSecretKeySetLittleEndianMod(blsSecretKey *sec, const void *buf, mclSize bufSize)
{
    bool b;
    cast(&sec->v)->setArray(&b, (const char *)buf, bufSize, mcl::fp::Mod);
    return b ? 0 : -1;
}

void blsGetPublicKey(blsPublicKey *pub, const blsSecretKey *sec)
{
    Gmul(*cast(&pub->v), getBasePoint(), *cast(&sec->v));
}

void blsSign(blsSignature *sig, const blsSecretKey *sec, const void *m, mclSize size)
{
#ifdef BLS_SWAP_G
    G2 Hm;
    hashAndMapToG2(Hm, m, size);
#else
    G1 Hm;
    hashAndMapToG1(Hm, m, size);
#endif
    GmulCT(*cast(&sig->v), Hm, *cast(&sec->v));
}

#ifdef BLS_SWAP_G
/*
    e(P, sHm) == e(sP, Hm)
    <=> finalExp(ML(P, sHm) * e(-sP, Hm)) == 1
*/
bool isEqualTwoPairings(const G2& sHm, const G1& sP, const G2& Hm)
{
    GT e1, e2;
    millerLoop(e1, getBasePoint(), sHm);
    G1 neg_sP;
    G1::neg(neg_sP, sP);
    millerLoop(e2, neg_sP, Hm);
    e1 *= e2;
    finalExp(e1, e1);
    return e1.isOne();
}
#else
/*
    e(P1, Q1) == e(P2, Q2)
    <=> finalExp(ML(P1, Q1)) == finalExp(ML(P2, Q2))
    <=> finalExp(ML(P1, Q1) / ML(P2, Q2)) == 1
    <=> finalExp(ML(P1, Q1) * ML(-P2, Q2)) == 1
    Q1 is precomputed
*/
bool isEqualTwoPairings(const G1& P1, const Fp6* Q1coeff, const G1& P2, const G2& Q2)
{
    GT e;
    precomputedMillerLoop2mixed(e, P2, Q2, -P1, Q1coeff);
    finalExp(e, e);
    return e.isOne();
}
#endif

int blsVerify(const blsSignature *sig, const blsPublicKey *pub, const void *m, mclSize size)
{
#ifdef BLS_SWAP_G
    G2 Hm;
    hashAndMapToG2(Hm, m, size);
    return isEqualTwoPairings(*cast(&sig->v), *cast(&pub->v), Hm);
#else
    G1 Hm;
    hashAndMapToG1(Hm, m, size);
    /*
        e(sHm, Q) = e(Hm, sQ)
        e(sig, Q) = e(Hm, pub)
    */
    return isEqualTwoPairings(*cast(&sig->v), getQcoeff().data(), Hm, *cast(&pub->v));
#endif
}

mclSize blsIdSerialize(void *buf, mclSize maxBufSize, const blsId *id)
{
    return cast(&id->v)->serialize(buf, maxBufSize);
}

mclSize blsSecretKeySerialize(void *buf, mclSize maxBufSize, const blsSecretKey *sec)
{
    return cast(&sec->v)->serialize(buf, maxBufSize);
}

mclSize blsPublicKeySerialize(void *buf, mclSize maxBufSize, const blsPublicKey *pub)
{
    return cast(&pub->v)->serialize(buf, maxBufSize);
}

mclSize blsSignatureSerialize(void *buf, mclSize maxBufSize, const blsSignature *sig)
{
    return cast(&sig->v)->serialize(buf, maxBufSize);
}

mclSize blsIdDeserialize(blsId *id, const void *buf, mclSize bufSize)
{
    return cast(&id->v)->deserialize(buf, bufSize);
}

mclSize blsSecretKeyDeserialize(blsSecretKey *sec, const void *buf, mclSize bufSize)
{
    return cast(&sec->v)->deserialize(buf, bufSize);
}

mclSize blsPublicKeyDeserialize(blsPublicKey *pub, const void *buf, mclSize bufSize)
{
    return cast(&pub->v)->deserialize(buf, bufSize);
}

mclSize blsSignatureDeserialize(blsSignature *sig, const void *buf, mclSize bufSize)
{
    return cast(&sig->v)->deserialize(buf, bufSize);
}

int blsIdIsEqual(const blsId *lhs, const blsId *rhs)
{
    return *cast(&lhs->v) == *cast(&rhs->v);
}

int blsSecretKeyIsEqual(const blsSecretKey *lhs, const blsSecretKey *rhs)
{
    return *cast(&lhs->v) == *cast(&rhs->v);
}

int blsPublicKeyIsEqual(const blsPublicKey *lhs, const blsPublicKey *rhs)
{
    return *cast(&lhs->v) == *cast(&rhs->v);
}

int blsSignatureIsEqual(const blsSignature *lhs, const blsSignature *rhs)
{
    return *cast(&lhs->v) == *cast(&rhs->v);
}

int blsSecretKeyShare(blsSecretKey *sec, const blsSecretKey* msk, mclSize k, const blsId *id)
{
    bool b;
    mcl::evaluatePolynomial(&b, *cast(&sec->v), cast(&msk->v), k, *cast(&id->v));
    return b ? 0 : -1;
}

int blsPublicKeyShare(blsPublicKey *pub, const blsPublicKey *mpk, mclSize k, const blsId *id)
{
    bool b;
    mcl::evaluatePolynomial(&b, *cast(&pub->v), cast(&mpk->v), k, *cast(&id->v));
    return b ? 0 : -1;
}

int blsSecretKeyRecover(blsSecretKey *sec, const blsSecretKey *secVec, const blsId *idVec, mclSize n)
{
    bool b;
    mcl::LagrangeInterpolation(&b, *cast(&sec->v), cast(&idVec->v), cast(&secVec->v), n);
    return b ? 0 : -1;
}

int blsPublicKeyRecover(blsPublicKey *pub, const blsPublicKey *pubVec, const blsId *idVec, mclSize n)
{
    bool b;
    mcl::LagrangeInterpolation(&b, *cast(&pub->v), cast(&idVec->v), cast(&pubVec->v), n);
    return b ? 0 : -1;
}

int blsSignatureRecover(blsSignature *sig, const blsSignature *sigVec, const blsId *idVec, mclSize n)
{
    bool b;
    mcl::LagrangeInterpolation(&b, *cast(&sig->v), cast(&idVec->v), cast(&sigVec->v), n);
    return b ? 0 : -1;
}

void blsSecretKeyAdd(blsSecretKey *sec, const blsSecretKey *rhs)
{
    *cast(&sec->v) += *cast(&rhs->v);
}

void blsPublicKeyAdd(blsPublicKey *pub, const blsPublicKey *rhs)
{
    *cast(&pub->v) += *cast(&rhs->v);
}

void blsSignatureAdd(blsSignature *sig, const blsSignature *rhs)
{
    *cast(&sig->v) += *cast(&rhs->v);
}

void blsSignatureVerifyOrder(int doVerify)
{
#ifdef BLS_SWAP_G
    verifyOrderG2(doVerify != 0);
#else
    verifyOrderG1(doVerify != 0);
#endif
}
void blsPublicKeyVerifyOrder(int doVerify)
{
#ifdef BLS_SWAP_G
    verifyOrderG1(doVerify != 0);
#else
    verifyOrderG2(doVerify != 0);
#endif
}
int blsSignatureIsValidOrder(const blsSignature *sig)
{
    return cast(&sig->v)->isValidOrder();
}
int blsPublicKeyIsValidOrder(const blsPublicKey *pub)
{
    return cast(&pub->v)->isValidOrder();
}

#ifndef BLS_MINIMUM_API
template<class G>
inline bool toG(G& Hm, const void *h, mclSize size)
{
    Fp t;
    t.setArrayMask((const char *)h, size);
    bool b;
#ifdef BLS_SWAP_G
    BN::mapToG2(&b, Hm, Fp2(t, 0));
#else
    BN::mapToG1(&b, Hm, t);
#endif
    return b;
}

int blsVerifyAggregatedHashes(const blsSignature *aggSig, const blsPublicKey *pubVec, const void *hVec, size_t sizeofHash, mclSize n)
{
    if (n == 0) return 0;
    GT e1, e2;
    const char *ph = (const char*)hVec;
#ifdef BLS_SWAP_G
    millerLoop(e1, getBasePoint(), -*cast(&aggSig->v));
    G2 h;
    if (!toG(h, &ph[0], sizeofHash)) return 0;
    BN::millerLoop(e2, *cast(&pubVec[0].v), h);
    e1 *= e2;
    for (size_t i = 1; i < n; i++) {
        if (!toG(h, &ph[i * sizeofHash], sizeofHash)) return 0;
        millerLoop(e2, *cast(&pubVec[i].v), h);
        e1 *= e2;
    }
#else
    /*
        e(aggSig, Q) = prod_i e(hVec[i], pubVec[i])
        <=> finalExp(ML(-aggSig, Q) * prod_i ML(hVec[i], pubVec[i])) == 1
    */
    BN::precomputedMillerLoop(e1, -*cast(&aggSig->v), g_Qcoeff.data());
    G1 h;
    if (!toG(h, &ph[0], sizeofHash)) return 0;
    BN::millerLoop(e2, h, *cast(&pubVec[0].v));
    e1 *= e2;
    for (size_t i = 1; i < n; i++) {
        if (!toG(h, &ph[i * sizeofHash], sizeofHash)) return 0;
        BN::millerLoop(e2, h, *cast(&pubVec[i].v));
        e1 *= e2;
    }
#endif
    BN::finalExp(e1, e1);
    return e1.isOne();
}

int blsSignHash(blsSignature *sig, const blsSecretKey *sec, const void *h, mclSize size)
{
#ifdef BLS_SWAP_G
    G2 Hm;
#else
    G1 Hm;
#endif
    if (!toG(Hm, h, size)) return -1;
    GmulCT(*cast(&sig->v), Hm, *cast(&sec->v));
    return 0;
}

int blsVerifyPairing(const blsSignature *X, const blsSignature *Y, const blsPublicKey *pub)
{
#ifdef BLS_SWAP_G
    return isEqualTwoPairings(*cast(&X->v), *cast(&pub->v), *cast(&Y->v));
#else
    return isEqualTwoPairings(*cast(&X->v), getQcoeff().data(), *cast(&Y->v), *cast(&pub->v));
#endif
}

int blsVerifyHash(const blsSignature *sig, const blsPublicKey *pub, const void *h, mclSize size)
{
    blsSignature Hm;
    if (!toG(*cast(&Hm.v), h, size)) return 0;
    return blsVerifyPairing(sig, &Hm, pub);
}

void blsSecretKeySub(blsSecretKey *sec, const blsSecretKey *rhs)
{
    *cast(&sec->v) -= *cast(&rhs->v);
}

void blsPublicKeySub(blsPublicKey *pub, const blsPublicKey *rhs)
{
    *cast(&pub->v) -= *cast(&rhs->v);
}

void blsSignatureSub(blsSignature *sig, const blsSignature *rhs)
{
    *cast(&sig->v) -= *cast(&rhs->v);
}

mclSize blsGetOpUnitSize() // FpUint64Size
{
    return Fp::getUnitSize() * sizeof(mcl::fp::Unit) / sizeof(uint64_t);
}

int blsGetCurveOrder(char *buf, mclSize maxBufSize)
{
    return (int)Fr::getModulo(buf, maxBufSize);
}

int blsGetFieldOrder(char *buf, mclSize maxBufSize)
{
    return (int)Fp::getModulo(buf, maxBufSize);
}

int blsGetG1ByteSize()
{
    return (int)Fp::getByteSize();
}

int blsGetFrByteSize()
{
    return (int)Fr::getByteSize();
}

#ifdef BLS_SWAP_G
void blsGetGeneratorOfG1(blsPublicKey *pub)
{
    *cast(&pub->v) = getBasePoint();
}
#else
void blsGetGeneratorOfG2(blsPublicKey *pub)
{
    *cast(&pub->v) = getBasePoint();
}
#endif

int blsIdSetDecStr(blsId *id, const char *buf, mclSize bufSize)
{
    return cast(&id->v)->deserialize(buf, bufSize, 10) > 0 ? 0 : -1;
}
int blsIdSetHexStr(blsId *id, const char *buf, mclSize bufSize)
{
    return cast(&id->v)->deserialize(buf, bufSize, 16) > 0 ? 0 : -1;
}

int blsIdSetLittleEndian(blsId *id, const void *buf, mclSize bufSize)
{
    cast(&id->v)->setArrayMask((const char *)buf, bufSize);
    return 0;
}

mclSize blsIdGetDecStr(char *buf, mclSize maxBufSize, const blsId *id)
{
    return cast(&id->v)->getStr(buf, maxBufSize, 10);
}

mclSize blsIdGetHexStr(char *buf, mclSize maxBufSize, const blsId *id)
{
    return cast(&id->v)->getStr(buf, maxBufSize, 16);
}

int blsHashToSecretKey(blsSecretKey *sec, const void *buf, mclSize bufSize)
{
    cast(&sec->v)->setHashOf(buf, bufSize);
    return 0;
}

#ifndef MCL_DONT_USE_CSPRNG
int blsSecretKeySetByCSPRNG(blsSecretKey *sec)
{
    bool b;
    cast(&sec->v)->setByCSPRNG(&b);
    return b ? 0 : -1;
}
void blsSetRandFunc(void *self, unsigned int (*readFunc)(void *self, void *buf, unsigned int bufSize))
{
    mcl::fp::RandGen::setRandFunc(self, readFunc);
}
#endif

void blsGetPop(blsSignature *sig, const blsSecretKey *sec)
{
    blsPublicKey pub;
    blsGetPublicKey(&pub, sec);
    char buf[1024];
    mclSize n = cast(&pub.v)->serialize(buf, sizeof(buf));
    assert(n);
    blsSign(sig, sec, buf, n);
}

int blsVerifyPop(const blsSignature *sig, const blsPublicKey *pub)
{
    char buf[1024];
    mclSize n = cast(&pub->v)->serialize(buf, sizeof(buf));
    if (n == 0) return 0;
    return blsVerify(sig, pub, buf, n);
}

mclSize blsIdGetLittleEndian(void *buf, mclSize maxBufSize, const blsId *id)
{
    return cast(&id->v)->serialize(buf, maxBufSize);
}
int blsSecretKeySetDecStr(blsSecretKey *sec, const char *buf, mclSize bufSize)
{
    return cast(&sec->v)->deserialize(buf, bufSize, 10) > 0 ? 0 : -1;
}
int blsSecretKeySetHexStr(blsSecretKey *sec, const char *buf, mclSize bufSize)
{
    return cast(&sec->v)->deserialize(buf, bufSize, 16) > 0 ? 0 : -1;
}
mclSize blsSecretKeyGetLittleEndian(void *buf, mclSize maxBufSize, const blsSecretKey *sec)
{
    return cast(&sec->v)->serialize(buf, maxBufSize);
}
mclSize blsSecretKeyGetDecStr(char *buf, mclSize maxBufSize, const blsSecretKey *sec)
{
    return cast(&sec->v)->getStr(buf, maxBufSize, 10);
}
mclSize blsSecretKeyGetHexStr(char *buf, mclSize maxBufSize, const blsSecretKey *sec)
{
    return cast(&sec->v)->getStr(buf, maxBufSize, 16);
}
int blsPublicKeySetHexStr(blsPublicKey *pub, const char *buf, mclSize bufSize)
{
    return cast(&pub->v)->deserialize(buf, bufSize, 16) > 0 ? 0 : -1;
}
mclSize blsPublicKeyGetHexStr(char *buf, mclSize maxBufSize, const blsPublicKey *pub)
{
    return cast(&pub->v)->getStr(buf, maxBufSize, 16);
}
int blsSignatureSetHexStr(blsSignature *sig, const char *buf, mclSize bufSize)
{
    return cast(&sig->v)->deserialize(buf, bufSize, 16) > 0 ? 0 : -1;
}
mclSize blsSignatureGetHexStr(char *buf, mclSize maxBufSize, const blsSignature *sig)
{
    return cast(&sig->v)->getStr(buf, maxBufSize, 16);
}
void blsDHKeyExchange(blsPublicKey *out, const blsSecretKey *sec, const blsPublicKey *pub)
{
    GmulCT(*cast(&out->v), *cast(&pub->v), *cast(&sec->v));
}

#endif