/usr/include/roboptim/core/numeric-linear-function.hxx is in libroboptim-core-dev 2.0-6.
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 | // Copyright (C) 2009 by Thomas Moulard, AIST, CNRS, INRIA.
//
// This file is part of the roboptim.
//
// roboptim 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 3 of the License, or
// (at your option) any later version.
//
// roboptim 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 roboptim. If not, see <http://www.gnu.org/licenses/>.
#ifndef ROBOPTIM_CORE_NUMERIC_LINEAR_FUNCTION_HXX
# define ROBOPTIM_CORE_NUMERIC_LINEAR_FUNCTION_HXX
# include "debug.hh"
# include <roboptim/core/indent.hh>
# include <roboptim/core/numeric-linear-function.hh>
# include <roboptim/core/util.hh>
namespace roboptim
{
template <typename T>
GenericNumericLinearFunction<T>::GenericNumericLinearFunction
(const matrix_t& a, const vector_t& b)
throw ()
: GenericLinearFunction<T>
(a.cols (), a.rows (), "numeric linear function"),
a_ (a),
b_ (b)
{
assert (b.size () == this->outputSize ());
}
template <typename T>
GenericNumericLinearFunction<T>::~GenericNumericLinearFunction () throw ()
{
}
// A * x + b
template <typename T>
void
GenericNumericLinearFunction<T>::impl_compute (result_t& result,
const argument_t& argument)
const throw ()
{
result.noalias () = a_* argument;
result += b_;
}
// A
template <typename T>
void
GenericNumericLinearFunction<T>::impl_jacobian
(jacobian_t& jacobian, const argument_t&) const throw ()
{
jacobian = this->a_;
}
// A(i)
template <typename T>
void
GenericNumericLinearFunction<T>::impl_gradient (gradient_t& gradient,
const argument_t&,
size_type idFunction) const throw ()
{
for (size_type j = 0; j < this->inputSize (); ++j)
gradient[j] = a_ (idFunction, j);
}
template <typename T>
std::ostream&
GenericNumericLinearFunction<T>::print (std::ostream& o) const throw ()
{
return o << "Numeric linear function" << incindent << iendl
<< "A = " << this->a_ << iendl
<< "B = " << this->b_
<< decindent;
}
} // end of namespace roboptim
#endif //! ROBOPTIM_CORE_QUADRATIC_FUNCTION_HXX
|