Boost C++ Libraries: Ticket #8423: miller_rabin_test core dump

anonymous — Wed, 10 Apr 2013 07:19:51 GMT

Here's the properly formatted code. Please ignore the malformed snippet above.

#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 for 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;
}

status changed; resolution set

John Maddock — Wed, 10 Apr 2013 12:20:11 GMT

status new → closed
resolution → worksforme

Confirmed the issue was in 1.53, also that it's already been fixed for the next release.

stephan@… — Fri, 02 Jan 2015 17:08:32 GMT

On my 64-bit machine, it crashed for 2⁽²6)+1. I can confirm it's fixed in 1.54. Thanks for documenting the issue (saved me some time) and for fixing it!