Function complex_horner_eval

src/transformations/horner.cpp:225–244

}

inline static PJ_UV complex_horner_eval(uint32_t order, const double *c,
                                        PJ_UV en, uint32_t order_offset = 0) {
    // the coefficients are ordered like this:
    // (Cn0+i*Ce0, Cn1+i*Ce1, ...)
    const uint32_t sz = horner_number_of_complex_coefficients(order);
    const double e = en.u;
    const double n = en.v;
    const double *cbeg = c + order_offset * 2;
    c += sz;

    double E = *--c;
    double N = *--c;
    double w;
    while (c > cbeg) {
        w = n * E + e * N + *--c;
        N = n * N - e * E + *--c;
        E = w;
    }
    return {E, N};
}

inline static PJ_UV generate_error_coords() {
    PJ_UV uv_error;