#ifndef FOSSIL_MATH_TRIG_H
#define FOSSIL_MATH_TRIG_H

#include "math.h"

#ifdef __cplusplus
extern "C"
{
#endif

// *****************************************************************************
// Function prototypes
// *****************************************************************************

// ======================================================
// Conversion
// ======================================================

/**
 * @brief Converts degrees to radians.
 * 
 * @param degrees Angle in degrees.
 * @return Angle in radians.
 */
double fossil_math_trig_deg_to_rad(double degrees);

/**
 * @brief Converts radians to degrees.
 * 
 * @param radians Angle in radians.
 * @return Angle in degrees.
 */
double fossil_math_trig_rad_to_deg(double radians);

// ======================================================
// Basic trig functions
// ======================================================

/**
 * @brief Computes the sine of an angle (in radians).
 * 
 * @param x Angle in radians.
 * @return Sine of the angle.
 */
double fossil_math_trig_sin(double x);

/**
 * @brief Computes the cosine of an angle (in radians).
 * 
 * @param x Angle in radians.
 * @return Cosine of the angle.
 */
double fossil_math_trig_cos(double x);

/**
 * @brief Computes the tangent of an angle (in radians).
 * 
 * @param x Angle in radians.
 * @return Tangent of the angle.
 */
double fossil_math_trig_tan(double x);

// Inverse

/**
 * @brief Computes the arcsine (inverse sine) of a value.
 * 
 * @param x Value whose arcsine is to be computed.
 * @return Angle in radians.
 */
double fossil_math_trig_asin(double x);

/**
 * @brief Computes the arccosine (inverse cosine) of a value.
 * 
 * @param x Value whose arccosine is to be computed.
 * @return Angle in radians.
 */
double fossil_math_trig_acos(double x);

/**
 * @brief Computes the arctangent (inverse tangent) of a value.
 * 
 * @param x Value whose arctangent is to be computed.
 * @return Angle in radians.
 */
double fossil_math_trig_atan(double x);

/**
 * @brief Computes the arctangent of y/x using the signs of the arguments to determine the correct quadrant.
 * 
 * @param y Ordinate value.
 * @param x Abscissa value.
 * @return Angle in radians.
 */
double fossil_math_trig_atan2(double y, double x);

// ======================================================
// Hyperbolic
// ======================================================

/**
 * @brief Computes the hyperbolic sine of a value.
 * 
 * @param x Value whose hyperbolic sine is to be computed.
 * @return Hyperbolic sine of the value.
 */
double fossil_math_trig_sinh(double x);

/**
 * @brief Computes the hyperbolic cosine of a value.
 * 
 * @param x Value whose hyperbolic cosine is to be computed.
 * @return Hyperbolic cosine of the value.
 */
double fossil_math_trig_cosh(double x);

/**
 * @brief Computes the hyperbolic tangent of a value.
 * 
 * @param x Value whose hyperbolic tangent is to be computed.
 * @return Hyperbolic tangent of the value.
 */
double fossil_math_trig_tanh(double x);

// Inverse hyperbolic

/**
 * @brief Computes the inverse hyperbolic sine of a value.
 * 
 * @param x Value whose inverse hyperbolic sine is to be computed.
 * @return Inverse hyperbolic sine of the value.
 */
double fossil_math_trig_asinh(double x);

/**
 * @brief Computes the inverse hyperbolic cosine of a value.
 * 
 * @param x Value whose inverse hyperbolic cosine is to be computed.
 * @return Inverse hyperbolic cosine of the value.
 */
double fossil_math_trig_acosh(double x);

/**
 * @brief Computes the inverse hyperbolic tangent of a value.
 * 
 * @param x Value whose inverse hyperbolic tangent is to be computed.
 * @return Inverse hyperbolic tangent of the value.
 */
double fossil_math_trig_atanh(double x);

#ifdef __cplusplus
}
#include <stdexcept>
#include <vector>
#include <string>

namespace fossil {

namespace math {

    /**
     * @class Trigonometry
     * @brief Provides static methods for trigonometric and hyperbolic operations.
     *
     * This class offers a C++ interface to the underlying Fossil Logic C trigonometric API.
     * All methods are static and directly wrap the corresponding C functions.
     */
    class Trigonometry {
    public:
        // ======================================================
        // Conversion
        // ======================================================

        /**
         * @brief Converts degrees to radians.
         * @param degrees Angle in degrees.
         * @return Angle in radians.
         */
        static double deg_to_rad(double degrees) {
            return fossil_math_trig_deg_to_rad(degrees);
        }

        /**
         * @brief Converts radians to degrees.
         * @param radians Angle in radians.
         * @return Angle in degrees.
         */
        static double rad_to_deg(double radians) {
            return fossil_math_trig_rad_to_deg(radians);
        }

        // ======================================================
        // Basic trig functions
        // ======================================================

        /**
         * @brief Computes the sine of an angle (in radians).
         * @param x Angle in radians.
         * @return Sine of the angle.
         */
        static double sin(double x) {
            return fossil_math_trig_sin(x);
        }

