This file is indexed.

/usr/include/NTL/mat_RR.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
#ifndef NTL_mat_RR__H
#define NTL_mat_RR__H

#include <NTL/matrix.h>
#include <NTL/vec_vec_RR.h>

NTL_OPEN_NNS

typedef Mat<RR> mat_RR;

void add(mat_RR& X, const mat_RR& A, const mat_RR& B); 
void sub(mat_RR& X, const mat_RR& A, const mat_RR& B); 
void negate(mat_RR& X, const mat_RR& A);
void mul(mat_RR& X, const mat_RR& A, const mat_RR& B); 
void mul(vec_RR& x, const mat_RR& A, const vec_RR& b); 
void mul(vec_RR& x, const vec_RR& a, const mat_RR& B); 

void mul(mat_RR& X, const mat_RR& A, const RR& b);
void mul(mat_RR& X, const mat_RR& A, double b);

inline void mul(mat_RR& X, const RR& a, const mat_RR& B)
   { mul(X, B, a); }

inline void mul(mat_RR& X, double a, const mat_RR& B)
   { mul(X, B, a); }

void ident(mat_RR& X, long n); 
inline mat_RR ident_mat_RR(long n)
   { mat_RR X; ident(X, n); NTL_OPT_RETURN(mat_RR, X); }

void determinant(RR& d, const mat_RR& A);
long IsIdent(const mat_RR& A, long n);
void transpose(mat_RR& X, const mat_RR& A);
void solve(RR& d, vec_RR& X,
           const mat_RR& A, const vec_RR& b);
void inv(RR& d, mat_RR& X, const mat_RR& A);

inline void sqr(mat_RR& X, const mat_RR& A)
   { mul(X, A, A); }

inline mat_RR sqr(const mat_RR& A)
   { mat_RR X; sqr(X, A); NTL_OPT_RETURN(mat_RR, X); }

void inv(mat_RR& X, const mat_RR& A);

inline mat_RR inv(const mat_RR& A)
   { mat_RR X; inv(X, A); NTL_OPT_RETURN(mat_RR, X); }

void power(mat_RR& X, const mat_RR& A, const ZZ& e);
inline mat_RR power(const mat_RR& A, const ZZ& e)
   { mat_RR X; power(X, A, e); NTL_OPT_RETURN(mat_RR, X); }

inline void power(mat_RR& X, const mat_RR& A, long e)
   { power(X, A, ZZ_expo(e)); }
inline mat_RR power(const mat_RR& A, long e)
   { mat_RR X; power(X, A, e); NTL_OPT_RETURN(mat_RR, X); }



void diag(mat_RR& X, long n, const RR& d);
inline mat_RR diag(long n, const RR& d)
   { mat_RR X; diag(X, n, d); NTL_OPT_RETURN(mat_RR, X); }

long IsDiag(const mat_RR& A, long n, const RR& d);


// miscellaneous:

RR determinant(const mat_RR& a);
// functional variant of determinant

inline mat_RR transpose(const mat_RR & a)
   { mat_RR x; transpose(x, a); NTL_OPT_RETURN(mat_RR, x); }


void clear(mat_RR& a);
// x = 0 (dimension unchanged)

long IsZero(const mat_RR& a);
// test if a is the zero matrix (any dimension)


// operator notation:

mat_RR operator+(const mat_RR& a, const mat_RR& b);
mat_RR operator-(const mat_RR& a, const mat_RR& b);
mat_RR operator*(const mat_RR& a, const mat_RR& b);

mat_RR operator-(const mat_RR& a);


// matrix/vector multiplication:

vec_RR operator*(const mat_RR& a, const vec_RR& b);

vec_RR operator*(const vec_RR& a, const mat_RR& b);



// matrix/scalar multiplication:

inline mat_RR operator*(const mat_RR& a, const RR& b)
   { mat_RR x; mul(x, a, b); NTL_OPT_RETURN(mat_RR, x); }

inline mat_RR operator*(const mat_RR& a, double b)
   { mat_RR x; mul(x, a, b); NTL_OPT_RETURN(mat_RR, x); }

inline mat_RR operator*(const RR& a, const mat_RR& b)
   { mat_RR x; mul(x, a, b); NTL_OPT_RETURN(mat_RR, x); }

inline mat_RR operator*(double a, const mat_RR& b)
   { mat_RR x; mul(x, a, b); NTL_OPT_RETURN(mat_RR, x); }




// assignment operator notation:

inline mat_RR& operator+=(mat_RR& x, const mat_RR& a)
{
   add(x, x, a);
   return x;
}   

inline mat_RR& operator-=(mat_RR& x, const mat_RR& a)
{
   sub(x, x, a);
   return x;
}   


inline mat_RR& operator*=(mat_RR& x, const mat_RR& a)
{
   mul(x, x, a);
   return x;
}   

inline mat_RR& operator*=(mat_RR& x, const RR& a)
{
   mul(x, x, a);
   return x;
}   

inline mat_RR& operator*=(mat_RR& x, double a)
{
   mul(x, x, a);
   return x;
}   
   

inline vec_RR& operator*=(vec_RR& x, const mat_RR& a)
{
   mul(x, x, a);
   return x;
}   


NTL_CLOSE_NNS


#endif