compf.h

Go to the documentation of this file.

//
//  compf.h
//
//  Authors: Nicolae Mihalache & François Vigneron
//
//  This software is released under BSD licence, with an attribution clause (see Licence file).
//  Please cite the reference below if you use or distribute this software.
//
//  • [1] R. Anton, N. Mihalache & F. Vigneron. Fast evaluation of real and complex polynomials. 2022.
//        https://hal.archives-ouvertes.fr/hal-03820369
//
//  Copyright 2022 Univ. Paris-Est Créteil & Univ. de Reims Champagne-Ardenne.
//
 
#ifndef compf_h
#define compf_h
 
#include "ntypes.h"
 
// MARK: data types
 
typedef struct {
    coeff_t x;  
    coeff_t y;  
} compf_struct;
 
typedef compf_struct compf[1];
 
typedef compf_struct *compf_ptr;
 
// MARK: macros & functions prototypes
 
#define compf_set(d, a)         (d)->x = (a)->x; \
                                (d)->y = (a)->y;
 
#define compf_setr(d, a)        (d)->x = (a); \
                                (d)->y = 0;
 
#define compf_neg(d, a)         (d)->x = -(a)->x; \
                                (d)->y = -(a)->y;
 
#define compf_add(d, a, b)      (d)->x = (a)->x + (b)->x; \
                                (d)->y = (a)->y + (b)->y;
 
#define compf_addr(d, a, r)     (d)->x = (a)->x + (r); \
                                (d)->y = (a)->y;
 
#define compf_sub(d, a, b)      (d)->x = (a)->x - (b)->x; \
                                (d)->y = (a)->y - (b)->y;
 
#define compf_subr(d, a, r)     (d)->x = (a)->x - (r); \
                                (d)->y = (a)->y;
 
#define compf_mul(d, a, b)      coeff_t px = (a)->x * (b)->x - (a)->y * (b)->y; \
                                (d)->y = (a)->x * (b)->y + (a)->y * (b)->x; \
                                (d)->x = px;
 
#define compf_mulr(d, a, r)     (d)->x = (a)->x * (r); \
                                (d)->y = (a)->y * (r);
 
#define compf_amr(d, a, r)      (d)->x += (a)->x * (r); \
                                (d)->y += (a)->y * (r);
 
#define compf_sqr(d, a)         coeff_t px = (a)->x * (a)->x - (a)->y * (a)->y; \
                                (d)->y = 2 * (a)->x * (a)->y; \
                                (d)->x = px;
 
#define compf_div(d, a, b)      { \
                                coeff_t m2 = compf_mod2(b); \
                                coeff_t px = ((a)->x * (b)->x + (a)->y * (b)->y) / m2; \
                                (d)->y = ((a)->y * (b)->x - (a)->x * (b)->y) / m2; \
                                (d)->x = px;\
                                }
 
#define compf_mod(a)            fhypot((a)->x, (a)->y)
 
#define compf_dist(a, b)        fhypot((a)->x - (b)->x, (a)->y - (b)->y)
 
#define compf_mod2(a)           ((a)->x * (a)->x + (a)->y * (a)->y)
 
#define compf_log2(a)           (plog2(compf_mod(a)))
 
#define coeff_log2(a)           (plog2((a) < 0 ? -(a) : (a)))
 
#define compf_s(a)              (pfloor(plog2(compf_mod(a))) + 1)
 
#define coeff_s(a)              (pfloor(coeff_log2(a)) + 1)
 
#endif /* compf_h */