// nextprobprime().

// General includes.
#include "cl_sysdep.h"

// Specification.
#include "cln/numtheory.h"


// Implementation.

#include "cln/real.h"

namespace cln {

const cl_I nextprobprime (const cl_R& x)
{
	if (minusp(x))
		return 2;
	// Now x >= 0.
	var cl_I n = ceiling1(x);
	if (n <= 2)
		return 2;
	// Now n>=3.
	if (evenp(n))
		n = n+1;
	// Now n>=3 odd.
	loop {
		if (isprobprime(n))
			return n;
		n = n+2;
	}
}

}  // namespace cln


syntax highlighted by Code2HTML, v. 0.9.1