123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102 |
- ul=(unsigned long)a*b; \
- if (ul != 0) \
- { \
- r=(ul&0xffff)-(ul>>16); \
- r-=((r)>>16); \
- } \
- else \
- r=(-(int)a-b+1);
- c+=n; \
- l1=l2=0; \
- switch (n) { \
- case 8: l2 =((unsigned long)(*(--(c)))) ; \
- \
- case 7: l2|=((unsigned long)(*(--(c))))<< 8; \
- \
- case 6: l2|=((unsigned long)(*(--(c))))<<16; \
- \
- case 5: l2|=((unsigned long)(*(--(c))))<<24; \
- \
- case 4: l1 =((unsigned long)(*(--(c)))) ; \
- \
- case 3: l1|=((unsigned long)(*(--(c))))<< 8; \
- \
- case 2: l1|=((unsigned long)(*(--(c))))<<16; \
- \
- case 1: l1|=((unsigned long)(*(--(c))))<<24; \
- } \
- }
- c+=n; \
- switch (n) { \
- case 8: *(--(c))=(unsigned char)(((l2) )&0xff); \
- \
- case 7: *(--(c))=(unsigned char)(((l2)>> 8)&0xff); \
- \
- case 6: *(--(c))=(unsigned char)(((l2)>>16)&0xff); \
- \
- case 5: *(--(c))=(unsigned char)(((l2)>>24)&0xff); \
- \
- case 4: *(--(c))=(unsigned char)(((l1) )&0xff); \
- \
- case 3: *(--(c))=(unsigned char)(((l1)>> 8)&0xff); \
- \
- case 2: *(--(c))=(unsigned char)(((l1)>>16)&0xff); \
- \
- case 1: *(--(c))=(unsigned char)(((l1)>>24)&0xff); \
- } \
- }
- l|=((unsigned long)(*((c)++)))<<16L, \
- l|=((unsigned long)(*((c)++)))<< 8L, \
- l|=((unsigned long)(*((c)++))))
- *((c)++)=(unsigned char)(((l)>>16L)&0xff), \
- *((c)++)=(unsigned char)(((l)>> 8L)&0xff), \
- *((c)++)=(unsigned char)(((l) )&0xff))
- *((c)++)=(unsigned char)(((l)>> 8L)&0xff))
- l|=((IDEA_INT)(*((c)++))) )
- x1&=0xffff; \
- idea_mul(x1,x1,*p,ul); p++; \
- x2+= *(p++); \
- x3+= *(p++); \
- x4&=0xffff; \
- idea_mul(x4,x4,*p,ul); p++; \
- t0=(x1^x3)&0xffff; \
- idea_mul(t0,t0,*p,ul); p++; \
- t1=(t0+(x2^x4))&0xffff; \
- idea_mul(t1,t1,*p,ul); p++; \
- t0+=t1; \
- x1^=t1; \
- x4^=t0; \
- ul=x2^t0; \
- x2=x3^t1; \
- x3=ul;
|