evalf.h

Go to the documentation of this file.

//
//  evalf.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 evalf_h
#define evalf_h
 
#include "ntypes.h"
#include "compf.h"
#include "concave.h"
#include "polyf.h"
#include "polyfr.h"
#include "powsf.h"
#include "powsfr.h"
 
// MARK: data types
 
typedef struct {
    polyf P;        
    polyfr Q;       
    bool real;      
    concave f;      
    powsf zn;       
    powsfr xn;      
    real_t valErr;  
    real_t derErr;  
    real_t ntErr;   
    deg_t terms;    
} evalf_struct;
 
typedef evalf_struct evalf_t[1];
 
typedef evalf_struct *evalf;
 
// MARK: macros & functions prototypes
 
#define evalf_lastValError(ev) (ev->valErr)
 
#define evalf_lastDerError(ev) (ev->derErr)
 
evalf evalf_new(polyf P);
 
evalf evalf_new_r(polyfr Q);
 
bool evalf_free(evalf ev);
 
bool evalf_val(compf v, evalf ev, compf z);
 
bool evalf_der(compf d, evalf ev, compf z);
 
bool evalf_val_der(compf v, compf d, evalf ev, compf z);
 
bool evalf_newton(compf nt, evalf ev, compf z);
 
// MARK: Bare-metal variants
 
bool evalf_analyse(evalf ev);
 
bool evalf_analyse_r(evalf ev);
 
bool evalf_val_cc(compf v, evalf ev, compf z);
 
bool evalf_der_cc(compf d, evalf ev, compf z);
 
bool evalf_val_der_cc(compf v, compf d, evalf ev, compf z);
 
bool evalf_newton_cc(compf nt, evalf ev, compf z);
 
bool evalf_val_cr(compf v, evalf ev, coeff_t x);
 
bool evalf_der_cr(compf d, evalf ev, coeff_t x);
 
bool evalf_val_der_cr(compf v, compf d, evalf ev, coeff_t x);
 
bool evalf_newton_cr(compf nt, evalf ev, coeff_t x);
 
bool evalf_val_rc(compf v, evalf ev, compf z);
 
bool evalf_der_rc(compf d, evalf ev, compf z);
 
bool evalf_val_der_rc(compf v, compf d, evalf ev, compf z);
 
bool evalf_newton_rc(compf nt, evalf ev, compf z);
 
coeff_t evalf_val_rr(evalf ev, coeff_t x);
 
coeff_t evalf_der_rr(evalf ev, coeff_t x);
 
bool evalf_val_der_rr(coeff_t *v, coeff_t *d, evalf ev, coeff_t x);
 
coeff_t evalf_newton_rr(evalf ev, coeff_t x);
 
#endif /* evalf_h */