g1p_compare.cu

Go to the documentation of this file.

// bls12_381: Arithmetic for BLS12-381
// Copyright 2022-2023 Dag Arne Osvik
// Copyright 2022-2023 Luan Cardoso dos Santos
 
#include <stdio.h>
 
#include "fp.cuh"
#include "fr.cuh"
#include "g1.cuh"
 
__device__ bool g1p_eq(const g1p_t &p, const g1p_t &q) {
    fp_t px, py, qx, qy;
 
#ifndef NDEBUG
    if (!g1p_isPoint(p) || !(g1p_isPoint(q))) {
       // printf("ERROR in g1p_eq(): Invalid point(s)\n");
       // g1p_print("p: ", p);
       // g1p_print("q: ", q);
 
        return false;
    }
#endif
 
    // (X1/Z1 == X2/Z2) && (Y1/Z1 == Y2/Z2)
 
    fp_cpy(px, p.x);
    fp_cpy(py, p.y);
 
    fp_cpy(qx, q.x);
    fp_cpy(qy, q.y);
 
    fp_mul(px, px, q.z);    // X1*Z2
    fp_mul(qx, qx, p.z);    // X2*Z1
 
    if (fp_neq(px, qx))
        return false;
 
    fp_mul(py, py, q.z);    // Y1*Z2
    fp_mul(qy, qy, p.z);    // Y2*Z1
 
    return fp_eq(py, qy);
}
 
__device__ bool g1p_neq(const g1p_t &p, const g1p_t &q) {
// observation: This function could be inlined, but
#ifndef NDEBUG
    if (!g1p_isPoint(p) || !(g1p_isPoint(q))) {
        printf("ERROR in g1p_neq(): Invalid point(s)\n");
        g1p_print("p: ", p);
        g1p_print("q: ", q);
 
        return true;
    }
#endif
 
    return !g1p_eq(p, q);
}
 
// vim: ts=4 et sw=4 si