eval.h

Go to the documentation of this file.

//
//  eval.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 eval_h
#define eval_h
 
#include <mpfr.h>
 
#include "ntypes.h"
#include "comp.h"
#include "concave.h"
#include "poly.h"
#include "polyr.h"
#include "pows.h"
#include "powsr.h"
 
// MARK: constants &data types
 
#define MIN_EVAL_PREC      8
#define MIN_EVAL_PREC_STR "8"
 
typedef struct {
    poly P;         
    polyr Q;        
    bool real;      
    prec_t prec;    
    concave f;      
    pows zn;        
    powsr xn;       
    real_t valErr;  
    real_t derErr;  
    real_t ntErr;   
    deg_t terms;    
    comp buf;       
    mpfr_t br;      
} eval_struct;
 
typedef eval_struct eval_t[1];
 
typedef eval_struct *eval;
 
// MARK: macros & functions prototypes
 
#define eval_lastValError(ev) (ev->valErr)
 
#define eval_lastDerError(ev) (ev->derErr)
 
eval eval_new(poly P, prec_t prec);
 
eval eval_new_r(polyr Q, prec_t prec);
 
bool eval_free(eval ev);
 
bool eval_val(comp v, eval ev, comp z);
 
bool eval_der(comp d, eval ev, comp z);
 
bool eval_val_der(comp v, comp d, eval ev, comp z);
 
bool eval_newton(comp nt, eval ev, comp z);
 
// MARK: Bare-metal variants
 
bool eval_analyse(eval ev);
 
bool eval_analyse_r(eval ev);
 
bool eval_val_cc(comp v, eval ev, comp z);
 
bool eval_der_cc(comp d, eval ev, comp z);
 
bool eval_val_der_cc(comp v, comp d, eval ev, comp z);
 
bool eval_newton_cc(comp nt, eval ev, comp z);
 
bool eval_val_cr(comp v, eval ev, mpfr_t x);
 
bool eval_der_cr(comp d, eval ev, mpfr_t x);
 
bool eval_val_der_cr(comp v, comp d, eval ev, mpfr_t x);
 
bool eval_newton_cr(comp nt, eval ev, mpfr_t x);
 
bool eval_val_rc(comp v, eval ev, comp z);
 
bool eval_der_rc(comp d, eval ev, comp z);
 
bool eval_val_der_rc(comp v, comp d, eval ev, comp z);
 
bool eval_newton_rc(comp nt, eval ev, comp z);
 
bool eval_val_rr(mpfr_t v, eval ev, mpfr_t x);
 
bool eval_der_rr(mpfr_t d, eval ev, mpfr_t x);
 
bool eval_val_der_rr(mpfr_t v, mpfr_t d, eval ev, mpfr_t x);
 
bool eval_newton_rr(mpfr_t nt, eval ev, mpfr_t x);
 
#endif /* eval_h */