This file is indexed.

/usr/include/rdkit/Numerics/EigenSolvers/PowerEigenSolver.h is in librdkit-dev 201309-1.

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
//
//  Copyright (C) 2004-2006 Rational Discovery LLC
//
//   @@ All Rights Reserved @@
//  This file is part of the RDKit.
//  The contents are covered by the terms of the BSD license
//  which is included in the file license.txt, found at the root
//  of the RDKit source tree.
//

#ifndef _RD_POWER_EIGENSOLVER_H
#define _RD_POWER_EIGENSOLVER_H

#include <Numerics/Vector.h>
#include <Numerics/Matrix.h>
#include <Numerics/SymmMatrix.h>

namespace RDNumeric {
  namespace EigenSolvers {
    //! Compute the \c numEig largest eigenvalues and, optionally,  the corresponding
    //! eigenvectors. 
    /*!
      
    \param numEig       the number of eigenvalues we are interested in
    \param mat          symmetric input matrix of dimension N*N
    \param eigenValues  Vector used to return the eigenvalues (size = numEig)
    \param eigenVectors Optional matrix used to return the eigenvectors (size = N*numEig)
    \param seed         Optional values to seed the random value generator used to 
                        initialize the eigen vectors
    \return a boolean indicating whether or not the calculation converged.
    
    <b>Notes:</b>
    - The matrix, \c mat, is changed in this function
    
    <b>Algorithm:</b>
    
    We use the iterative power method, which works like this:

    \verbatim
     u = arbitrary unit vector
     tol = 0.001
     currEigVal = 0.0;
     prevEigVal = -1.0e100
     while (abs(currEigVal - prevEigVal) > tol) :
         v = Au
         prevEigVal = currEigVal
         currEigVal = v[i] // where i is the id os the largest absolute component
         u = c*v
    \endverbatim

      
    */
    bool powerEigenSolver(unsigned int numEig, DoubleSymmMatrix &mat,
                          DoubleVector &eigenValues,
                          DoubleMatrix *eigenVectors=0, 
                          int seed=-1);
    //! \overload
    static bool powerEigenSolver(unsigned int numEig, DoubleSymmMatrix &mat,
                          DoubleVector &eigenValues,
                          DoubleMatrix &eigenVectors, 
                          int seed=-1) {
      return powerEigenSolver(numEig,mat,eigenValues,&eigenVectors,seed);
    }
  };
};

#endif