/usr/include/linbox/solutions/getentry.inl is in liblinbox-dev 1.4.2-5build1.
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 | /* linbox/solutions/getentry.inl
* Copyright(C) LinBox
* Evolved from an earlier one by Bradford Hovinen <hovinen@cis.udel.edu>
* -bds
*
*
* ========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_getentry_INL
#define __LINBOX_getentry_INL
#include <vector>
#include "linbox/util/debug.h"
#include "linbox/vector/vector-domain.h"
#include "linbox/matrix/dense-matrix.h"
#include "linbox/matrix/sparse-matrix.h"
#include "linbox/matrix/sparse-matrix.h"
#include "linbox/blackbox/scalar-matrix.h"
#include "linbox/blackbox/compose.h"
#include "linbox/blackbox/diagonal.h"
#include "linbox/solutions/methods.h"
#include "linbox/solutions/solution-tags.h"
#include "linbox/vector/blas-vector.h"
namespace LinBox
{
template <class BB>
typename BB::Field::Element& getEntry(typename BB::Field::Element& x, const BB& A, const size_t i, const size_t j)
{
typename GetEntryCategory<BB>::Tag t;
return getEntry(x, A, i, j, t);
}
// General case category of BBs requires use of apply.
template<class BB> struct GetEntryCategory{ typedef SolutionTags::Generic Tag; };
// Look for local specializations of this in blackbox/*
// To ignore methods
template <class BB, class Method>
typename BB::Field::Element& getEntry(typename BB::Field::Element& x, const BB& A, const size_t i, const size_t j, Method & m)
{
return getEntry(x, A, i, j);
}
// Generic BBs require use of apply.
template <class BB>
typename BB::Field::Element& getEntry(typename BB::Field::Element& x, const BB& A, const size_t i, const size_t j, SolutionTags::Generic t)
{
typedef typename BB::Field Field;
typedef BlasVector<Field> Vector ;
const Field& F = A.field();
Vector v(F,A.coldim(), F.zero), w(F,A.rowdim(), F.zero);
F.assign(v[j],F.one); // use standard basis ?
A.apply (w, v);
return F.assign(x, w[i]);
}
// some BBs have their own.
template <class BB>
typename BB::Field::Element& getEntry(typename BB::Field::Element& x, const BB& A, const size_t i, const size_t j, SolutionTags::Local t )
{
return A.getEntry(x, i, j);
}
// Compose< Diagonal, BB > specialization
template <class Field, class Trait, class BB>
typename Field::Element& getEntry(typename Field::Element& x, const Compose<Diagonal<Field, Trait>, BB>& A, const size_t i, const size_t j)
{
typename Field::Element y;
getEntry(y, *(A.getLeftPtr()), i, i);
getEntry(x, *(A.getRightPtr()), i, j);
return A.field().mulin(x, y);
}
// Compose< BB, Diagonal > specialization
template <class BB, class Field, class Trait>
typename Field::Element& getEntry(typename Field::Element& x, const Compose<BB, Diagonal<Field, Trait> >& A, const size_t i, const size_t j)
{
typename Field::Element y;
getEntry(y, *(A.getLeftPtr()), i, j);
getEntry(x, *(A.getRightPtr()), j, j);
return A.field().mulin(x, y);
}
// Compose< Diagonal, Diagonal > specialization
template <class Field, class T1, class T2>
typename Field::Element& getEntry(typename Field::Element& x, const Compose<Diagonal<Field,T1>, Diagonal<Field, T2> >& A, const size_t i, const size_t j)
{
if (i != j) return A.field().assign(x, A.field().zero);
typename Field::Element y;
getEntry(y, *(A.getLeftPtr()), i, i);
getEntry(x, *(A.getRightPtr()), j, j);
return A.field().mulin(x, y);
}
}
#endif // __LINBOX_getentry_INL
// Local Variables:
// mode: C++
// tab-width: 8
// indent-tabs-mode: nil
// c-basic-offset: 8
// End:
// vim:sts=8:sw=8:ts=8:noet:sr:cino=>s,f0,{0,g0,(0,\:0,t0,+0,=s
|