/usr/include/linbox/blackbox/jit-matrix.h is in liblinbox-dev 1.3.2-1.1build2.
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 | /* linbox/blackbox/jit-matrix.h
* Copyright (c) LinBox
*
* bds, jpm
*
*
* ========LICENCE========
* This file is part of the library LinBox.
*
* LinBox is free software: you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with this library; if not, write to the Free Software
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
* ========LICENCE========
*.
*/
#ifndef __LINBOX_jitmatrix_H
#define __LINBOX_jitmatrix_H
#include "linbox/blackbox/blackbox-interface.h"
namespace LinBox
{
/**
\brief Example of a blackbox that is space efficient, though not time efficient.
\ingroup blackbox
Just In Time Matrix.
The matrix itself is not stored in memory. Rather, an EntryGenerator
function is called to provide the entries. The entry generator is
called once for each entry during an apply or applyTranspose
operation.
An toy example of its use is the Hilbert matrix, whose \f$i,j\f$ entry is
generated by the formula \f$1/(i+j+2)\f$ in zero based indexing.
The motivating examples were matrices also defined by formula, the Paley type matrices.
\sa{MSW07}% ISSAC 07 paper
In that context block structured turned out to be essential and the
JIT_Matrix class is primarily intended for block structured matrices,
the JIT entries being matrix blocks.
@param _Field only need provide the \c init() and \c axpyin() functions.
@param JIT_EntryGenerator \c gen() is a function object defining the
matrix by providing <code>gen(e, i, j)</code> which sets field element e to the \c i,j entry
of the matrix. Indexing is zero based.
*/
template <class _Field, class JIT_EntryGenerator>
class JIT_Matrix {
public:
typedef _Field Field;
typedef typename Field::Element Element;
typedef MatrixCategories::BlackboxTag MatrixCategory;
/**
* m by n matrix is constructed.
* JIT(Field::Element& e, size_t i, size_t j) is a function object which
* assigns the i,j entry to e (and returns a reference to e)
* and must be valid for 0 <= i < m, 0 <= j < n.
**/
JIT_Matrix (_Field& F, const size_t m, const size_t n,
const JIT_EntryGenerator& JIT) :
_field(F), _m(m), _n(n), _gen(JIT)
{};
template<class OutVector, class InVector>
OutVector& apply (OutVector& y, const InVector& x) ;
//OutVector& apply (OutVector& y, const InVector& x) const;
template<class OutVector, class InVector>
OutVector& applyTranspose (OutVector& y, const InVector& x);
//OutVector& applyTranspose (OutVector& y, const InVector& x) const;
size_t rowdim (void) const { return _m; }
size_t coldim (void) const { return _n; }
const Field& field() const { return _field; }
protected:
// Field for arithmetic
Field _field;
// Number of rows and columns of matrix.
size_t _m;
size_t _n;
// STL vector of field elements used in applying matrix.
JIT_EntryGenerator _gen;
}; // class JIT_Matrix
// Method implementations
template <class Field, class JIT_EntryGenerator>
template <class OutVector, class InVector>
inline OutVector& JIT_Matrix<Field, JIT_EntryGenerator>::apply (OutVector& y, const InVector& x)
{
Element entry;
_field.init(entry,0);
for (size_t i = 0; i < _m; ++i) {
_field.init(y[i], 0);
for (size_t j = 0; j < _n; ++j) {
_gen(entry, i, j);
_field.axpyin (y[i], entry, x[j]);
}
}
return y;
} //apply
template <class Field, class JIT_EntryGenerator>
template <class OutVector, class InVector>
inline OutVector& JIT_Matrix<Field, JIT_EntryGenerator>::applyTranspose (OutVector& y, const InVector& x)
{
Element entry;
_field.init(entry,0);
for (size_t i = 0; i < _m; ++i) {
_field.init(y[i], 0);
for (size_t j = 0; j < _n; ++j) {
_field.axpyin ( y[i], x[j], _gen(entry, j, i) );
}
}
return y;
} // applyTranspose
// Example: Generator to create psuedo-random entries
// !WARNING! repeated calls will give different values for the same entry
template < class Field >
class JIT_RandomEntryGenerator {
typename Field::RandIter _r;
size_t _b;
public:
JIT_RandomEntryGenerator(Field& F, size_t b) :
_r(F), _b(b)
{}
typename Field::Element& operator()(typename Field::Element& e,
size_t k, size_t l)
{
return _r.random(e);
}
};
} // namespace LinBox
#endif // __LINBOX_jitmatrix_H
// vim:sts=8:sw=8:ts=8:noet:sr:cino=>s,f0,{0,g0,(0,:0,t0,+0,=s
// Local Variables:
// mode: C++
// tab-width: 8
// indent-tabs-mode: nil
// c-basic-offset: 8
// End:
|