/usr/include/NTL/mat_lzz_p.h is in libntl-dev 9.9.1-3.
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 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 | #ifndef NTL_mat_zz_p__H
#define NTL_mat_zz_p__H
#include <NTL/matrix.h>
#include <NTL/vec_vec_lzz_p.h>
NTL_OPEN_NNS
typedef Mat<zz_p> mat_zz_p;
void add(mat_zz_p& X, const mat_zz_p& A, const mat_zz_p& B);
void sub(mat_zz_p& X, const mat_zz_p& A, const mat_zz_p& B);
void negate(mat_zz_p& X, const mat_zz_p& A);
void mul(mat_zz_p& X, const mat_zz_p& A, const mat_zz_p& B);
void mul(vec_zz_p& x, const mat_zz_p& A, const vec_zz_p& b);
void mul(vec_zz_p& x, const vec_zz_p& a, const mat_zz_p& B);
void mul(mat_zz_p& X, const mat_zz_p& A, zz_p b);
void mul(mat_zz_p& X, const mat_zz_p& A, long b);
inline void mul(mat_zz_p& X, zz_p a, const mat_zz_p& B)
{ mul(X, B, a); }
inline void mul(mat_zz_p& X, long a, const mat_zz_p& B)
{ mul(X, B, a); }
void ident(mat_zz_p& X, long n);
inline mat_zz_p ident_mat_zz_p(long n)
{ mat_zz_p X; ident(X, n); NTL_OPT_RETURN(mat_zz_p, X); }
long IsIdent(const mat_zz_p& A, long n);
void transpose(mat_zz_p& X, const mat_zz_p& A);
// ************************
void relaxed_solve(zz_p& d, vec_zz_p& x, const mat_zz_p& A, const vec_zz_p& b, bool relax=true);
void relaxed_solve(zz_p& d, const mat_zz_p& A, vec_zz_p& x, const vec_zz_p& b, bool relax=true);
void relaxed_inv(zz_p& d, mat_zz_p& X, const mat_zz_p& A, bool relax=true);
inline void relaxed_inv(mat_zz_p& X, const mat_zz_p& A, bool relax=true)
{ zz_p d; relaxed_inv(d, X, A, relax); if (d == 0) ArithmeticError("inv: non-invertible matrix"); }
inline mat_zz_p relaxed_inv(const mat_zz_p& A, bool relax=true)
{ mat_zz_p X; relaxed_inv(X, A, relax); NTL_OPT_RETURN(mat_zz_p, X); }
void relaxed_determinant(zz_p& d, const mat_zz_p& A, bool relax=true);
inline zz_p relaxed_determinant(const mat_zz_p& a, bool relax=true)
{ zz_p x; relaxed_determinant(x, a, relax); return x; }
void relaxed_power(mat_zz_p& X, const mat_zz_p& A, const ZZ& e, bool relax=true);
inline mat_zz_p relaxed_power(const mat_zz_p& A, const ZZ& e, bool relax=true)
{ mat_zz_p X; relaxed_power(X, A, e, relax); NTL_OPT_RETURN(mat_zz_p, X); }
inline void relaxed_power(mat_zz_p& X, const mat_zz_p& A, long e, bool relax=true)
{ relaxed_power(X, A, ZZ_expo(e), relax); }
inline mat_zz_p relaxed_power(const mat_zz_p& A, long e, bool relax=true)
{ mat_zz_p X; relaxed_power(X, A, e, relax); NTL_OPT_RETURN(mat_zz_p, X); }
// ***********************
inline void solve(zz_p& d, vec_zz_p& x, const mat_zz_p& A, const vec_zz_p& b)
{ relaxed_solve(d, x, A, b, false); }
inline void solve(zz_p& d, const mat_zz_p& A, vec_zz_p& x, const vec_zz_p& b)
{ relaxed_solve(d, A, x, b, false); }
inline void inv(zz_p& d, mat_zz_p& X, const mat_zz_p& A)
{ relaxed_inv(d, X, A, false); }
inline void inv(mat_zz_p& X, const mat_zz_p& A)
{ relaxed_inv(X, A, false); }
inline mat_zz_p inv(const mat_zz_p& A)
{ return relaxed_inv(A, false); }
inline void determinant(zz_p& d, const mat_zz_p& A)
{ relaxed_determinant(d, A, false); }
inline zz_p determinant(const mat_zz_p& a)
{ return relaxed_determinant(a, false); }
inline void power(mat_zz_p& X, const mat_zz_p& A, const ZZ& e)
{ relaxed_power(X, A, e, false); }
inline mat_zz_p power(const mat_zz_p& A, const ZZ& e)
{ return relaxed_power(A, e, false); }
inline void power(mat_zz_p& X, const mat_zz_p& A, long e)
{ relaxed_power(X, A, e, false); }
inline mat_zz_p power(const mat_zz_p& A, long e)
{ return relaxed_power(A, e, false); }
// ************************
inline void sqr(mat_zz_p& X, const mat_zz_p& A)
{ mul(X, A, A); }
inline mat_zz_p sqr(const mat_zz_p& A)
{ mat_zz_p X; sqr(X, A); NTL_OPT_RETURN(mat_zz_p, X); }
void diag(mat_zz_p& X, long n, zz_p d);
inline mat_zz_p diag(long n, zz_p d)
{ mat_zz_p X; diag(X, n, d); NTL_OPT_RETURN(mat_zz_p, X); }
long IsDiag(const mat_zz_p& A, long n, zz_p d);
long gauss(mat_zz_p& M);
long gauss(mat_zz_p& M, long w);
void image(mat_zz_p& X, const mat_zz_p& A);
void kernel(mat_zz_p& X, const mat_zz_p& A);
// miscellaneous:
inline mat_zz_p transpose(const mat_zz_p& a)
{ mat_zz_p x; transpose(x, a); NTL_OPT_RETURN(mat_zz_p, x); }
void clear(mat_zz_p& a);
// x = 0 (dimension unchanged)
long IsZero(const mat_zz_p& a);
// test if a is the zero matrix (any dimension)
// operator notation:
mat_zz_p operator+(const mat_zz_p& a, const mat_zz_p& b);
mat_zz_p operator-(const mat_zz_p& a, const mat_zz_p& b);
mat_zz_p operator*(const mat_zz_p& a, const mat_zz_p& b);
mat_zz_p operator-(const mat_zz_p& a);
// matrix/scalar multiplication:
inline mat_zz_p operator*(const mat_zz_p& a, zz_p b)
{ mat_zz_p x; mul(x, a, b); NTL_OPT_RETURN(mat_zz_p, x); }
inline mat_zz_p operator*(const mat_zz_p& a, long b)
{ mat_zz_p x; mul(x, a, b); NTL_OPT_RETURN(mat_zz_p, x); }
inline mat_zz_p operator*(zz_p a, const mat_zz_p& b)
{ mat_zz_p x; mul(x, a, b); NTL_OPT_RETURN(mat_zz_p, x); }
inline mat_zz_p operator*(long a, const mat_zz_p& b)
{ mat_zz_p x; mul(x, a, b); NTL_OPT_RETURN(mat_zz_p, x); }
// matrix/vector multiplication:
vec_zz_p operator*(const mat_zz_p& a, const vec_zz_p& b);
vec_zz_p operator*(const vec_zz_p& a, const mat_zz_p& b);
// assignment operator notation:
inline mat_zz_p& operator+=(mat_zz_p& x, const mat_zz_p& a)
{
add(x, x, a);
return x;
}
inline mat_zz_p& operator-=(mat_zz_p& x, const mat_zz_p& a)
{
sub(x, x, a);
return x;
}
inline mat_zz_p& operator*=(mat_zz_p& x, const mat_zz_p& a)
{
mul(x, x, a);
return x;
}
inline mat_zz_p& operator*=(mat_zz_p& x, zz_p a)
{
mul(x, x, a);
return x;
}
inline mat_zz_p& operator*=(mat_zz_p& x, long a)
{
mul(x, x, a);
return x;
}
inline vec_zz_p& operator*=(vec_zz_p& x, const mat_zz_p& a)
{
mul(x, x, a);
return x;
}
NTL_CLOSE_NNS
#endif
|