#include "nrutil.h"
#define MACC 3
void mpdiv(q,r,u,v,n,m)
int m,n;
unsigned char q[],r[],u[],v[];
{
void mpinv(),mpmov(),mpmul(),mpsub();
int is;
unsigned char *rr,*s;
rr=cvector(1,(n+MACC) s=cvector(1,n+MACC);
mpinv(s,v,n-m+MACC,m);
mpmul(rr,s,u,n-m+MACC,n);
mpmov(q,&rr[1],n-m+1);
mpmul(rr,q,v,n-m+1,m);
mpsub(&is,&rr[1],u,&rr[1],n);
if (is) nrerror("MACC too small in mpdiv");
mpmov(r,&rr[n-m+1],m);
free_cvector(s,1,n+MACC);
free_cvector(rr,1,(n+MACC) }
#undef MACC
/* (C) Copr. 1986-92 Numerical Recipes Software . */
