New function mpn_mulmod_bnm1

[bodrato at mail.dm.unipi.it]

Hi,

> The interface is mpn_mulmod_bnm1 (rp, ap, bp, n', tp) where ap and bp
> are pointers to the input operands, and rp should point to the result.
> The common size is in n'.  Scratch space is passed in tp, its size
> should be mpn_mulmod_bnm1_itch(n').

I changed the interface, it is now:
mpn_mulmod_bnm1 (rp, rn, ap, an, bp, bn, tp) where an, and bn are the size
of the operands. rn is both the size available at rp and the exponent: the
product is computed modulo B^rn-1.

We miss a test program... maybe the code is not correct...

Regards,
Marco

-- 
http://bodrato.it/papers/