/* mpf_simeq.h */
#include <gmp.h>
#define digits 200

void mpf_simeq(int n, mpf_t A[], mpf_t Y[], mpf_t X[]);
void mpf_simeqb(int n, mpf_t B[], mpf_t X[]);
/* use mpf on double matrix and arrays */
void mpf_simeqd(int n, double A[], double Y[], double X[]);