/usr/include/linbox/field/ntl-pid-lzz_p.h is in liblinbox-dev 1.1.6~rc0-4.1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 | /* -*- mode: C++; tab-width: 8; indent-tabs-mode: t; c-basic-offset: 8 -*- */
/* linbox/field/ntl-pid-lzz_p.h
* written by bds
*
*/
#ifndef __NTL_PID_zz_p_H
#define __NTL_PID_zz_p_H
#include "linbox/field/ntl-lzz_p.h"
#include "linbox/util/debug.h"
#include "linbox/linbox-config.h"
#include <NTL/ZZ.h>
#include <linbox/field/field-traits.h>
// Namespace in which all LinBox library code resides
namespace LinBox
{
template <class Ring>
struct ClassifyRing;
class NTL_PID_zz_p;
template<>
struct ClassifyRing<NTL_PID_zz_p> {
typedef RingCategories::ModularTag categoryTag;
};
/** \brief extend Wrapper of zz_p from NTL. Add PID functions
\ingroup field
*/
struct NTL_PID_zz_p: public NTL_zz_p
{
protected: long _modulus;
public:
NTL_PID_zz_p(long pp, int exp = 1)
: NTL_zz_p(pp), _modulus(pp) {
if( exp != 1 ) throw PreconditionFailed(__FUNCTION__,__LINE__,"exponent must be 1");
}
Element& gcd(Element& g, const Element& a, const Element& b) const
{ g = NTL::GCD(NTL::rep(a), NTL::rep(b));
g = NTL::GCD(NTL::rep(g), _modulus);
return g;
}
Element& gcdin(Element& a, const Element& b) const
{ return gcd(a,a, b); }
bool isUnit(const Element& a) const
{ return 1 == NTL::GCD(NTL::rep(a), _modulus); }
Element& div(Element& c, const Element& a, const Element& b) const
{ return c = NTL::rep(a)/NTL::GCD(NTL::rep(a),NTL::rep(b)); }
Element& divin(Element& a, const Element& b) const
{ return div(a, a, b); }
static inline double getMaxModulus() { return (double)NTL_SP_BOUND; }
};
} // namespace LinBox
#endif // __NTL_PID_zz_p_H
|