examples/FiniteField/domain_to_operatorstyle.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.
 
#include <iostream>
#include <givaro/gfq.h>
#include <givaro/StaticElement.h>
 
using namespace Givaro;
 
namespace Givaro
{
 
    // Finite Field with Domain style
    typedef GFqDom<int64_t> Field;
 
    // Wrapper to give an operator style to the elements
    typedef StaticElement< Field > Element;
 
    // Mandatory declaration (because of static template)
    // and an actual constructed field is mandatory (the "(2)") for g++ 3.4
 
    template<>
    Field Element::_domain(2);
}
 
int main(int argc, char ** argv) {
    uint64_t P = (argc>1 ? (uint64_t)atoi(argv[1]) : 5009U);
 
    // Initialization of static member
    Element::setDomain( Field(P) );
 
    // Initialisations of elements
    Element a(2),b(-29.8),c(33),d(Integer("123456789012345678901234567890"));
 
 
    // Computations
    a = b;      std::cerr << a << " = " << b << " mod " << P << ";" << std::endl;
    a = b + c;  std::cerr << a << " = " << b << " + " << c << " mod " << P << ";" << std::endl;
    a = b - c;  std::cerr << a << " = " << b << " - " << c << " mod " << P << ";" << std::endl;
    a = b * c;  std::cerr << a << " = " << b << " * " << c << " mod " << P << ";" << std::endl;
    a = b / c;  std::cerr << a << " = " << b << " / " << c << " mod " << P << ";" << std::endl;
 
    // Computations in place
    std::cerr << d << " + " << a << " mod " << P << " = ";
    d += a; std::cerr << d << ";" << std::endl;
 
    std::cerr << d << " - " << a << " mod " << P << " = ";
    d -= a; std::cerr << d << ";" << std::endl;
 
    std::cerr << d << " * " << a << " mod " << P << " = ";
    d *= a; std::cerr << d << ";" << std::endl;
 
    std::cerr << d << " / " << a << " mod " << P << " = ";
    d /= a; std::cerr << d << ";" << std::endl;
 
    // Tests
 
    const Element zero(0);
 
    std::cerr << a << " is non zero is " << (a != zero ) << std::endl;
    std::cerr << a << " is non zero is " << (! Element::isZero(a)) << std::endl;
    std::cerr << a << " is non zero is " << (! a.isZero()) << std::endl;
 
    a = 0; std::cerr << a << " is zero is " << (a == zero ) << std::endl;
    std::cerr << a << " is zero is " << Element::isZero(a) << std::endl;
    std::cerr << a << " is zero is " << a.isZero() << std::endl;
 
    // Access to Field object
    const Field & F = Element::getDomain();
    F.write( std::cerr << "Test: within ") << 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