util.hpp

#pragma once
#include "algo/ad/dual.hpp"
#include <boost/numeric/ublas/vector.hpp>

namespace algo { namespace ad {
    template <typename V>
    boost::numeric::ublas::vector<V>
    apply_element(
        const std::vector<std::function<
            V (const boost::numeric::ublas::vector<V>&)> >& functions,
        const boost::numeric::ublas::vector<V>& x)
    {
        namespace ublas = boost::numeric::ublas;
        typedef std::function<V (const ublas::vector<V>&)> function_type;

        ublas::vector<V> result(functions.size());
        auto transformer = [&x](const function_type& f) {
            return f(x);
        };
        std::transform(
            functions.begin(), functions.end(), result.begin(), transformer);
        return result;
    }
    dual<boost::numeric::ublas::vector<double> >
    make_dual(
        const double value,
        const boost::numeric::ublas::vector<double>& infinitesimal);
    dual<boost::numeric::ublas::vector<double> >
    make_unit_dual(
        const double value,
        const std::size_t size,
        const std::size_t index);
    boost::numeric::ublas::vector<dual<boost::numeric::ublas::vector<double> > >
    make_vector_dual(const boost::numeric::ublas::vector<double>& value);
} } // namespace algo { namespace ad {