powsr.h

Go to the documentation of this file.

//
//  powsr.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 powsr_h
#define powsr_h
 
#include <mpfr.h>
 
#include "ntypes.h"
 
// MARK: data types
 
// comment this line to use repeated squares method and buffers to compute powers of @c x
// leave uncommented tu use a direct computation with @c mpfr_pow_ui()
#define POWSR_DIRECT_POWERS
 
#ifdef POWSR_DIRECT_POWERS
 
typedef struct {
    prec_t prec;    
    bool inited;    
    mpfr_t x;       
    mpfr_t res;     
} powsr_struct;
 
#else
 
typedef struct {
    prec_t prec;    
    deg_t size;     
    byte tps;       
    bool *computed; 
    mpfr_ptr xn;    
    mpfr_t pth;     
    mpfr_t res;     
} powsr_struct;
 
#endif
 
typedef powsr_struct powsr_t[1];
 
typedef powsr_struct *powsr;
 
// MARK: macros & functions prototypes
 
powsr powsr_new(prec_t prec, deg_t size);
 
bool powsr_free(powsr xn);
 
bool powsr_set(powsr xn, mpfr_t x);
 
mpfr_ptr powsr_pow(powsr xn, deg_t pow);
 
mpfr_ptr powsr_pow_once(powsr xn, deg_t pow);
 
#endif /* powsr_h */