bravais 0.1.3 → 0.2.0

package/dist/bravais.js

@@ -1,15 +1,15 @@
-import { ACESFilmicToneMapping as e, AlwaysStencilFunc as t, AmbientLight as n, Box3 as r, BoxGeometry as i, BufferGeometry as a, CanvasTexture as o, Color as s, ConeGeometry as c, Controls as l, CylinderGeometry as u, DataTexture as d, DirectionalLight as f, DoubleSide as p, EqualStencilFunc as m, Float32BufferAttribute as h, FloatType as g, Fog as _, Group as v, InstancedMesh as y, LineBasicMaterial as ee, LineSegments as b, LinearFilter as te, LinearMipmapLinearFilter as ne, MOUSE as x, MathUtils as S, Matrix4 as C, Mesh as w, MeshStandardMaterial as T, NearestFilter as E, NotEqualStencilFunc as re, Object3D as ie, OrthographicCamera as ae, PlaneGeometry as D, Quaternion as O, RGBAFormat as oe, ReplaceStencilOp as se, SRGBColorSpace as ce, Scene as k, SphereGeometry as le, Sprite as ue, SpriteMaterial as de, Vector2 as A, Vector3 as j, WebGLRenderer as fe } from "three";
+import { ACESFilmicToneMapping as e, AlwaysStencilFunc as t, AmbientLight as n, Box3 as r, BoxGeometry as i, BufferGeometry as a, CanvasTexture as o, Color as s, ConeGeometry as c, Controls as l, CylinderGeometry as u, DataTexture as d, DirectionalLight as f, DoubleSide as p, EqualStencilFunc as m, Float32BufferAttribute as h, FloatType as g, Fog as ee, Group as _, InstancedMesh as v, LineBasicMaterial as y, LineSegments as b, LinearFilter as x, LinearMipmapLinearFilter as te, MOUSE as S, MathUtils as ne, Matrix4 as C, Mesh as w, MeshStandardMaterial as T, NearestFilter as E, NotEqualStencilFunc as D, Object3D as O, OrthographicCamera as k, PlaneGeometry as re, Quaternion as A, RGBAFormat as ie, ReplaceStencilOp as ae, SRGBColorSpace as oe, Scene as se, SphereGeometry as ce, Sprite as le, SpriteMaterial as ue, Vector2 as j, Vector3 as M, WebGLRenderer as N } from "three";
 //#region src/cif/tokenizer.ts
