dual_binary_operator_helper.hpp

Go to the documentation of this file.

#pragma once

namespace algo { namespace ad { namespace detail {
    /*--------------------------------------------------------------------------
     * dual_add_function
     *------------------------------------------------------------------------*/
    template<typename E1, typename E2>
    auto dual_add_function(const dual_expression<E1>& e1, const dual_expression<E2>& e2)
    -> decltype(e1().getDerivative() + e2().getDerivative())
    {
        return e1().getDerivative() + e2().getDerivative();
    }
    template<typename E2>
    auto dual_add_function(const double e1, const dual_expression<E2>& e2)
    -> decltype(e2().getDerivative())
    {
        return e2().getDerivative();
    }
    template<typename E1>
    auto dual_add_function(const dual_expression<E1>& e1, const double e2)
    -> decltype(e1().getDerivative())
    {
        return e1().getDerivative();
    }
    /*--------------------------------------------------------------------------
     * dual_minus_function
     *------------------------------------------------------------------------*/
    template<typename E1, typename E2>
    auto dual_minus_function(const dual_expression<E1>& e1, const dual_expression<E2>& e2)
    -> decltype(e1().getDerivative() - e2().getDerivative())
    {
        return e1().getDerivative() - e2().getDerivative();
    }
    template<typename E2>
    auto dual_minus_function(const double e1, const dual_expression<E2>& e2)
    -> decltype(e2().getDerivative())
    {
        return e2().getDerivative();
    }
    template<typename E1>
    auto dual_minus_function(const dual_expression<E1>& e1, const double e2)
    -> decltype(e1().getDerivative())
    {
        return e1().getDerivative();
    }
    /*--------------------------------------------------------------------------
     * dual_multiplies_function
     *------------------------------------------------------------------------*/
    template<typename E1, typename E2>
    auto dual_multiplies_function(
        const dual_expression<E1>& e1, 
        const dual_expression<E2>& e2)
    -> decltype(e1().getDerivative() * e2().getValue() + e1().getValue() * e2().getDerivative())
    {
        return e1().getDerivative() * e2().getValue() + e1().getValue() * e2().getDerivative();
    }
    template<typename E2>
    auto dual_multiplies_function(const double e1, const dual_expression<E2>& e2)
    -> decltype(e2().getDerivative())
    {
        return e2().getDerivative();
    }
    template<typename E1>
    auto dual_multiplies_function(const dual_expression<E1>& e1, const double e2)
    -> decltype(e1().getDerivative())
    {
        return e1().getDerivative();
    }
    /*--------------------------------------------------------------------------
     * dual_divide_function
     *------------------------------------------------------------------------*/
    template<typename E1, typename E2>
    auto dual_divide_function(
        const dual_expression<E1>& e1, 
        const dual_expression<E2>& e2)
    -> decltype((e1().getDerivative() * e2().getValue() - e1().getValue() * e2().getDerivative())
            / (e1().getValue() * e1().getValue()))
    {
        return (e1().getDerivative() * e2().getValue() - e1().getValue() * e2().getDerivative()) 
            / (e1().getValue() * e1().getValue());
    }
    template<typename E2>
    auto dual_divide_function(const double e1, const dual_expression<E2>& e2)
    -> decltype(e2().getDerivative())
    {
        return e2().getDerivative();
    }
    template<typename E1>
    auto dual_divide_function(const dual_expression<E1>& e1, const double e2)
    -> decltype(e1().getDerivative())
    {
        return e1().getDerivative();
    }
} } } // namespace algo { namespace ad { namespace detail {