This file is indexed.

/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