/*
    Copyright (C) 2012 Andres Goens
    Copyright (C) 2013 Mike Hansen

    This file is part of FLINT.

    FLINT is free software: you can redistribute it and/or modify it under
    the terms of the GNU Lesser General Public License (LGPL) as published
    by the Free Software Foundation; either version 2.1 of the License, or
    (at your option) any later version.  See <https://www.gnu.org/licenses/>.
*/

#include "flint.h"
#include "fq.h"
#include <stdio.h>
#include "profiler.h"

#ifndef REPS
#define REPS 1000000
#endif

int
main()
{
    timeit_t t0;

    slong i, d, cpu, wall;
    fmpz_t p;
    fq_ctx_t ctx;
    fq_t a,b;

    FLINT_TEST_INIT(state);
    
    fmpz_init(p);
    fmpz_set_ui(p, n_randprime(state, 2+ n_randint(state,3),1));
    d = n_randint(state,10)+1;
    fq_ctx_init_conway(ctx,p,d,"a");

    fq_init(a, ctx);
    fq_init(b, ctx);

    fq_randtest_not_zero(a,state,ctx);

    flint_printf("INV benchmark:single repeated inversion: \n");
    timeit_start(t0);
    for(i=0;i<REPS;i++) fq_inv(b,a,ctx);
    timeit_stop(t0);
    flint_printf ( " cpu = %wd ms, wall = %wd ms \n " , t0->cpu , t0->wall );

    flint_printf("random inversions: \n");

    wall = 0;
    cpu = 0;
    for(i=0;i<REPS;i++)
    {
    fq_randtest_not_zero(a,state,ctx);
    timeit_start(t0);
    fq_inv(b,a,ctx);
    timeit_stop(t0);
    cpu = cpu + t0->cpu;
    wall = wall + t0->wall;
    }

    flint_printf ( " cpu = %wd ms, wall = %wd ms \n " , cpu , wall );

    FLINT_TEST_CLEANUP(state);
    
    return 0;

}