matterviz 0.3.1 → 0.3.3

package/dist/brillouin/BrillouinZoneScene.svelte

@@ -1,59 +1,150 @@
-<script lang="ts">import { AXIS_COLORS, NEG_AXIS_COLORS } from '../colors';
-import * as math from '../math';
-import { DEFAULTS } from '../settings';
-import Arrow from '../structure/Arrow.svelte';
-import Cylinder from '../structure/Cylinder.svelte';
-import { T, useThrelte } from '@threlte/core';
-import * as extras from '@threlte/extras';
-import { BufferAttribute, BufferGeometry, Vector3 } from 'three';
-let { bz_data = $bindable(), camera_position = $bindable(), camera_projection = $bindable(`perspective`), surface_color = $bindable(`#4488ff`), surface_opacity = $bindable(0.3), edge_color = $bindable(`#000000`), edge_width = $bindable(0.05), show_vectors = $bindable(true), vector_scale = $bindable(1.0), 
-// Irreducible BZ options
-show_ibz = false, ibz_data = null, ibz_color = `#ff8844`, ibz_opacity = 0.5, rotation_damping = DEFAULTS.structure.rotation_damping, max_zoom = DEFAULTS.structure.max_zoom, min_zoom = DEFAULTS.structure.min_zoom, rotate_speed = DEFAULTS.structure.rotate_speed, zoom_speed = DEFAULTS.structure.zoom_speed, pan_speed = DEFAULTS.structure.pan_speed, zoom_to_cursor = DEFAULTS.structure.zoom_to_cursor, fov = DEFAULTS.structure.fov, initial_zoom = DEFAULTS.structure.initial_zoom, ambient_light = DEFAULTS.structure.ambient_light, directional_light = DEFAULTS.structure.directional_light, gizmo = DEFAULTS.structure.show_gizmo, auto_rotate = DEFAULTS.structure.auto_rotate, camera_is_moving = $bindable(false), scene = $bindable(), camera = $bindable(), k_path_points = [], k_path_labels = [], hovered_k_point = null, hovered_qpoint_index = null, hover_data = $bindable(null), } = $props();
-const threlte = useThrelte();
-$effect(() => {
-    scene = threlte.scene;
-    camera = threlte.camera.current;
+<script lang="ts">
+  import { AXIS_COLORS, NEG_AXIS_COLORS } from '../colors'
+  import type { Vec3 } from '../math'
+  import * as math from '../math'
+  import { type CameraProjection, DEFAULTS } from '../settings'
+  import Arrow from '../structure/Arrow.svelte'
+  import Cylinder from '../structure/Cylinder.svelte'
+  import { T, useThrelte } from '@threlte/core'
+  import * as extras from '@threlte/extras'
+  import type { ComponentProps } from 'svelte'
+  import type { Camera, Scene } from 'three'
+  import { BufferAttribute, BufferGeometry, Vector3 } from 'three'
+  import type { BrillouinZoneData, BZHoverData, IrreducibleBZData } from './types'
+
+  // Threlte pointer event type for mesh interactions
+  type ThreltePointerEvent = { point: Vector3; nativeEvent: PointerEvent }
+
+  let {
+    bz_data = $bindable(),
+    camera_position = $bindable(),
+    camera_projection = $bindable(`perspective`),
+    surface_color = $bindable(`#4488ff`),
+    surface_opacity = $bindable(0.3),
+    edge_color = $bindable(`#000000`),
+    edge_width = $bindable(0.05),
+    show_vectors = $bindable(true),
+    vector_scale = $bindable(1.0),
+    // Irreducible BZ options
+    show_ibz = false,
+    ibz_data = null as IrreducibleBZData | null,
+    ibz_color = `#ff8844`,
+    ibz_opacity = 0.5,
+    rotation_damping = DEFAULTS.structure.rotation_damping,
+    max_zoom = DEFAULTS.structure.max_zoom,
+    min_zoom = DEFAULTS.structure.min_zoom,
+    rotate_speed = DEFAULTS.structure.rotate_speed,
+    zoom_speed = DEFAULTS.structure.zoom_speed,
+    pan_speed = DEFAULTS.structure.pan_speed,
+    zoom_to_cursor = DEFAULTS.structure.zoom_to_cursor,
+    fov = DEFAULTS.structure.fov,
+    initial_zoom = DEFAULTS.structure.initial_zoom,
+    ambient_light = DEFAULTS.structure.ambient_light,
+    directional_light = DEFAULTS.structure.directional_light,
+    gizmo = DEFAULTS.structure.show_gizmo,
+    auto_rotate = DEFAULTS.structure.auto_rotate,
+    camera_is_moving = $bindable(false),
+    scene = $bindable(),
+    camera = $bindable(),
+    k_path_points = [],
+    k_path_labels = [],
+    hovered_k_point = null,
+    hovered_qpoint_index = null,
+    hover_data = $bindable<BZHoverData | null>(null),
+  }: {
+    bz_data?: BrillouinZoneData
+    camera_position?: Vec3 | undefined
+    camera_projection?: CameraProjection
+    surface_color?: string
+    surface_opacity?: number
+    edge_color?: string
+    edge_width?: number
+    show_vectors?: boolean
+    vector_scale?: number
+    // Irreducible BZ options
+    show_ibz?: boolean
+    ibz_data?: IrreducibleBZData | null
+    ibz_color?: string
+    ibz_opacity?: number
+    rotation_damping?: number
+    max_zoom?: number
+    min_zoom?: number
+    rotate_speed?: number
+    zoom_speed?: number
+    pan_speed?: number
+    zoom_to_cursor?: boolean
+    fov?: number
+    initial_zoom?: number
+    ambient_light?: number
+    directional_light?: number
+    gizmo?: boolean | ComponentProps<typeof extras.Gizmo>
+    auto_rotate?: number
+    camera_is_moving?: boolean
+    scene?: Scene
+    camera?: Camera
+    k_path_points?: Vec3[]
+    k_path_labels?: { position: Vec3; label: string | null }[]
+    hovered_k_point?: Vec3 | null
+    hovered_qpoint_index?: number | null
+    hover_data?: BZHoverData | null
+  } = $props()
+
+  const threlte = useThrelte()
+  $effect(() => {
+    scene = threlte.scene
+    camera = threlte.camera.current
     if (threlte.renderer) {
-        Object.assign(threlte.renderer.domElement, { __renderer: threlte.renderer });
+      Object.assign(threlte.renderer.domElement, { __renderer: threlte.renderer })
     }
-});
-extras.interactivity();
-// Compute centroid of BZ vertices for proper rotation center
-const rotation_target = $derived.by(() => {
-    if (!bz_data?.vertices || bz_data.vertices.length === 0)
-        return [0, 0, 0];
-    const sum = bz_data.vertices.reduce((acc, v) => math.add(acc, v), [0, 0, 0]);
-    return math.scale(sum, 1 / bz_data.vertices.length);
-});
-// BZ size for camera positioning: average magnitude of k-vectors
-const bz_size = $derived.by(() => {
-    if (!bz_data?.k_lattice)
-        return 10;
-    const mags = bz_data.k_lattice.map((vec) => Math.hypot(...vec));
-    return mags.reduce((sum, mag) => sum + mag, 0) / 3;
-});
-const computed_camera_position = $derived.by(() => camera_position || [10, 3, 8].map((x) => x * Math.max(1, bz_size)));
-const gizmo_props = $derived({
+  })
+
+  extras.interactivity()
+
+  // Compute centroid of BZ vertices for proper rotation center
+  const rotation_target = $derived.by((): Vec3 => {
+    if (!bz_data?.vertices || bz_data.vertices.length === 0) return [0, 0, 0]
+    const sum = bz_data.vertices.reduce(
+      (acc, v) => math.add(acc, v),
+      [0, 0, 0] as Vec3,
+    )
+    return math.scale(sum, 1 / bz_data.vertices.length)
+  })
+
+  // BZ size for camera positioning: average magnitude of k-vectors
+  const bz_size = $derived.by(() => {
+    if (!bz_data?.k_lattice) return 10
+    const mags = bz_data.k_lattice.map((vec) => Math.hypot(...vec))
+    return mags.reduce((sum, mag) => sum + mag, 0) / 3
+  })
+
+  const computed_camera_position = $derived.by(() =>
+    camera_position || ([10, 3, 8].map((x) => x * Math.max(1, bz_size)) as Vec3)
+  )
+
+  const gizmo_props = $derived({
     background: { enabled: false },
     className: `responsive-gizmo`,
-    ...Object.fromEntries([...AXIS_COLORS, ...NEG_AXIS_COLORS].map(([axis, color, hover]) => [
+    ...Object.fromEntries(
+      [...AXIS_COLORS, ...NEG_AXIS_COLORS].map(([axis, color, hover]) => [
         axis,
         {
-            color,
-            labelColor: `#111`,
-            opacity: axis.startsWith(`n`) ? 0.9 : 0.8,
-            hover: {
-                color: hover,
-                labelColor: `#222`,
-                opacity: axis.startsWith(`n`) ? 1 : 0.9,
-            },
+          color,
+          labelColor: `#111`,
+          opacity: axis.startsWith(`n`) ? 0.9 : 0.8,
+          hover: {
+            color: hover,
+            labelColor: `#222`,
+            opacity: axis.startsWith(`n`) ? 1 : 0.9,
+          },
         },
-    ])),
+      ]),
+    ),
     ...(typeof gizmo === `object` ? gizmo : {}),
     offset: { left: 5, bottom: 5 },
-});
-const is_ortho = $derived(camera_projection === `orthographic`);
-const orbit_controls_props = $derived({
+  })
+
+  const is_ortho = $derived(camera_projection === `orthographic`)
+  const orbit_controls_props = $derived({
     position: [0, 0, 0],
     target: rotation_target,
     enableRotate: rotate_speed > 0,
@@ -71,134 +162,155 @@ const orbit_controls_props = $derived({
     dampingFactor: rotation_damping,
     onstart: () => (camera_is_moving = true),
     onend: () => (camera_is_moving = false),
-});
-const vector_colors = [`red`, `green`, `blue`];
-const vector_labels = [`b₁`, `b₂`, `b₃`];
-// Create mesh geometry from faces with fan triangulation
-function create_mesh_geometry(vertices, faces) {
-    if (faces.length === 0)
-        return null;
-    const positions = [];
-    const normals = [];
+  })
+
+  const vector_colors = [`red`, `green`, `blue`]
+  const vector_labels = [`b₁`, `b₂`, `b₃`]
+
+  // Create mesh geometry from faces with fan triangulation
+  function create_mesh_geometry(
+    vertices: Vec3[],
+    faces: number[][],
+  ): BufferGeometry | null {
+    if (faces.length === 0) return null
+
+    const positions: number[] = []
+    const normals: number[] = []
+
     for (const face of faces) {
-        if (face.length < 3)
-            continue;
-        for (let face_idx = 1; face_idx < face.length - 1; face_idx++) {
-            const indices = [face[0], face[face_idx], face[face_idx + 1]];
-            if (indices.some((idx) => idx < 0 || idx >= vertices.length))
-                continue;
-            const [v0, v1, v2] = indices.map((idx) => vertices[idx]);
-            positions.push(...v0, ...v1, ...v2);
-            const e1 = math.subtract(v1, v0);
-            const e2 = math.subtract(v2, v0);
-            const normal_vec = math.cross_3d(e1, e2);
-            const len = Math.hypot(...normal_vec);
-            const norm = len > 1e-10 ? normal_vec.map((x) => x / len) : [0, 0, 0];
-            normals.push(...norm, ...norm, ...norm);
-        }
+      if (face.length < 3) continue
+      for (let face_idx = 1; face_idx < face.length - 1; face_idx++) {
+        const indices = [face[0], face[face_idx], face[face_idx + 1]]
+        if (indices.some((idx) => idx < 0 || idx >= vertices.length)) continue
+        const [v0, v1, v2] = indices.map((idx) => vertices[idx])
+        positions.push(...v0, ...v1, ...v2)
+
+        const e1: Vec3 = math.subtract(v1, v0)
+        const e2: Vec3 = math.subtract(v2, v0)
+        const normal_vec = math.cross_3d(e1, e2)
+        const len = Math.hypot(...normal_vec)
+        const norm = len > 1e-10 ? normal_vec.map((x) => x / len) : [0, 0, 0]
+        normals.push(...norm, ...norm, ...norm)
+      }
     }
-    const geometry = new BufferGeometry();
-    geometry.setAttribute(`position`, new BufferAttribute(new Float32Array(positions), 3));
-    geometry.setAttribute(`normal`, new BufferAttribute(new Float32Array(normals), 3));
-    geometry.computeBoundingSphere();
-    return geometry;
-}
-const bz_geometry = $derived(bz_data ? create_mesh_geometry(bz_data.vertices, bz_data.faces) : null);
-const ibz_geometry = $derived(show_ibz && ibz_data
-    ? create_mesh_geometry(ibz_data.vertices, ibz_data.faces)
-    : null);
-// Separate effects to avoid disposing one geometry when only the other changes
-$effect(() => {
-    const prev = bz_geometry;
-    return () => prev?.dispose();
-});
-$effect(() => {
-    const prev = ibz_geometry;
-    return () => prev?.dispose();
-});
-// Compute inverse of k_lattice for Cartesian->fractional conversion
-const k_lattice_inv = $derived.by(() => {
-    if (!bz_data?.k_lattice)
-        return null;
+
+    const geometry = new BufferGeometry()
+    geometry.setAttribute(
+      `position`,
+      new BufferAttribute(new Float32Array(positions), 3),
+    )
+    geometry.setAttribute(`normal`, new BufferAttribute(new Float32Array(normals), 3))
+    geometry.computeBoundingSphere()
+    return geometry
+  }
+
+  const bz_geometry = $derived(
+    bz_data ? create_mesh_geometry(bz_data.vertices, bz_data.faces) : null,
+  )
+  const ibz_geometry = $derived(
+    show_ibz && ibz_data
+      ? create_mesh_geometry(ibz_data.vertices, ibz_data.faces)
+      : null,
+  )
+
+  // Separate effects to avoid disposing one geometry when only the other changes
+  $effect(() => {
+    const prev = bz_geometry
+    return () => prev?.dispose()
+  })
+  $effect(() => {
+    const prev = ibz_geometry
+    return () => prev?.dispose()
+  })
+
+  // Compute inverse of k_lattice for Cartesian->fractional conversion
+  const k_lattice_inv = $derived.by(() => {
+    if (!bz_data?.k_lattice) return null
     try {
-        return math.matrix_inverse_3x3(bz_data.k_lattice);
-    }
-    catch {
-        return null;
+      return math.matrix_inverse_3x3(bz_data.k_lattice)
+    } catch {
+      return null
     }
-});
-// Convert Cartesian k-coordinates to fractional (reciprocal lattice units)
-function cartesian_to_fractional(cart) {
-    if (!k_lattice_inv)
-        return null;
-    return math.mat3x3_vec3_multiply(k_lattice_inv, cart);
-}
-// Throttle state for pointer move events
-let last_hover_time = 0;
-let last_hover_mesh = null;
-const HOVER_THROTTLE_MS = 16; // ~60fps
-// Reset throttle when bz_data changes to ensure immediate response
-$effect(() => {
+  })
+
+  // Convert Cartesian k-coordinates to fractional (reciprocal lattice units)
+  function cartesian_to_fractional(cart: Vec3): Vec3 | null {
+    if (!k_lattice_inv) return null
+    return math.mat3x3_vec3_multiply(k_lattice_inv, cart)
+  }
+
+  // Throttle state for pointer move events
+  let last_hover_time = 0
+  let last_hover_mesh: `bz` | `ibz` | null = null
+  const HOVER_THROTTLE_MS = 16 // ~60fps
+
+  // Reset throttle when bz_data changes to ensure immediate response
+  $effect(() => {
     if (bz_data) {
-        last_hover_time = 0;
-        last_hover_mesh = null;
+      last_hover_time = 0
+      last_hover_mesh = null
     }
-});
-// Track IBZ hover state - IBZ takes priority over BZ
-let ibz_hovered = false;
-$effect(() => {
+  })
+
+  // Track IBZ hover state - IBZ takes priority over BZ
+  let ibz_hovered = false
+  $effect(() => {
     if (!show_ibz) {
-        ibz_hovered = false;
-        // Clear hover tooltip if it was showing IBZ data
-        if (hover_data?.is_ibz)
-            hover_data = null;
+      ibz_hovered = false
+      // Clear hover tooltip if it was showing IBZ data
+      if (hover_data?.is_ibz) hover_data = null
     }
-});
-// Create hover data from pointer event
-function create_hover_data(event, is_ibz) {
-    if (!bz_data)
-        return null;
-    const position_cartesian = [event.point.x, event.point.y, event.point.z];
-    const position_fractional = cartesian_to_fractional(position_cartesian);
-    const { clientX, clientY } = event.nativeEvent;
-    const ibz_vol = ibz_data?.volume ?? null;
+  })
+
+  // Create hover data from pointer event
+  function create_hover_data(
+    event: ThreltePointerEvent,
+    is_ibz: boolean,
+  ): BZHoverData | null {
+    if (!bz_data) return null
+
+    const position_cartesian: Vec3 = [event.point.x, event.point.y, event.point.z]
+    const position_fractional = cartesian_to_fractional(position_cartesian)
+
+    const { clientX, clientY } = event.nativeEvent
+    const ibz_vol = ibz_data?.volume ?? null
     // Round to nearest integer since symmetry multiplicity is the point group order
     const symmetry_multiplicity = ibz_vol != null && ibz_vol > 0
-        ? Math.round(bz_data.volume / ibz_vol)
-        : null;
+      ? Math.round(bz_data.volume / ibz_vol)
+      : null
+
     return {
-        position_cartesian,
-        position_fractional,
-        screen_position: { x: clientX, y: clientY },
-        is_ibz,
-        bz_order: bz_data.order,
-        bz_volume: bz_data.volume,
-        ibz_volume: ibz_vol,
-        symmetry_multiplicity,
-    };
-}
-// Throttled hover handler - IBZ takes priority over BZ
-function handle_hover(event, is_ibz) {
-    if (is_ibz)
-        ibz_hovered = true;
-    else if (ibz_hovered)
-        return; // BZ defers to IBZ
-    const mesh = is_ibz ? `ibz` : `bz`;
-    const now = performance.now();
+      position_cartesian,
+      position_fractional,
+      screen_position: { x: clientX, y: clientY },
+      is_ibz,
+      bz_order: bz_data.order,
+      bz_volume: bz_data.volume,
+      ibz_volume: ibz_vol,
+      symmetry_multiplicity,
+    }
+  }
+
+  // Throttled hover handler - IBZ takes priority over BZ
+  function handle_hover(event: ThreltePointerEvent, is_ibz: boolean): void {
+    if (is_ibz) ibz_hovered = true
+    else if (ibz_hovered) return // BZ defers to IBZ
+
+    const mesh = is_ibz ? `ibz` : `bz`
+    const now = performance.now()
     // Bypass throttle when switching meshes for responsive transitions
-    if (last_hover_mesh === mesh && now - last_hover_time < HOVER_THROTTLE_MS)
-        return;
-    last_hover_time = now;
-    last_hover_mesh = mesh;
-    hover_data = create_hover_data(event, is_ibz);
-}
-// Leave handler - IBZ clears state, BZ only clears if IBZ not hovered
-function handle_leave(is_ibz) {
-    if (is_ibz)
-        ibz_hovered = false;
-    if (is_ibz || !ibz_hovered)
-        hover_data = null;
-}
+    if (last_hover_mesh === mesh && now - last_hover_time < HOVER_THROTTLE_MS) return
+
+    last_hover_time = now
+    last_hover_mesh = mesh
+    hover_data = create_hover_data(event, is_ibz)
+  }
+
+  // Leave handler - IBZ clears state, BZ only clears if IBZ not hovered
+  function handle_leave(is_ibz: boolean): void {
+    if (is_ibz) ibz_hovered = false
+    if (is_ibz || !ibz_hovered) hover_data = null
+  }
 </script>
 
 {#if camera_projection === `perspective`}
@@ -243,7 +355,7 @@ function handle_leave(is_ibz) {
     {/if}
 
     <!-- BZ edges -->
-    {#each bz_data.edges as edge_segment (JSON.stringify(edge_segment))}
+    {#each bz_data.edges as edge_segment, edge_idx (`bz-edge-${edge_idx}`)}
       {@const [from, to] = edge_segment}
       <Cylinder {from} {to} thickness={edge_width} color={edge_color} />
     {/each}
@@ -267,7 +379,7 @@ function handle_leave(is_ibz) {
 
     <!-- IBZ edges -->
     {#if show_ibz && ibz_data}
-      {#each ibz_data.edges as edge_segment (JSON.stringify(edge_segment))}
+      {#each ibz_data.edges as edge_segment, edge_idx (`ibz-edge-${edge_idx}`)}
         {@const [from, to] = edge_segment}
         <Cylinder {from} {to} thickness={edge_width * 1.5} color={ibz_color} />
       {/each}

package/dist/brillouin/BrillouinZoneScene.svelte.d.ts

@@ -43,6 +43,6 @@ type $$ComponentProps = {
     hovered_qpoint_index?: number | null;
     hover_data?: BZHoverData | null;
 };
-declare const BrillouinZoneScene: import("svelte").Component<$$ComponentProps, {}, "camera_position" | "camera_projection" | "camera_is_moving" | "scene" | "camera" | "surface_color" | "surface_opacity" | "edge_color" | "edge_width" | "show_vectors" | "bz_data" | "vector_scale" | "hover_data">;
+declare const BrillouinZoneScene: import("svelte").Component<$$ComponentProps, {}, "camera_position" | "camera_projection" | "vector_scale" | "camera" | "scene" | "camera_is_moving" | "surface_color" | "surface_opacity" | "edge_color" | "edge_width" | "show_vectors" | "bz_data" | "hover_data">;
 type BrillouinZoneScene = ReturnType<typeof BrillouinZoneScene>;
 export default BrillouinZoneScene;

package/dist/brillouin/BrillouinZoneTooltip.svelte

@@ -1,10 +1,20 @@
-<script lang="ts">// Tooltip component for Brillouin zone hover information
-// Displays k-coordinates, BZ order, volume, and IBZ-specific info
-import { format_num, format_vec3 } from '../labels';
-import { TooltipContent } from '../tooltip';
-let { hover_data, tooltip, } = $props();
-// Ordinal for BZ order (only 1-3 in practice)
-const ordinal = (num) => `${num}${[`th`, `st`, `nd`, `rd`][num] ?? `th`}`;
+<script lang="ts">
+  // Tooltip component for Brillouin zone hover information
+  // Displays k-coordinates, BZ order, volume, and IBZ-specific info
+  import { format_num, format_vec3 } from '../labels'
+  import { TooltipContent } from '../tooltip'
+  import type { BZHoverData, BZTooltipProp } from './types'
+
+  let {
+    hover_data,
+    tooltip,
+  }: {
+    hover_data: BZHoverData
+    tooltip?: BZTooltipProp
+  } = $props()
+
+  // Ordinal for BZ order (only 1-3 in practice)
+  const ordinal = (num: number) => `${num}${[`th`, `st`, `nd`, `rd`][num] ?? `th`}`
 </script>
 
 <TooltipContent data={hover_data} snippet_arg={{ hover_data }} {tooltip}>

package/dist/brillouin/compute.js

@@ -29,11 +29,14 @@ export function extract_point_group_from_operations(operations) {
 }
 // Multiply two 3x3 matrices: C = A · B
 function mat3x3_multiply(A, B) {
-    const result = [[0, 0, 0], [0, 0, 0], [0, 0, 0]];
+    const result = [
+        [0, 0, 0],
+        [0, 0, 0],
+        [0, 0, 0],
+    ];
     for (let row = 0; row < 3; row++) {
         for (let col = 0; col < 3; col++) {
-            result[row][col] = A[row][0] * B[0][col] + A[row][1] * B[1][col] +
-                A[row][2] * B[2][col];
+            result[row][col] = A[row][0] * B[0][col] + A[row][1] * B[1][col] + A[row][2] * B[2][col];
         }
     }
     return result;
@@ -164,7 +167,8 @@ max_planes_by_order = { 1: 26, 2: 80, 3: 150 }) {
                 // Count how many planes this vertex is beyond (with early termination)
                 let beyond_count = 0;
                 for (let p_idx = 0; p_idx < normals.length; p_idx++) {
-                    const dot = vertex[0] * normals[p_idx][0] + vertex[1] * normals[p_idx][1] +
+                    const dot = vertex[0] * normals[p_idx][0] +
+                        vertex[1] * normals[p_idx][1] +
                         vertex[2] * normals[p_idx][2];
                     if (dot > distances[p_idx] + TOL) {
                         beyond_count++;
@@ -207,9 +211,10 @@ export function compute_convex_hull(vertices, edge_sharp_angle_deg = 5) {
     const vert_map = new Map();
     for (let idx_vertex = 0; idx_vertex < pos.count; idx_vertex++) {
         const vert = [pos.getX(idx_vertex), pos.getY(idx_vertex), pos.getZ(idx_vertex)];
-        const existing_idx = unique_verts.findIndex((u) => Math.abs(u[0] - vert[0]) < TOL && Math.abs(u[1] - vert[1]) < TOL &&
+        const existing_idx = unique_verts.findIndex((u) => Math.abs(u[0] - vert[0]) < TOL &&
+            Math.abs(u[1] - vert[1]) < TOL &&
             Math.abs(u[2] - vert[2]) < TOL);
-        vert_map.set(idx_vertex, existing_idx === -1 ? (unique_verts.push(vert) - 1) : existing_idx);
+        vert_map.set(idx_vertex, existing_idx === -1 ? unique_verts.push(vert) - 1 : existing_idx);
     }
     // Build faces with deduplicated vertex indices
     const faces = [];