Add primality testing and next_prime()/prev_prime()

1 files changed, 29 insertions, 0 deletions

diff --git a/test/integers.tm b/test/integers.tm
index 11603b22..5743cbb2 100644
--- a/test/integers.tm
+++ b/test/integers.tm
@@ -82,3 +82,32 @@ func main():
 			//! n=$n, d=$d:
 			>> (n/d)*d + (n mod d) == n
 			= yes
+
+	>> 0:next_prime()
+	= 2
+	>> 7:next_prime()
+	= 11
+	>> 11:prev_prime()
+	= 7
+	>> (and) p:is_prime() for p in [
+		2, 3, 5, 7,
+		137372146048179869781170214707,
+		811418847921670560768224995279,
+		292590241572454328697048860273,
+		754893741683930091960170890717,
+		319651808258437169510475301537,
+		323890224935694708770556249787,
+		507626552342376235511933571091,
+		548605069630614185274710840981,
+		121475876690852432982324195553,
+		771958616175795150904761471637,
+	]
+	= yes
+
+	>> (or) p:is_prime() for p in [
+		-1, 0, 1, 4, 6,
+		137372146048179869781170214707*2,
+		811418847921670560768224995279*3,
+		292590241572454328697048860273*754893741683930091960170890717,
+	]
+	= no