/usr/include/roboptim/core/differentiable-function.hxx is in libroboptim-core-dev 2.0-7.
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 | // 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_DIFFERENTIABLE_FUNCTION_HXX
# define ROBOPTIM_CORE_DIFFERENTIABLE_FUNCTION_HXX
# include <boost/algorithm/string/replace.hpp>
# include "roboptim/core/indent.hh"
# include "roboptim/core/util.hh"
namespace roboptim
{
template <typename T>
GenericDifferentiableFunction<T>::GenericDifferentiableFunction
(size_type inputSize, size_type outputSize, std::string name)
throw ()
: GenericFunction<T> (inputSize, outputSize, name)
{
}
template <>
inline void
GenericDifferentiableFunction<EigenMatrixSparse>::impl_jacobian
(jacobian_t& jacobian, const argument_t& argument)
const throw ()
{
#ifndef ROBOPTIM_DO_NOT_CHECK_ALLOCATION
Eigen::internal::set_is_malloc_allowed (true);
#endif //! ROBOPTIM_DO_NOT_CHECK_ALLOCATION
typedef Eigen::Triplet<double> triplet_t;
std::vector<triplet_t> coefficients;
for (jacobian_t::Index i = 0; i < this->outputSize (); ++i)
{
gradient_t grad = gradient (argument, i);
const unsigned int i_ = static_cast<const unsigned int> (i);
for (gradient_t::InnerIterator it (grad); it; ++it)
{
const unsigned int idx = static_cast<const unsigned int> (it.index ());
coefficients.push_back
(triplet_t (i_, idx, it.value ()));
}
}
jacobian.setFromTriplets (coefficients.begin (), coefficients.end ());
}
template <typename T>
void
GenericDifferentiableFunction<T>::impl_jacobian (jacobian_t& jacobian,
const argument_t& argument)
const throw ()
{
#ifndef ROBOPTIM_DO_NOT_CHECK_ALLOCATION
Eigen::internal::set_is_malloc_allowed (true);
#endif //! ROBOPTIM_DO_NOT_CHECK_ALLOCATION
for (typename jacobian_t::Index i = 0; i < this->outputSize (); ++i)
jacobian.row (i) = gradient (argument, i);
}
template <typename T>
std::ostream&
GenericDifferentiableFunction<T>::print (std::ostream& o) const throw ()
{
if (this->getName ().empty ())
return o << "Differentiable function";
std::stringstream ss;
ss << std::endl;
char fill = o.fill (' ');
ss << std::setw ((int)indent (o))
<< ""
<< std::setfill (fill);
std::string name = this->getName ();
boost::algorithm::replace_all (name, "\n", ss.str ());
return o << name << " (differentiable function)";
}
} // end of namespace roboptim
#endif //! ROBOPTIM_CORE_DIFFERENTIABLE_FUNCTION
|