///////////////////////////////////////////////////////////////////////////////////
/// OpenGL Mathematics (glm.g-truc.net)
///
/// Copyright (c) 2005 - 2015 G-Truc Creation (www.g-truc.net)
/// Permission is hereby granted, free of charge, to any person obtaining a copy
/// of this software and associated documentation files (the "Software"), to deal
/// in the Software without restriction, including without limitation the rights
/// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
/// copies of the Software, and to permit persons to whom the Software is
/// furnished to do so, subject to the following conditions:
/// 
/// The above copyright notice and this permission notice shall be included in
/// all copies or substantial portions of the Software.
/// 
/// Restrictions:
///             By making use of the Software for military purposes, you choose to make
///             a Bunny unhappy.
/// 
/// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
/// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
/// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
/// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
/// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
/// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
/// THE SOFTWARE.
///
/// @ref gtc_reciprocal
/// @file glm/gtc/reciprocal.inl
/// @date 2008-10-09 / 2012-04-07
/// @author Christophe Riccio
///////////////////////////////////////////////////////////////////////////////////

#include "../trigonometric.hpp"
#include <limits>

namespace glm
{
        // sec
        template <typename genType>
        GLM_FUNC_QUALIFIER genType sec(genType angle)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'sec' only accept floating-point values");
                return genType(1) / glm::cos(angle);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> sec(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'sec' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(sec, x);
        }

        // csc
        template <typename genType>
        GLM_FUNC_QUALIFIER genType csc(genType angle)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'csc' only accept floating-point values");
                return genType(1) / glm::sin(angle);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> csc(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'csc' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(csc, x);
        }

        // cot
        template <typename genType>
        GLM_FUNC_QUALIFIER genType cot(genType angle)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'cot' only accept floating-point values");
        
                genType const pi_over_2 = genType(3.1415926535897932384626433832795 / 2.0);
                return glm::tan(pi_over_2 - angle);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> cot(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'cot' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(cot, x);
        }

        // asec
        template <typename genType>
        GLM_FUNC_QUALIFIER genType asec(genType x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'asec' only accept floating-point values");
                return acos(genType(1) / x);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> asec(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'asec' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(asec, x);
        }

        // acsc
        template <typename genType>
        GLM_FUNC_QUALIFIER genType acsc(genType x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'acsc' only accept floating-point values");
                return asin(genType(1) / x);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> acsc(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'acsc' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(acsc, x);
        }

        // acot
        template <typename genType>
        GLM_FUNC_QUALIFIER genType acot(genType x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'acot' only accept floating-point values");

                genType const pi_over_2 = genType(3.1415926535897932384626433832795 / 2.0);
                return pi_over_2 - atan(x);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> acot(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'acot' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(acot, x);
        }

        // sech
        template <typename genType>
        GLM_FUNC_QUALIFIER genType sech(genType angle)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'sech' only accept floating-point values");
                return genType(1) / glm::cosh(angle);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> sech(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'sech' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(sech, x);
        }

        // csch
        template <typename genType>
        GLM_FUNC_QUALIFIER genType csch(genType angle)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'csch' only accept floating-point values");
                return genType(1) / glm::sinh(angle);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> csch(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'csch' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(csch, x);
        }

        // coth
        template <typename genType>
        GLM_FUNC_QUALIFIER genType coth(genType angle)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'coth' only accept floating-point values");
                return glm::cosh(angle) / glm::sinh(angle);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> coth(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'coth' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(coth, x);
        }

        // asech
        template <typename genType>
        GLM_FUNC_QUALIFIER genType asech(genType x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'asech' only accept floating-point values");
                return acosh(genType(1) / x);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> asech(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'asech' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(asech, x);
        }

        // acsch
        template <typename genType>
        GLM_FUNC_QUALIFIER genType acsch(genType x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'acsch' only accept floating-point values");
                return acsch(genType(1) / x);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> acsch(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'acsch' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(acsch, x);
        }

        // acoth
        template <typename genType>
        GLM_FUNC_QUALIFIER genType acoth(genType x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_iec559, "'acoth' only accept floating-point values");
                return atanh(genType(1) / x);
        }

        template <typename T, precision P, template <typename, precision> class vecType>
        GLM_FUNC_QUALIFIER vecType<T, P> acoth(vecType<T, P> const & x)
        {
                GLM_STATIC_ASSERT(std::numeric_limits<T>::is_iec559, "'acoth' only accept floating-point inputs");
                return detail::functor1<T, T, P, vecType>::call(acoth, x);
        }
}//namespace glm