dual_unary_operator_helper.hpp

Go to the documentation of this file.

#pragma once
#include "algo/ad/ublas_unary_operator.hpp"

namespace algo { namespace ad { namespace detail {
    /*--------------------------------------------------------------------------
     * dual_negate_helper
     *------------------------------------------------------------------------*/
    template<typename E>
    auto dual_negate_helper(const dual_expression<E>& e)
    -> decltype(-e().getDerivative())
    {
        return -e().getDerivative();
    }
    double dual_negate_helper(const double e)
    {
        return -e;
    }
    /*--------------------------------------------------------------------------
     * dual_exp_helper
     *------------------------------------------------------------------------*/
    template<typename E>
    auto dual_exp_helper(const dual_expression<E>& e)
    -> decltype(boost::numeric::ublas::exp(e().getDerivative()))
    {
        return boost::numeric::ublas::exp(e().getDerivative());
    }
    double dual_exp_helper(const double e)
    {
        return std::exp(e);
    }
    /*--------------------------------------------------------------------------
     * dual_log_helper
     *------------------------------------------------------------------------*/
    template<typename E>
    auto dual_log_helper(const dual_expression<E>& e)
    -> decltype(boost::numeric::ublas::element_inverse(e().getDerivative()))
    {
        return boost::numeric::ublas::element_inverse(e().getDerivative());
    }
    double dual_log_helper(const double e)
    {
        return 1.0 / e;
    }
    /*--------------------------------------------------------------------------
     * dual_sin_helper
     *------------------------------------------------------------------------*/
    template<typename E>
    auto dual_sin_helper(const dual_expression<E>& e)
    -> decltype(boost::numeric::ublas::cos(e().getDerivative()))
    {
        return boost::numeric::ublas::cos(e().getDerivative());
    }
    double dual_sin_helper(const double e)
    {
        return std::cos(e);
    }
    /*--------------------------------------------------------------------------
     * dual_cos_helper
     *------------------------------------------------------------------------*/
    template<typename E>
    auto dual_cos_helper(const dual_expression<E>& e)
    -> decltype(-boost::numeric::ublas::sin(e().getDerivative()))
    {
        return -boost::numeric::ublas::sin(e().getDerivative());
    }
    double dual_cos_helper(const double e)
    {
        return -std::sin(e);
    }
} } } // namespace algo { namespace ad { namespace detail {