libntl-dev 10.5.0-2 / usr / include / NTL / MatPrime.h

This file is indexed.

/usr/include/NTL/MatPrime.h is in libntl-dev 10.5.0-2.

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
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
#ifndef NTL_MatPrime__H
#define NTL_MatPrime__H

#include <NTL/ZZ.h>
#include <NTL/ZZVec.h>
#include <NTL/lzz_p.h>
#include <NTL/vector.h>
#include <NTL/SmartPtr.h>
#include <NTL/LazyTable.h>

NTL_OPEN_NNS

#define NTL_MatPrimeFudge (3)
// similar to the FFTPrime strategy...ensures
// we can use floating point to approximate a quotient

#define NTL_MatPrimeLimit (1L << 20)
// Limit on dimension for matrix mul



#ifdef NTL_HAVE_AVX
#define NTL_MatPrime_NBITS (23)
#else
#define NTL_MatPrime_NBITS NTL_SP_NBITS
#endif

#if (NTL_MatPrime_NBITS > NTL_SP_NBITS)
// This is mainly academic
#define NTL_MatPrime_NBITS NTL_SP_NBITS
#endif

#if (NTL_MatPrime_NBITS < NTL_BITS_PER_INT)
typedef int MatPrime_residue_t;
#else
typedef long MatPrime_residue_t;
#endif

#if (2*NTL_MatPrime_NBITS+1 <= NTL_SP_NBITS)
#define NTL_MatPrime_HALF_SIZE_STRATEGY
#endif


struct MatPrimeInfo {
   long q;   
   zz_pContext context; 
};

void InitMatPrimeInfo(MatPrimeInfo& info, long q, long w);


#define NTL_MAX_MATPRIMES (20000)


typedef LazyTable<MatPrimeInfo, NTL_MAX_MATPRIMES> MatPrimeTablesType;

extern MatPrimeTablesType MatPrimeTables;
// a truly GLOBAL variable, shared among all threads


inline 
long GetMatPrime(long i)
{
   return MatPrimeTables[i]->q;
}

inline
void RestoreMatPrime(long i)
{
   MatPrimeTables[i]->context.restore();
}

void UseMatPrime(long index);
// allocates and initializes information for Mat prime



#ifndef NTL_MatPrime_HALF_SIZE_STRATEGY


struct MatPrime_crt_helper_scratch {
   ZZ t;
};

struct MatPrime_crt_helper {

   long NumPrimes;
   long sz;
   ZZ MinusMModP;  //  -M mod p, M = product of primes

   // the following arrays are indexed 0..NumPrimes-1
   // q[i] = MatPrime[i]
   Vec<long> prime;  // prime[i] = q[i]
   Vec<double> prime_recip;  // prime_recip[i] = 1/double(q[i])
   Vec<long> u;  // u[i] = (M/q[i])^{-1} mod q[i]
   Vec<mulmod_precon_t> uqinv;
   Vec<const sp_ZZ_reduce_struct*> ZZ_red_struct;

   ZZVec coeff;

   ZZ_ReduceStructAdapter montgomery_struct;

   long GetNumPrimes() const { return NumPrimes; }

   double cost;
   double GetCost() const { return cost; }


};


#else



struct MatPrime_crt_helper_scratch {
   ZZ t;
};

struct MatPrime_crt_helper {

   long NumPrimes;
   long sz;
   ZZ MinusMModP;  //  -M mod p, M = product of primes

   // the following arrays are indexed 0..NumPrimes-1
   // q[i] = MatPrime[i]
   Vec<long> prime;  // prime[i] = q[i]
   Vec<double> prime_recip;  // prime_recip[i] = 1/double(q[i])
   Vec<long> u;  // u[i] = (M/q[i])^{-1} mod q[i]
   Vec<mulmod_precon_t> uqinv;
   Vec<sp_reduce_struct> red_struct;

   // Indexed 0..ceil(NumPrimes/2)-1
   Vec<sp_ZZ_reduce_struct> ZZ_red_struct;

   // Indexed 0..ceil(NumPrimes/2)-1
   ZZVec coeff;

   ZZ_ReduceStructAdapter montgomery_struct;

   long GetNumPrimes() const { return NumPrimes; }

   double cost;
   double GetCost() const { return cost; }


};

#endif

void build(MatPrime_crt_helper& H, const ZZ& P);

void init_scratch(const MatPrime_crt_helper& H, MatPrime_crt_helper_scratch& scratch);

void reduce(const MatPrime_crt_helper& H, const ZZ& value, MatPrime_residue_t *remainders, 
	    MatPrime_crt_helper_scratch& scratch);

void reconstruct(const MatPrime_crt_helper& H, ZZ& value, const MatPrime_residue_t *remainders, 
		 MatPrime_crt_helper_scratch& scratch);







NTL_CLOSE_NNS

#endif

APT Browse - Built by Thomas Orozco.