        /**
         * @brief Computes the cosine of an angle (in radians).
         * @param x Angle in radians.
         * @return Cosine of the angle.
         */
        static double cos(double x) {
            return fossil_math_trig_cos(x);
        }

        /**
         * @brief Computes the tangent of an angle (in radians).
         * @param x Angle in radians.
         * @return Tangent of the angle.
         */
        static double tan(double x) {
            return fossil_math_trig_tan(x);
        }

        // Inverse

        /**
         * @brief Computes the arcsine (inverse sine) of a value.
         * @param x Value whose arcsine is to be computed.
         * @return Angle in radians.
         */
        static double asin(double x) {
            return fossil_math_trig_asin(x);
        }

        /**
         * @brief Computes the arccosine (inverse cosine) of a value.
         * @param x Value whose arccosine is to be computed.
         * @return Angle in radians.
         */
        static double acos(double x) {
            return fossil_math_trig_acos(x);
        }

        /**
         * @brief Computes the arctangent (inverse tangent) of a value.
         * @param x Value whose arctangent is to be computed.
         * @return Angle in radians.
         */
        static double atan(double x) {
            return fossil_math_trig_atan(x);
        }

        /**
         * @brief Computes the arctangent of y/x using the signs of the arguments to determine the correct quadrant.
         * @param y Ordinate value.
         * @param x Abscissa value.
         * @return Angle in radians.
         */
        static double atan2(double y, double x) {
            return fossil_math_trig_atan2(y, x);
        }

        // ======================================================
        // Hyperbolic
        // ======================================================

        /**
         * @brief Computes the hyperbolic sine of a value.
         * @param x Value whose hyperbolic sine is to be computed.
         * @return Hyperbolic sine of the value.
         */
        static double sinh(double x) {
            return fossil_math_trig_sinh(x);
        }

        /**
         * @brief Computes the hyperbolic cosine of a value.
         * @param x Value whose hyperbolic cosine is to be computed.
         * @return Hyperbolic cosine of the value.
         */
        static double cosh(double x) {
            return fossil_math_trig_cosh(x);
        }

        /**
         * @brief Computes the hyperbolic tangent of a value.
         * @param x Value whose hyperbolic tangent is to be computed.
         * @return Hyperbolic tangent of the value.
         */
        static double tanh(double x) {
            return fossil_math_trig_tanh(x);
        }

        // Inverse hyperbolic

        /**
         * @brief Computes the inverse hyperbolic sine of a value.
         * @param x Value whose inverse hyperbolic sine is to be computed.
         * @return Inverse hyperbolic sine of the value.
         */
        static double asinh(double x) {
            return fossil_math_trig_asinh(x);
        }

        /**
         * @brief Computes the inverse hyperbolic cosine of a value.
         * @param x Value whose inverse hyperbolic cosine is to be computed.
         * @return Inverse hyperbolic cosine of the value.
         */
        static double acosh(double x) {
            return fossil_math_trig_acosh(x);
        }

        /**
         * @brief Computes the inverse hyperbolic tangent of a value.
         * @param x Value whose inverse hyperbolic tangent is to be computed.
         * @return Inverse hyperbolic tangent of the value.
         */
        static double atanh(double x) {
            return fossil_math_trig_atanh(x);
        }
    };

} // namespace math

} // namespace fossil

#endif

#endif /* FOSSIL_MATH_TRIG_H */

#include <fossil/math/trig.h>
#include <stdio.h>

int main() {
    double deg = 45.0;
    double rad = fossil_math_trig_deg_to_rad(deg);

    printf("%f degrees = %f radians\n", deg, rad);

    double s = fossil_math_trig_sin(rad);
    double c = fossil_math_trig_cos(rad);
    double t = fossil_math_trig_tan(rad);

    printf("sin: %f, cos: %f, tan: %f\n", s, c, t);

    double angle = fossil_math_trig_atan2(1.0, 1.0);
    printf("atan2(1,1) = %f radians\n", angle);

    double hs = fossil_math_trig_sinh(1.0);
    double hc = fossil_math_trig_cosh(1.0);
    double ht = fossil_math_trig_tanh(1.0);

    printf("sinh: %f, cosh: %f, tanh: %f\n", hs, hc, ht);

    return 0;
}

#include <fossil/math/trig.h>
#include <iostream>

int main() {
    using namespace fossil::math;

    double deg = 60.0;
    double rad = Trigonometry::deg_to_rad(deg);
    std::cout << deg << " degrees = " << rad << " radians\n";

    std::cout << "sin: " << Trigonometry::sin(rad)
              << ", cos: " << Trigonometry::cos(rad)
              << ", tan: " << Trigonometry::tan(rad) << "\n";

    double angle = Trigonometry::atan2(1.0, 1.0);
    std::cout << "atan2(1,1) = " << angle << " radians\n";

    std::cout << "sinh: " << Trigonometry::sinh(1.0)
              << ", cosh: " << Trigonometry::cosh(1.0)
              << ", tanh: " << Trigonometry::tanh(1.0) << "\n";

    double asin_val = Trigonometry::asin(0.5);
    std::cout << "asin(0.5) = " << asin_val << " radians\n";

    return 0;
}