/usr/include/givaro/givrnscstor.inl is in libgivaro-dev 4.0.2-8ubuntu1.
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 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 | // ==========================================================================
// $Source: /var/lib/cvs/Givaro/src/kernel/zpz/givrnscstor.inl,v $
// 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.
// Authors: T. Gaut%ier
// $Id: givrnscstor.inl,v 1.4 2011-02-02 16:23:56 briceboyer Exp $
// ==========================================================================
// Description:
#ifndef __GIVARO_rns_cstor_INL
#define __GIVARO_rns_cstor_INL
namespace Givaro {
// -- free memory allocated in array !
template<class RING, class Domain>
RNSsystem<RING,Domain>::~RNSsystem()
{
}
template<class RING, class Domain>
RNSsystem<RING,Domain>::RNSsystem ()
: _primes(0), _ck(0)
{}
template<class RING, class Domain>
RNSsystem<RING,Domain>::RNSsystem (const Self_t & R)
: _primes(R._primes, givWithCopy()),
_ck(R._ck, givWithCopy())
{}
template<class RING, class Domain>
void RNSsystem<RING,Domain>::setPrimes (const domains& inprimes)
{
_primes.allocate(0);
_primes.copy( inprimes );
_ck.reallocate(0);
}
// -- Array of primes are given
template<class RING, class Domain>
RNSsystem<RING,Domain>::RNSsystem( const domains& inprimes)
: _primes(inprimes, givWithCopy()),
_ck(0)
{
GIVARO_ASSERT( inprimes.size()>0, "[RNSsystem<RING,Domain>::RNSsystem] bad size of array");
}
// -- Computes Ck , Ck = (\prod_{i=0}^{k-1} primes[i])^(-1) % primes[k],
// for k=1..(_sz-1)
template<class RING, class Domain>
void RNSsystem<RING,Domain>::ComputeCk()
{
if (_ck.size() !=0) return; // -- already computed
// - reallocation of a new array :
int Size = (int) _primes.size();
_ck.reallocate((size_t)Size);
// _ck[0] = Neutral::zero; // -- undefined and never used
for (int k=1; k < Size; ++k)
{
modulo prod, tmp;
_primes[k].init(prod, _primes[0].characteristic());
for (int i= 1; i < k; ++i)
_primes[k].mulin(prod, _primes[k].init(tmp,_primes[i].characteristic()));
_primes[k].inv(_ck[k], prod);
}
}
template<class RING, class Domain>
const typename RNSsystem<RING,Domain>::domains& RNSsystem<RING,Domain>::Primes() const
{
return _primes;
}
template<class RING, class Domain>
const Domain RNSsystem<RING,Domain>::ith(const size_t i) const
{
return _primes[i];
}
template<class RING, class Domain>
const typename RNSsystem<RING,Domain>::array& RNSsystem<RING,Domain>::Reciprocals() const
{
if (_ck.size() ==0) ((RNSsystem<RING,Domain>*)this)->ComputeCk();
return _ck;
}
template<class RING, class Domain>
const typename RNSsystem<RING,Domain>::modulo RNSsystem<RING,Domain>::reciprocal(const size_t i) const
{
if (_ck.size() ==0) ((RNSsystem<RING,Domain>*)this)->ComputeCk();
return _ck[i];
}
} // namespace Givaro
#endif // __GIVARO_rns_cstor_INL
// vim:sts=8:sw=8:ts=8:noet:sr:cino=>s,f0,{0,g0,(0,\:0,t0,+0,=s:syntax=cpp.doxygen
|