acosh.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: BSD-2-Clause
 
#ifndef pRC_CORE_COMPLEX_FUNCTIONS_ACOSH_H
#define pRC_CORE_COMPLEX_FUNCTIONS_ACOSH_H
 
#include <prc/core/complex/complex.hpp>
 
namespace pRC
{
    template<IsComplex T>
    static inline constexpr auto acosh(T const &a)
    {
        using R = typename T::Type;
 
        return identity<R>(2) *
            log(sqrt(identity<R>(0.5) * (a + unit<R>())) +
                sqrt(identity<R>(0.5) * (a - unit<R>())));
    }
}
#endif // pRC_CORE_COMPLEX_FUNCTIONS_ACOSH_H