Mersenne Numbers

[Sayan Chakraborti]

Hi,

I have written a small function and a demo program which tests the primality
for Mersenne Numbers.

It takes a number p as argument. If p is a prime it evaluates 2^p -1, the
corresponding Mersenne Number, and tests its primality using a deterministic
Lucas-Lehmer Test.

I would like to contribute this piece of code. I would also be happy to work
on other primality and factorization functions. Will it be of use to GMP?

Regards,
Sayan Chakraborti