/****************************************************************************** * * * 64x32 MULTIPLY (UNSIGNED) * * 11/14/06 (dkc) * * * ******************************************************************************/ unsigned int carry(unsigned int a, unsigned int b, unsigned int sum); void bigprod(unsigned int a0, unsigned int a2, unsigned int m0, unsigned int *product) { unsigned int a1,a3,m1,temp; unsigned int p0,p1,p2,p3,p4,p5,p6,p7; unsigned int s1,s2,s3; unsigned int c2,c3; a1=a0&0xffff; a0=a0>>16; a3=a2&0xffff; a2=a2>>16; m1=m0&0xffff; m0=m0>>16; p0=a0*m0; p1=a0*m1; p2=a1*m0; p3=a1*m1; p4=a2*m0; p5=a2*m1; p6=a3*m0; p7=a3*m1; s3=p7+(p6<<16); c3=carry(p7,(p6<<16),s3); temp=s3+(p5<<16); c3+=carry(s3,(p5<<16),temp); s3=temp; s2=p4+(p6>>16); c2=carry(p4,(p6>>16),s2); temp=s2+(p5>>16); c2+=carry(s2,(p5>>16),temp); s2=temp; temp=s2+p3; c2+=carry(s2,p3,temp); s2=temp; temp=s2+(p2<<16); c2+=carry(s2,(p2<<16),temp); s2=temp; temp=s2+(p1<<16); c2+=carry(s2,(p1<<16),temp); s2=temp; s1=p0+(p2>>16); s1=s1+(p1>>16); temp=s2+c3; c2+=carry(s2,c3,temp); s2=temp; s1=s1+c2; *product=s1; *(product+1)=s2; *(product+2)=s3; return; }