-var M = class extends Error {
+var P = class extends Error {
 	line;
 	constructor(e, t) {
 		super(`CIF line ${t}: ${e}`), this.name = "CifSyntaxError", this.line = t;
 	}
 };
-function pe(e) {
+function F(e) {
 	return e === " " || e === "	" || e === "\r" || e === "\n";
 }
-function me(e, t) {
+function de(e, t) {
 	let n = e.toLowerCase();
 	return e.startsWith("_") ? {
 		kind: "tag",
@@ -37,7 +37,7 @@ function me(e, t) {
 		line: t
 	};
 }
-function he(e) {
+function fe(e) {
 	let t = [], n = e.split(/\r\n|\r|\n/), r = 0;
 	for (; r < n.length;) {
 		let e = n[r] ?? "", i = r + 1;
@@ -54,7 +54,7 @@ function he(e) {
 				}
 				a.push(e), r += 1;
 			}
-			if (!o) throw new M("unterminated text field", i);
+			if (!o) throw new P("unterminated text field", i);
 			t.push({
 				kind: "value",
 				text: a.join("\n"),
@@ -65,7 +65,7 @@ function he(e) {
 		let a = 0;
 		for (; a < e.length;) {
 			let n = e.charAt(a);
-			if (pe(n)) {
+			if (F(n)) {
 				a += 1;
 				continue;
 			}
@@ -74,15 +74,15 @@ function he(e) {
 				let r = a + 1;
 				for (; r < e.length;) {
 					let t = e.indexOf(n, r);
-					if (t === -1) throw new M("unterminated quoted string", i);
+					if (t === -1) throw new P("unterminated quoted string", i);
 					let a = e.charAt(t + 1);
-					if (t + 1 >= e.length || pe(a)) {
+					if (t + 1 >= e.length || F(a)) {
 						r = t;
 						break;
 					}
 					r = t + 1;
 				}
-				if (r >= e.length) throw new M("unterminated quoted string", i);
+				if (r >= e.length) throw new P("unterminated quoted string", i);
 				t.push({
 					kind: "value",
 					text: e.slice(a + 1, r),
@@ -91,8 +91,8 @@ function he(e) {
 				continue;
 			}
 			let r = a;
-			for (; r < e.length && !pe(e.charAt(r));) r += 1;
-			t.push(me(e.slice(a, r), i)), a = r;
+			for (; r < e.length && !F(e.charAt(r));) r += 1;
+			t.push(de(e.slice(a, r), i)), a = r;
 		}
 		r += 1;
 	}
@@ -100,10 +100,10 @@ function he(e) {
 }
 //#endregion
 //#region src/cif/parser.ts
-function ge(e) {
+function pe(e) {
 	return typeof e == "string" ? [e] : e;
 }
-var _e = class {
+var me = class {
 	name;
 	items = /* @__PURE__ */ new Map();
 	loops = [];
@@ -111,41 +111,41 @@ var _e = class {
 		this.name = e;
 	}
 	get(e) {
-		for (let t of ge(e)) {
+		for (let t of pe(e)) {
 			let e = this.items.get(t.toLowerCase());
 			if (e !== void 0) return e;
 		}
 		return null;
 	}
 	findLoop(e) {
-		for (let t of ge(e)) {
+		for (let t of pe(e)) {
 			let e = t.toLowerCase(), n = this.loops.find((t) => t.tags.includes(e));
 			if (n) return n;
 		}
 		return null;
 	}
 };
-function N(e, t) {
-	for (let n of ge(t)) {
+function I(e, t) {
+	for (let n of pe(t)) {
 		let t = e.tags.indexOf(n.toLowerCase());
 		if (t !== -1) return e.rows.map((e) => e[t] ?? "?");
 	}
 	return null;
 }
-function P(e) {
+function L(e) {
 	return e === "?" || e === ".";
 }
-function F(e) {
-	if (e == null || P(e)) return null;
+function R(e) {
+	if (e == null || L(e)) return null;
 	let t = e.replace(/\(\d+\)$/, "").trim();
 	return /^[+-]?(\d+\.?\d*|\.\d+)([eE][+-]?\d+)?$/.test(t) ? Number(t) : null;
 }
-function ve(e) {
+function he(e) {
 	return e === void 0 || e.kind !== "value";
 }
-function ye(e) {
-	let t = he(e), n = [], r = null, i = 0, a = 0, o = (e) => {
-		if (r || (r = new _e(""), n.push(r)), a > 0) throw new M(`unexpected ${e.kind} inside a save frame`, e.line);
+function ge(e) {
+	let t = fe(e), n = [], r = null, i = 0, a = 0, o = (e) => {
+		if (r || (r = new me(""), n.push(r)), a > 0) throw new P(`unexpected ${e.kind} inside a save frame`, e.line);
 		return r;
 	};
 	for (; i < t.length;) {
@@ -153,7 +153,7 @@ function ye(e) {
 		if (!e) break;
 		switch (e.kind) {
 			case "data":
-				r = new _e(e.text), n.push(r), a = 0, i += 1;
+				r = new me(e.text), n.push(r), a = 0, i += 1;
 				break;
 			case "save":
 				a = e.text === "" ? Math.max(0, a - 1) : a + 1, i += 1;
@@ -167,18 +167,18 @@ function ye(e) {
 					break;
 				}
 				let n = o(e), r = t[i + 1];
-				if (!r || r.kind !== "value") throw new M(`tag ${e.text} has no value`, e.line);
+				if (!r || r.kind !== "value") throw new P(`tag ${e.text} has no value`, e.line);
 				n.items.set(e.text.toLowerCase(), r.text), i += 2;
 				break;
 			}
 			case "loop": {
 				let n = [];
 				for (i += 1; t[i]?.kind === "tag";) n.push(t[i].text.toLowerCase()), i += 1;
-				if (n.length === 0) throw new M("loop_ without tags", e.line);
+				if (n.length === 0) throw new P("loop_ without tags", e.line);
 				let r = [];
-				for (; !ve(t[i]);) r.push(t[i].text), i += 1;
+				for (; !he(t[i]);) r.push(t[i].text), i += 1;
 				if (a > 0) break;
-				if (r.length % n.length !== 0) throw new M(`loop with ${n.length} tags has ${r.length} values`, e.line);
+				if (r.length % n.length !== 0) throw new P(`loop with ${n.length} tags has ${r.length} values`, e.line);
 				let s = [];
 				for (let e = 0; e < r.length; e += n.length) s.push(r.slice(e, e + n.length));
 				o(e).loops.push({
@@ -187,60 +187,60 @@ function ye(e) {
 				});
 				break;
 			}
-			case "value": throw new M(`unexpected value ${JSON.stringify(e.text)}`, e.line);
+			case "value": throw new P(`unexpected value ${JSON.stringify(e.text)}`, e.line);
 		}
 	}
 	return { blocks: n };
 }
 //#endregion
 //#region src/structure/elements.ts
-var I = /* @__PURE__ */ "H.He.Li.Be.B.C.N.O.F.Ne.Na.Mg.Al.Si.P.S.Cl.Ar.K.Ca.Sc.Ti.V.Cr.Mn.Fe.Co.Ni.Cu.Zn.Ga.Ge.As.Se.Br.Kr.Rb.Sr.Y.Zr.Nb.Mo.Tc.Ru.Rh.Pd.Ag.Cd.In.Sn.Sb.Te.I.Xe.Cs.Ba.La.Ce.Pr.Nd.Pm.Sm.Eu.Gd.Tb.Dy.Ho.Er.Tm.Yb.Lu.Hf.Ta.W.Re.Os.Ir.Pt.Au.Hg.Tl.Pb.Bi.Po.At.Rn.Fr.Ra.Ac.Th.Pa.U.Np.Pu.Am.Cm.Bk.Cf.Es.Fm.Md.No.Lr.Rf.Db.Sg.Bh.Hs.Mt.Ds.Rg.Cn.Nh.Fl.Mc.Lv.Ts.Og".split("."), be = new Set(I);
-function xe(e) {
-	return be.has(e) || e === "D";
+var _e = /* @__PURE__ */ "H.He.Li.Be.B.C.N.O.F.Ne.Na.Mg.Al.Si.P.S.Cl.Ar.K.Ca.Sc.Ti.V.Cr.Mn.Fe.Co.Ni.Cu.Zn.Ga.Ge.As.Se.Br.Kr.Rb.Sr.Y.Zr.Nb.Mo.Tc.Ru.Rh.Pd.Ag.Cd.In.Sn.Sb.Te.I.Xe.Cs.Ba.La.Ce.Pr.Nd.Pm.Sm.Eu.Gd.Tb.Dy.Ho.Er.Tm.Yb.Lu.Hf.Ta.W.Re.Os.Ir.Pt.Au.Hg.Tl.Pb.Bi.Po.At.Rn.Fr.Ra.Ac.Th.Pa.U.Np.Pu.Am.Cm.Bk.Cf.Es.Fm.Md.No.Lr.Rf.Db.Sg.Bh.Hs.Mt.Ds.Rg.Cn.Nh.Fl.Mc.Lv.Ts.Og".split("."), ve = new Set(_e);
+function ye(e) {
+	return ve.has(e) || e === "D";
 }
-function Se(e) {
+function be(e) {
 	return e.charAt(0).toUpperCase() + e.slice(1).toLowerCase();
 }
-function Ce(e) {
+function xe(e) {
 	let t = e.match(/^[A-Za-z]{1,2}/)?.[0];
 	if (!t) return null;
-	let n = Se(t);
-	if (t.length === 2 && xe(n)) return n;
+	let n = be(t);
+	if (t.length === 2 && ye(n)) return n;
 	let r = t.charAt(0).toUpperCase();
-	return xe(r) ? r : null;
+	return ye(r) ? r : null;
 }
-function we(e) {
+function Se(e) {
 	return e === "D" ? "H" : e;
 }
 //#endregion
 //#region src/structure/fractional.ts
-var Te = .001, Ee = 1e-9;
-function L(e, t = Ee) {
+var Ce = .001, we = 1e-9;
+function z(e, t = we) {
 	let n = e - Math.floor(e);
 	return n >= 1 - t ? 0 : n;
 }
-function De(e) {
+function Te(e) {
 	return [
-		L(e[0]),
-		L(e[1]),
-		L(e[2])
+		z(e[0]),
+		z(e[1]),
+		z(e[2])
 	];
 }
-function Oe(e, t, n) {
+function Ee(e, t, n) {
 	for (let r = 0; r < 3; r += 1) {
 		let i = Math.abs(e[r] - t[r]);
 		if (i -= Math.round(i), Math.abs(i) > n) return !1;
 	}
 	return !0;
 }
-var R = 100;
-function z(e) {
-	return Math.floor(L(e) * R) % R;
+var B = 100;
+function De(e) {
+	return Math.floor(z(e) * B) % B;
 }
-function ke(e, t, n) {
-	return (e * R + t) * R + n;
+function Oe(e, t, n) {
+	return (e * B + t) * B + n;
 }
-var Ae = class {
+var ke = class {
 	tolerance;
 	buckets = /* @__PURE__ */ new Map();
 	positions = [];
@@ -250,25 +250,25 @@ var Ae = class {
 	add(e) {
 		let t = this.positions.length;
 		this.positions.push(e);
-		let n = ke(z(e[0]), z(e[1]), z(e[2])), r = this.buckets.get(n);
+		let n = Oe(De(e[0]), De(e[1]), De(e[2])), r = this.buckets.get(n);
 		return r ? r.push(t) : this.buckets.set(n, [t]), t;
 	}
 	near(e) {
-		let t = De(e), [n, r, i] = t.map(z), a = [];
+		let t = Te(e), [n, r, i] = t.map(De), a = [];
 		for (let e = -1; e <= 1; e += 1) {
-			let o = (n + e + R) % R;
+			let o = (n + e + B) % B;
 			for (let e = -1; e <= 1; e += 1) {
-				let n = (r + e + R) % R;
+				let n = (r + e + B) % B;
 				for (let e = -1; e <= 1; e += 1) {
-					let r = (i + e + R) % R;
-					for (let e of this.buckets.get(ke(o, n, r)) ?? []) Oe(this.positions[e], t, this.tolerance) && a.push(e);
+					let r = (i + e + B) % B;
+					for (let e of this.buckets.get(Oe(o, n, r)) ?? []) Ee(this.positions[e], t, this.tolerance) && a.push(e);
 				}
 			}
 		}
 		return a.sort((e, t) => e - t);
 	}
 };
-function je(e, t) {
+function Ae(e, t) {
 	let n = e.rotation, [r, i, a] = t, o = e.translation;
 	return [
 		(n[0] ?? 0) * r + (n[1] ?? 0) * i + (n[2] ?? 0) * a + o[0],
@@ -278,12 +278,12 @@ function je(e, t) {
 }
 //#endregion
 //#region src/structure/lattice.ts
-var B = Math.PI / 180;
-function Me(e, t, n) {
+var V = Math.PI / 180;
+function je(e, t, n) {
 	return Math.min(n, Math.max(t, e));
 }
-function Ne(e) {
-	let { a: t, b: n, c: r } = e, i = e.alpha * B, a = e.beta * B, o = e.gamma * B, s = Math.cos(i), c = Math.cos(a), l = Math.cos(o), u = Math.sin(i), d = Math.sin(a), f = Math.acos(Me((s * c - l) / (u * d), -1, 1)), p = [
+function Me(e) {
+	let { a: t, b: n, c: r } = e, i = e.alpha * V, a = e.beta * V, o = e.gamma * V, s = Math.cos(i), c = Math.cos(a), l = Math.cos(o), u = Math.sin(i), d = Math.sin(a), f = Math.acos(je((s * c - l) / (u * d), -1, 1)), p = [
 		t * d,
 		0,
 		t * c
@@ -305,8 +305,8 @@ function Ne(e) {
 		]
 	};
 }
-function Pe(e) {
-	let [t, n, r] = e, i = (e) => Math.hypot(e[0], e[1], e[2]), a = (e, t) => Math.acos(Me((e[0] * t[0] + e[1] * t[1] + e[2] * t[2]) / (i(e) * i(t)), -1, 1)) / B;
+function Ne(e) {
+	let [t, n, r] = e, i = (e) => Math.hypot(e[0], e[1], e[2]), a = (e, t) => Math.acos(je((e[0] * t[0] + e[1] * t[1] + e[2] * t[2]) / (i(e) * i(t)), -1, 1)) / V;
 	return {
 		a: i(t),
 		b: i(n),
@@ -317,7 +317,7 @@ function Pe(e) {
 		vectors: e
 	};
 }
-function Fe(e, t) {
+function Pe(e, t) {
 	let [n, r, i] = e.vectors, [a, o, s] = t;
 	return [
 		a * n[0] + o * r[0] + s * i[0],
@@ -325,49 +325,49 @@ function Fe(e, t) {
 		a * n[2] + o * r[2] + s * i[2]
 	];
 }
-function Ie(e, t) {
+function Fe(e, t) {
 	let n = e[0] - t[0], r = e[1] - t[1], i = e[2] - t[2];
 	return Math.sqrt(n * n + r * r + i * i);
 }
 //#endregion
 //#region src/cif/errors.ts
-var V = class extends Error {
+var H = class extends Error {
 	constructor(e) {
 		super(e), this.name = "CifStructureError";
 	}
-}, Le = [
+}, Ie = [
 	"x",
 	"y",
 	"z"
 ];
-function Re(e) {
+function Le(e) {
 	if (e === "" || e === "+") return 1;
 	if (e === "-") return -1;
 	let t = e.match(/^([+-]?\d+)\/(\d+)$/), n = t ? Number(t[1]) / Number(t[2]) : Number(e);
-	if (!Number.isFinite(n)) throw new V(`Cannot read the number ${JSON.stringify(e)} in a symmetry operator.`);
+	if (!Number.isFinite(n)) throw new H(`Cannot read the number ${JSON.stringify(e)} in a symmetry operator.`);
 	return n;
 }
-function ze(e) {
+function Re(e) {
 	let t = e.split(",").map((e) => e.toLowerCase().replace(/[^xyz0-9./+-]/g, ""));
-	if (t.length !== 3) throw new V(`A symmetry operator needs three components: ${JSON.stringify(e)}.`);
+	if (t.length !== 3) throw new H(`A symmetry operator needs three components: ${JSON.stringify(e)}.`);
 	let n = [], r = [];
 	for (let i of t) {
 		let t = i.match(/[+-]?[^+-]+/g);
-		if (!t || t.join("") !== i) throw new V(`Cannot read the symmetry operator ${JSON.stringify(e)}.`);
+		if (!t || t.join("") !== i) throw new H(`Cannot read the symmetry operator ${JSON.stringify(e)}.`);
 		let a = [
 			0,
 			0,
 			0
 		], o = 0;
 		for (let e of t) {
-			let t = Le.findIndex((t) => e.includes(t));
-			if (t === -1) o += Re(e);
+			let t = Ie.findIndex((t) => e.includes(t));
+			if (t === -1) o += Le(e);
 			else {
-				let n = Re(e.replace(Le[t], "").replace(/\*$/, ""));
+				let n = Le(e.replace(Ie[t], "").replace(/\*$/, ""));
 				a[t] = (a[t] ?? 0) + n;
 			}
 		}
-		n.push(...a), r.push(L(o));
+		n.push(...a), r.push(z(o));
 	}
 	return {
 		rotation: n,
@@ -376,99 +376,99 @@ function ze(e) {
 }
 //#endregion
 //#region src/cif/structure.ts
-var Be = 1e-4, Ve = ["_space_group_symop_operation_xyz", "_symmetry_equiv_pos_as_xyz"], He = ["_space_group_it_number", "_symmetry_int_tables_number"], Ue = [
+var ze = 1e-4, Be = ["_space_group_symop_operation_xyz", "_symmetry_equiv_pos_as_xyz"], Ve = ["_space_group_it_number", "_symmetry_int_tables_number"], He = [
 	"_space_group_name_h-m_alt",
 	"_symmetry_space_group_name_h-m",
 	"_space_group_name_h-m_full"
-], We = ["_space_group_name_hall", "_symmetry_space_group_name_hall"];
-function H(e, t) {
-	let n = F(e.get(t));
-	if (n === null) throw new V(`${e.name || "The CIF"} has no usable ${t}.`);
+], Ue = ["_space_group_name_hall", "_symmetry_space_group_name_hall"];
+function U(e, t) {
+	let n = R(e.get(t));
+	if (n === null) throw new H(`${e.name || "The CIF"} has no usable ${t}.`);
 	return n;
 }
-function Ge(e, t, n) {
+function We(e, t, n) {
 	let r = e.find((e) => e.element === t);
 	r ? r.occupancy += n : e.push({
 		element: t,
 		occupancy: n
 	}), e.sort((e, t) => t.occupancy - e.occupancy);
 }
-function Ke(e) {
+function Ge(e) {
 	for (let t of e) {
 		let e = t.species.reduce((e, t) => e + t.occupancy, 0);
 		if (e > 1.000001) for (let n of t.species) n.occupancy /= e;
 	}
 }
-function qe(e) {
+function Ke(e) {
 	let t = e.findLoop("_atom_site_fract_x");
-	if (!t) throw new V(`${e.name || "The CIF"} has no atom sites.`);
-	let n = N(t, "_atom_site_fract_x") ?? [], r = N(t, "_atom_site_fract_y") ?? [], i = N(t, "_atom_site_fract_z") ?? [], a = N(t, "_atom_site_label"), o = N(t, "_atom_site_type_symbol"), s = N(t, "_atom_site_occupancy"), c = [], l = new Ae(Be);
+	if (!t) throw new H(`${e.name || "The CIF"} has no atom sites.`);
+	let n = I(t, "_atom_site_fract_x") ?? [], r = I(t, "_atom_site_fract_y") ?? [], i = I(t, "_atom_site_fract_z") ?? [], a = I(t, "_atom_site_label"), o = I(t, "_atom_site_type_symbol"), s = I(t, "_atom_site_occupancy"), c = [], l = new ke(ze);
 	for (let e = 0; e < t.rows.length; e += 1) {
-		let t = F(n[e]), u = F(r[e]), d = F(i[e]);
+		let t = R(n[e]), u = R(r[e]), d = R(i[e]);
 		if (t === null || u === null || d === null) continue;
-		let f = a?.[e] ?? o?.[e] ?? `site${e + 1}`, p = o?.[e], m = (p && !P(p) ? Ce(p) : null) ?? Ce(f);
-		if (!m) throw new V(`Atom site ${JSON.stringify(f)} names no element (type symbol ${JSON.stringify(p ?? "")}).`);
-		let h = F(s?.[e]) ?? 1, g = [
-			L(t),
-			L(u),
-			L(d)
-		], [_] = l.near(g);
-		_ === void 0 ? (l.add(g), c.push({
+		let f = a?.[e] ?? o?.[e] ?? `site${e + 1}`, p = o?.[e], m = (p && !L(p) ? xe(p) : null) ?? xe(f);
+		if (!m) throw new H(`Atom site ${JSON.stringify(f)} names no element (type symbol ${JSON.stringify(p ?? "")}).`);
+		let h = R(s?.[e]) ?? 1, g = [
+			z(t),
+			z(u),
+			z(d)
+		], [ee] = l.near(g);
+		ee === void 0 ? (l.add(g), c.push({
 			label: f,
 			species: [{
 				element: m,
 				occupancy: h
 			}],
 			fractional: g
-		})) : Ge(c[_].species, m, h);
+		})) : We(c[ee].species, m, h);
 	}
-	if (c.length === 0) throw new V(`${e.name || "The CIF"} has no atom sites with coordinates.`);
-	return Ke(c), c;
+	if (c.length === 0) throw new H(`${e.name || "The CIF"} has no atom sites with coordinates.`);
+	return Ge(c), c;
 }
-function Je(e) {
-	let t = e.findLoop(Ve);
+function qe(e) {
+	let t = e.findLoop(Be);
 	if (!t) {
-		let t = e.get(Ve);
-		return t === null ? null : [ze(t)];
+		let t = e.get(Be);
+		return t === null ? null : [Re(t)];
 	}
-	let n = (N(t, Ve) ?? []).filter((e) => !P(e)).map(ze);
+	let n = (I(t, Be) ?? []).filter((e) => !L(e)).map(Re);
 	return n.length > 0 ? n : null;
 }
-function Ye(e) {
-	let t = e.get(He), n = t === null ? null : F(t), r = e.get(Ue), i = e.get(We);
+function Je(e) {
+	let t = e.get(Ve), n = t === null ? null : R(t), r = e.get(He), i = e.get(Ue);
 	return {
 		number: n !== null && Number.isInteger(n) && n >= 1 && n <= 230 ? n : null,
-		symbol: r !== null && !P(r) ? r.trim() : null,
-		hallSymbol: i !== null && !P(i) ? i.trim() : null
+		symbol: r !== null && !L(r) ? r.trim() : null,
+		hallSymbol: i !== null && !L(i) ? i.trim() : null
 	};
 }
-function Xe(e) {
+function Ye(e) {
 	return e.findLoop("_atom_site_fract_x") !== null;
 }
-function Ze(e) {
+function Xe(e) {
 	return {
-		lattice: Ne({
-			a: H(e, "_cell_length_a"),
-			b: H(e, "_cell_length_b"),
-			c: H(e, "_cell_length_c"),
-			alpha: H(e, "_cell_angle_alpha"),
-			beta: H(e, "_cell_angle_beta"),
-			gamma: H(e, "_cell_angle_gamma")
+		lattice: Me({
+			a: U(e, "_cell_length_a"),
+			b: U(e, "_cell_length_b"),
+			c: U(e, "_cell_length_c"),
+			alpha: U(e, "_cell_angle_alpha"),
+			beta: U(e, "_cell_angle_beta"),
+			gamma: U(e, "_cell_angle_gamma")
 		}),
-		sites: qe(e),
-		operations: Je(e),
-		spaceGroup: Ye(e)
+		sites: Ke(e),
+		operations: qe(e),
+		spaceGroup: Je(e)
 	};
 }
 //#endregion
 //#region src/structure/closure.ts
-function Qe(e, t) {
-	let n = new Ae(Te);
+function Ze(e, t) {
+	let n = new ke(Ce);
 	for (let t of e) n.add(t.fractional);
 	let r = !1, i = !0;
 	return e.forEach((e, a) => {
 		for (let o of t) {
-			let t = n.near(je(o, e.fractional));
+			let t = n.near(Ae(o, e.fractional));
 			t.some((e) => e !== a) && (r = !0), t.length === 0 && (i = !1);
 		}
 	}), {
@@ -478,7 +478,7 @@ function Qe(e, t) {
 }
 //#endregion
 //#region src/structure/types.ts
-var $e = {
+var Qe = {
 	rotation: [
 		1,
 		0,
@@ -495,29 +495,29 @@ var $e = {
 		0,
 		0
 	]
-}, et = null, tt = null;
-function nt() {
-	return tt ||= import("./moyo-wasm.js").then(async (e) => {
-		await e.instantiate(), et = e;
+}, $e = null, et = null;
+function tt() {
+	return et ||= import("./moyo-wasm.js").then(async (e) => {
+		await e.instantiate(), $e = e;
 	}).catch((e) => {
-		throw tt = null, e;
-	}), tt;
+		throw et = null, e;
+	}), et;
 }
-function U() {
-	if (!et) throw Error("moyo is not loaded yet; await ensureMoyo() first.");
-	return et;
+function W() {
+	if (!$e) throw Error("moyo is not loaded yet; await ensureMoyo() first.");
+	return $e;
 }
 //#endregion
 //#region src/symmetry/moyo.ts
-var rt = class extends Error {
+var nt = class extends Error {
 	constructor(e) {
 		super(e), this.name = "SpaceGroupLookupError";
 	}
-}, W = null;
-function it() {
-	return nt();
+}, G = null;
+function rt() {
+	return tt();
 }
-function at(e) {
+function it(e) {
 	let t = e.rotation;
 	return {
 		rotation: [
@@ -538,119 +538,119 @@ function at(e) {
 		]
 	};
 }
-function ot() {
-	if (!W) {
-		W = [];
-		for (let e = 1; e <= 230; e += 1) W.push(...U().hall_symbol_entries_from_number(e));
+function at() {
+	if (!G) {
+		G = [];
+		for (let e = 1; e <= 230; e += 1) G.push(...W().hall_symbol_entries_from_number(e));
 	}
-	return W;
+	return G;
 }
-function G(e) {
+function K(e) {
 	return e.replace(/[\s_]/g, "").toLowerCase();
 }
-function st(e) {
+function ot(e) {
 	let t = e.trim(), n = t.match(/^(.*?)\s*:\s*(\S+)$/);
 	if (n && n[1] && n[2]) return {
-		base: G(n[1]),
+		base: K(n[1]),
 		setting: n[2].toLowerCase()
 	};
 	let r = t.match(/^(.*\S)\s+([HRhr])$/);
 	return r && r[1] && r[2] ? {
-		base: G(r[1]),
+		base: K(r[1]),
 		setting: r[2].toLowerCase()
 	} : {
-		base: G(t),
+		base: K(t),
 		setting: null
 	};
 }
-function ct(e) {
+function st(e) {
 	let t = e[0];
 	if (!t) return e;
-	let n = U().space_group_type(t.number).hall_number;
+	let n = W().space_group_type(t.number).hall_number;
 	return [...e].sort((e, t) => Number(t.hall_number === n) - Number(e.hall_number === n));
 }
-function lt(e) {
-	let t = ot().filter((t) => e.number === null || t.number === e.number);
+function ct(e) {
+	let t = at().filter((t) => e.number === null || t.number === e.number);
 	if (e.hallSymbol) {
-		let n = G(e.hallSymbol), r = t.filter((e) => G(e.hall_symbol) === n);
+		let n = K(e.hallSymbol), r = t.filter((e) => K(e.hall_symbol) === n);
 		if (r.length > 0) return r;
 	}
 	if (e.symbol) {
-		let { base: n, setting: r } = st(e.symbol), i = t.filter((e) => G(e.hm_short) === n || G(e.hm_full) === n);
+		let { base: n, setting: r } = ot(e.symbol), i = t.filter((e) => K(e.hm_short) === n || K(e.hm_full) === n);
 		if (r) {
 			let e = i.filter((e) => e.setting.toLowerCase() === r);
 			if (e.length > 0) return e;
 		}
-		return ct(i);
+		return st(i);
 	}
 	return [];
 }
-function ut(e) {
+function lt(e) {
 	return [
 		e.hallSymbol,
 		e.symbol,
 		e.number
 	].filter((e) => e !== null).join(" / ");
 }
-async function dt(e) {
-	await it();
-	let { operations_from_number: t } = U();
+async function ut(e) {
+	await rt();
+	let { operations_from_number: t } = W();
 	if (e.hallSymbol || e.symbol) {
-		let n = lt(e);
-		if (n.length === 0) throw new rt(`Space group ${ut(e)} is not in the symmetry table.`);
+		let n = ct(e);
+		if (n.length === 0) throw new nt(`Space group ${lt(e)} is not in the symmetry table.`);
 		return n.map((e) => ({
 			setting: e.setting,
 			operations: t(e.number, {
 				type: "HallNumber",
 				value: e.hall_number
-			}, !1).map(at)
+			}, !1).map(it)
 		}));
 	}
-	if (e.number === null) throw new rt("The CIF names no space group.");
+	if (e.number === null) throw new nt("The CIF names no space group.");
 	return [{
 		setting: "",
-		operations: t(e.number, { type: "Spglib" }, !1).map(at)
+		operations: t(e.number, { type: "Spglib" }, !1).map(it)
 	}];
 }
-async function ft(e) {
-	let [t] = await dt(e);
+async function dt(e) {
+	let [t] = await ut(e);
 	return t.operations;
 }
 //#endregion
 //#region src/symmetry/resolve.ts
-var pt = .001, mt = .01;
-function ht() {
-	return new V("The symmetry does not fit the listed atoms: they are already expanded, yet the operations generate positions the file does not list.");
+var ft = .001, pt = .01;
+function mt() {
+	return new H("The symmetry does not fit the listed atoms: they are already expanded, yet the operations generate positions the file does not list.");
 }
-function gt(e, t) {
-	let n = Qe(e, t);
-	if (n.preExpanded && !n.closed) throw ht();
+function ht(e, t) {
+	let n = Ze(e, t);
+	if (n.preExpanded && !n.closed) throw mt();
 }
-function _t({ a: e, b: t, c: n, alpha: r, beta: i, gamma: a }) {
-	let o = (e, t) => Math.abs(e - t) <= pt * Math.max(e, t), s = (e, t) => Math.abs(e - t) <= mt;
+function gt({ a: e, b: t, c: n, alpha: r, beta: i, gamma: a }) {
+	let o = (e, t) => Math.abs(e - t) <= ft * Math.max(e, t), s = (e, t) => Math.abs(e - t) <= pt;
 	return o(e, t) && o(t, n) && s(r, i) && s(i, a) && !s(r, 90) ? "R" : o(e, t) && s(r, 90) && s(i, 90) && s(a, 120) ? "H" : null;
 }
-async function vt(e) {
-	if (e.operations) return gt(e.sites, e.operations), e.operations;
+async function _t(e) {
+	if (e.operations) return ht(e.sites, e.operations), e.operations;
 	let { spaceGroup: t } = e;
-	if (t.hallSymbol === null && t.symbol === null && t.number === null) return [$e];
-	let n = await dt(t);
+	if (t.hallSymbol === null && t.symbol === null && t.number === null) return [Qe];
+	let n = await ut(t);
 	if (n.length === 1) {
 		let [t] = n;
-		return gt(e.sites, t.operations), t.operations;
+		return ht(e.sites, t.operations), t.operations;
 	}
-	let r = n.map((t) => Qe(e.sites, t.operations)), i = n.find((e, t) => r[t].preExpanded && r[t].closed);
+	let r = n.map((t) => Ze(e.sites, t.operations)), i = n.find((e, t) => r[t].preExpanded && r[t].closed);
 	if (i) return i.operations;
-	if (r.some((e) => e.preExpanded)) throw ht();
-	let a = _t(e.lattice);
+	if (r.some((e) => e.preExpanded)) throw mt();
+	let a = gt(e.lattice);
 	return ((a === null ? void 0 : n.find((e) => e.setting.toUpperCase() === a)) ?? n[0]).operations;
 }
 //#endregion
 //#region src/structure/expand.ts
-function yt(e) {
-	let t = [], n = new Ae(Te);
+function vt(e) {
+	let t = [], n = new ke(Ce);
 	for (let r of e.sites) for (let i of e.operations) {
-		let e = De(je(i, r.fractional));
+		let e = Te(Ae(i, r.fractional));
 		n.near(e).length === 0 && (n.add(e), t.push({
 			label: r.label,
 			species: r.species,
@@ -660,17 +660,17 @@ function yt(e) {
 	return {
 		...e,
 		sites: t,
-		operations: [$e]
+		operations: [Qe]
 	};
 }
 //#endregion
 //#region src/symmetry/dataset.ts
-var bt = .01, xt = class extends Error {
+var yt = .01, bt = class extends Error {
 	constructor(e) {
 		super(e), this.name = "SymmetryAnalysisError";
 	}
 };
-function St(e) {
+function xt(e) {
 	let t = /* @__PURE__ */ new Map(), n = /* @__PURE__ */ new Map();
 	return {
 		numbers: e.map((e) => {
@@ -680,7 +680,7 @@ function St(e) {
 		speciesByLabel: n
 	};
 }
-function Ct(e, t, n) {
+function St(e, t, n) {
 	let r = e.std_cell.lattice.basis, i = [
 		[
 			r[0],
@@ -702,13 +702,13 @@ function Ct(e, t, n) {
 		return a.set(s, c), {
 			label: `${s}${c}`,
 			species: o,
-			fractional: De(n)
+			fractional: Te(n)
 		};
 	});
 	return {
-		lattice: Pe(i),
+		lattice: Ne(i),
 		sites: o,
-		operations: [$e],
+		operations: [Qe],
 		spaceGroup: {
 			number: e.number,
 			symbol: e.hm_symbol,
@@ -716,9 +716,9 @@ function Ct(e, t, n) {
 		}
 	};
 }
-async function wt(e, t = {}) {
-	await it();
-	let { analyze_cell: n, hall_symbol_entry: r, space_group_type: i } = U(), a = yt(e), { numbers: o, speciesByLabel: s } = St(a.sites), [c, l, u] = a.lattice.vectors, d = {
+async function Ct(e, t = {}) {
+	await rt();
+	let { analyze_cell: n, hall_symbol_entry: r, space_group_type: i } = W(), a = vt(e), { numbers: o, speciesByLabel: s } = xt(a.sites), [c, l, u] = a.lattice.vectors, d = {
 		lattice: { basis: [
 			...c,
 			...l,
@@ -728,11 +728,11 @@ async function wt(e, t = {}) {
 		numbers: o
 	}, f;
 	try {
-		f = n(JSON.stringify(d), t.symprec ?? bt, "Standard");
+		f = n(JSON.stringify(d), t.symprec ?? yt, "Standard");
 	} catch (e) {
-		throw new xt(`moyo could not analyse the structure: ${e instanceof Error ? e.message : String(e)}`);
+		throw new bt(`moyo could not analyse the structure: ${e instanceof Error ? e.message : String(e)}`);
 	}
-	let p = r(f.hall_number), m = i(f.number), h = Ct(f, s, p);
+	let p = r(f.hall_number), m = i(f.number), h = St(f, s, p);
 	return {
 		number: f.number,
 		symbol: f.hm_symbol,
@@ -741,7 +741,7 @@ async function wt(e, t = {}) {
 		hallNumber: f.hall_number,
 		setting: p.setting,
 		crystalSystem: m.crystal_system,
-		operations: f.operations.map(at),
+		operations: f.operations.map(it),
 		wyckoffs: f.wyckoffs,
 		siteSymmetries: f.site_symmetry_symbols,
 		structure: a,
@@ -750,34 +750,34 @@ async function wt(e, t = {}) {
 }
 //#endregion
 //#region src/symmetry/symbol.ts
-var Tt = "₀₁₂₃₄₅₆₇₈₉", Et = "̅";
-function Dt(e) {
-	return Tt[Number(e)] ?? e;
+var wt = "₀₁₂₃₄₅₆₇₈₉", Tt = "̅";
+function Et(e) {
+	return wt[Number(e)] ?? e;
 }
-function Ot(e) {
+function Dt(e) {
 	let t = e.trim(), n = t.indexOf(":"), r = n === -1 ? t : t.slice(0, n).trim(), i = n === -1 ? "" : ` ${t.slice(n)}`, a = r.split(/\s+/);
-	return a.map((e) => a.length > 1 ? e.replace(/^([2346])([1-5])(?=$|\/)/, (e, t, n) => Number(n) < Number(t) ? t + Dt(n) : e) : e).join("").replace(/(\d)_(\d)/g, (e, t, n) => t + Dt(n)).replace(/-(\d)/g, (e, t) => t + Et) + i;
+	return a.map((e) => a.length > 1 ? e.replace(/^([2346])([1-5])(?=$|\/)/, (e, t, n) => Number(n) < Number(t) ? t + Et(n) : e) : e).join("").replace(/(\d)_(\d)/g, (e, t, n) => t + Et(n)).replace(/-(\d)/g, (e, t) => t + Tt) + i;
 }
 //#endregion
 //#region src/structure/bonds.ts
-var kt = {
+var Ot = {
 	source: "VESTA bond-search cutoffs (K. Momma and F. Izumi), as distributed in pymatgen's vesta_cutoffs.yaml (MIT). Where VESTA lists a pair twice (H-O, D-O), the shorter cutoff is kept.",
 	pairs: /* @__PURE__ */ JSON.parse("{\"Ac-Br\":3.22726,\"Ac-Cl\":3.08646,\"Ac-F\":2.58646,\"Ac-O\":2.7326,\"Ag-As\":2.45642,\"Ag-Br\":2.37642,\"Ag-Cl\":3.05939,\"Ag-F\":2.76939,\"Ag-H\":1.65642,\"Ag-I\":2.53642,\"Ag-N\":2.00642,\"Ag-O\":2.81139,\"Ag-P\":2.37642,\"Ag-S\":3.08839,\"Ag-Se\":2.41642,\"Ag-Te\":2.66642,\"Al-As\":2.75646,\"Al-Br\":2.65646,\"Al-Cl\":2.48846,\"Al-F\":2.00146,\"Al-H\":1.90646,\"Al-I\":2.86646,\"Al-N\":2.24646,\"Al-O\":2.1074,\"Al-P\":2.69646,\"Al-S\":2.66646,\"Al-Se\":2.72646,\"Al-Te\":2.93646,\"Am-Br\":3.22944,\"Am-Cl\":3.08945,\"Am-F\":2.61945,\"Am-O\":2.71649,\"As-Au\":2.26998,\"As-B\":2.42646,\"As-Ba\":3.69,\"As-Be\":2.60649,\"As-Bi\":2.87642,\"As-Br\":2.80646,\"As-C\":2.38646,\"As-Ca\":3.48295,\"As-Cd\":3.29295,\"As-Ce\":3.08644,\"As-Cl\":2.61646,\"As-Co\":2.43642,\"As-Cr\":2.49642,\"As-Cs\":4.15942,\"As-Cu\":2.71895,\"As-Dy\":3.14,\"As-Er\":3.18,\"As-Eu\":2.93898,\"As-F\":2.15646,\"As-Fe\":2.50642,\"As-Ga\":2.38998,\"As-Gd\":3.18,\"As-Ge\":2.47998,\"As-Hf\":2.71642,\"As-Hg\":2.65642,\"As-Ho\":2.87898,\"As-I\":3.03646,\"As-In\":2.66642,\"As-Ir\":2.58998,\"As-K\":3.94942,\"As-La\":3.51293,\"As-Lu\":3.19649,\"As-Mg\":3.09293,\"As-Mn\":2.51642,\"As-Mo\":2.62701,\"As-Na\":3.64942,\"As-Nb\":2.69642,\"As-Ni\":2.28998,\"As-O\":2.24546,\"As-P\":2.29998,\"As-Pb\":3.02644,\"As-Pd\":2.34998,\"As-Pr\":3.35649,\"As-Pt\":2.30998,\"As-Rb\":4.04645,\"As-Re\":2.61998,\"As-Rh\":2.41998,\"As-Ru\":2.40998,\"As-S\":2.84649,\"As-Sb\":2.64998,\"As-Sc\":3.04291,\"As-Se\":2.98649,\"As-Si\":2.66002,\"As-Sm\":3.30649,\"As-Sn\":2.77642,\"As-Sr\":3.62295,\"As-Ta\":2.70642,\"As-Tb\":3.26649,\"As-Te\":3.10646,\"As-Th\":3.40649,\"As-Ti\":2.87642,\"As-Tl\":3.09644,\"As-Tm\":3.22649,\"As-U\":3.02644,\"As-V\":2.54642,\"As-W\":2.58998,\"As-Y\":3.24649,\"As-Yb\":3.19649,\"As-Zn\":2.95293,\"As-Zr\":3.07004,\"Au-Br\":2.77646,\"Au-Cl\":2.88295,\"Au-F\":2.34646,\"Au-H\":1.41998,\"Au-I\":3.21295,\"Au-N\":2.3826,\"Au-O\":2.34646,\"Au-P\":2.18998,\"Au-S\":2.8326,\"Au-Se\":2.22998,\"Au-Te\":2.45998,\"B-B\":1.85846,\"B-Br\":2.33646,\"B-Cl\":2.19646,\"B-F\":1.76646,\"B-H\":1.59646,\"B-I\":2.55646,\"B-N\":1.93846,\"B-O\":1.82746,\"B-P\":2.37646,\"B-S\":2.27646,\"B-Se\":2.40646,\"B-Te\":2.65646,\"Ba-Br\":3.74295,\"Ba-Cl\":3.55295,\"Ba-F\":3.05095,\"Ba-H\":3.08295,\"Ba-I\":3.99295,\"Ba-N\":3.33295,\"Ba-O\":3.14795,\"Ba-P\":3.48295,\"Ba-S\":3.66195,\"Ba-Se\":3.74295,\"Ba-Te\":3.94295,\"Be-Br\":2.50649,\"Be-Cl\":2.36649,\"Be-F\":1.88749,\"Be-H\":1.71649,\"Be-I\":2.70649,\"Be-N\":2.10649,\"Be-O\":1.98749,\"Be-P\":2.55649,\"Be-S\":2.43649,\"Be-Se\":2.57649,\"Be-Te\":2.81649,\"Bi-Br\":3.17993,\"Bi-Cl\":3.04291,\"Bi-F\":2.55291,\"Bi-H\":2.12642,\"Bi-I\":3.38291,\"Bi-N\":2.56329,\"Bi-O\":2.6608,\"Bi-P\":2.78642,\"Bi-S\":3.13291,\"Bi-Se\":3.24329,\"Bi-Te\":3.02642,\"Bk-Br\":3.15233,\"Bk-Cl\":3.02291,\"Bk-F\":2.54233,\"Bk-O\":2.64329,\"Br-C\":2.26002,\"Br-Ca\":3.36995,\"Br-Cd\":3.21295,\"Br-Ce\":3.40452,\"Br-Cf\":3.14233,\"Br-Cl\":2.33296,\"Br-Co\":2.33642,\"Br-Cr\":2.97293,\"Br-Cs\":4.06942,\"Br-Cu\":2.89295,\"Br-Dy\":3.16945,\"Br-Er\":3.14945,\"Br-Eu\":3.46549,\"Br-F\":2.20646,\"Br-Fe\":2.8952,\"Br-Ga\":2.6426,\"Br-Gd\":3.19945,\"Br-Ge\":2.34998,\"Br-Hf\":2.62642,\"Br-Hg\":3.09293,\"Br-Ho\":3.20159,\"Br-In\":3.05329,\"Br-Ir\":2.49998,\"Br-K\":3.8513,\"Br-La\":3.43293,\"Br-Li\":3.11654,\"Br-Lu\":3.11945,\"Br-Mg\":2.99293,\"Br-Mn\":3.05293,\"Br-Mo\":2.8535,\"Br-Na\":3.57715,\"Br-Nb\":3.07646,\"Br-Nd\":3.37293,\"Br-Ni\":2.80649,\"Br-Np\":3.21233,\"Br-O\":2.35646,\"Br-Os\":2.72002,\"Br-P\":2.44293,\"Br-Pa\":3.18437,\"Br-Pb\":3.62451,\"Br-Pd\":2.80649,\"Br-Pm\":3.18233,\"Br-Pr\":3.27649,\"Br-Pt\":2.94191,\"Br-Pu\":3.19233,\"Br-Rb\":4.05498,\"Br-Re\":2.70002,\"Br-Rh\":2.7126,\"Br-Ru\":2.30998,\"Br-S\":2.21998,\"Br-Sb\":2.96646,\"Br-Sc\":2.94291,\"Br-Se\":2.78002,\"Br-Si\":2.55002,\"Br-Sm\":3.26649,\"Br-Sn\":3.2152,\"Br-Sr\":3.54295,\"Br-Ta\":2.60642,\"Br-Tb\":3.18649,\"Br-Te\":2.90002,\"Br-Th\":3.31945,\"Br-Ti\":3.20293,\"Br-Tl\":3.80942,\"Br-Tm\":3.13945,\"Br-U\":3.39452,\"Br-V\":2.87329,\"Br-W\":2.49998,\"Br-Y\":3.15649,\"Br-Yb\":3.12945,\"Br-Zn\":2.86293,\"Br-Zr\":2.98004,\"C-C\":1.89002,\"C-Cl\":2.11002,\"C-Co\":2.19691,\"C-Cu\":2.32649,\"C-F\":1.76002,\"C-Fe\":2.25191,\"C-H\":1.2,\"C-I\":2.16998,\"C-N\":1.79202,\"C-O\":1.97249,\"C-P\":1.93998,\"C-Pd\":2.33649,\"C-Pt\":2.36649,\"C-S\":2.15002,\"C-Se\":2.01998,\"C-Si\":2.23302,\"C-Te\":2.25998,\"Ca-Cl\":3.23295,\"Ca-F\":2.70495,\"Ca-H\":2.69295,\"Ca-I\":3.58295,\"Ca-N\":3.00295,\"Ca-O\":2.83062,\"Ca-P\":3.41295,\"Ca-S\":3.31295,\"Ca-Se\":3.42295,\"Ca-Te\":3.62295,\"Cd-Cl\":3.09295,\"Cd-F\":2.67395,\"Cd-H\":2.52295,\"Cd-I\":3.46295,\"Cd-N\":2.82295,\"Cd-O\":2.76695,\"Cd-P\":3.20295,\"Cd-S\":3.16695,\"Cd-Se\":3.26295,\"Cd-Te\":3.45295,\"Ce-Cl\":3.25293,\"Ce-F\":2.75452,\"Ce-H\":2.34644,\"Ce-I\":3.62452,\"Ce-N\":2.78549,\"Ce-O\":2.86393,\"Ce-P\":3.00644,\"Ce-S\":3.36293,\"Ce-Se\":3.04644,\"Ce-Te\":3.22644,\"Cf-Cl\":3.01291,\"Cf-F\":2.53233,\"Cf-O\":2.63291,\"Cl-Cl\":2.14296,\"Cl-Cm\":3.18291,\"Cl-Co\":2.74593,\"Cl-Cr\":2.80293,\"Cl-Cs\":3.91042,\"Cl-Cu\":2.75295,\"Cl-D\":1.5,\"Cl-Dy\":3.01945,\"Cl-Er\":2.99944,\"Cl-Eu\":3.32549,\"Cl-F\":2.14646,\"Cl-Fe\":2.86293,\"Cl-Ga\":2.52646,\"Cl-Gd\":3.06349,\"Cl-Ge\":2.49002,\"Cl-H\":1.5,\"Cl-Hf\":2.75646,\"Cl-Hg\":3.24939,\"Cl-Ho\":3.05159,\"Cl-I\":3.33295,\"Cl-In\":3.52939,\"Cl-Ir\":2.56746,\"Cl-K\":3.65976,\"Cl-La\":3.33452,\"Cl-Li\":2.91814,\"Cl-Lu\":2.96944,\"Cl-Mg\":2.79293,\"Cl-Mn\":2.84593,\"Cl-Mo\":2.80447,\"Cl-N\":2.20646,\"Cl-NH\":3.48195,\"Cl-Na\":3.39412,\"Cl-Nb\":2.76291,\"Cl-Nd\":3.22493,\"Cl-Ni\":2.62649,\"Cl-Np\":3.07233,\"Cl-O\":2.16646,\"Cl-Os\":2.54002,\"Cl-P\":2.28746,\"Cl-Pa\":3.01437,\"Cl-Pb\":3.39295,\"Cl-Pd\":2.65649,\"Cl-Pm\":3.41233,\"Cl-Pr\":3.12749,\"Cl-Pt\":2.75646,\"Cl-Pu\":3.05233,\"Cl-Rb\":3.86664,\"Cl-Re\":3.44712,\"Cl-Rh\":2.62646,\"Cl-Ru\":2.70646,\"Cl-S\":2.37002,\"Cl-Sb\":2.80646,\"Cl-Sc\":2.92291,\"Cl-Se\":2.57002,\"Cl-Si\":2.38002,\"Cl-Sm\":3.08749,\"Cl-Sn\":3.13111,\"Cl-Sr\":3.37295,\"Cl-Ta\":2.6739,\"Cl-Tb\":3.04349,\"Cl-Tc\":2.56002,\"Cl-Te\":2.73906,\"Cl-Th\":3.15945,\"Cl-Ti\":3.02293,\"Cl-Tl\":3.72942,\"Cl-Tm\":2.98944,\"Cl-U\":3.24452,\"Cl-V\":3.15293,\"Cl-W\":2.47,\"Cl-Y\":3.00649,\"Cl-Yb\":2.98249,\"Cl-Zn\":2.72293,\"Cl-Zr\":3.33651,\"Cm-F\":2.68291,\"Cm-O\":2.79291,\"Co-F\":2.35293,\"Co-H\":1.9278,\"Co-I\":2.52878,\"Co-In\":3.13629,\"Co-N\":2.36293,\"Co-O\":2.40493,\"Co-P\":2.36642,\"Co-S\":2.65293,\"Co-Se\":2.39642,\"Co-Te\":2.61642,\"Cr-F\":2.45293,\"Cr-H\":1.67642,\"Cr-I\":3.19293,\"Cr-N\":2.5152,\"Cr-O\":2.44293,\"Cr-P\":2.42642,\"Cr-S\":2.72491,\"Cr-Se\":2.44642,\"Cr-Te\":2.67642,\"Cs-F\":3.49942,\"Cs-H\":3.55942,\"Cs-I\":4.40942,\"Cs-N\":3.94942,\"Cs-O\":3.53642,\"Cs-P\":3.64194,\"Cs-S\":4.24942,\"Cs-Se\":4.21655,\"Cs-Te\":4.4631,\"Cu-F\":2.46295,\"Cu-H\":1.81649,\"Cu-I\":3.01795,\"Cu-N\":2.49295,\"Cu-O\":2.47295,\"Cu-P\":2.65649,\"Cu-S\":2.76095,\"Cu-Se\":2.76295,\"Cu-Te\":2.87649,\"D-F\":1.1,\"D-N\":1.2,\"D-O\":1.2,\"Dy-F\":2.52944,\"Dy-H\":2.09,\"Dy-I\":3.39945,\"Dy-N\":2.38,\"Dy-O\":2.65651,\"Dy-P\":2.77,\"Dy-S\":2.823,\"Dy-Se\":2.81,\"Dy-Te\":3.22,\"Er-F\":2.51049,\"Er-H\":2.06,\"Er-I\":3.38945,\"Er-N\":2.36,\"Er-O\":2.63651,\"Er-P\":2.75,\"Er-S\":3.27651,\"Er-Se\":3.18649,\"Er-Te\":3.22,\"Es-O\":2.70139,\"Eu-F\":2.83549,\"Eu-H\":2.18898,\"Eu-I\":3.69549,\"Eu-N\":2.95549,\"Eu-O\":2.94249,\"Eu-P\":2.85898,\"Eu-S\":3.37949,\"Eu-Se\":2.89898,\"Eu-Te\":3.08898,\"F-Fe\":2.36293,\"F-Ga\":2.14646,\"F-Gd\":3.15651,\"F-Ge\":2.01002,\"F-H\":1.1,\"F-Hf\":3.18291,\"F-Hg\":2.88293,\"F-Ho\":2.56159,\"F-I\":3.18295,\"F-In\":2.35491,\"F-Ir\":2.15002,\"F-K\":3.11142,\"F-Kr\":2.67549,\"F-La\":2.79293,\"F-Li\":2.34276,\"F-Lu\":2.48249,\"F-Mg\":2.29093,\"F-Mn\":2.41093,\"F-Mo\":2.2998,\"F-N\":1.82646,\"F-NH\":2.99195,\"F-Na\":2.80398,\"F-Nb\":2.35646,\"F-Nd\":2.73452,\"F-Ni\":2.20249,\"F-Np\":2.59233,\"F-Os\":2.07746,\"F-P\":2.01002,\"F-Pa\":2.54437,\"F-Pb\":2.92045,\"F-Pd\":2.34649,\"F-Pm\":2.55233,\"F-Po\":2.83646,\"F-Pr\":2.62945,\"F-Pt\":2.54002,\"F-Pu\":2.58233,\"F-Rb\":3.37645,\"F-Re\":2.16002,\"F-Rh\":2.16646,\"F-Ru\":2.57646,\"F-S\":1.95002,\"F-Sb\":2.35646,\"F-Sc\":2.32291,\"F-Se\":2.08002,\"F-Si\":1.93002,\"F-Sm\":2.60649,\"F-Sn\":2.63793,\"F-Sr\":2.88195,\"F-Ta\":2.2539,\"F-Tb\":2.54249,\"F-Tc\":2.24219,\"F-Te\":2.22002,\"F-Th\":2.68945,\"F-Ti\":2.86293,\"F-Tl\":3.26942,\"F-Tm\":2.51649,\"F-U\":2.80293,\"F-V\":2.87293,\"F-W\":2.03,\"F-Xe\":2.62649,\"F-Y\":2.51049,\"F-Yb\":2.50649,\"F-Zn\":2.38293,\"F-Zr\":2.99651,\"Fe-H\":1.68642,\"Fe-I\":3.1552,\"Fe-N\":2.48193,\"Fe-O\":2.44693,\"Fe-P\":2.42642,\"Fe-S\":2.83793,\"Fe-Se\":2.43642,\"Fe-Te\":2.68642,\"Ga-H\":1.55998,\"Ga-I\":2.91646,\"Ga-N\":1.96,\"Ga-O\":2.18646,\"Ga-P\":2.46998,\"Ga-S\":2.61946,\"Ga-Se\":3.41295,\"Ga-Te\":2.58998,\"Gd-H\":2.13,\"Gd-I\":3.41945,\"Gd-N\":2.42,\"Gd-O\":2.76651,\"Gd-P\":2.81,\"Gd-S\":3.13649,\"Gd-Se\":2.85,\"Gd-Te\":3.24,\"Ge-Ge\":2.6,\"Ge-H\":1.59998,\"Ge-I\":2.54998,\"Ge-N\":1.92998,\"Ge-O\":2.09802,\"Ge-P\":2.36998,\"Ge-S\":2.56702,\"Ge-Se\":2.70002,\"Ge-Te\":2.60998,\"H-Hf\":1.93642,\"H-Hg\":1.86642,\"H-Ho\":2.11898,\"H-In\":1.87642,\"H-Ir\":1.80998,\"H-K\":3.21942,\"H-La\":2.77293,\"H-Lu\":2.42649,\"H-Mg\":2.24293,\"H-Mn\":1.70642,\"H-Mo\":1.83701,\"H-N\":1.2,\"H-Na\":2.79942,\"H-Nb\":1.90642,\"H-Ni\":1.44998,\"H-O\":1.2,\"H-P\":1.45998,\"H-Pb\":2.27644,\"H-Pd\":1.51998,\"H-Pr\":2.62649,\"H-Pt\":1.44998,\"H-Rb\":3.43645,\"H-Re\":1.79998,\"H-Rh\":1.59998,\"H-Ru\":1.65998,\"H-S\":1.42998,\"H-Sb\":2.81998,\"H-Sc\":2.24291,\"H-Se\":1.58998,\"H-Si\":1.82002,\"H-Sm\":2.56649,\"H-Sn\":2.00642,\"H-Sr\":2.87295,\"H-Ta\":1.91642,\"H-Tb\":2.51649,\"H-Te\":1.87998,\"H-Th\":2.67649,\"H-Ti\":1.76642,\"H-Tl\":2.35644,\"H-Tm\":2.45649,\"H-U\":2.27644,\"H-V\":1.73642,\"H-W\":1.80998,\"H-Y\":2.46649,\"H-Yb\":2.42649,\"H-Zn\":2.13293,\"H-Zr\":2.29004,\"Hf-I\":2.83642,\"Hf-N\":2.24642,\"Hf-O\":2.37946,\"Hf-P\":2.63642,\"Hf-S\":2.64642,\"Hf-Se\":2.67642,\"Hf-Te\":2.87642,\"Hg-Hg\":3.1952,\"Hg-I\":3.33293,\"Hg-N\":2.17642,\"Hg-O\":2.86939,\"Hg-P\":2.57642,\"Hg-S\":3.42093,\"Hg-Se\":2.82642,\"Hg-Te\":2.76642,\"Ho-I\":3.44159,\"Ho-N\":2.41898,\"Ho-O\":2.67047,\"Ho-P\":2.79898,\"Ho-S\":3.13547,\"Ho-Se\":2.84898,\"Ho-Te\":3.03898,\"I-I\":2.4,\"I-In\":3.19291,\"I-Ir\":2.70998,\"I-K\":4.11717,\"I-La\":3.64293,\"I-Li\":3.37676,\"I-Lu\":3.36945,\"I-Mg\":3.17293,\"I-Mn\":3.23293,\"I-Mo\":2.74701,\"I-Na\":3.88251,\"I-Nb\":3.1439,\"I-Nd\":3.59452,\"I-Ni\":3.00649,\"I-Np\":3.44233,\"I-O\":2.47646,\"I-P\":2.44998,\"I-Pb\":3.69562,\"I-Pd\":2.96649,\"I-Pr\":3.49649,\"I-Pt\":2.41998,\"I-Pu\":3.43233,\"I-Rb\":4.33462,\"I-Re\":2.65998,\"I-Rh\":2.52998,\"I-Ru\":2.52998,\"I-S\":2.40998,\"I-Sb\":3.21646,\"I-Sc\":3.15291,\"I-Se\":2.58998,\"I-Si\":2.76002,\"I-Sm\":3.44649,\"I-Sn\":3.52293,\"I-Sr\":3.74295,\"I-Ta\":2.81642,\"I-Tb\":3.40945,\"I-Te\":3.13702,\"I-Th\":3.56649,\"I-Ti\":3.08291,\"I-Tl\":3.94142,\"I-Tm\":3.37945,\"I-U\":3.62452,\"I-V\":2.66642,\"I-W\":2.70998,\"I-Y\":3.37649,\"I-Yb\":3.37945,\"I-Zn\":3.07293,\"I-Zr\":3.19004,\"In-Mn\":3.14729,\"In-N\":2.18642,\"In-O\":2.46491,\"In-P\":2.68642,\"In-S\":2.93291,\"In-Se\":2.62642,\"In-Te\":2.84642,\"Ir-N\":2.10998,\"Ir-O\":2.27746,\"Ir-P\":2.50998,\"Ir-S\":2.42998,\"Ir-Se\":2.55998,\"Ir-Te\":2.75998,\"K-N\":3.41942,\"K-O\":3.25142,\"K-P\":3.44942,\"K-S\":3.79285,\"K-Se\":4.00186,\"K-Te\":4.23284,\"La-N\":3.05293,\"La-O\":2.90983,\"La-P\":3.32293,\"La-S\":3.35593,\"La-Se\":3.45293,\"La-Te\":3.65293,\"Li-N\":2.66213,\"Li-O\":2.60087,\"Li-S\":3.02481,\"Li-Se\":3.2433,\"Li-Te\":3.42496,\"Lu-N\":2.71649,\"Lu-O\":2.57749,\"Lu-P\":3.11649,\"Lu-S\":3.03649,\"Lu-Se\":3.16649,\"Lu-Te\":3.35649,\"Mg-N\":2.56293,\"Mg-O\":2.41824,\"Mg-P\":3.00293,\"Mg-S\":2.89293,\"Mg-Se\":3.03293,\"Mg-Te\":3.24293,\"Mn-N\":2.56193,\"Mn-O\":2.51652,\"Mn-P\":2.39642,\"Mn-S\":2.93293,\"Mn-Se\":2.47642,\"Mn-Te\":2.70642,\"Mo-N\":2.4735,\"Mo-O\":2.3475,\"Mo-P\":2.54701,\"Mo-S\":2.80067,\"Mo-Se\":2.59701,\"Mo-Te\":2.79701,\"N-N\":1.8826,\"N-Na\":3.12942,\"N-Nb\":2.46046,\"N-Nd\":3.01293,\"N-Ni\":2.30649,\"N-O\":1.81746,\"N-P\":1.97146,\"N-Pb\":3.0967,\"N-Pd\":2.40451,\"N-Pr\":2.90649,\"N-Pt\":2.41649,\"N-Rb\":3.79645,\"N-Re\":2.10998,\"N-Rh\":2.2626,\"N-Ru\":2.2626,\"N-S\":2.28849,\"N-Sb\":2.56446,\"N-Sc\":2.54291,\"N-Se\":2.1,\"N-Si\":2.12002,\"N-Sm\":3.02351,\"N-Sn\":2.7152,\"N-Sr\":3.09295,\"N-Ta\":2.16642,\"N-Tb\":2.80649,\"N-Te\":2.16998,\"N-Th\":2.94649,\"N-Ti\":2.08642,\"N-Tl\":2.59644,\"N-Tm\":2.74649,\"N-U\":2.78649,\"N-V\":2.38329,\"N-W\":2.10998,\"N-Y\":2.77649,\"N-Yb\":2.84851,\"N-Zn\":2.43293,\"N-Zr\":2.65004,\"NH-O\":3.08895,\"Na-O\":2.95693,\"Na-P\":3.47942,\"Na-S\":3.57685,\"Na-Se\":3.71593,\"Na-Te\":3.95459,\"Nb-O\":2.45329,\"Nb-P\":2.61642,\"Nb-S\":2.74,\"Nb-Se\":2.66642,\"Nb-Te\":2.85642,\"Nd-O\":2.8587,\"Nd-S\":3.42712,\"Nd-Se\":3.42293,\"Nd-Te\":3.60293,\"Ni-O\":2.28149,\"Ni-P\":2.31998,\"Ni-S\":2.58649,\"Ni-Se\":2.18998,\"Ni-Te\":2.47998,\"Np-O\":2.63646,\"Np-S\":3.2,\"O-O\":1.7,\"O-Os\":2.23,\"O-P\":2.08646,\"O-Pa\":2.63383,\"O-Pb\":3.04096,\"O-Pd\":2.39849,\"O-Po\":2.64646,\"O-Pr\":2.74449,\"O-Pt\":2.40649,\"O-Pu\":2.68329,\"O-Rb\":3.43945,\"O-Re\":2.3426,\"O-Rh\":2.24946,\"O-Ru\":2.22646,\"O-S\":1.8,\"O-Sb\":2.45237,\"O-Sc\":2.42029,\"O-Se\":2.16102,\"O-Si\":1.99002,\"O-Sm\":2.88251,\"O-Sn\":2.82146,\"O-Sr\":2.98095,\"O-Ta\":2.74646,\"O-Tb\":2.65549,\"O-Tc\":2.22446,\"O-Te\":2.3334,\"O-Th\":2.77349,\"O-Ti\":2.35391,\"O-Tl\":3.36945,\"O-Tm\":2.60649,\"O-U\":2.94295,\"O-V\":2.84939,\"O-W\":2.20102,\"O-Xe\":2.63451,\"O-Y\":2.62549,\"O-Yb\":2.74551,\"O-Zn\":2.41693,\"O-Zr\":3.09651,\"Os-S\":2.56002,\"P-P\":2.48381,\"P-Pb\":2.94644,\"P-Pd\":2.46998,\"P-Pr\":3.28649,\"P-Pt\":2.23998,\"P-Rb\":3.50645,\"P-Re\":2.50998,\"P-Rh\":2.43998,\"P-Ru\":2.33998,\"P-S\":2.57646,\"P-Sb\":2.56998,\"P-Sc\":2.96291,\"P-Se\":2.69646,\"P-Si\":2.58002,\"P-Sm\":3.23649,\"P-Sn\":2.60642,\"P-Sr\":3.43295,\"P-Ta\":2.62642,\"P-Tb\":3.19649,\"P-Te\":2.56998,\"P-Th\":3.33649,\"P-Ti\":2.51642,\"P-Tl\":3.01644,\"P-Tm\":3.13649,\"P-U\":2.94644,\"P-V\":2.46642,\"P-W\":2.50998,\"P-Y\":3.17649,\"P-Yb\":3.13649,\"P-Zn\":2.86293,\"P-Zr\":3.02004,\"Pb-S\":3.40395,\"Pb-Se\":3.55295,\"Pb-Te\":3.14644,\"Pd-S\":2.69649,\"Pd-Se\":2.26998,\"Pd-Te\":2.52998,\"Pr-S\":3.20649,\"Pr-Se\":3.32649,\"Pr-Te\":3.50649,\"Pt-S\":2.76649,\"Pt-Se\":2.23998,\"Pt-Te\":2.49998,\"Pu-S\":3.3,\"Rb-S\":3.97645,\"Rb-Se\":4.13773,\"Rb-Te\":4.28802,\"Re-S\":2.61998,\"Re-Se\":2.54998,\"Re-Te\":2.74998,\"Rh-S\":2.19998,\"Rh-Se\":2.37998,\"Rh-Te\":2.59998,\"Ru-S\":2.6426,\"Ru-Se\":2.69451,\"Ru-Te\":2.58998,\"S-S\":2.2,\"S-Sb\":3.05,\"S-Sc\":2.88391,\"S-Se\":2.81649,\"S-Si\":2.47602,\"S-Sm\":3.15649,\"S-Sn\":3.15293,\"S-Sr\":3.51295,\"S-Ta\":2.8439,\"S-Tb\":3.11649,\"S-Te\":2.79002,\"S-Th\":3.24649,\"S-Ti\":2.74646,\"S-Tl\":3.66442,\"S-Tm\":3.05649,\"S-U\":3.25293,\"S-V\":2.82293,\"S-W\":2.43998,\"S-Y\":3.08649,\"S-Yb\":3.03649,\"S-Zn\":2.80293,\"S-Zr\":2.91004,\"Sb-Se\":3.05646,\"Sb-Te\":2.82998,\"Sc-Se\":3.00291,\"Sc-Te\":3.20291,\"Se-Se\":2.93649,\"Se-Si\":2.61002,\"Se-Sm\":3.27649,\"Se-Sn\":2.96646,\"Se-Sr\":3.58295,\"Se-Ta\":2.66642,\"Se-Tb\":3.23649,\"Se-Te\":2.57998,\"Se-Th\":3.36649,\"Se-Ti\":2.53642,\"Se-Tl\":3.00644,\"Se-Tm\":3.18649,\"Se-U\":3.00644,\"Se-V\":2.48642,\"Se-W\":2.55998,\"Se-Y\":3.21649,\"Se-Yb\":3.16649,\"Se-Zn\":2.93293,\"Se-Zr\":3.03004,\"Si-Si\":2.6,\"Si-Te\":2.84002,\"Sm-Te\":3.46649,\"Sn-Te\":2.91642,\"Sr-Te\":3.73295,\"Ta-Te\":2.85642,\"Tb-Te\":3.42649,\"Te-Te\":2.80998,\"Te-Th\":3.54649,\"Te-Ti\":2.95642,\"Te-Tl\":3.23644,\"Te-Tm\":3.37649,\"Te-U\":3.16644,\"Te-V\":2.72642,\"Te-W\":2.75998,\"Te-Y\":3.40649,\"Te-Yb\":3.36649,\"Te-Zn\":3.16293,\"Te-Zr\":3.17004}")
-}.pairs, At = .4;
-function jt(e, t) {
-	return kt[e < t ? `${e}-${t}` : `${t}-${e}`] ?? null;
+}.pairs, kt = .4;
+function At(e, t) {
+	return Ot[e < t ? `${e}-${t}` : `${t}-${e}`] ?? null;
 }
-function Mt(e) {
+function jt(e) {
 	let t = /* @__PURE__ */ new Map(), n = e.map((e) => {
 		let n = t.get(e.element);
 		return n === void 0 && (n = t.size, t.set(e.element, n)), n;
 	}), r = [...t.keys()], i = new Float64Array(r.length * r.length), a = 0;
 	if (r.forEach((e, t) => {
 		r.forEach((n, o) => {
-			let s = jt(e, n) ?? 0;
+			let s = At(e, n) ?? 0;
 			i[t * r.length + o] = s, a = Math.max(a, s);
 		});
-	}), a < At) return [];
+	}), a < kt) return [];
 	let o = e.map((e) => e.position.map((e) => Math.floor(e / a) | 0)), s = [
 		Infinity,
 		Infinity,
@@ -799,9 +799,9 @@ function Mt(e) {
 		m.length = 0;
 		for (let r = -1; r <= 1; r += 1) for (let o = -1; o <= 1; o += 1) for (let p = -1; p <= 1; p += 1) for (let h of f.get(d(s + r, c + o, l + p)) ?? []) {
 			let r = i[u + n[h]];
-			if (h <= a || r < At) continue;
-			let o = Ie(t.position, e[h].position);
-			o >= At && o <= r && m.push({
+			if (h <= a || r < kt) continue;
+			let o = Fe(t.position, e[h].position);
+			o >= kt && o <= r && m.push({
 				start: a,
 				end: h,
 				length: o
@@ -812,11 +812,11 @@ function Mt(e) {
 }
 //#endregion
 //#region src/structure/images.ts
-var Nt = 1e-4;
-function Pt(e, t) {
+var Mt = 1e-4;
+function Nt(e, t) {
 	let n = [];
 	return e.forEach((e, r) => {
-		let i = e.fractional.map((e) => L(e, Nt));
+		let i = e.fractional.map((e) => z(e, Mt));
 		if (n.push({
 			siteIndex: r,
 			fractional: i,
@@ -826,7 +826,7 @@ function Pt(e, t) {
 			0,
 			1,
 			2
-		].filter((e) => i[e] < Nt), o = 1 << a.length;
+		].filter((e) => i[e] < Mt), o = 1 << a.length;
 		for (let e = 1; e < o; e += 1) {
 			let t = [
 				i[0],
@@ -845,16 +845,16 @@ function Pt(e, t) {
 }
 //#endregion
 //#region src/scene/build.ts
-function Ft(e) {
+function Pt(e) {
 	let t = [];
 	for (let n of e.sites) for (let e of n.species) {
-		let n = we(e.element);
+		let n = Se(e.element);
 		t.includes(n) || t.push(n);
 	}
 	return t;
 }
-function It(e, t = {}) {
-	let n = Pt(e.sites, t.imageAtoms ?? !0).map((t) => {
+function Ft(e, t = {}) {
+	let n = Nt(e.sites, t.imageAtoms ?? !0).map((t) => {
 		let n = e.sites[t.siteIndex];
 		if (!n) throw Error(`Site ${t.siteIndex} is missing.`);
 		let r = n.species[0];
@@ -863,7 +863,7 @@ function It(e, t = {}) {
 			siteIndex: t.siteIndex,
 			element: r.element,
 			species: n.species,
-			position: Fe(e.lattice, t.fractional),
+			position: Pe(e.lattice, t.fractional),
 			fractional: t.fractional,
 			isImage: t.isImage
 		};
@@ -871,20 +871,20 @@ function It(e, t = {}) {
 	return {
 		cell: e.lattice.vectors,
 		atoms: n,
-		bonds: Mt(n),
-		elements: Ft(e)
+		bonds: jt(n),
+		elements: Pt(e)
 	};
 }
 //#endregion
 //#region node_modules/three/examples/jsm/controls/TrackballControls.js
-var Lt = { type: "change" }, Rt = { type: "start" }, zt = { type: "end" }, K = 1e-6, q = {
+var It = { type: "change" }, Lt = { type: "start" }, Rt = { type: "end" }, zt = 1e-6, q = {
 	NONE: -1,
 	ROTATE: 0,
 	ZOOM: 1,
 	PAN: 2,
 	TOUCH_ROTATE: 3,
 	TOUCH_ZOOM_PAN: 4
-}, Bt = new A(), J = new A(), Vt = new j(), Ht = new j(), Ut = new j(), Y = new O(), Wt = new j(), Gt = new j(), Kt = new j(), qt = new j(), Jt = class extends l {
+}, Bt = new j(), J = new j(), Vt = new M(), Ht = new M(), Ut = new M(), Y = new A(), Wt = new M(), Gt = new M(), Kt = new M(), qt = new M(), Jt = class extends l {
 	constructor(e, t = null) {
 		super(e, t), this.screen = {
 			left: 0,
@@ -896,10 +896,10 @@ var Lt = { type: "change" }, Rt = { type: "start" }, zt = { type: "end" }, K = 1
 			"KeyS",
 			"KeyD"
 		], this.mouseButtons = {
-			LEFT: x.ROTATE,
-			MIDDLE: x.DOLLY,
-			RIGHT: x.PAN
-		}, this.target = new j(), this.state = q.NONE, this.keyState = q.NONE, this._lastPosition = new j(), this._lastUp = new j(), this._lastZoom = 1, this._touchZoomDistanceStart = 0, this._touchZoomDistanceEnd = 0, this._touchRollAngle = 0, this._touchRollPointerIds = "", this._lastAngle = 0, this._lastRollAngle = 0, this._eye = new j(), this._movePrev = new A(), this._moveCurr = new A(), this._lastAxis = new j(), this._zoomStart = new A(), this._zoomEnd = new A(), this._panStart = new A(), this._panEnd = new A(), this._pointers = [], this._pointerPositions = {}, this._onPointerMove = Xt.bind(this), this._onPointerDown = Yt.bind(this), this._onPointerUp = Zt.bind(this), this._onPointerCancel = Qt.bind(this), this._onContextMenu = on.bind(this), this._onMouseWheel = an.bind(this), this._onKeyDown = en.bind(this), this._onKeyUp = $t.bind(this), this._onTouchStart = sn.bind(this), this._onTouchMove = cn.bind(this), this._onTouchEnd = ln.bind(this), this._onMouseDown = tn.bind(this), this._onMouseMove = nn.bind(this), this._onMouseUp = rn.bind(this), this._target0 = this.target.clone(), this._position0 = this.object.position.clone(), this._up0 = this.object.up.clone(), this._zoom0 = this.object.zoom, t !== null && (this.connect(t), this.handleResize()), this.update();
+			LEFT: S.ROTATE,
+			MIDDLE: S.DOLLY,
+			RIGHT: S.PAN
+		}, this.target = new M(), this.state = q.NONE, this.keyState = q.NONE, this._lastPosition = new M(), this._lastUp = new M(), this._lastZoom = 1, this._touchZoomDistanceStart = 0, this._touchZoomDistanceEnd = 0, this._touchRollAngle = 0, this._touchRollPointerIds = "", this._lastAngle = 0, this._lastRollAngle = 0, this._eye = new M(), this._movePrev = new j(), this._moveCurr = new j(), this._lastAxis = new M(), this._zoomStart = new j(), this._zoomEnd = new j(), this._panStart = new j(), this._panEnd = new j(), this._pointers = [], this._pointerPositions = {}, this._onPointerMove = Xt.bind(this), this._onPointerDown = Yt.bind(this), this._onPointerUp = Zt.bind(this), this._onPointerCancel = Qt.bind(this), this._onContextMenu = on.bind(this), this._onMouseWheel = an.bind(this), this._onKeyDown = en.bind(this), this._onKeyUp = $t.bind(this), this._onTouchStart = sn.bind(this), this._onTouchMove = cn.bind(this), this._onTouchEnd = ln.bind(this), this._onMouseDown = tn.bind(this), this._onMouseMove = nn.bind(this), this._onMouseUp = rn.bind(this), this._target0 = this.target.clone(), this._position0 = this.object.position.clone(), this._up0 = this.object.up.clone(), this._zoom0 = this.object.zoom, t !== null && (this.connect(t), this.handleResize()), this.update();
 	}
 	connect(e) {
 		super.connect(e), window.addEventListener("keydown", this._onKeyDown), window.addEventListener("keyup", this._onKeyUp), this.domElement.addEventListener("pointerdown", this._onPointerDown), this.domElement.addEventListener("pointercancel", this._onPointerCancel), this.domElement.addEventListener("wheel", this._onMouseWheel, { passive: !1 }), this.domElement.addEventListener("contextmenu", this._onContextMenu), this.domElement.style.touchAction = "none";
@@ -917,16 +917,16 @@ var Lt = { type: "change" }, Rt = { type: "start" }, zt = { type: "end" }, K = 1
 	update() {
 		if (this._eye.subVectors(this.object.position, this.target), this.noRotate || (this._rotateCamera(), this.multiTouchRoll === !0 && this._rollCamera()), this.noZoom || this._zoomCamera(), this.noPan || this._panCamera(), this.object.position.addVectors(this.target, this._eye), this.object.isPerspectiveCamera) {
 			this._checkDistances(), this.object.lookAt(this.target);
-			let e = this._lastPosition.distanceToSquared(this.object.position) > K, t = this._lastUp.distanceToSquared(this.object.up) > K;
-			(e || t) && (this.dispatchEvent(Lt), this._lastPosition.copy(this.object.position), this._lastUp.copy(this.object.up));
+			let e = this._lastPosition.distanceToSquared(this.object.position) > zt, t = this._lastUp.distanceToSquared(this.object.up) > zt;
+			(e || t) && (this.dispatchEvent(It), this._lastPosition.copy(this.object.position), this._lastUp.copy(this.object.up));
 		} else if (this.object.isOrthographicCamera) {
 			this.object.lookAt(this.target);
-			let e = this._lastPosition.distanceToSquared(this.object.position) > K, t = this._lastUp.distanceToSquared(this.object.up) > K, n = this._lastZoom !== this.object.zoom;
-			(e || t || n) && (this.dispatchEvent(Lt), this._lastPosition.copy(this.object.position), this._lastUp.copy(this.object.up), this._lastZoom = this.object.zoom);
+			let e = this._lastPosition.distanceToSquared(this.object.position) > zt, t = this._lastUp.distanceToSquared(this.object.up) > zt, n = this._lastZoom !== this.object.zoom;
+			(e || t || n) && (this.dispatchEvent(It), this._lastPosition.copy(this.object.position), this._lastUp.copy(this.object.up), this._lastZoom = this.object.zoom);
 		} else console.warn("THREE.TrackballControls: Unsupported camera type.");
 	}
 	reset() {
-		this.state = q.NONE, this.keyState = q.NONE, this.target.copy(this._target0), this.object.position.copy(this._position0), this.object.up.copy(this._up0), this.object.zoom = this._zoom0, this.object.updateProjectionMatrix(), this._eye.subVectors(this.object.position, this.target), this.object.lookAt(this.target), this.dispatchEvent(Lt), this._lastPosition.copy(this.object.position), this._lastUp.copy(this.object.up), this._lastZoom = this.object.zoom;
+		this.state = q.NONE, this.keyState = q.NONE, this.target.copy(this._target0), this.object.position.copy(this._position0), this.object.up.copy(this._up0), this.object.zoom = this._zoom0, this.object.updateProjectionMatrix(), this._eye.subVectors(this.object.position, this.target), this.object.lookAt(this.target), this.dispatchEvent(It), this._lastPosition.copy(this.object.position), this._lastUp.copy(this.object.up), this._lastZoom = this.object.zoom;
 	}
 	_panCamera() {
 		if (J.copy(this._panEnd).sub(this._panStart), J.lengthSq()) {
@@ -959,7 +959,7 @@ var Lt = { type: "change" }, Rt = { type: "start" }, zt = { type: "end" }, K = 1
 	}
 	_zoomCamera() {
 		let e;
-		this.state === q.TOUCH_ZOOM_PAN ? (e = this._touchZoomDistanceStart / this._touchZoomDistanceEnd, this._touchZoomDistanceStart = this._touchZoomDistanceEnd, this.object.isPerspectiveCamera ? this._eye.multiplyScalar(e) : this.object.isOrthographicCamera ? (this.object.zoom = S.clamp(this.object.zoom / e, this.minZoom, this.maxZoom), this._lastZoom !== this.object.zoom && this.object.updateProjectionMatrix()) : console.warn("THREE.TrackballControls: Unsupported camera type")) : (e = 1 + (this._zoomEnd.y - this._zoomStart.y) * this.zoomSpeed, e !== 1 && e > 0 && (this.object.isPerspectiveCamera ? this._eye.multiplyScalar(e) : this.object.isOrthographicCamera ? (this.object.zoom = S.clamp(this.object.zoom / e, this.minZoom, this.maxZoom), this._lastZoom !== this.object.zoom && this.object.updateProjectionMatrix()) : console.warn("THREE.TrackballControls: Unsupported camera type")), this.staticMoving ? this._zoomStart.copy(this._zoomEnd) : this._zoomStart.y += (this._zoomEnd.y - this._zoomStart.y) * this.dynamicDampingFactor);
+		this.state === q.TOUCH_ZOOM_PAN ? (e = this._touchZoomDistanceStart / this._touchZoomDistanceEnd, this._touchZoomDistanceStart = this._touchZoomDistanceEnd, this.object.isPerspectiveCamera ? this._eye.multiplyScalar(e) : this.object.isOrthographicCamera ? (this.object.zoom = ne.clamp(this.object.zoom / e, this.minZoom, this.maxZoom), this._lastZoom !== this.object.zoom && this.object.updateProjectionMatrix()) : console.warn("THREE.TrackballControls: Unsupported camera type")) : (e = 1 + (this._zoomEnd.y - this._zoomStart.y) * this.zoomSpeed, e !== 1 && e > 0 && (this.object.isPerspectiveCamera ? this._eye.multiplyScalar(e) : this.object.isOrthographicCamera ? (this.object.zoom = ne.clamp(this.object.zoom / e, this.minZoom, this.maxZoom), this._lastZoom !== this.object.zoom && this.object.updateProjectionMatrix()) : console.warn("THREE.TrackballControls: Unsupported camera type")), this.staticMoving ? this._zoomStart.copy(this._zoomEnd) : this._zoomStart.y += (this._zoomEnd.y - this._zoomStart.y) * this.dynamicDampingFactor);
 	}
 	_getMouseOnScreen(e, t) {
 		return Bt.set((e - this.screen.left) / this.screen.width, (t - this.screen.top) / this.screen.height), Bt;
@@ -979,7 +979,7 @@ var Lt = { type: "change" }, Rt = { type: "start" }, zt = { type: "end" }, K = 1
 	}
 	_trackPointer(e) {
 		let t = this._pointerPositions[e.pointerId];
-		t === void 0 && (t = new A(), this._pointerPositions[e.pointerId] = t), t.set(e.pageX, e.pageY);
+		t === void 0 && (t = new j(), this._pointerPositions[e.pointerId] = t), t.set(e.pageX, e.pageY);
 	}
 	_getSecondPointerPosition(e) {
 		let t = e.pointerId === this._pointers[0].pointerId ? this._pointers[1] : this._pointers[0];
@@ -1022,26 +1022,26 @@ function tn(e) {
 		default: t = -1;
 	}
 	switch (t) {
-		case x.DOLLY:
+		case S.DOLLY:
 			this.state = q.ZOOM;
 			break;
-		case x.ROTATE:
+		case S.ROTATE:
 			this.state = q.ROTATE;
 			break;
-		case x.PAN:
+		case S.PAN:
 			this.state = q.PAN;
 			break;
 		default: this.state = q.NONE;
 	}
 	let n = this.keyState === q.NONE ? this.state : this.keyState;
-	n === q.ROTATE && !this.noRotate ? (this._moveCurr.copy(this._getMouseOnCircle(e.pageX, e.pageY)), this._movePrev.copy(this._moveCurr)) : n === q.ZOOM && !this.noZoom ? (this._zoomStart.copy(this._getMouseOnScreen(e.pageX, e.pageY)), this._zoomEnd.copy(this._zoomStart)) : n === q.PAN && !this.noPan && (this._panStart.copy(this._getMouseOnScreen(e.pageX, e.pageY)), this._panEnd.copy(this._panStart)), this.dispatchEvent(Rt);
+	n === q.ROTATE && !this.noRotate ? (this._moveCurr.copy(this._getMouseOnCircle(e.pageX, e.pageY)), this._movePrev.copy(this._moveCurr)) : n === q.ZOOM && !this.noZoom ? (this._zoomStart.copy(this._getMouseOnScreen(e.pageX, e.pageY)), this._zoomEnd.copy(this._zoomStart)) : n === q.PAN && !this.noPan && (this._panStart.copy(this._getMouseOnScreen(e.pageX, e.pageY)), this._panEnd.copy(this._panStart)), this.dispatchEvent(Lt);
 }
 function nn(e) {
 	let t = this.keyState === q.NONE ? this.state : this.keyState;
 	t === q.ROTATE && !this.noRotate ? this._moveCurr.copy(this._getMouseOnCircle(e.pageX, e.pageY)) : t === q.ZOOM && !this.noZoom ? this._zoomEnd.copy(this._getMouseOnScreen(e.pageX, e.pageY)) : t === q.PAN && !this.noPan && this._panEnd.copy(this._getMouseOnScreen(e.pageX, e.pageY));
 }
 function rn() {
-	this.state = q.NONE, this.dispatchEvent(zt);
+	this.state = q.NONE, this.dispatchEvent(Rt);
 }
 function an(e) {
 	if (this.enabled !== !1 && this.noZoom !== !0) {
@@ -1054,7 +1054,7 @@ function an(e) {
 				break;
 			default: this._zoomStart.y -= e.deltaY * 25e-5;
 		}
-		this.dispatchEvent(Rt), this.dispatchEvent(zt);
+		this.dispatchEvent(Lt), this.dispatchEvent(Rt);
 	}
 }
 function on(e) {
@@ -1072,7 +1072,7 @@ function sn(e) {
 			let n = (this._pointers[0].pageX + this._pointers[1].pageX) / 2, r = (this._pointers[0].pageY + this._pointers[1].pageY) / 2;
 			this._panStart.copy(this._getMouseOnScreen(n, r)), this._panEnd.copy(this._panStart);
 	}
-	this.dispatchEvent(Rt);
+	this.dispatchEvent(Lt);
 }
 function cn(e) {
 	switch (this._trackPointer(e), this._pointers.length) {
@@ -1102,28 +1102,28 @@ function ln(e) {
 				break;
 			}
 	}
-	this.dispatchEvent(zt);
+	this.dispatchEvent(Rt);
 }
 //#endregion
 //#region src/render/camera.ts
-var un = 1 / 3, dn = 1 / 6, fn = 2, pn = 1e-12, mn = new j(0, 0, -1), hn = new j(0, 1, 0);
+var un = 1 / 3, dn = 1 / 6, fn = 2, pn = 1e-12, mn = new M(0, 0, -1), hn = new M(0, 1, 0);
 function gn(e, t) {
 	let n = e.clone().sub(t.clone().multiplyScalar(e.dot(t)));
 	return n.lengthSq() < pn ? null : n.normalize();
 }
 function _n(e) {
 	for (let t of [
-		new j(0, 0, 1),
-		new j(0, 1, 0),
-		new j(1, 0, 0)
+		new M(0, 0, 1),
+		new M(0, 1, 0),
+		new M(1, 0, 0)
 	]) {
 		let n = gn(t, e);
 		if (n) return n;
 	}
-	return new j(0, 1, 0);
+	return new M(0, 1, 0);
 }
 function vn(e) {
-	let t = new j(...e[0]), n = new j(...e[1]), r = new j(...e[2]).normalize(), i = gn(t, r) ?? gn(n, r) ?? _n(r), a = r.clone().cross(i).normalize(), o = i.clone().add(a.multiplyScalar(un)).add(r.clone().multiplyScalar(dn)).normalize();
+	let t = new M(...e[0]), n = new M(...e[1]), r = new M(...e[2]).normalize(), i = gn(t, r) ?? gn(n, r) ?? _n(r), a = r.clone().cross(i).normalize(), o = i.clone().add(a.multiplyScalar(un)).add(r.clone().multiplyScalar(dn)).normalize();
 	return {
 		outward: o,
 		up: gn(r, o) ?? _n(o)
@@ -1134,7 +1134,7 @@ function yn(e) {
 }
 function bn(e) {
 	let { outward: t, up: n } = vn(e), r = n.clone().cross(t).normalize(), i = t.clone().cross(r).normalize(), a = new C().makeBasis(r, i, t);
-	return new O().setFromRotationMatrix(a).normalize();
+	return new A().setFromRotationMatrix(a).normalize();
 }
 function xn(e) {
 	return {
@@ -1148,7 +1148,7 @@ function xn(e) {
 	};
 }
 function Sn(e, t, n, r) {
-	let i = new O(...t.quaternion).normalize(), a = new j(...t.target), o = mn.clone().applyQuaternion(i);
+	let i = new A(...t.quaternion).normalize(), a = new M(...t.target), o = mn.clone().applyQuaternion(i);
 	e.position.copy(a).sub(o.multiplyScalar(n)), e.up.copy(hn.clone().applyQuaternion(i)), e.quaternion.copy(i), e.near = .01, e.far = Math.max(1e3, n + r * 8), e.updateProjectionMatrix();
 }
 function Cn(e) {
@@ -1160,9 +1160,9 @@ function Cn(e) {
 	];
 }
 function wn(e, t) {
-	let n = new j(1, 0, 0).applyQuaternion(t), r = new j(0, 1, 0).applyQuaternion(t), i = Infinity, a = -Infinity, o = Infinity, s = -Infinity;
+	let n = new M(1, 0, 0).applyQuaternion(t), r = new M(0, 1, 0).applyQuaternion(t), i = Infinity, a = -Infinity, o = Infinity, s = -Infinity;
 	for (let t of [e.min.x, e.max.x]) for (let c of [e.min.y, e.max.y]) for (let l of [e.min.z, e.max.z]) {
-		let e = new j(t, c, l), u = e.dot(n), d = e.dot(r);
+		let e = new M(t, c, l), u = e.dot(n), d = e.dot(r);
 		i = Math.min(i, u), a = Math.max(a, u), o = Math.min(o, d), s = Math.max(s, d);
 	}
 	return {
@@ -1184,17 +1184,17 @@ var Dn = .1, On = .6, kn = .8, An = 2.2, jn = [
 	.2,
 	0
 ], Mn = "#fafafa", Nn = "#444444", Pn = "#f2f2f2", Fn = .3, In = "6.28318530718";
-function X(e, t) {
+function Ln(e, t) {
 	return t.reduce((e, [t, n]) => {
 		if (!e.includes(t)) throw Error(`bravais cannot draw with this three.js: its shader has no ${t}.`);
 		return e.replace(t, n);
 	}, e);
 }
-function Ln(e) {
-	return X(e, [["#include <common>", "#include <common>\nflat varying int vSectorSite;\nflat varying vec3 vSphereCenter;\nflat varying float vSphereRadius;\nflat varying vec2 vDepthCoefficients;\nvarying vec3 vSpherePoint;"], ["#include <project_vertex>", "mat4 sphereMatrix = modelViewMatrix * instanceMatrix;\nvSphereCenter = sphereMatrix[3].xyz;\nvSphereRadius = length(sphereMatrix[0].xyz);\nvSpherePoint = vSphereCenter + vec3(position.xy * vSphereRadius, 0.0);\nvDepthCoefficients = vec2(projectionMatrix[2][2], projectionMatrix[3][2]);\nvSectorSite = gl_InstanceID;\nvec4 mvPosition = vec4(vSpherePoint, 1.0);\ngl_Position = projectionMatrix * mvPosition;"]]);
-}
 function Rn(e) {
-	return X(e, [
+	return Ln(e, [["#include <common>", "#include <common>\nflat varying int vSectorSite;\nflat varying vec3 vSphereCenter;\nflat varying float vSphereRadius;\nflat varying vec2 vDepthCoefficients;\nvarying vec3 vSpherePoint;"], ["#include <project_vertex>", "mat4 sphereMatrix = modelViewMatrix * instanceMatrix;\nvSphereCenter = sphereMatrix[3].xyz;\nvSphereRadius = length(sphereMatrix[0].xyz);\nvSpherePoint = vSphereCenter + vec3(position.xy * vSphereRadius, 0.0);\nvDepthCoefficients = vec2(projectionMatrix[2][2], projectionMatrix[3][2]);\nvSectorSite = gl_InstanceID;\nvec4 mvPosition = vec4(vSpherePoint, 1.0);\ngl_Position = projectionMatrix * mvPosition;"]]);
+}
+function zn(e) {
+	return Ln(e, [
 		["#include <common>", "#include <common>\nuniform sampler2D sectorTexture;\nuniform float vacancyOpacity;\nflat varying int vSectorSite;\nflat varying vec3 vSphereCenter;\nflat varying float vSphereRadius;\nflat varying vec2 vDepthCoefficients;\nvarying vec3 vSpherePoint;"],
 		["#include <clipping_planes_fragment>", "// Where the view ray meets the sphere. The outline is smoothed over one\n// pixel through the alpha: sample coverage, or blending in a see-through layer.\nvec2 sphereOffset = vSpherePoint.xy - vSphereCenter.xy;\nfloat sphereDistance = length(sphereOffset);\nfloat sphereCoverage = clamp(\n  (vSphereRadius - sphereDistance) / max(fwidth(sphereDistance), 1e-6) + 0.5,\n  0.0,\n  1.0\n);\nif (sphereCoverage <= 0.0) {\n  discard;\n}\nfloat sphereHeight = sqrt(\n  max(0.0, vSphereRadius * vSphereRadius - dot(sphereOffset, sphereOffset))\n);\nvec3 sphereNormal = vec3(sphereOffset, sphereHeight) / vSphereRadius;\nfloat sphereZ = vSphereCenter.z + sphereHeight;\ngl_FragDepth = (vDepthCoefficients.x * sphereZ + vDepthCoefficients.y) * 0.5 + 0.5;\n#include <clipping_planes_fragment>"],
 		["#include <color_fragment>", `#include <color_fragment>
@@ -1222,7 +1222,7 @@ function Rn(e) {
 		["#include <normal_fragment_begin>", "#include <normal_fragment_begin>\nnormal = sphereNormal;\nnonPerturbedNormal = normal;"]
 	]);
 }
-function zn(e, t, n) {
+function Bn(e, t, n) {
 	let r = e < 1 || n, i = new T({
 		color: "#ffffff",
 		metalness: Dn,
@@ -1232,21 +1232,21 @@ function zn(e, t, n) {
 		alphaToCoverage: !r
 	});
 	return i.onBeforeCompile = (e) => {
-		e.uniforms.sectorTexture = { value: t }, e.uniforms.vacancyOpacity = { value: Fn }, e.vertexShader = Ln(e.vertexShader), e.fragmentShader = Rn(e.fragmentShader);
+		e.uniforms.sectorTexture = { value: t }, e.uniforms.vacancyOpacity = { value: Fn }, e.vertexShader = Rn(e.vertexShader), e.fragmentShader = zn(e.fragmentShader);
 	}, i.customProgramCacheKey = () => "bravais-atom-spheres", i;
 }
-function Bn(e) {
-	return X(e, [["#include <common>", "#include <common>\nuniform float capRadius;\nflat varying vec3 vTubeCenter;\nflat varying vec3 vTubeAxis;\nflat varying vec3 vAtomCenter;\nflat varying float vTubeRadius;\nflat varying float vTubeHalfLength;\nflat varying vec2 vDepthCoefficients;\nvarying vec3 vBoxPoint;"], ["#include <project_vertex>", "mat4 tubeMatrix = modelViewMatrix * instanceMatrix;\nfloat tubeLength = max(length(tubeMatrix[1].xyz), 1e-6);\nvTubeRadius = length(tubeMatrix[0].xyz);\nvTubeHalfLength = tubeLength * 0.5;\nvTubeCenter = tubeMatrix[3].xyz;\nvTubeAxis = tubeMatrix[1].xyz / tubeLength;\nfloat rim = sqrt(max(0.0, capRadius * capRadius - vTubeRadius * vTubeRadius));\nvAtomCenter = vTubeCenter - vTubeAxis * (rim + vTubeHalfLength);\nvDepthCoefficients = vec2(projectionMatrix[2][2], projectionMatrix[3][2]);\nvec4 mvPosition = tubeMatrix * vec4(position, 1.0);\nvBoxPoint = mvPosition.xyz;\ngl_Position = projectionMatrix * mvPosition;"]]);
-}
 function Vn(e) {
-	return X(e, [
+	return Ln(e, [["#include <common>", "#include <common>\nuniform float capRadius;\nflat varying vec3 vTubeCenter;\nflat varying vec3 vTubeAxis;\nflat varying vec3 vAtomCenter;\nflat varying float vTubeRadius;\nflat varying float vTubeHalfLength;\nflat varying vec2 vDepthCoefficients;\nvarying vec3 vBoxPoint;"], ["#include <project_vertex>", "mat4 tubeMatrix = modelViewMatrix * instanceMatrix;\nfloat tubeLength = max(length(tubeMatrix[1].xyz), 1e-6);\nvTubeRadius = length(tubeMatrix[0].xyz);\nvTubeHalfLength = tubeLength * 0.5;\nvTubeCenter = tubeMatrix[3].xyz;\nvTubeAxis = tubeMatrix[1].xyz / tubeLength;\nfloat rim = sqrt(max(0.0, capRadius * capRadius - vTubeRadius * vTubeRadius));\nvAtomCenter = vTubeCenter - vTubeAxis * (rim + vTubeHalfLength);\nvDepthCoefficients = vec2(projectionMatrix[2][2], projectionMatrix[3][2]);\nvec4 mvPosition = tubeMatrix * vec4(position, 1.0);\nvBoxPoint = mvPosition.xyz;\ngl_Position = projectionMatrix * mvPosition;"]]);
+}
+function Hn(e) {
+	return Ln(e, [
 		["#include <common>", "#include <common>\nuniform float capRadius;\nflat varying vec3 vTubeCenter;\nflat varying vec3 vTubeAxis;\nflat varying vec3 vAtomCenter;\nflat varying float vTubeRadius;\nflat varying float vTubeHalfLength;\nflat varying vec2 vDepthCoefficients;\nvarying vec3 vBoxPoint;"],
 		["#include <clipping_planes_fragment>", "// The view ray through this pixel runs down view z. The tube's surface is\n// the nearest of: the wall outside the atom's sphere, the sphere's far side\n// inside the tube (the cut face) and the midpoint disc. Along the ray,\n// q(z) = p(z) - centre has q.a = k + (z - cz) a.z and |q|^2 = u + (z - cz)^2.\nvec2 tubePixel = vBoxPoint.xy;\nvec2 tubeOffset = tubePixel - vTubeCenter.xy;\nfloat tubeK = dot(tubeOffset, vTubeAxis.xy);\nfloat tubeU = dot(tubeOffset, tubeOffset);\nfloat axisZ = vTubeAxis.z;\nfloat quadA = 1.0 - axisZ * axisZ;\nfloat quadB = -2.0 * tubeK * axisZ;\nfloat quadC = tubeU - tubeK * tubeK - vTubeRadius * vTubeRadius;\nfloat tubeRadius2 = vTubeRadius * vTubeRadius;\n// The ray's closest approach to the axis, for the outline's coverage.\nfloat axisDistance2 = tubeU - tubeK * tubeK;\nif (quadA > 1e-6) {\n  axisDistance2 -= quadB * quadB / (4.0 * quadA);\n}\nfloat axisDistance = sqrt(max(0.0, axisDistance2));\nfloat tubeCoverage = clamp(\n  (vTubeRadius - axisDistance) / max(fwidth(axisDistance), 1e-6) + 0.5,\n  0.0,\n  1.0\n);\nif (tubeCoverage <= 0.0) {\n  discard;\n}\nfloat tubeHitZ = -1e30;\nvec3 tubeNormal = vec3(0.0, 0.0, 1.0);\nif (quadA > 1e-6) {\n  float discriminant = quadB * quadB - 4.0 * quadA * quadC;\n  if (discriminant >= 0.0) {\n    float z = vTubeCenter.z + (-quadB + sqrt(discriminant)) / (2.0 * quadA);\n    vec3 d = vec3(tubePixel, z) - vTubeCenter;\n    float s = dot(d, vTubeAxis);\n    if (abs(s) <= vTubeHalfLength && distance(vec3(tubePixel, z), vAtomCenter) >= capRadius) {\n      tubeHitZ = z;\n      tubeNormal = (d - s * vTubeAxis) / vTubeRadius;\n    }\n  }\n}\nif (axisZ > 1e-6) {\n  float z = vTubeCenter.z + (vTubeHalfLength - tubeK) / axisZ;\n  vec3 d = vec3(tubePixel, z) - vTubeCenter;\n  if (dot(d, d) - vTubeHalfLength * vTubeHalfLength <= tubeRadius2 && z > tubeHitZ &&\n      distance(vec3(tubePixel, z), vAtomCenter) >= capRadius) {\n    tubeHitZ = z;\n    tubeNormal = vTubeAxis;\n  }\n}\n{\n  vec2 m = tubePixel - vAtomCenter.xy;\n  float inside = capRadius * capRadius - dot(m, m);\n  if (inside >= 0.0) {\n    float z = vAtomCenter.z - sqrt(inside);\n    vec3 d = vec3(tubePixel, z) - vTubeCenter;\n    float s = dot(d, vTubeAxis);\n    if (abs(s) <= vTubeHalfLength && dot(d, d) - s * s <= tubeRadius2 && z > tubeHitZ) {\n      tubeHitZ = z;\n      tubeNormal = (vAtomCenter - vec3(tubePixel, z)) / capRadius;\n    }\n  }\n}\nif (tubeHitZ < -1e29) {\n  discard;\n}\ngl_FragDepth = (vDepthCoefficients.x * tubeHitZ + vDepthCoefficients.y) * 0.5 + 0.5;\n#include <clipping_planes_fragment>"],
 		["#include <color_fragment>", "#include <color_fragment>\ndiffuseColor.a *= tubeCoverage;"],
 		["#include <normal_fragment_begin>", "#include <normal_fragment_begin>\nnormal = tubeNormal;\nnonPerturbedNormal = normal;"]
 	]);
 }
-function Hn(e, t) {
+function Un(e, t) {
 	let n = e < 1, r = new T({
 		color: "#ffffff",
 		metalness: Dn,
@@ -1256,18 +1256,18 @@ function Hn(e, t) {
 		alphaToCoverage: !n
 	});
 	return r.onBeforeCompile = (e) => {
-		e.uniforms.capRadius = { value: t }, e.vertexShader = Bn(e.vertexShader), e.fragmentShader = Vn(e.fragmentShader);
+		e.uniforms.capRadius = { value: t }, e.vertexShader = Vn(e.vertexShader), e.fragmentShader = Hn(e.fragmentShader);
 	}, r.customProgramCacheKey = () => "bravais-bond-tubes", r;
 }
-function Un(e) {
-	return new ee({
+function Wn(e) {
+	return new y({
 		color: Nn,
 		opacity: e,
 		transparent: !0,
 		fog: !1
 	});
 }
-function Wn(e, n) {
+function Gn(e, n) {
 	return new T({
 		color: e,
 		metalness: 0,
@@ -1276,19 +1276,16 @@ function Wn(e, n) {
 		stencilWrite: !0,
 		stencilRef: n,
 		stencilFunc: t,
-		stencilZPass: se
+		stencilZPass: ae
 	});
 }
-function Gn(e) {
-	for (let t of Array.isArray(e) ? e : [e]) t.dispose();
-}
 //#endregion
 //#region src/render/gizmo.ts
 var Kn = {
 	a: "#d27686",
 	b: "#80c393",
 	c: "#75a4dd"
-}, qn = .82, Jn = .055, Yn = .24, Xn = .13, Zn = .13, Qn = 1.3, $n = .5, er = 128, tr = "#2a2a2a", nr = "#ffffff", rr = 3.2, ir = new j(0, 1, 0);
+}, qn = .82, Jn = .055, Yn = .24, Xn = .13, Zn = .13, Qn = 1.3, $n = .5, er = 128, tr = "#2a2a2a", nr = "#ffffff", rr = 3.2, ir = new M(0, 1, 0);
 function ar(e, t) {
 	let n = document.createElement("canvas");
 	n.width = er, n.height = er;
@@ -1298,19 +1295,17 @@ function ar(e, t) {
 		r.clearRect(0, 0, n.width, n.height), r.font = `italic 600 ${i}px ${t}`, r.textAlign = "center", r.textBaseline = "middle", r.fillStyle = "#ffffff", r.fillText(e, n.width / 2, n.height / 2);
 	}
 	let i = new o(n);
-	return i.colorSpace = ce, i.premultiplyAlpha = !0, i.magFilter = te, i.minFilter = ne, i;
-}
-function or(e, t, n, r) {
-	let i = new v();
-	i.quaternion.setFromUnitVectors(ir, e.clone().normalize());
-	let a = Wn(t, 1), o = new w(new u(Jn, Jn, qn, 24), a);
-	o.position.y = qn / 2;
-	let s = new w(new c(Xn, Yn, 32), a);
-	s.position.y = .94, i.add(o, s);
-	let l = ar(n, r), d = [[tr, re], [nr, m]];
-	for (let [e, t] of d) {
-		let n = new ue(new de({
-			map: l,
+	return i.colorSpace = oe, i.premultiplyAlpha = !0, i.magFilter = x, i.minFilter = te, i;
+}
+function or(e, t, n) {
+	let r = new _(), i = Gn(e, 1), a = new w(new u(Jn, Jn, qn, 24), i);
+	a.position.y = qn / 2;
+	let o = new w(new c(Xn, Yn, 32), i);
+	o.position.y = .94, r.add(a, o);
+	let s = ar(t, n), l = [[tr, D], [nr, m]];
+	for (let [e, t] of l) {
+		let n = new le(new ue({
+			map: s,
 			color: e,
 			depthTest: !1,
 			depthWrite: !1,
@@ -1319,30 +1314,31 @@ function or(e, t, n, r) {
 			stencilRef: 1,
 			stencilFunc: t
 		}));
-		n.scale.set($n, $n, 1), n.position.y = Qn, n.renderOrder = 10, i.add(n);
+		n.scale.set($n, $n, 1), n.position.y = Qn, n.renderOrder = 10, r.add(n);
 	}
-	return i;
+	return r;
 }
 var sr = class {
-	fontFamily;
-	scene = new k();
-	camera = new ae(-1, 1, 1, -1, .1, 20);
-	axes = new v();
+	scene = new se();
+	camera = new k(-1, 1, 1, -1, .1, 20);
+	axes = new _();
+	arrows;
 	headlight = new f("#ffffff", 1.78);
-	inverse = new O();
+	inverse = new A();
+	direction = new M();
 	constructor(e) {
-		this.fontFamily = e, this.camera.position.set(0, 0, 5), this.camera.lookAt(0, 0, 0), this.scene.add(new n("#ffffff", .68)), this.headlight.position.set(1.6, 1.1, 5), this.scene.add(this.headlight), this.scene.add(this.axes);
+		this.camera.position.set(0, 0, 5), this.camera.lookAt(0, 0, 0), this.scene.add(new n("#ffffff", .68)), this.headlight.position.set(1.6, 1.1, 5), this.scene.add(this.headlight), this.arrows = [
+			"a",
+			"b",
+			"c"
+		].map((t) => or(Kn[t], t, e));
+		let t = new w(new ce(Zn, 32, 16), Gn("#f3f2ee", 0));
+		this.axes.add(...this.arrows, t), this.scene.add(this.axes);
 	}
 	setCell(e) {
-		this.axes.clear();
-		let t = [
-			["a", new j(...e[0])],
-			["b", new j(...e[1])],
-			["c", new j(...e[2])]
-		];
-		for (let [e, n] of t) this.axes.add(or(n, Kn[e], e, this.fontFamily));
-		let n = new w(new le(Zn, 32, 16), Wn("#f3f2ee", 0));
-		this.axes.add(n);
+		this.arrows.forEach((t, n) => {
+			t.quaternion.setFromUnitVectors(ir, this.direction.set(...e[n]).normalize());
+		});
 	}
 	render(e, t, n) {
 		this.axes.quaternion.copy(this.inverse.copy(t).invert());
@@ -1353,7 +1349,7 @@ var sr = class {
 	}
 	dispose() {
 		this.axes.traverse((e) => {
-			if (e instanceof w || e instanceof ue) {
+			if (e instanceof w || e instanceof le) {
 				e.geometry.dispose();
 				let t = e.material;
 				"map" in t && t.map && t.map.dispose(), t.dispose();
@@ -1532,7 +1528,7 @@ var vr = cr.elements, yr = {
 	Cr: "#6d87c9",
 	V: "#bf4b62",
 	F: "#9bd7d0"
-}, br = /* @__PURE__ */ "O.Li.Na.K.Mg.Ca.Sr.Ba.Al.Si.Ti.Mn.Fe.Co.Ni.Cu.Zn.V.Cr.Ga.Ge.Y.Zr.Nb.Mo.Ru.Ag.Cd.In.Sn.Sb.La.Ce.Nd.Sm.Gd.Ta.W.Pt.Au.Pb.Bi.H.C.N.F.P.S.Cl.Se.Te.I.Br.B.Sc.Rb.Cs.Hf.Re.Ir.Pd.Rh.Os.Pr.Eu.Tb.Dy.Ho.Er.Tm.Yb.Lu.Th.U".split("."), xr = new Map(br.map((e, t) => [e, t])), Sr = .1, Cr = .7, Z = {
+}, br = /* @__PURE__ */ "O.Li.Na.K.Mg.Ca.Sr.Ba.Al.Si.Ti.Mn.Fe.Co.Ni.Cu.Zn.V.Cr.Ga.Ge.Y.Zr.Nb.Mo.Ru.Ag.Cd.In.Sn.Sb.La.Ce.Nd.Sm.Gd.Ta.W.Pt.Au.Pb.Bi.H.C.N.F.P.S.Cl.Se.Te.I.Br.B.Sc.Rb.Cs.Hf.Re.Ir.Pd.Rh.Os.Pr.Eu.Tb.Dy.Ho.Er.Tm.Yb.Lu.Th.U".split("."), xr = new Map(br.map((e, t) => [e, t])), Sr = .1, Cr = .7, X = {
 	minL: .45,
 	maxL: .88,
 	maxC: .2
@@ -1601,7 +1597,7 @@ function Rr(e) {
 function zr(e) {
 	let t = Rr(e);
 	if (t) return hr(dr(t));
-	let n = I.indexOf(e) + 1;
+	let n = _e.indexOf(e) + 1;
 	return {
 		L: .68,
 		C: .12,
@@ -1611,13 +1607,13 @@ function zr(e) {
 function Br(e) {
 	let t = xr.get(e);
 	if (t !== void 0) return t;
-	let n = I.indexOf(e);
-	return br.length + (n >= 0 ? n : I.length);
+	let n = _e.indexOf(e);
+	return br.length + (n >= 0 ? n : _e.length);
 }
 function Vr(e) {
 	return _r({
-		L: Math.min(Z.maxL, Math.max(Z.minL, e.L)),
-		C: Math.min(Z.maxC, Math.max(0, e.C)),
+		L: Math.min(X.maxL, Math.max(X.minL, e.L)),
+		C: Math.min(X.maxC, Math.max(0, e.C)),
 		h: (e.h % 360 + 360) % 360
 	});
 }
@@ -1655,7 +1651,7 @@ function Wr(e) {
 	}
 	let n = e.C / Math.max(1e-6, Hr(e.L, e.h));
 	for (let r of Dr) for (let i of Or) for (let a of kr) {
-		let o = Math.min(Z.maxL, Math.max(Z.minL, e.L + i)), s = e.h + r, c = Math.min(1, Math.max(.1, n + a)) * Hr(o, s), l = Math.abs(r) / jr + Math.abs(i) / Mr + Math.abs(a) / Nr;
+		let o = Math.min(X.maxL, Math.max(X.minL, e.L + i)), s = e.h + r, c = Math.min(1, Math.max(.1, n + a)) * Hr(o, s), l = Math.abs(r) / jr + Math.abs(i) / Mr + Math.abs(a) / Nr;
 		t.push(Ur({
 			L: o,
 			C: c,
@@ -1664,9 +1660,9 @@ function Wr(e) {
 	}
 	return t.sort((e, t) => e.cost - t.cost);
 }
-var Q = /* @__PURE__ */ new Map();
+var Z = /* @__PURE__ */ new Map();
 function Gr(e, t) {
-	let n = [...Q.values()].map((e) => e.lab), r = (e) => n.reduce((t, n) => Math.min(t, Ir(e, n)), Infinity), i = gr(t);
+	let n = [...Z.values()].map((e) => e.lab), r = (e) => n.reduce((t, n) => Math.min(t, Ir(e, n)), Infinity), i = gr(t);
 	if (r(i) >= Sr) return {
 		hex: e ?? mr(i),
 		lab: i
@@ -1689,18 +1685,18 @@ function Gr(e, t) {
 	};
 }
 function Kr(e) {
-	let t = [...new Set(e)].filter((e) => !Q.has(e)).sort((e, t) => Br(e) - Br(t));
-	for (let e of t) Q.set(e, Gr(Rr(e), zr(e)));
+	let t = [...new Set(e)].filter((e) => !Z.has(e)).sort((e, t) => Br(e) - Br(t));
+	for (let e of t) Z.set(e, Gr(Rr(e), zr(e)));
 }
 function qr(e) {
-	let t = Q.get(e);
-	return t ? t.hex : (Kr([e]), Q.get(e).hex);
+	let t = Z.get(e);
+	return t ? t.hex : (Kr([e]), Z.get(e).hex);
 }
 function Jr() {
-	return [...Q].map(([e, t]) => [e, t.hex]);
+	return [...Z].map(([e, t]) => [e, t.hex]);
 }
 function Yr() {
-	Q.clear();
+	Z.clear();
 }
 //#endregion
 //#region src/render/legend.ts
@@ -1730,11 +1726,11 @@ var ei = class {
 		}
 		this.element.hidden = e.length === 0;
 	}
-}, $ = .4, ti = .1, ni = .398, ri = Math.sqrt(ni * ni - ti * ti), ii = new j(0, 1, 0);
-function ai(e) {
-	let [t, n, r] = e.map((e) => new j(...e));
+}, Q = .4, $ = .1, ti = .398, ni = Math.sqrt(ti * ti - $ * $), ri = new M(0, 1, 0);
+function ii(e) {
+	let [t, n, r] = e.map((e) => new M(...e));
 	return [
-		new j(),
+		new M(),
 		t.clone(),
 		n.clone(),
 		r.clone(),
@@ -1744,7 +1740,7 @@ function ai(e) {
 		t.clone().add(n).add(r)
 	];
 }
-function oi(e) {
+function ai(e) {
 	let [t, n, r] = e;
 	return [
 		(t[0] + n[0] + r[0]) / 2,
@@ -1752,7 +1748,7 @@ function oi(e) {
 		(t[2] + n[2] + r[2]) / 2
 	];
 }
-var si = [
+var oi = [
 	[0, 1],
 	[0, 2],
 	[0, 3],
@@ -1766,85 +1762,85 @@ var si = [
 	[5, 7],
 	[6, 7]
 ];
-function ci(e, t) {
-	let n = ai(e), r = [];
-	for (let [e, t] of si) {
+function si(e, t) {
+	let n = ii(e), r = [];
+	for (let [e, t] of oi) {
 		let i = n[e], a = n[t];
 		r.push(i.x, i.y, i.z, a.x, a.y, a.z);
 	}
 	let i = new a();
-	return i.setAttribute("position", new h(r, 3)), new b(i, Un(t));
+	return i.setAttribute("position", new h(r, 3)), new b(i, Wn(t));
 }
-function li(e) {
+function ci(e) {
 	let t = 2048, n = Math.max(1, Math.ceil(e.length / 256)), r = new Float32Array(t * n * 4).fill(1), i = new s(Pn), a = new s(), o = !1;
 	e.forEach((e, n) => {
 		let s = Math.floor(n / 256), c = n % 256 * 8, l = (s * t + c) * 4, u = 0, d = Math.min(1, e.species[0]?.occupancy ?? 1);
 		for (let t of e.species.slice(1)) u < 7 && (a.set(qr(t.element)).toArray(r, l), r[l + 3] = d, l += 4, u += 1), d = Math.min(1, d + t.occupancy);
 		d < 1 - 1e-6 && (i.toArray(r, l), r[l + 3] = d + 2, o = !0);
 	});
-	let c = new d(r, t, n, oe, g);
+	let c = new d(r, t, n, ie, g);
 	return c.magFilter = E, c.minFilter = E, c.generateMipmaps = !1, c.needsUpdate = !0, {
 		texture: c,
 		hasVacancies: o
 	};
 }
-function ui() {
+function li() {
 	let e = /* @__PURE__ */ new Map();
 	return (t) => {
 		let n = e.get(t);
 		return n || (n = new s(qr(t)), e.set(t, n)), n;
 	};
 }
-function di(e, t) {
-	let n = new D(2, 2), r = li(e), i = new y(n, zn(t, r.texture, r.hasVacancies), e.length), a = new C(), o = ui();
+function ui(e, t, n = li()) {
+	let r = new re(2, 2), i = ci(e), a = new v(r, Bn(t, i.texture, i.hasVacancies), e.length), o = new C();
 	return e.forEach((e, t) => {
-		a.makeScale($, $, $), a.setPosition(e.position[0], e.position[1], e.position[2]), i.setMatrixAt(t, a), i.setColorAt(t, o(e.element));
-	}), i.instanceMatrix.needsUpdate = !0, i.instanceColor && (i.instanceColor.needsUpdate = !0), i.userData.sectorTexture = r.texture, i;
+		o.makeScale(Q, Q, Q), o.setPosition(e.position[0], e.position[1], e.position[2]), a.setMatrixAt(t, o), a.setColorAt(t, n(e.element));
+	}), a.instanceMatrix.needsUpdate = !0, a.instanceColor && (a.instanceColor.needsUpdate = !0), a.userData.sectorTexture = i.texture, a;
 }
-function fi(e, t) {
+function di(e, t, n = li()) {
 	if (e.bonds.length === 0) return null;
-	let n = new i(2, 1, 2), r = new y(n, Hn(t, ni), e.bonds.length * 2), a = new ie(), o = new j(), s = new j(), c = new j(), l = new O(), u = ui();
-	return e.bonds.forEach((t, n) => {
+	let r = new i(2, 1, 2), a = new v(r, Un(t, ti), e.bonds.length * 2), o = new O(), s = new M(), c = new M(), l = new M(), u = new A();
+	return e.bonds.forEach((t, r) => {
 		let i = e.atoms[t.start], d = e.atoms[t.end];
-		o.set(...i.position), s.set(...d.position), c.subVectors(s, o);
-		let f = c.length() / 2;
-		c.normalize(), l.setFromUnitVectors(ii, c);
+		s.set(...i.position), c.set(...d.position), l.subVectors(c, s);
+		let f = l.length() / 2;
+		l.normalize(), u.setFromUnitVectors(ri, l);
 		let p = [[
 			0,
 			i,
-			ri,
+			ni,
 			f
 		], [
 			1,
 			d,
 			f,
-			2 * f - ri
+			2 * f - ni
 		]];
-		for (let [e, t, i, s] of p) {
-			let d = Math.max(0, s - i);
-			a.position.copy(o).addScaledVector(c, i + d / 2), a.quaternion.copy(l), a.scale.set(ti, d, ti), a.updateMatrix(), r.setMatrixAt(n * 2 + e, a.matrix), r.setColorAt(n * 2 + e, u(t.element));
+		for (let [e, t, i, c] of p) {
+			let d = Math.max(0, c - i);
+			o.position.copy(s).addScaledVector(l, i + d / 2), o.quaternion.copy(u), o.scale.set($, d, $), o.updateMatrix(), a.setMatrixAt(r * 2 + e, o.matrix), a.setColorAt(r * 2 + e, n(t.element));
 		}
-	}), r.instanceMatrix.needsUpdate = !0, r.instanceColor && (r.instanceColor.needsUpdate = !0), r;
+	}), a.instanceMatrix.needsUpdate = !0, a.instanceColor && (a.instanceColor.needsUpdate = !0), a;
 }
-function pi(e) {
+function fi(e) {
 	let t = e[0];
 	if (!t) throw Error("A scene model needs at least one layer.");
 	Kr(e.flatMap((e) => [...e.scene.elements, ...e.scene.atoms.map((e) => e.element)]));
-	let n = oi(t.scene.cell), i = new v();
+	let n = ai(t.scene.cell), i = new _();
 	i.position.set(-n[0], -n[1], -n[2]);
-	let a = new r(), o = [], s = 1, c = [], l = new j(...n), u = new j(), f = new j();
+	let a = new r(), o = [], s = 1, c = [], l = new M(...n), u = new M(), d = new M();
 	for (let t of e) {
-		for (let e of ai(t.scene.cell)) a.expandByPoint(e.sub(l));
-		for (let e of t.scene.atoms) u.set(...e.position).sub(l), a.expandByPoint(f.copy(u).addScalar($)), a.expandByPoint(f.copy(u).subScalar($));
+		for (let e of ii(t.scene.cell)) a.expandByPoint(e.sub(l));
+		for (let e of t.scene.atoms) u.set(...e.position).sub(l), a.expandByPoint(d.copy(u).addScalar(Q)), a.expandByPoint(d.copy(u).subScalar(Q));
 		for (let e of t.scene.cell) s = Math.max(s, Math.hypot(...e));
 		let e = t.opacity ?? 1;
 		if (e <= 0) continue;
-		let n = di(t.scene.atoms, e);
+		let n = ui(t.scene.atoms, e);
 		i.add(n), c.push(n);
-		let r = fi(t.scene, e);
+		let r = di(t.scene, e);
 		r && (i.add(r), c.push(r));
-		let d = ci(t.scene.cell, e);
-		i.add(d), c.push(d);
+		let f = si(t.scene.cell, e);
+		i.add(f), c.push(f);
 		for (let e of t.scene.elements) o.includes(e) || o.push(e);
 	}
 	return {
@@ -1855,27 +1851,83 @@ function pi(e) {
 		span: s,
 		elements: o,
 		dispose() {
-			for (let e of c) {
-				e.geometry.dispose(), Gn(e.material);
-				let t = e.userData.sectorTexture;
-				t instanceof d && t.dispose();
-			}
-			i.clear();
+			pi(c), i.clear();
 		}
 	};
 }
+function pi(e) {
+	for (let t of e) {
+		t.geometry.dispose();
+		let e = t.userData.sectorTexture;
+		e instanceof d && e.dispose();
+	}
+}
+var mi = {
+	cell: [
+		[
+			1,
+			0,
+			0
+		],
+		[
+			0,
+			1,
+			0
+		],
+		[
+			0,
+			0,
+			1
+		]
+	],
+	atoms: [0, 2].map((e, t) => ({
+		siteIndex: t,
+		element: "X",
+		species: [{
+			element: "X",
+			occupancy: 1
+		}],
+		position: [
+			0,
+			0,
+			e
+		],
+		fractional: [
+			0,
+			0,
+			e
+		],
+		isImage: !1
+	})),
+	bonds: [{
+		start: 0,
+		end: 1,
+		length: 2
+	}],
+	elements: []
+};
+function hi() {
+	let e = new s("#ffffff"), t = [];
+	for (let n of [1, .5]) t.push(ui(mi.atoms, n, () => e), di(mi, n, () => e));
+	t.push(si(mi.cell, 1));
+	let n = new _();
+	return n.add(...t), {
+		group: n,
+		dispose: () => pi(t)
+	};
+}
 //#endregion
 //#region src/render/stage.ts
-function mi(t = {}) {
-	let r = new fe({
+function gi(t = {}) {
+	let r = new N({
 		antialias: !0,
 		alpha: !0,
 		stencil: t.stencil ?? !1
 	});
 	r.setClearColor(Mn, 1), r.toneMapping = e;
-	let i = new k(), a = new ae(-1, 1, 1, -1, .01, 1e3), o = new f("#ffffff", An);
+	let i = new se(), a = new k(-1, 1, 1, -1, .01, 1e3), o = new f("#ffffff", An);
 	i.add(new n("#ffffff", kn), o, o.target);
-	let s = new j();
+	let s = new M();
 	return {
 		renderer: r,
 		scene: i,
@@ -1888,76 +1940,86 @@ function mi(t = {}) {
 }
 //#endregion
 //#region src/render/viewer.ts
-var hi = .4, gi = .4, _i = .4, vi = 2, yi = 1.2, bi = .2, xi = 1500, Si = .2, Ci = 5, wi = 3, Ti = 1e-5, Ei = 4, Di = .35, Oi = 110, ki = 180;
-function Ai(e, t, n) {
-	let r = new j(0, 0, 1).applyQuaternion(n), i = 0, a = 0, o = new j();
+var _i = .4, vi = .4, yi = .4, bi = 2, xi = 1.2, Si = .2, Ci = 1500, wi = .2, Ti = 5, Ei = 3, Di = 1e-5, Oi = 4, ki = .35, Ai = 110, ji = 180;
+function Mi(e, t, n) {
+	let r = new M(0, 0, 1).applyQuaternion(n), i = 0, a = 0, o = new M();
 	for (let t of [e.bounds.min.x, e.bounds.max.x]) for (let n of [e.bounds.min.y, e.bounds.max.y]) for (let s of [e.bounds.min.z, e.bounds.max.z]) {
 		let e = o.set(t, n, s).dot(r);
 		i = Math.max(i, e), a = Math.min(a, e);
 	}
-	let s = Math.max(1, e.span), c = Math.max(.01 * s, -a), l = Math.min(c, -i), u = s * _i, d = l - u, f = t + (d + (Math.max(d, c - u) - d) * gi), p = f + (t + c - f) / hi;
-	return new _(Mn, f, p);
+	let s = Math.max(1, e.span), c = Math.max(.01 * s, -a), l = Math.min(c, -i), u = s * yi, d = l - u, f = t + (d + (Math.max(d, c - u) - d) * vi), p = f + (t + c - f) / _i;
+	return new ee(Mn, f, p);
 }
-function ji(e, t = {}) {
+function Ni(e, t = {}) {
 	let n = e.ownerDocument;
 	Qr(n), e.classList.add("bravais-viewer");
-	let { renderer: r, scene: i, camera: a, updateHeadlight: o } = mi({ stencil: !0 });
+	let { renderer: r, scene: i, camera: a, updateHeadlight: o } = gi({ stencil: !0 });
 	r.setPixelRatio(n.defaultView?.devicePixelRatio ?? 1), r.domElement.className = "bravais-canvas", e.appendChild(r.domElement);
 	let s = t.legend === !1 ? null : new ei(n);
 	s && e.appendChild(s.element);
 	let c = n.defaultView?.getComputedStyle(e).fontFamily ?? "sans-serif", l = t.gizmo === !1 ? null : new sr(c), u = new Jt(a, r.domElement);
-	u.rotateSpeed = vi, u.zoomSpeed = yi, u.dynamicDampingFactor = bi, u.staticMoving = !1, u.noPan = !0, u.mouseButtons.LEFT = x.ROTATE, u.mouseButtons.MIDDLE = x.DOLLY, u.mouseButtons.RIGHT = null;
-	let d = null, f = 4, p = 1, m = 1, h = 1, g = null, _ = !1, v = 0, y = !1, ee = 0, b = !1, te = new O(), ne = 0, S = new j(), C = /* @__PURE__ */ new Set();
-	function w() {
-		ee = 0, b = !1;
-		for (let e of C) e();
-	}
-	function T() {
-		if (o(), r.setViewport(0, 0, m, h), r.render(i, a), l && d) {
-			let e = Math.round(Math.min(ki, Math.max(Oi, Math.min(m, h) * Di)));
+	u.rotateSpeed = bi, u.zoomSpeed = xi, u.dynamicDampingFactor = Si, u.staticMoving = !1, u.noPan = !0, u.mouseButtons.LEFT = S.ROTATE, u.mouseButtons.MIDDLE = S.DOLLY, u.mouseButtons.RIGHT = null;
+	let d = null, f = null, p = 4, m = 1, h = 1, g = 1, _ = null, v = !1, y = 0, b = !1, x = 0, te = !1, ne = new j(), C = new A(), w = 0, T = new M(), E = /* @__PURE__ */ new Set();
+	function D() {
+		x = 0, te = !1;
+		for (let e of E) e();
+	}
+	function O() {
+		if (!f && (o(), r.setViewport(0, 0, h, g), r.render(i, a), l && d)) {
+			let e = Math.round(Math.min(ji, Math.max(Ai, Math.min(h, g) * ki)));
 			l.render(r, a.quaternion, e);
 		}
 	}
-	function E() {
+	function k() {
 		let t = e.getBoundingClientRect();
-		m = Math.max(1, Math.round(t.width)), h = Math.max(1, Math.round(t.height)), r.setSize(m, h, !1), En(a, m, h), u.handleResize(), u.rotateSpeed = m / xi * vi;
+		h = Math.max(1, Math.round(t.width)), g = Math.max(1, Math.round(t.height));
+		let n = r.getSize(ne);
+		(n.x !== h || n.y !== g) && r.setSize(h, g, !1), En(a, h, g), u.handleResize(), u.rotateSpeed = h / Ci * bi;
 	}
 	function re(e) {
-		d && (p = Tn(d.bounds, e, m, h), u.minZoom = p * Si, u.maxZoom = p * Ci);
+		d && (m = Tn(d.bounds, e, h, g), u.minZoom = m * wi, u.maxZoom = m * Ti);
 	}
 	function ie() {
-		let e = te.angleTo(a.quaternion) > Ti || Math.abs(ne - a.zoom) > Ti * p || S.distanceToSquared(u.target) > Ti;
-		return te.copy(a.quaternion), ne = a.zoom, S.copy(u.target), e;
+		let e = C.angleTo(a.quaternion) > Di || Math.abs(w - a.zoom) > Di * m || T.distanceToSquared(u.target) > Di;
+		return C.copy(a.quaternion), w = a.zoom, T.copy(u.target), e;
 	}
 	function ae() {
-		g = null, u.update(), ie() ? (v = 0, T(), ee += 1, b = !0, ee >= Ei && w()) : v += 1, y || v < wi ? g = requestAnimationFrame(ae) : (_ = !1, b && w());
-	}
-	function D() {
-		v = 0, _ || (_ = !0, g = requestAnimationFrame(ae));
+		_ = null, u.update(), ie() ? (y = 0, O(), x += 1, te = !0, x >= Oi && D()) : y += 1, b || y < Ei ? _ = requestAnimationFrame(ae) : (v = !1, te && D());
 	}
 	function oe() {
-		g !== null && (cancelAnimationFrame(g), g = null), _ = !1;
+		y = 0, v || (v = !0, _ = requestAnimationFrame(ae));
 	}
-	u.addEventListener("start", () => {
-		y = !0, D();
-	}), u.addEventListener("end", () => {
-		y = !1, D();
-	}), u.addEventListener("change", D);
 	function se() {
-		let e = p > 0 ? a.zoom / p : 1;
-		E(), re(a.quaternion), a.zoom = p * e, a.updateProjectionMatrix(), T();
+		_ !== null && (cancelAnimationFrame(_), _ = null), v = !1;
 	}
-	let ce = typeof ResizeObserver > "u" ? null : new ResizeObserver(se);
-	ce?.observe(e), E();
-	function k(e) {
-		d && (Sn(a, e, f, d.span), u.target.set(...e.target), re(a.quaternion), a.zoom = p * e.scale, a.updateProjectionMatrix(), u.update(), ie());
+	u.addEventListener("start", () => {
+		b = !0, oe();
+	}), u.addEventListener("end", () => {
+		b = !1, oe();
+	}), u.addEventListener("change", oe);
+	function ce() {
+		let e = m > 0 ? a.zoom / m : 1;
+		k(), re(a.quaternion), a.zoom = m * e, a.updateProjectionMatrix(), O();
+	}
+	let le = typeof ResizeObserver > "u" ? null : new ResizeObserver(ce);
+	le?.observe(e), k();
+	let ue = hi();
+	i.add(ue.group), i.fog = new ee(Mn, 1, 2);
+	let N = Promise.all([r.compileAsync(i, a), l && r.compileAsync(l.scene, l.camera)]);
+	i.remove(ue.group), i.fog = null, ue.dispose(), f = N;
+	let P = () => {
+		f === N && (f = null, O());
+	};
+	N.then(P, P);
+	function F(e) {
+		d && (Sn(a, e, p, d.span), u.target.set(...e.target), re(a.quaternion), a.zoom = m * e.scale, a.updateProjectionMatrix(), u.update(), ie());
 	}
-	let le = {
+	let de = {
 		setLayers(e, t = null) {
-			le.clear(), e.length !== 0 && (d = pi(e), i.add(d.group), f = yn(d.span), i.fog = Ai(d, f, bn(d.cell)), s?.setElements(d.elements), l?.setCell(d.cell), E(), k(t ?? xn(d.cell)), T());
+			de.clear(), e.length !== 0 && (d = fi(e), i.add(d.group), p = yn(d.span), i.fog = Mi(d, p, bn(d.cell)), s?.setElements(d.elements), l?.setCell(d.cell), k(), F(t ?? xn(d.cell)), O());
 		},
 		clear() {
-			oe(), d &&= (i.remove(d.group), d.dispose(), null), i.fog = null, s?.setElements([]), r.setViewport(0, 0, m, h), r.clear();
+			se(), d &&= (i.remove(d.group), d.dispose(), null), i.fog = null, s?.setElements([]), r.setViewport(0, 0, h, g), r.clear();
 		},
 		getView() {
 			return d ? {
@@ -1967,48 +2029,49 @@ function ji(e, t = {}) {
 					u.target.y,
 					u.target.z
 				],
-				scale: p > 0 ? a.zoom / p : 1
+				scale: m > 0 ? a.zoom / m : 1
 			} : null;
 		},
 		setView(e) {
-			k(e), T();
+			F(e), O();
 		},
 		resetView() {
-			d && k(xn(d.cell)), T(), w();
+			d && F(xn(d.cell)), O(), D();
 		},
-		resize: se,
-		render: T,
+		resize: ce,
+		render: O,
 		onInteraction(e) {
-			return C.add(e), () => C.delete(e);
+			return E.add(e), () => E.delete(e);
 		},
 		dispose() {
-			le.clear(), y = !1, ce?.disconnect(), u.dispose(), l?.dispose(), r.dispose(), r.domElement.remove(), s?.element.remove(), e.classList.remove("bravais-viewer");
+			de.clear(), f = null, b = !1, le?.disconnect(), u.dispose(), l?.dispose(), r.dispose(), r.domElement.remove(), s?.element.remove(), e.classList.remove("bravais-viewer");
 		}
 	};
-	return le;
+	return de;
 }
 //#endregion
 //#region src/render/thumbnails.ts
-function Mi() {
-	let { renderer: e, scene: t, camera: n, updateHeadlight: r } = mi(), i = /* @__PURE__ */ new Map();
-	function a(e, t) {
+function Pi() {
+	let { renderer: e, scene: t, camera: n, updateHeadlight: r } = gi(), i = /* @__PURE__ */ new Map(), a = hi();
+	t.add(a.group), e.compile(t, n), t.remove(a.group), a.dispose();
+	function o(e, t) {
 		let n = i.get(e);
 		if (n && n.layers === t) return n.model;
 		n?.model.dispose();
-		let r = pi(t);
+		let r = fi(t);
 		return i.set(e, {
 			layers: t,
 			model: r
 		}), r;
 	}
 	return {
-		render(i, o, s, c) {
+		render(i, a, s, c) {
 			let l = c.getContext("2d");
-			if (!l || o.length === 0) return;
-			let u = a(i, o), d = Math.max(1, c.width), f = Math.max(1, c.height);
+			if (!l || a.length === 0) return;
+			let u = o(i, a), d = Math.max(1, c.width), f = Math.max(1, c.height);
 			(e.domElement.width !== d || e.domElement.height !== f) && e.setSize(d, f, !1), En(n, d, f);
 			let p = yn(u.span), m = s ?? xn(u.cell);
-			Sn(n, m, p, u.span), n.zoom = Tn(u.bounds, new O(...m.quaternion), d, f) * m.scale, n.updateProjectionMatrix(), r(), t.add(u.group), e.render(t, n), t.remove(u.group), l.clearRect(0, 0, d, f), l.drawImage(e.domElement, 0, 0, d, f);
+			Sn(n, m, p, u.span), n.zoom = Tn(u.bounds, new A(...m.quaternion), d, f) * m.scale, n.updateProjectionMatrix(), r(), t.add(u.group), e.render(t, n), t.remove(u.group), l.clearRect(0, 0, d, f), l.drawImage(e.domElement, 0, 0, d, f);
 		},
 		release(e) {
 			i.get(e)?.model.dispose(), i.delete(e);
@@ -2021,18 +2084,18 @@ function Mi() {
 }
 //#endregion
 //#region src/index.ts
-async function Ni(e) {
-	let t = ye(e), n = t.blocks.find(Xe) ?? t.blocks[0];
-	if (!n) throw new V("The CIF has no data block.");
-	let r = Ze(n);
+async function Fi(e) {
+	let t = ge(e), n = t.blocks.find(Ye) ?? t.blocks[0];
+	if (!n) throw new H("The CIF has no data block.");
+	let r = Xe(n);
 	return {
 		lattice: r.lattice,
 		sites: r.sites,
-		operations: await vt(r),
+		operations: await _t(r),
 		spaceGroup: r.spaceGroup
 	};
 }
 //#endregion
-export { V as CifStructureError, M as CifSyntaxError, rt as SpaceGroupLookupError, xt as SymmetryAnalysisError, It as buildScene, Mi as createThumbnailRenderer, ji as createViewer, qr as elementColor, Jr as elementColorAssignments, Kr as ensureElementColors, it as ensureMoyo, yt as expandStructure, Ne as latticeFromParameters, Pe as latticeFromVectors, ft as operationsForSpaceGroup, ye as parseCif, Ze as parseCifStructure, ze as parseSymmetryOperator, Ni as readCif, Yr as resetElementColors, Ft as structureElements, wt as symmetryDataset, Ot as typesetSymbol };
+export { H as CifStructureError, P as CifSyntaxError, nt as SpaceGroupLookupError, bt as SymmetryAnalysisError, Ft as buildScene, Pi as createThumbnailRenderer, Ni as createViewer, qr as elementColor, Jr as elementColorAssignments, Kr as ensureElementColors, rt as ensureMoyo, vt as expandStructure, Me as latticeFromParameters, Ne as latticeFromVectors, dt as operationsForSpaceGroup, ge as parseCif, Xe as parseCifStructure, Re as parseSymmetryOperator, Fi as readCif, Yr as resetElementColors, Pt as structureElements, Ct as symmetryDataset, Dt as typesetSymbol };
 
 //# sourceMappingURL=bravais.js.map