/******************************************************************
 *  Interface to translate Fortran calls to MPUTIL routines
 *    into calls to GMP
 *****************************************************************/

#include <stdio.h>
#include <string.h>

#ifdef _CRAY


#include "gmp.h"
#include <fortran.h>

/* CVS info   */
/* $Date: 2005/01/10 20:46:01 $     */
/* $Revision: 1.2 $ */
/* $RCSfile: wrap.c,v $  */
/* $Name: rel_5 $     */

static char cvs_info[] = "BMARKGRP $Date: 2005/01/10 20:46:01 $ $Revision: 1.2 $ $RCSfile: wrap.c,v $ $Name: rel_5 $";

MPADD ( MP_INT *sum, MP_INT *u, MP_INT *v )
{
	mpz_add ( sum, u, v );
}

MPMULT ( MP_INT *prod, MP_INT *u, MP_INT *v )
{
	mpz_mul ( prod, u, v );
}

MPMOD ( MP_INT *rem, MP_INT *dividend, MP_INT *modulus )
{
	mpz_mod ( rem, dividend, modulus );
}

MPSET ( MP_INT *dest, MP_INT *source )
{
	mpz_set ( dest, source );
}
 
MPFROMSTRING ( MP_INT *X, _fcd fstring )
{
	char *cstring;
        int i, length;

	cstring = _fcdtocp(fstring);
        length = strlen(cstring);
        for (i=0;i<length;i++)
                if (cstring[i]==' ') { cstring[i]='\0'; length=i; break;}
 	mpz_set_str ( X, cstring, 10 );
        return;
}

MPREAD ( MP_INT *X )
{
	mpz_inp_str ( X, stdin, 10 );
        return;
}

MPWRITE ( MP_INT *X )
{
	mpz_out_str ( stdout, 10, X );
}
#endif

#ifdef __alpha

/*** Alpha version *********************************/
#include "gmp.h"

mpadd_ ( MP_INT *sum, MP_INT *u, MP_INT *v )
{
	mpz_add ( sum, u, v );
}

mpmult_ ( MP_INT *prod, MP_INT *u, MP_INT *v )
{
	mpz_mul ( prod, u, v );
}

mpmod_ ( MP_INT *rem, MP_INT *dividend, MP_INT *modulus )
{
	mpz_mod ( rem, dividend, modulus );
}

mpset_ ( MP_INT *dest, MP_INT *source )
{
	mpz_set ( dest, source );
}
 
mpfromstring_ ( MP_INT *X, char *fstring )
{
	int i, length;
	length = strlen(fstring);
	for (i=0; i<length; i++)
		if (fstring[i]==' ') { fstring[i]='\0'; length=i; break;}
	mpz_set_str ( X, fstring, 10 );
	return length;
}

mpread_ ( MP_INT *X )
{
	mpz_inp_str ( X, stdin, 10 );
}

mpwrite_ ( MP_INT *X )
{
#if 1
	mpz_out_str ( stdout, 10, X );
	printf("\n");
	fflush( stdout );
#else
	printf("%d\n", X[3]);
	printf("%016X  %016X  %016X  %016X \n",
		X[0], X[1], X[2], X[3]);
#endif
}

#endif