examples/FiniteField/GFirreducible.C

NO DOC

// ========================================================== //
// Copyright(c)'1994-2009 by The Givaro group
// This file is part of Givaro.
// Givaro is governed by the CeCILL-B license under French law
// and abiding by the rules of distribution of free software.
// see the COPYRIGHT file for more details.
// Time-stamp: <12 Jun 15 18:47:22 Jean-Guillaume.Dumas@imag.fr>
// Check of irreducible polynomial and generators of GFq
// ========================================================== //
//
#include <givaro/gfq.h>
#include <givaro/givpower.h>
#include <givaro/givtimer.h>
 
using namespace Givaro;
 
 
 
int main(int argc, char** argv)
{
 
    int64_t p = (argc>1?atoi(argv[1]):5);
    int64_t e = (argc>2?atoi(argv[2]):3);
    typedef GFqDom<int64_t>::Residu_t TT ;
    GFqDom<int64_t> GFq((TT)p, (TT)e);
    GFqDom<int64_t> PrimeField((TT)p,1);
    std::cout << "Working in GF(" << p << '^' << e << ')' << std::endl;
    std::cout << "Elements are polynomials in X modulo " << p << std::endl;
 
    Poly1Dom< GFqDom<int64_t>, Dense > Pdom( PrimeField, Indeter("X") );
 
    // First get the irreducible polynomial irred
    // via irred = X^e - mQ
    // and X^e = X^(e-1) * X
    GFqDom<int64_t>::Element temo, t, tmp;
    Poly1Dom< GFqDom<int64_t>, Dense >::Element G, H, J, mQ, irred, modP;
 
    Pdom.init(G, Degree((int64_t)GFq.exponent()-1)); // X^(e-1)
    GFq.init(temo,G); // internal representation
 
    Pdom.init(H, Degree(1) ); // X
    GFq.init(t,H); // internal representation
 
    GFq.init(tmp);
 
    GFq.mul(tmp, temo, t); // internal representation of X^e
    GFq.negin(tmp);        // internal representation of -X^e, i.e. of the irreducible polynomial without the largest monomial, X^e
 
    int64_t lowerpart;
    // p-adic value of the lower part of the irreducible polynomial
    GFq.convert(lowerpart, tmp);
    std::cout << ' ' << p << "-adic value of the lower part of the irreducible : " << lowerpart << std::endl;
 
    int64_t ptoe = power(p,e); // p-adic value of X^e
    std::cout << ' ' << p << '^' << e << " is : " << ptoe << std::endl;
 
    std::cout << " --> Computed irreducible: " << ptoe+lowerpart << std::endl;
    std::cout << "Stored        irreducible: " << GFq.irreducible() << std::endl;
 
    Poly1PadicDom< GFqDom<int64_t>, Dense > PAD(Pdom);
    Poly1PadicDom< GFqDom<int64_t>, Dense >::Element Polynomial;
 
    PAD.radix(Polynomial, GFq.irreducible());
 
    std::cout << "Irreducible polynomial coefficients: ";
    for(Poly1PadicDom< GFqDom<int64_t>, Dense >::Element::iterator it = Polynomial.begin(); it != Polynomial.end(); ++it)
        PrimeField.write(std::cout << ' ', *it);
    std::cout << std::endl;
 
 
    PAD.write(std::cout << "The latter " << GFq.irreducible() << " represents: ", Polynomial)
    << " in " << p << "-adic"
    << std::endl;
 
 
 
    return 0;
 
}
/* -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */
// vim:sts=4:sw=4:ts=4:et:sr:cino=>s,f0,{0,g0,(0,\:0,t0,+0,=s