summaryrefslogtreecommitdiff
path: root/mpfi.h
blob: 64ca685184530ca0c5cb52c2d6b67b138043054b (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
/* mpfi.h -- Include file for mpfi.
 
Copyright 1999 PolKA project, Inria Lorraine and Loria

This file is part of the MPFR Library.
 
The MPFR Library is free software; you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as published by
the Free Software Foundation; either version 2.1 of the License, or (at your
option) any later version.
 
The MPFR Library is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
License for more details.
 
You should have received a copy of the GNU Lesser General Public License
along with the MPFR Library; see the file COPYING.LIB.  If not, write to
the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
MA 02111-1307, USA. */

#ifndef _mpfi_h
#define _mpfi_h

#define MPFI_RNDD GMP_RNDD
#define MPFI_RNDU GMP_RNDU

#include <math.h>

/* Machine dependant implementations */

#include "gmp.h"
#include "mpfr.h"

typedef struct {
  __mpfr_struct left;
  __mpfr_struct right;
}__mpfi_struct;

typedef __mpfi_struct mpfi_t[1];
typedef __mpfi_struct *mpfi_ptr;
typedef __gmp_const __mpfi_struct *mpfi_srcptr;

/* Gestion de la precision */

void mpfi_set_prec(const unsigned int);

unsigned int mpfi_get_prec();

/* Gestion de memoire */

void mpfi_clear(mpfi_ptr);

/* initializations */

void   mpfi_init    (mpfi_ptr);
void   mpfi_set     (mpfi_ptr, mpfi_srcptr);
void   mpfi_set_si  (mpfi_ptr, const long);
void   mpfi_set_ui  (mpfi_ptr, const unsigned long);
void   mpfi_set_d   (mpfi_ptr, const double);


/* void   mpfi_set_str (mpfi_ptr, const char *); */

/* arithmetic operations */

void   mpfi_add     (mpfi_ptr, mpfi_srcptr, mpfi_srcptr);
void   mpfi_sub     (mpfi_ptr, mpfi_srcptr, mpfi_srcptr);
void   mpfi_mul     (mpfi_ptr, mpfi_srcptr, mpfi_srcptr);
void   mpfi_div     (mpfi_ptr, mpfi_srcptr, mpfi_srcptr);

void   mpfi_add_d   (mpfi_ptr, mpfi_srcptr, const double);
void   mpfi_sub_d   (mpfi_ptr, mpfi_srcptr, const double);
void   mpfi_mul_d   (mpfi_ptr, mpfi_srcptr, const double);
void   mpfi_div_d   (mpfi_ptr, mpfi_srcptr, const double);

#define mpfi_add_si(a,b,c) mpfi_add_d(a,b,(double) c);
#define mpfi_sub_si(a,b,c) mpfi_sub_d(a,b,(double) c);
#define mpfi_mul_si(a,b,c) mpfi_mul_d(a,b,(double) c);
#define mpfi_div_si(a,b,c) mpfi_div_d(a,b,(double) c);

/* destructive arithmetic operations */

void   mpfi_dadd    (mpfi_ptr, mpfi_srcptr);
void   mpfi_dsub    (mpfi_ptr, mpfi_srcptr);
void   mpfi_dmul    (mpfi_ptr, mpfi_srcptr);
void   mpfi_ddiv    (mpfi_ptr, mpfi_srcptr);

void   mpfi_dadd_d  (mpfi_ptr, const double);
void   mpfi_dsub_d  (mpfi_ptr, const double);
void   mpfi_dmul_d  (mpfi_ptr, const double);
void   mpfi_ddiv_d  (mpfi_ptr, const double);

#define mpfi_dadd_si(a,c) mpfi_dadd_d(a,(double) c);
#define mpfi_dsub_si(a,c) mpfi_dsub_d(a,(double) c);
#define mpfi_dmul_si(a,c) mpfi_dmul_d(a,(double) c);
#define mpfi_ddiv_si(a,c) mpfi_ddiv_d(a,(double) c);

/* destructive composite arithmetic operations */


/* various operations */

void mpfi_neg(mpfi_ptr, mpfi_srcptr);
void mpfi_inv(mpfi_ptr, mpfi_srcptr);
void mpfi_mul_2exp(mpfi_ptr, mpfi_srcptr,unsigned long);

/* interval manipulation */

void   mpfi_set_left  (mpfi_ptr, mpfr_srcptr);
void   mpfi_set_right (mpfi_ptr, mpfr_srcptr);
void   mpfi_get_left  (mpfi_srcptr, mpfr_ptr);
void   mpfi_get_right (mpfi_srcptr, mpfr_ptr);


/* customizable comparison functions */

int    (*mpfi_cmp)  (mpfi_srcptr,mpfi_srcptr);
int    (*mpfi_cmp_d)   (mpfi_srcptr,const double);
int    (*mpfi_cmp_ui)  (mpfi_srcptr,const unsigned long);
int    (*mpfi_cmp_si)  (mpfi_srcptr,const long);

int    (*mpfi_is_pos)  (mpfi_srcptr);
int    (*mpfi_is_neg)  (mpfi_srcptr);
int    (*mpfi_is_zero) (mpfi_srcptr);


/* Customizable error handling */

#ifdef MPFI_USE_CHECK_ERROR
int    reset_error();
int    get_error();
char * get_error_str();
#endif

#endif