/* mpf_tsimeq.h   declare 'digits' externally before use */
#include <gmp.h>
void mpf_tsimeq(int np, int n, mpf_t A[], mpf_t Y[], mpf_t X[]);
void mpf_tsimeqb(int np, int n, mpf_t B[], mpf_t X[]);
/* use mpf on double matrix and arrays */
void mpf_tsimeqd(int np, int n, double A[], double Y[], double X[]);