sin_cos_series

gistools::proj::geodesic

Function sin_cos_series 

Source 
pub fn sin_cos_series(
    sinp: bool,
    sinx: f64,
    cosx: f64,
    c: &[f64],
    n: usize,
) -> f64
Expand description

Evaluate

y = sinp
    ? sum(c[i] * sin( 2*i    * x), i, 1, n)
    : sum(c[i] * cos((2*i+1) * x), i, 0, n-1)

using Clenshaw summation. N.B. c[0] is unused for sin series Approx operation count = (n + 5) mult and (2 * n + 2) add