f256 0.11.2

Octuple-precision floating-point arithmetic.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

// ---------------------------------------------------------------------------
// Copyright:   (c) 2024 ff. Michael Amrhein (michael@adrhinum.de)
// License:     This program is part of a larger application. For license
//              details please read the file LICENSE.TXT provided together
//              with the application.
// ---------------------------------------------------------------------------
// $Source$
// $Revision$

use super::{
    approx_cos::{approx_cos, SMALL_CUT_OFF},
    approx_sin::approx_sin,
    FP492,
};

pub(crate) fn approx_sin_cos(x: &FP492) -> (FP492, FP492) {
    let mut x_abs = *x;
    x_abs.iabs();
    // If x is zero or very small, cosine x == 1 and sine x == x.
    if x_abs <= SMALL_CUT_OFF {
        return (*x, FP492::ONE);
    };
    (approx_sin(x), approx_cos(x))
}