/* cddmp.h (cddlib arithmetic operations using gmp)
Copyright: Komei Fukuda 2000, fukuda@ifor.math.ethz.ch
Version 0.94f, February 7, 2008
*/
/* This program is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
This program 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 General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program; if not, write to the Free Software
Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
*/
#ifndef __CDDMP_H
#define __CDDMP_H
#endif /* __CDDMP_H */
/**********************************/
/* MACROS */
/* dependent on mp implementation */
/**********************************/
#if defined GMPRATIONAL
#include "gmp.h"
#define dd_ARITHMETIC "GMP rational"
#define dd_init(a) mpq_init(a)
#define dd_clear(a) mpq_clear(a)
#define dd_set(a, b) mpq_set(a,b)
#define dd_set_si(a, b) ddd_mpq_set_si(a,b) /* defined in cddgmp.c */
#define dd_set_si2(a, b, c) mpq_set_si(a,b,c) /* gmp 3.1 or higher */
#define dd_add(a, b, c) mpq_add(a,b,c)
#define dd_sub(a, b, c) mpq_sub(a,b,c)
#define dd_mul(a, b, c) mpq_mul(a,b,c)
#define dd_div(a, b, c) mpq_div(a,b,c)
#define dd_neg(a, b) mpq_neg(a,b)
#define dd_inv(a, b) mpq_inv(a,b)
#define dd_cmp(a, b) mpq_cmp(a,b)
/* returns pos if a>b, 0 if a=b, negative if a<b */
#define dd_sgn(a) mpq_sgn(a)
/* returns nonzero if equal. much faster than mpq_cmp. */
#define dd_get_d(a) mpq_get_d(a)
#elif defined GMPFLOAT
#include "gmp.h"
#define dd_ARITHMETIC "GMP float"
#define dd_init(a) mpf_init(a)
#define dd_clear(a) mpf_clear(a)
#define dd_set(a, b) mpf_set(a,b)
#define dd_set_d(a, b) mpf_set_d(a,b)
#define dd_set_si(a, b) mpf_set_si(a,b)
#define dd_set_si2(a, b, c) mpf_set_si(a,b,c) /* gmp 3.1 or higher */
#define dd_add(a, b, c) mpf_add(a,b,c)
#define dd_sub(a, b, c) mpf_sub(a,b,c)
#define dd_mul(a, b, c) mpf_mul(a,b,c)
#define dd_div(a, b, c) mpf_div(a,b,c)
#define dd_neg(a, b) mpf_neg(a,b)
#define dd_inv(a, b) mpf_inv(a,b)
#define dd_cmp(a, b) mpf_cmp(a,b)
/* returns pos if a>b, 0 if a=b, negative if a<b */
#define dd_sgn(a) mpf_sgn(a)
#define dd_get_d(a) mpf_get_d(a)
#else /* built-in C double */
#define dd_ARITHMETIC "C double"
#define dd_CDOUBLE
#define dd_init(a) ddd_init(a)
#define dd_clear(a) ddd_clear(a)
#define dd_set(a, b) ddd_set(a,b)
#define dd_set_si(a, b) ddd_set_si(a,b)
#define dd_set_si2(a, b, c) ddd_set_si2(a,b,c)
#define dd_set_d(a, b) ddd_set_d(a,b)
#define dd_add(a, b, c) ddd_add(a,b,c)
#define dd_sub(a, b, c) ddd_sub(a,b,c)
#define dd_mul(a, b, c) ddd_mul(a,b,c)
#define dd_div(a, b, c) ddd_div(a,b,c)
#define dd_neg(a, b) ddd_neg(a,b)
#define dd_inv(a, b) ddd_inv(a,b)
#define dd_cmp(a, b) ddd_cmp(a,b)
/* returns pos if a>b, 0 if a=b, negative if a<b */
#define dd_sgn(a) ddd_sgn(a)
#define dd_get_d(a) ddd_get_d(a)
#endif
#if defined GMPRATIONAL
typedef mpq_t mytype;
#elif defined GMPFLOAT
typedef mpf_t mytype;
#else /* built-in C double */
typedef double mytype[1];
#endif
#if defined(__cplusplus)
extern "C" {
#endif
void ddd_mpq_set_si(mytype,signed long);
void ddd_init(mytype);
void ddd_clear(mytype);
void ddd_set(mytype,mytype);
void ddd_set_d(mytype,double);
void ddd_set_si(mytype,signed long);
void ddd_set_si2(mytype,signed long, unsigned long);
void ddd_add(mytype,mytype,mytype);
void ddd_sub(mytype,mytype,mytype);
void ddd_mul(mytype,mytype,mytype);
void ddd_div(mytype,mytype,mytype);
void ddd_neg(mytype,mytype);
void ddd_inv(mytype,mytype);
int ddd_cmp(mytype,mytype);
int ddd_sgn(mytype);
double ddd_get_d(mytype);
void ddd_mpq_set_si(mytype,signed long);
void dd_set_global_constants(void);
void dd_free_global_constants(void); /* 094d */
#if defined(__cplusplus)
}
#endif
/* end of cddmp.h */
