From c0f59edac697e7e6b33ac16076ba4e46d176e12a Mon Sep 17 00:00:00 2001
From: Mike Lynch <mike@mikelynch.org>
Date: Fri, 28 Jul 2023 07:18:27 +1000
Subject: [PATCH] Added 600-cell

---
 main.js      |  7 ++--
 polytopes.js | 96 +++++++++++++++++++++++++++++++++-------------------
 2 files changed, 64 insertions(+), 39 deletions(-)

diff --git a/main.js b/main.js
index a7c9015..14f46cc 100644
--- a/main.js
+++ b/main.js
@@ -157,10 +157,10 @@ const STRUCTURES = {
 	'16-cell': POLYTOPES.cell16(),
 	'tesseract': POLYTOPES.tesseract(),
 	'24-cell': POLYTOPES.cell24(),
-	'120-cell': POLYTOPES.cell120()
+	'120-cell': POLYTOPES.cell120(),
+	'600-cell': POLYTOPES.cell600()
 };
 
-console.log(STRUCTURES);
 
 let shape = false;
 
@@ -168,7 +168,6 @@ function createShape(name) {
 	if( shape ) {
 		scene.remove(shape);
 	}
-	console.log(STRUCTURES[name]);
 	shape = new FourDShape(node_ms, link_ms, STRUCTURES[name]);
 	scene.add(shape);
 
@@ -220,7 +219,7 @@ const gui_params = {
 };
 
 gui.add(gui_params, 'shape',
-	[ '5-cell', '16-cell', 'tesseract', '24-cell', '120-cell' ]
+	[ '5-cell', '16-cell', 'tesseract', '24-cell', '120-cell', '600-cell' ]
 	).onChange(createShape)
 
 gui.add(gui_params, 'hyperplane', 1.5, 4);
diff --git a/polytopes.js b/polytopes.js
index a2d3899..3cd91e4 100644
--- a/polytopes.js
+++ b/polytopes.js
@@ -244,6 +244,41 @@ export const cell24 = () => {
 }
 
 
+// above ^^ fix me convert to new style
+
+function dist2(n1, n2) {
+	return (n1.x - n2.x) ** 2 + (n1.y - n2.y) ** 2 + (n1.z - n2.z) ** 2 + (n1.w - n2.w) ** 2;
+}
+
+
+
+function auto_detect_edges(nodes, neighbours) {
+	const seen = {};
+	const nnodes = nodes.length;
+	const links = [];
+	let id = 1;
+	for( const n1 of nodes ) {
+		const d2 = [];
+		for( const n2 of nodes ) {
+			d2.push({ d2: dist2(n1, n2), id: n2.id });
+		}
+		d2.sort((a, b) => a.d2 - b.d2);
+		const closest = d2.slice(1, neighbours + 1);
+		for( const e of closest ) {
+			const ids = [ n1.id, e.id ];
+			ids.sort();
+			const fp = ids.join(',');
+			if( !seen[fp] ) {
+				seen[fp] = true;
+				links.push({ id: id, label: 0, source: n1.id, target: e.id });
+				id++;
+			}
+		}
+	}
+	return links;
+}
+
+
 function make_120cell_vertices() {
 	const phi = 0.5 * (1 + Math.sqrt(5));  
 	const r5 = Math.sqrt(5);   
@@ -265,42 +300,9 @@ function make_120cell_vertices() {
 }
 
 
-function dist2(n1, n2) {
-	return (n1.x - n2.x) ** 2 + (n1.y - n2.y) ** 2 + (n1.z - n2.z) ** 2 + (n1.w - n2.w) ** 2;
-}
-
-
-function make_120cell_edges(nodes) {
-	const seen = {};
-	const nnodes = nodes.length;
-	const links = [];
-	let id = 1;
-	for( const n1 of nodes ) {
-		console.log(n1);
-		const d2 = [];
-		for( const n2 of nodes ) {
-			d2.push({ d2: dist2(n1, n2), id: n2.id });
-		}
-		d2.sort((a, b) => a.d2 - b.d2);
-		const closest = d2.slice(1, 5);
-		for( const e of closest ) {
-			const ids = [ n1.id, e.id ];
-			ids.sort();
-			const fp = ids.join(',');
-			if( !seen[fp] ) {
-				seen[fp] = true;
-				links.push({ id: id, label: 0, source: n1.id, target: e.id });
-				id++;
-			}
-		}
-	}
-	return links;
-}
-
-
 export const cell120 = () => {
 	const nodes  = make_120cell_vertices();
-	const links = make_120cell_edges(nodes);
+	const links = auto_detect_edges(nodes, 4);
 	return {
 		nodes: nodes,
 		links: links
@@ -308,4 +310,28 @@ export const cell120 = () => {
 }
 
 
-const nodes = cell120();
\ No newline at end of file
+function make_600cell_vertices() {
+	const phi = 0.5 * (1 + Math.sqrt(5));  
+
+	const nodes = [
+		PERMUTE.coordinates([0, 0, 0, 2],  0),
+		PERMUTE.coordinates([1, 1, 1, 1], 1),
+
+		PERMUTE.coordinates([phi, 1, 1 / phi, 0], 1, true)
+		].flat();
+	return scale_and_index(nodes, 0.75);
+}
+
+
+export const cell600 = () => {
+	const nodes  = make_600cell_vertices();
+	const links = auto_detect_edges(nodes, 20);
+	return {
+		nodes: nodes,
+		links: links
+	}
+}
+
+
+
+