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,