id	summary	reporter	owner	description	type	status	milestone	component	version	severity	resolution	keywords	cc
8423	miller_rabin_test core dump	robbieye@…	John Maddock	"When running the following code in Ubuntu 12.04 (gcc 4.6) with Boost 1.53.0, the program had a core-dump in miller_rabin_test().

The program was doing primality test on 2^(2^5) + 1 when it crashed. Values such as 2^(2^4) + 1 or 2^(2^6) + 1 do not cause such problems.

Thanks in advance for looking into it.

-Rob

------------------------------------------------------
#include <iostream>

#include <boost/multiprecision/cpp_int.hpp>
#include <boost/multiprecision/miller_rabin.hpp>
#include <boost/multiprecision/random.hpp>

namespace mp = boost::multiprecision;

typedef mp::cpp_int BIGNUM;

int
main(int argc, char* argv[])
{
    // Test primality of 2^(2^5) + 1
    BIGNUM num = (BIGNUM(1) << (1u << 5)) + 1;
    boost::random::mt19937 prng(clock());

    std::cout << ""Running primality test for : "" << num; 
    bool isPrime = mp::miller_rabin_test(num, 25, prng);

    std::cout << "", result : ""
              << (isPrime ? """" : ""not "")
              << ""prime.""
              << std::endl;

    return 0;
}"	Bugs	closed	To Be Determined	multiprecision	Boost 1.53.0	Problem	worksforme	miller rabin primality