Ignore callgrind files.

Also updating dependencies in Cargo.lock

.gitignore

@@ -1,2 +1,3 @@
 /target
 /indivisible.1.gz
+callgrind.out.*

Cargo.lock

@@ -4,9 +4,9 @@ version = 4
 
 [[package]]
 name = "ansi_term"
-version = "0.11.0"
+version = "0.12.1"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "ee49baf6cb617b853aa8d93bf420db2383fab46d314482ca2803b40d5fde979b"
+checksum = "d52a9bb7ec0cf484c551830a7ce27bd20d67eac647e1befb56b0be4ee39a55d2"
 dependencies = [
  "winapi",
 ]
@@ -24,15 +24,15 @@ dependencies = [
 
 [[package]]
 name = "bitflags"
-version = "1.2.1"
+version = "1.3.2"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "cf1de2fe8c75bc145a2f577add951f8134889b4795d47466a54a5c846d691693"
+checksum = "bef38d45163c2f1dde094a7dfd33ccf595c92905c8f8f4fdc18d06fb1037718a"
 
 [[package]]
 name = "clap"
-version = "2.33.3"
+version = "2.34.0"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "37e58ac78573c40708d45522f0d80fa2f01cc4f9b4e2bf749807255454312002"
+checksum = "a0610544180c38b88101fecf2dd634b174a62eef6946f84dfc6a7127512b381c"
 dependencies = [
  "ansi_term",
  "atty",
@@ -45,18 +45,18 @@ dependencies = [
 
 [[package]]
 name = "heck"
-version = "0.3.2"
+version = "0.3.3"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "87cbf45460356b7deeb5e3415b5563308c0a9b057c85e12b06ad551f98d0a6ac"
+checksum = "6d621efb26863f0e9924c6ac577e8275e5e6b77455db64ffa6c65c904e9e132c"
 dependencies = [
  "unicode-segmentation",
 ]
 
 [[package]]
 name = "hermit-abi"
-version = "0.1.18"
+version = "0.1.19"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "322f4de77956e22ed0e5032c359a0f1273f1f7f0d79bfa3b8ffbc730d7fbcc5c"
+checksum = "62b467343b94ba476dcb2500d242dadbb39557df889310ac77c5d99100aaac33"
 dependencies = [
  "libc",
 ]
@@ -70,15 +70,15 @@ dependencies = [
 
 [[package]]
 name = "lazy_static"
-version = "1.4.0"
+version = "1.5.0"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "e2abad23fbc42b3700f2f279844dc832adb2b2eb069b2df918f455c4e18cc646"
+checksum = "bbd2bcb4c963f2ddae06a2efc7e9f3591312473c50c6685e1f298068316e66fe"
 
 [[package]]
 name = "libc"
-version = "0.2.87"
+version = "0.2.178"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "265d751d31d6780a3f956bb5b8022feba2d94eeee5a84ba64f4212eedca42213"
+checksum = "37c93d8daa9d8a012fd8ab92f088405fb202ea0b6ab73ee2482ae66af4f42091"
 
 [[package]]
 name = "proc-macro-error"
@@ -106,18 +106,18 @@ dependencies = [
 
 [[package]]
 name = "proc-macro2"
-version = "1.0.24"
+version = "1.0.103"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "1e0704ee1a7e00d7bb417d0770ea303c1bccbabf0ef1667dae92b5967f5f8a71"
+checksum = "5ee95bc4ef87b8d5ba32e8b7714ccc834865276eab0aed5c9958d00ec45f49e8"
 dependencies = [
- "unicode-xid",
+ "unicode-ident",
 ]
 
 [[package]]
 name = "quote"
-version = "1.0.9"
+version = "1.0.42"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "c3d0b9745dc2debf507c8422de05d7226cc1f0644216dfdfead988f9b1ab32a7"
+checksum = "a338cc41d27e6cc6dce6cefc13a0729dfbb81c262b1f519331575dd80ef3067f"
 dependencies = [
  "proc-macro2",
 ]
@@ -130,9 +130,9 @@ checksum = "8ea5119cdb4c55b55d432abb513a0429384878c15dde60cc77b1c99de1a95a6a"
 
 [[package]]
 name = "structopt"
-version = "0.3.21"
+version = "0.3.26"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "5277acd7ee46e63e5168a80734c9f6ee81b1367a7d8772a2d765df2a3705d28c"
+checksum = "0c6b5c64445ba8094a6ab0c3cd2ad323e07171012d9c98b0b15651daf1787a10"
 dependencies = [
  "clap",
  "lazy_static",
@@ -141,9 +141,9 @@ dependencies = [
 
 [[package]]
 name = "structopt-derive"
-version = "0.4.14"
+version = "0.4.18"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "5ba9cdfda491b814720b6b06e0cac513d922fc407582032e8706e9f137976f90"
+checksum = "dcb5ae327f9cc13b68763b5749770cb9e048a99bd9dfdfa58d0cf05d5f64afe0"
 dependencies = [
  "heck",
  "proc-macro-error",
@@ -154,13 +154,13 @@ dependencies = [
 
 [[package]]
 name = "syn"
-version = "1.0.60"
+version = "1.0.109"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "c700597eca8a5a762beb35753ef6b94df201c81cca676604f547495a0d7f0081"
+checksum = "72b64191b275b66ffe2469e8af2c1cfe3bafa67b529ead792a6d0160888b4237"
 dependencies = [
  "proc-macro2",
  "quote",
- "unicode-xid",
+ "unicode-ident",
 ]
 
 [[package]]
@@ -173,22 +173,22 @@ dependencies = [
 ]
 
 [[package]]
-name = "unicode-segmentation"
-version = "1.7.1"
+name = "unicode-ident"
+version = "1.0.22"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "bb0d2e7be6ae3a5fa87eed5fb451aff96f2573d2694942e40543ae0bbe19c796"
+checksum = "9312f7c4f6ff9069b165498234ce8be658059c6728633667c526e27dc2cf1df5"
+
+[[package]]
+name = "unicode-segmentation"
+version = "1.12.0"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "f6ccf251212114b54433ec949fd6a7841275f9ada20dddd2f29e9ceea4501493"
 
 [[package]]
 name = "unicode-width"
-version = "0.1.8"
+version = "0.1.14"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "9337591893a19b88d8d87f2cec1e73fad5cdfd10e5a6f349f498ad6ea2ffb1e3"
-
-[[package]]
-name = "unicode-xid"
-version = "0.2.1"
-source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "f7fe0bb3479651439c9112f72b6c505038574c9fbb575ed1bf3b797fa39dd564"
+checksum = "7dd6e30e90baa6f72411720665d41d89b9a3d039dc45b8faea1ddd07f617f6af"
 
 [[package]]
 name = "vec_map"
@@ -198,9 +198,9 @@ checksum = "f1bddf1187be692e79c5ffeab891132dfb0f236ed36a43c7ed39f1165ee20191"
 
 [[package]]
 name = "version_check"
-version = "0.9.2"
+version = "0.9.5"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "b5a972e5669d67ba988ce3dc826706fb0a8b01471c088cb0b6110b805cc36aed"
+checksum = "0b928f33d975fc6ad9f86c8f283853ad26bdd5b10b7f1542aa2fa15e2289105a"
 
 [[package]]
 name = "winapi"

Cargo.toml

@@ -4,11 +4,9 @@ version = "2.0.0"
 authors = ["Nicolás Ortega Froysa <nicolas@ortegas.org>"]
 description = "A prime number generator and tester"
 edition = "2018"
-readme = "README"
+readme = "README.md"
 license = "GPL-3.0"
 keywords = [ "primes", "math" ]
 
-# See more keys and their definitions at https://doc.rust-lang.org/cargo/reference/manifest.html
-
 [dependencies]
-structopt = "0.3.13"
+structopt = "0.3"

README.md

@@ -38,8 +38,7 @@ indivisible -i ./primes 600
 ## Legacy
 
 This project was originally written in C. It can be found in the
-[indivisible-legacy repository](https://gitlab.com/naortega/Indivisible-legacy)
-on my GitLab.
+[indivisible-legacy repository](https://code.ortegas.org/nortega/indivisible-legacy).
 
 ## License
 

benchmark.sh

@@ -0,0 +1,28 @@
+#!/bin/bash
+
+BIN="./target/release/indivisible"
+TRIALS=20
+
+if ! [ -f "$BIN" ]
+then
+	>&2 echo "Release build not available. Please run 'cargo build -r'."
+	exit 1
+fi
+
+if ! command -v calc &>/dev/null
+then
+	>&2 echo "Missing 'calc' program. Please install it for this script."
+	exit 1
+fi
+
+echo "Calculating primes up to 100,000,000"
+TOTAL="0"
+for _ in $(seq "$TRIALS")
+do
+	TIME=$(command time -f "%e" "$BIN" 100000000 2>&1 >/dev/null)
+	TOTAL=$(calc "$TOTAL + $TIME")
+done
+
+AVG=$(calc "$TOTAL / $TRIALS")
+
+echo "Average time: ${AVG}s"

src/candidate.rs

@@ -1,62 +0,0 @@
-/*
- * Copyright (C) 2025  Nicolás Ortega Froysa <nicolas@ortegas.org>
- * Author: Nicolás Ortega Froysa <nicolas@ortegas.org>
- *
- * This program is free software: you can redistribute it and/or modify
- * it under the terms of the GNU General Public License as published by
- * the Free Software Foundation, either version 3 of the License, or
- * (at your option) any later version.
- *
- * This program is distributed in the hope that it will be useful,
- * but WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
- * GNU General Public License for more details.
- *
- * You should have received a copy of the GNU General Public License
- * along with this program.  If not, see <http://www.gnu.org/licenses/>.
- */
-
-pub struct CandidateGenerator {
-	base:u64,
-	first_use:bool,
-}
-
-impl CandidateGenerator {
-	pub fn new() -> CandidateGenerator {
-		CandidateGenerator {
-			base: 0,
-			first_use: true,
-		}
-	}
-
-	pub fn next(&mut self) -> u64 {
-		/*
-		 * All primes, except 2 and 3, will be equal to (n * 6 ± 1). This avoids
-		 * multiples of three, optimizing our counting.
-		 */
-		let val;
-
-		if self.base != 0 {
-			if self.first_use {
-				val = self.base - 1;
-			} else {
-				val = self.base + 1;
-			}
-		} else {
-			if self.first_use {
-				val = 2;
-			} else {
-				val = 3;
-			}
-		}
-
-		if self.first_use {
-			self.first_use = false;
-		} else {
-			self.first_use = true;
-			self.base += 6;
-		}
-
-		val
-	}
-}

src/main.rs

@@ -16,17 +16,12 @@
  * along with this program.  If not, see <http://www.gnu.org/licenses/>.
  */
 
-use std::cell::RefCell;
-use std::collections::VecDeque;
-use std::fs::File;
-use std::io::{BufRead, BufReader};
-use std::path::PathBuf;
+//use std::fs::File;
+//use std::io::{BufRead, BufReader};
+//use std::path::PathBuf;
 use std::process;
-use std::rc::Rc;
 use structopt::StructOpt;
 
-mod candidate;
-use candidate::CandidateGenerator;
 mod worker;
 
 #[derive(StructOpt)]
@@ -34,74 +29,53 @@ mod worker;
 struct Opt {
 	#[structopt(short, long, help = "Print all found primes")]
 	verbose:bool,
-	#[structopt(short, long, name = "FILE", help = "Import prime numbers from FILE")]
-	import:Option<PathBuf>,
+	//#[structopt(short, long, name = "FILE", help = "Import prime numbers from FILE")]
+	//import:Option<PathBuf>,
 	#[structopt(short, long, help = "Test if num is prime instead of generation")]
 	test:bool,
-	#[structopt(help = "Ordinal of the prime to generate or number to test for primality")]
-	num:u64,
-	#[structopt(short, long, name = "n", help = "Number of threads to spawn")]
-	jobs:Option<u64>,
+	#[structopt(help = "Max of the prime to generate or number to test for primality")]
+	num:usize,
+	//#[structopt(short, long, name = "n", default_value = "1", help = "Number of threads to spawn")]
+	//jobs:u64,
 }
 
 fn main() {
 	let opts = Opt::from_args();
 
-	let primes_list = Rc::new(RefCell::new(VecDeque::<u64>::new()));
-
-	if opts.import.is_some() {
+	/*if opts.import.is_some() {
 		let in_file = File::open(opts.import.unwrap()).unwrap();
 		let reader = BufReader::new(in_file);
 		for p in reader.lines().into_iter() {
-			primes_list.borrow_mut().push_back(p.unwrap().parse().unwrap());
+			prime_list.push(p.unwrap().parse().unwrap());
 		}
-	}
+	}*/
 
-	let jobs = match opts.jobs {
-		Some(n) => n,
-		None => 1, // TODO: use number of CPUs
-	};
-
-	if opts.num == 0 {
+	if opts.num < 2 {
 		eprintln!("Invalid value for num: {}", opts.num);
 		process::exit(1);
 	}
 
-	if opts.test && *primes_list.borrow().back().unwrap_or(&0) >= opts.num {
-		for i in primes_list.borrow().iter() {
-			if *i == opts.num {
-				process::exit(0)
+	let prime_list = worker::work_segment(0, opts.num);
+
+	if opts.test {
+		if *prime_list.last().unwrap() == (opts.num as u64) {
+			if opts.verbose {
+				println!("{} is prime", opts.num);
 			}
+			process::exit(0);
+		} else {
+			if opts.verbose {
+				println!("{} is composite", opts.num);
+			}
+			process::exit(1);
 		}
-		process::exit(1)
-	} else if !opts.test && primes_list.borrow().len() >= opts.num as usize {
-		let res = *primes_list.borrow().get(opts.num as usize).unwrap();
-		println!("{}", res);
 	} else {
-		let mut cand_gen = CandidateGenerator::new();
-
-		while primes_list.borrow().len() < opts.num as usize {
-			let cand = cand_gen.next();
-			let mut is_prime = true;
-
-			for i in primes_list.borrow().iter() {
-				if cand % *i == 0 {
-					is_prime = false;
-					break;
-				}
-			}
-
-			if is_prime {
-				primes_list.borrow_mut().push_back(cand);
-				if opts.verbose {
-					println!("{}", cand);
-				}
-			}
-		}
-
 		if !opts.verbose {
-			let last_prime = *primes_list.borrow().back().unwrap();
-			println!("{}", last_prime);
+			println!("{}", prime_list.last().unwrap());
+		} else {
+			for p in prime_list {
+				println!("{}", p);
+			}
 		}
 	}
 }

src/worker.rs

@@ -16,21 +16,68 @@
  * along with this program.  If not, see <http://www.gnu.org/licenses/>.
  */
 
-use std::cell::RefCell;
-use std::collections::VecDeque;
-use std::rc::Rc;
+/**
+ * @brief Work on a segment.
+ *
+ * @param start:usize Beginning of the segment (inclusive).
+ * @param end:usize End of the segment (inclusive).
+ *
+ * @return List of primes found in segment.
+ */
+pub fn work_segment(start:usize, end:usize) -> Vec<u64> {
+	let mut found_primes = Vec::<u64>::new();
+	let mut arr = vec![false; end - start + 1];
 
-pub struct Worker {
-	primes_list:Rc<RefCell<VecDeque<u64>>>,
-}
+	if start < 2 && end > 2 {
+		arr[2] = true;
+	}
+	if start < 3 && end > 3 {
+		arr[3] = true;
+	}
 
-impl Worker {
-	pub fn new(primes_list:Rc<RefCell<VecDeque<u64>>>) -> Worker {
-		Worker {
-			primes_list,
+	let sqrt_of_num = f64::sqrt(end as f64) as usize;
+	for x in 1..=sqrt_of_num {
+		let xx4 = 4 * x * x;
+		let xx3 = 3 * x * x;
+		for y in 1..=sqrt_of_num {
+			let yy = y * y;
+
+			let n1 = xx4 + yy;
+			if n1 <= end && (n1 % 12 == 1 || n1 % 12 == 5) {
+				arr[n1] = !arr[n1];
+			}
+
+			let n2 = xx3 + yy;
+			if n2 <= end && n2 % 12 == 7 {
+				arr[n2] = !arr[n2];
+			}
+
+			if x > y {
+				let n3 = xx3 - yy;
+				if n3 <= end && n3 % 12 == 11 {
+					arr[n3] = !arr[n3];
+				}
+			}
 		}
 	}
 
-	pub fn run(&mut self) {
+	for i in 5..=sqrt_of_num {
+		if !arr[i] {
+			continue;
+		}
+
+		let mut j = i * i;
+		while j <= end {
+			arr[j] = false;
+			j += i * i;
+		}
 	}
+
+	for i in 2..=end {
+		if arr[i] {
+			found_primes.push(i as u64);
+		}
+	}
+
+	found_primes
 }

test.sh

@@ -1,28 +1,50 @@
 #!/bin/bash
 
-echo -n "Prime finding..."
-if [[ $(./target/release/indivisible 5) == 11 ]]
+DEBUG_BIN="./target/debug/indivisible"
+RELEASE_BIN="./target/release/indivisible"
+
+if [ -f "$DEBUG_BIN" ]
 then
-	echo " pass"
+	BINARY="$DEBUG_BIN"
+elif [ -f "$RELEASE_BIN" ]
+then
+	BINARY="$RELEASE_BIN"
 else
-	echo " FAIL"
+	>&2 echo "No valid binary found. Please compile the project."
 	exit 1
 fi
 
-echo -n "Positive prime test..."
-if ./target/release/indivisible -t 11
+tests=0
+passed=0
+
+((tests++))
+echo -n "${tests}: Find 5th prime number..."
+if [[ $("$BINARY" 5) == 11 ]]
 then
 	echo " pass"
+	((passed++))
 else
 	echo " FAIL"
-	exit 1
 fi
 
-echo -n "Negative prime test..."
-if ! ./target/release/indivisible -t 9
+((tests++))
+echo -n "${tests}: 11 is prime..."
+if "$BINARY" -t 11
 then
 	echo " pass"
+	((passed++))
 else
 	echo " FAIL"
-	exit 1
 fi
+
+((tests++))
+echo -n "${tests}: 9 is not prime..."
+if ! "$BINARY" -t 9
+then
+	echo " pass"
+	((passed++))
+else
+	echo " FAIL"
+fi
+
+echo "Results: $passed/$tests"