fk20test_fft_rand.cu

// bls12_381: Arithmetic for BLS12-381
// Copyright 2022-2023 Dag Arne Osvik
// Copyright 2022-2023 Luan Cardoso dos Santos

#warning "Functionality of this file has been replaced by fk20_512."

#include <stdio.h>
#include <unistd.h>

#include "g1.cuh"
#include "fk20test.cuh"
#include "fk20.cuh"


//static __managed__ uint8_t cmp[16*512];
//static __managed__ fr_t fr_tmp[16*512];
//static __managed__ g1p_t g1p_tmp[512];

//Global correct answers, read from the instrumented fk20 python implementation
static  __managed__ size_t elementsRead;
static  __managed__ fr_t  polynomial[4096];
static  __managed__ g1p_t setup[4097];
static  __managed__ g1p_t xext_fft[16][512];
static  __managed__ fr_t  toeplitz_coefficients[16][512];
static  __managed__ fr_t  toeplitz_coefficients_fft[16][512];
static  __managed__ g1p_t hext_fft[512];
static  __managed__ g1p_t h[512];
static  __managed__ g1p_t h_fft[512];

#define readAndCheck(var) do{                                                  \
        elementsRead = fread(var, 1, sizeof(var), inputStream);                \
        if(elementsRead != sizeof(var)){                                       \
            printf("NONFATAL: Unbale to read test vectors, aborting\n");         \
            return;                                                          \
        } else {                                                               \
            printf("    Read in " #var " %ld Bytes\n", elementsRead);          \
        }                                                                      \
    }while(0)

void FK20TestFFTRand(FILE *inputStream) {
    /*
     * Reads from inputStream values generated by the python implementation, and test against those.
     * if the inputstream is NULL, then read from stdin
     */

    if (inputStream == NULL){
        inputStream = stdin;
    }

    if (isatty(fileno(inputStream)))
        return;

    printf(">>>> Rand Tests\n");

    readAndCheck(polynomial);
    readAndCheck(setup);
    readAndCheck(xext_fft);
    readAndCheck(toeplitz_coefficients);
    readAndCheck(toeplitz_coefficients_fft);
    readAndCheck(hext_fft);
    readAndCheck(h);
    readAndCheck(h_fft);

    toeplitz_coefficients2toeplitz_coefficients_fft(toeplitz_coefficients, toeplitz_coefficients_fft);
    h2h_fft(h, h_fft);
    h_fft2h(h_fft, h);
    hext_fft2h(hext_fft, h);
    hext_fft2h_fft(hext_fft, h_fft);

    fk20_poly2toeplitz_coefficients_test(polynomial, toeplitz_coefficients);
    fk20_poly2hext_fft_test(polynomial, xext_fft, hext_fft);
    fk20_msmloop(hext_fft, toeplitz_coefficients_fft, xext_fft);
    fk20_poly2h_fft_test(polynomial, xext_fft, h_fft);


    //Null the pointer if we used STDIN
    if (inputStream == stdin){
        inputStream = NULL;
    }
}

// vim: ts=4 et sw=4 si