matterviz 0.4.1 → 0.4.2

package/dist/fermi-surface/FermiSurface.svelte

@@ -136,9 +136,7 @@
     spinner_props?: ComponentProps<typeof Spinner>
     loading?: boolean
     error_msg?: string
-    children?: Snippet<
-      [{ fermi_data?: FermiSurfaceData; bz_data?: BrillouinZoneData }]
-    >
+    children?: Snippet<[{ fermi_data?: FermiSurfaceData; bz_data?: BrillouinZoneData }]>
     on_file_drop?: (filename: string) => void
     on_file_load?: (data: FermiFileLoadData) => void
     on_error?: (data: FermiErrorData) => void
@@ -173,17 +171,12 @@
   // Yield to browser so spinner can render before heavy computation
   const tick = () =>
     new Promise<void>((resolve) =>
-      requestAnimationFrame(() => requestAnimationFrame(() => resolve()))
+      requestAnimationFrame(() => requestAnimationFrame(() => resolve())),
     )
 
   // Parse and load Fermi surface from content (async for UI responsiveness)
-  async function parse_fermi_content(
-    content: string | ArrayBuffer,
-    filename: string,
-  ) {
-    const text = content instanceof ArrayBuffer
-      ? new TextDecoder().decode(content)
-      : content
+  async function parse_fermi_content(content: string | ArrayBuffer, filename: string) {
+    const text = content instanceof ArrayBuffer ? new TextDecoder().decode(content) : content
 
     // parse_fermi_file throws a descriptive error when parsing fails
     const parsed = parse_fermi_file(text, filename)
@@ -208,9 +201,7 @@
     try {
       await parse_fermi_content(content, filename)
     } catch (err) {
-      error_msg = `Failed to parse ${filename}: ${
-        to_error(err).message
-      }`
+      error_msg = `Failed to parse ${filename}: ${to_error(err).message}`
       on_error?.({ error_msg, filename })
     }
   }
@@ -277,11 +268,10 @@
   // Compute BZ when structure or fermi_data changes
   $effect(() => {
     // Get k_lattice from available sources (priority order)
-    const k_lattice = fermi_data?.k_lattice ??
+    const k_lattice =
+      fermi_data?.k_lattice ??
       band_data?.k_lattice ??
-      (structure?.lattice?.matrix
-        ? reciprocal_lattice(structure.lattice.matrix)
-        : null)
+      (structure?.lattice?.matrix ? reciprocal_lattice(structure.lattice.matrix) : null)
 
     if (!k_lattice) {
       bz_data = undefined
@@ -361,8 +351,7 @@
     set_fullscreen: (val) => (fullscreen = val),
     bg_css_var: `--fermi-bg-fullscreen`,
     on_request_error: () => (fullscreen = false),
-    on_change: (val) =>
-      on_fullscreen_change?.({ fermi_data, bz_data, fullscreen: val }),
+    on_change: (val) => on_fullscreen_change?.({ fermi_data, bz_data, fullscreen: val }),
   })
 </script>
 
@@ -379,9 +368,12 @@
   onmouseenter={() => (hovered = true)}
   onmouseleave={() => (hovered = false)}
   ondrop={handle_file_drop}
-  {...drag_over_handlers({ allow: () => allow_file_drop, set_dragover: (over) => dragover = over })}
+  {...drag_over_handlers({
+    allow: () => allow_file_drop,
+    set_dragover: (over) => (dragover = over),
+  })}
   {...rest}
-  class="fermi-surface {rest.class ?? ``}"
+  class={[`fermi-surface`, rest.class]}
 >
   {@render children?.({ fermi_data, bz_data })}
   {#if loading}

package/dist/fermi-surface/FermiSurface.svelte.d.ts

@@ -52,12 +52,10 @@ type $$ComponentProps = {
     spinner_props?: ComponentProps<typeof Spinner>;
     loading?: boolean;
     error_msg?: string;
-    children?: Snippet<[
-        {
-            fermi_data?: FermiSurfaceData;
-            bz_data?: BrillouinZoneData;
-        }
-    ]>;
+    children?: Snippet<[{
+        fermi_data?: FermiSurfaceData;
+        bz_data?: BrillouinZoneData;
+    }]>;
     on_file_drop?: (filename: string) => void;
     on_file_load?: (data: FermiFileLoadData) => void;
     on_error?: (data: FermiErrorData) => void;
@@ -69,6 +67,6 @@ type $$ComponentProps = {
     on_point_hover?: (data: FermiHoverData | null) => void;
     on_hover?: (data: FermiHoverData | null) => void;
 } & HTMLAttributes<HTMLDivElement>;
-declare const FermiSurface: import("svelte").Component<$$ComponentProps, {}, "height" | "width" | "dragover" | "fullscreen" | "wrapper" | "hovered" | "controls_open" | "color_scale" | "loading" | "camera_projection" | "error_msg" | "surface_opacity" | "show_vectors" | "bz_data" | "mu" | "selected_bands" | "interpolation_factor" | "fermi_data" | "band_data" | "color_property" | "representation" | "show_bz" | "bz_opacity" | "tile_bz" | "clip_enabled" | "clip_axis" | "clip_position" | "clip_flip">;
+declare const FermiSurface: import("svelte").Component<$$ComponentProps, {}, "height" | "width" | "dragover" | "fullscreen" | "wrapper" | "hovered" | "controls_open" | "camera_projection" | "color_scale" | "loading" | "error_msg" | "surface_opacity" | "show_vectors" | "bz_data" | "mu" | "selected_bands" | "interpolation_factor" | "fermi_data" | "band_data" | "color_property" | "representation" | "show_bz" | "bz_opacity" | "tile_bz" | "clip_enabled" | "clip_axis" | "clip_position" | "clip_flip">;
 type FermiSurface = ReturnType<typeof FermiSurface>;
 export default FermiSurface;

package/dist/fermi-surface/FermiSurfaceControls.svelte

@@ -2,7 +2,7 @@
   import SettingsSection from '../layout/SettingsSection.svelte'
   import DraggablePane from '../overlays/DraggablePane.svelte'
   import type { CameraProjection } from '../settings'
-  import { make_change_detector } from '../utils'
+  import { make_change_detector, parse_num_token } from '../utils'
   import type { Snippet } from 'svelte'
   import type {
     BandGridData,
@@ -90,8 +90,8 @@
   let available_bands = $derived(
     fermi_data
       ? [...new Set(fermi_data.isosurfaces.map((iso) => iso.band_index))].toSorted(
-        (a, b) => a - b,
-      )
+          (a, b) => a - b,
+        )
       : [],
   )
   let available_bands_key = $derived(available_bands.join(`,`))
@@ -119,8 +119,7 @@
   }
 
   function handle_mu_change(event: Event & { currentTarget: HTMLInputElement }) {
-    const trimmed = event.currentTarget.value.trim()
-    const parsed = parseFloat(trimmed)
+    const parsed = parse_num_token(event.currentTarget.value)
     // Only update mu when input is valid; keep last valid value during transient
     // invalid states (e.g. empty string while user is typing a new value)
     if (Number.isFinite(parsed)) {
@@ -147,14 +146,7 @@
   >
     <label>
       <span>μ offset (eV):</span>
-      <input
-        type="range"
-        min="-1"
-        max="1"
-        step="0.01"
-        value={mu}
-        oninput={handle_mu_change}
-      />
+      <input type="range" min="-1" max="1" step="0.01" value={mu} oninput={handle_mu_change} />
       <input
         type="number"
         step="0.01"
@@ -300,13 +292,7 @@
       </label>
       <label>
         <span>Position:</span>
-        <input
-          type="range"
-          min="-1"
-          max="1"
-          step="0.01"
-          bind:value={clip_position}
-        />
+        <input type="range" min="-1" max="1" step="0.01" bind:value={clip_position} />
         <span class="value">{clip_position.toFixed(2)}</span>
       </label>
       <label>

package/dist/fermi-surface/FermiSurfaceControls.svelte.d.ts

@@ -30,6 +30,6 @@ type $$ComponentProps = {
         band_data?: BandGridData;
     }]>;
 };
-declare const FermiSurfaceControls: import("svelte").Component<$$ComponentProps, {}, "controls_open" | "color_scale" | "camera_projection" | "surface_opacity" | "show_vectors" | "mu" | "selected_bands" | "interpolation_factor" | "color_property" | "representation" | "show_bz" | "bz_opacity" | "tile_bz" | "clip_enabled" | "clip_axis" | "clip_position" | "clip_flip">;
+declare const FermiSurfaceControls: import("svelte").Component<$$ComponentProps, {}, "controls_open" | "camera_projection" | "color_scale" | "surface_opacity" | "show_vectors" | "mu" | "selected_bands" | "interpolation_factor" | "color_property" | "representation" | "show_bz" | "bz_opacity" | "tile_bz" | "clip_enabled" | "clip_axis" | "clip_position" | "clip_flip">;
 type FermiSurfaceControls = ReturnType<typeof FermiSurfaceControls>;
 export default FermiSurfaceControls;

package/dist/fermi-surface/FermiSurfaceScene.svelte

@@ -57,7 +57,7 @@
     bz_color = `#888888`,
     bz_opacity = 0.1,
     bz_edge_color = `#333333`,
-    bz_edge_width = 0.02,
+    bz_edge_width = 0.002,
     show_vectors = true,
     tile_bz = false,
     // Clipping plane
@@ -153,8 +153,8 @@
 
   // Filter surfaces based on selected bands
   let visible_surfaces = $derived(
-    fermi_data?.isosurfaces.filter((surface) =>
-      selected_bands === undefined || selected_bands.includes(surface.band_index)
+    fermi_data?.isosurfaces.filter(
+      (surface) => selected_bands === undefined || selected_bands.includes(surface.band_index),
     ) ?? [],
   )
 
@@ -208,9 +208,6 @@
 
   // Get color for a surface/vertex
   function get_surface_color(surface: Isosurface, vertex_idx?: number): string {
-    if (color_property === `band`) {
-      return constants.BAND_COLORS[surface.band_index % constants.BAND_COLORS.length]
-    }
     if (
       (color_property === `velocity` || color_property === `custom`) &&
       surface.properties &&
@@ -218,9 +215,7 @@
     ) {
       const prop = surface.properties[vertex_idx]
       const [min_val, max_val] = property_range
-      const normalized = max_val > min_val
-        ? (prop - min_val) / (max_val - min_val)
-        : 0.5
+      const normalized = max_val > min_val ? (prop - min_val) / (max_val - min_val) : 0.5
       return get_d3_interpolator(color_scale as D3InterpolateName)(normalized)
     }
     // Spin coloring
@@ -240,8 +235,7 @@
     const normals: number[] = []
     const colors: number[] = []
 
-    const use_vertex_colors = color_property === `velocity` ||
-      color_property === `custom`
+    const use_vertex_colors = color_property === `velocity` || color_property === `custom`
 
     const n_vertices = surface.vertices.length
 
@@ -257,9 +251,12 @@
 
         // Validate face indices are within bounds (protects against malformed JSON files)
         if (
-          idx0 < 0 || idx0 >= n_vertices ||
-          idx1 < 0 || idx1 >= n_vertices ||
-          idx2 < 0 || idx2 >= n_vertices
+          idx0 < 0 ||
+          idx0 >= n_vertices ||
+          idx1 < 0 ||
+          idx1 >= n_vertices ||
+          idx2 < 0 ||
+          idx2 >= n_vertices
         ) {
           continue
         }
@@ -297,10 +294,7 @@
     }
 
     const geometry = new BufferGeometry()
-    geometry.setAttribute(
-      `position`,
-      new BufferAttribute(new Float32Array(positions), 3),
-    )
+    geometry.setAttribute(`position`, new BufferAttribute(new Float32Array(positions), 3))
     geometry.setAttribute(`normal`, new BufferAttribute(new Float32Array(normals), 3))
 
     if (use_vertex_colors) {
@@ -358,18 +352,20 @@
     camera_position || default_camera_position(scene_size),
   )
 
-  const orbit_controls_props = $derived(build_orbit_props({
-    camera_projection,
-    target: rotation_target,
-    rotate_speed,
-    zoom_speed,
-    zoom_to_cursor,
-    pan_speed,
-    max_zoom,
-    min_zoom,
-    auto_rotate,
-    rotation_damping,
-  }))
+  const orbit_controls_props = $derived(
+    build_orbit_props({
+      camera_projection,
+      target: rotation_target,
+      rotate_speed,
+      zoom_speed,
+      zoom_to_cursor,
+      pan_speed,
+      max_zoom,
+      min_zoom,
+      auto_rotate,
+      rotation_damping,
+    }),
+  )
 
   // Create BZ geometry
   const bz_geometry = $derived(
@@ -423,11 +419,7 @@
     let [min_dist, nearest_idx] = [Infinity, 0]
     for (let idx = 0; idx < surface.vertices.length; idx++) {
       const vertex = surface.vertices[idx]
-      const dist = Math.hypot(
-        point[0] - vertex[0],
-        point[1] - vertex[1],
-        point[2] - vertex[2],
-      )
+      const dist = Math.hypot(point[0] - vertex[0], point[1] - vertex[1], point[2] - vertex[2])
       if (dist < min_dist) {
         min_dist = dist
         nearest_idx = idx
@@ -459,9 +451,8 @@
     const nearest_idx = find_nearest_vertex(surface, local_position)
     const property_value = surface.properties?.[nearest_idx]
     const has_velocities = fermi_data?.metadata?.has_velocities
-    const property_name = property_value != null
-      ? (has_velocities ? `velocity` : `custom`)
-      : undefined
+    const property_name =
+      property_value != null ? (has_velocities ? `velocity` : `custom`) : undefined
 
     const { clientX, clientY } = event.nativeEvent
     return {
@@ -533,19 +524,11 @@
 
   <!-- Reciprocal lattice vectors -->
   {#if show_vectors && fermi_data?.k_lattice}
-    <ReciprocalVectors
-      k_lattice={fermi_data.k_lattice}
-      {vector_scale}
-      size={scene_size}
-    />
+    <ReciprocalVectors k_lattice={fermi_data.k_lattice} {vector_scale} size={scene_size} />
   {/if}
 
   <!-- Fermi surfaces (with optional symmetry tiling) -->
-  {#each visible_surfaces as
-    surface,
-    surface_idx
-    (`surface-${surface.band_index}-${surface.spin}-${surface_idx}`)
-  }
+  {#each visible_surfaces as surface, surface_idx (`surface-${surface.band_index}-${surface.spin}-${surface_idx}`)}
     {@const geo_key = `${surface.band_index}-${surface.spin}-${surface_idx}`}
     {@const geo_data = geometry_cache.get(geo_key)}
     {@const surface_color = get_surface_color(surface)}
@@ -562,7 +545,7 @@
             matrixAutoUpdate={false}
             {renderOrder}
             onpointermove={(event: ThreltePointerEvent) =>
-            handle_pointer_move(event, surface, surface_color, sym_idx, sym_matrix)}
+              handle_pointer_move(event, surface, surface_color, sym_idx, sym_matrix)}
             onpointerleave={clear_hover}
           >
             <T.MeshBasicMaterial
@@ -583,7 +566,7 @@
             matrixAutoUpdate={false}
             renderOrder={renderOrder * 2}
             onpointermove={(event: ThreltePointerEvent) =>
-            handle_pointer_move(event, surface, surface_color, sym_idx, sym_matrix)}
+              handle_pointer_move(event, surface, surface_color, sym_idx, sym_matrix)}
             onpointerleave={clear_hover}
           >
             <T.MeshStandardMaterial
@@ -600,16 +583,11 @@
             matrixAutoUpdate={false}
             renderOrder={renderOrder * 2 + 1}
             onpointermove={(event: ThreltePointerEvent) =>
-            handle_pointer_move(event, surface, surface_color, sym_idx, sym_matrix)}
+              handle_pointer_move(event, surface, surface_color, sym_idx, sym_matrix)}
             onpointerleave={clear_hover}
           >
             <T.MeshStandardMaterial
-              {...get_material_props(
-                surface_color,
-                use_vertex_colors,
-                surface_idx,
-                `front`,
-              )}
+              {...get_material_props(surface_color, use_vertex_colors, surface_idx, `front`)}
               metalness={0.1}
               roughness={0.6}
               flatShading={false}
@@ -623,16 +601,11 @@
             matrixAutoUpdate={false}
             {renderOrder}
             onpointermove={(event: ThreltePointerEvent) =>
-            handle_pointer_move(event, surface, surface_color, sym_idx, sym_matrix)}
+              handle_pointer_move(event, surface, surface_color, sym_idx, sym_matrix)}
             onpointerleave={clear_hover}
           >
             <T.MeshStandardMaterial
-              {...get_material_props(
-                surface_color,
-                use_vertex_colors,
-                surface_idx,
-                `front`,
-              )}
+              {...get_material_props(surface_color, use_vertex_colors, surface_idx, `front`)}
               metalness={0.1}
               roughness={0.6}
               flatShading={false}

package/dist/fermi-surface/FermiSurfaceScene.svelte.d.ts

@@ -25,6 +25,6 @@ type $$ComponentProps = SceneControlProps & {
     vector_scale?: number;
     hover_data?: FermiHoverData | null;
 };
-declare const FermiSurfaceScene: import("svelte").Component<$$ComponentProps, {}, "scene" | "camera" | "camera_projection" | "camera_position" | "surface_opacity" | "bz_data" | "hover_data" | "fermi_data">;
+declare const FermiSurfaceScene: import("svelte").Component<$$ComponentProps, {}, "camera" | "camera_projection" | "camera_position" | "scene" | "surface_opacity" | "bz_data" | "hover_data" | "fermi_data">;
 type FermiSurfaceScene = ReturnType<typeof FermiSurfaceScene>;
 export default FermiSurfaceScene;

package/dist/fermi-surface/FermiSurfaceTooltip.svelte

@@ -5,7 +5,10 @@
   import { KCoords, TooltipContent } from '../tooltip'
   import type { FermiHoverData, FermiTooltipProp } from './types'
 
-  let { hover_data, tooltip }: {
+  let {
+    hover_data,
+    tooltip,
+  }: {
     hover_data: FermiHoverData
     tooltip?: FermiTooltipProp
   } = $props()
@@ -27,15 +30,15 @@
 
     {#if hover_data.property_value != null}
       <div class="property-row">
-        {hover_data.property_name || `Property`}: {
-          format_num(hover_data.property_value, `.4~`)
-        }
+        {hover_data.property_name || `Property`}: {format_num(
+          hover_data.property_value,
+          `.4~`,
+        )}
         <span class="nearest-note">(nearest)</span>
       </div>
     {/if}
 
-    {#if hover_data.is_tiled && hover_data.symmetry_index != null &&
-        hover_data.symmetry_index > 0}
+    {#if hover_data.is_tiled && hover_data.symmetry_index != null && hover_data.symmetry_index > 0}
       <div class="tiling-info">
         Symmetry copy #{hover_data.symmetry_index + 1}/48
       </div>

package/dist/fermi-surface/compute.js

@@ -3,6 +3,8 @@ import * as math from '../math';
 import { EPS } from '../math';
 import { CLOSED_CONTOUR_TOLERANCE, IRREDUCIBLE_BZ_MIN_VERTICES, IRREDUCIBLE_BZ_TOLERANCE, SPANNING_THRESHOLD, } from './constants';
 import { marching_cubes } from '../marching-cubes';
+// clamp01 guards float jitter at the cell boundary
+import { clamp01 } from '../utils';
 const safe_mod = (val, dim) => ((val % dim) + dim) % dim;
 // Precompute Catmull-Rom coefficients for a given t value
 // Returns [c0, c1, c2, c3] where result = c0*p0 + c1*p1 + c2*p2 + c3*p3
@@ -204,7 +206,7 @@ export function compute_surface_area(surface) {
         const cy = e1z * e2x - e1x * e2z;
         const cz = e1x * e2y - e1y * e2x;
         // Area is half the magnitude of cross product
-        total_area += Math.sqrt(cx * cx + cy * cy + cz * cz) * 0.5;
+        total_area += Math.hypot(cx, cy, cz) * 0.5;
     }
     return total_area;
 }
@@ -218,8 +220,6 @@ function compute_fermi_velocities(surface, velocity_grid, k_lattice, k_grid, per
     const k_inv = math.matrix_inverse_3x3(math.transpose_3x3_matrix(k_lattice));
     // Grid point i sits at frac i/(n−1) (endpoint-inclusive BXSF) or i/n (periodic)
     const [sx, sy, sz] = periodic ? [nx, ny, nz] : [nx - 1, ny - 1, nz - 1];
-    // Clamp guards float jitter at the cell boundary
-    const clamp01 = (val) => Math.min(Math.max(val, 0), 1);
     for (const vertex of surface.vertices) {
         const frac = math.mat3x3_vec3_multiply(k_inv, vertex); // centered, in [-0.5, 0.5]
         const gx = clamp01(frac[0] + 0.5) * sx;

package/dist/fermi-surface/parse.js

@@ -1,12 +1,13 @@
 import * as math from '../math';
-import { is_plain_object } from '../utils';
+import { is_plain_object, normalize_scientific_notation, parse_leading_num, to_error, } from '../utils';
 import * as constants from './constants';
 import { compute_vertex_normals } from '../marching-cubes';
 const parse_number_tokens = (line) => line.split(/\s+/).filter(Boolean);
-// Parse whitespace-separated floats from a line (optimized with unary +)
-const parse_floats = (line) => parse_number_tokens(line).map(Number);
+// Parse whitespace-separated floats from a line. Normalizes Fortran D-exponents
+// (`0.1234D+01`) which Fortran codes emit in BXSF/FRMSF and Number() rejects
+const parse_floats = (line) => parse_number_tokens(normalize_scientific_notation(line)).map(Number);
 // Parse whitespace-separated integers from a line
-const parse_ints = (line) => parse_number_tokens(line).map((part) => parseInt(part, 10));
+const parse_ints = (line) => parse_number_tokens(line).map((part) => Math.trunc(Number(part)));
 // Parse BXSF (Band-XSF) format used by XCrySDen, Quantum ESPRESSO, etc.
 // Format specification: http://www.xcrysden.org/doc/XSF.html
 function parse_bxsf(content) {
@@ -40,7 +41,7 @@ function parse_bxsf(content) {
         throw new Error(`Expected BANDGRID_3D header, got: ${bandgrid_line}`);
     }
     // Parse number of bands
-    const n_bands = parseInt(next_line(), 10);
+    const n_bands = Math.trunc(Number(next_line()));
     if (isNaN(n_bands) || n_bands <= 0) {
         throw new Error(`Invalid number of bands in BXSF file`);
     }
@@ -114,9 +115,9 @@ function parse_bxsf(content) {
         const lower = line.toLowerCase();
         if (lower.includes(`fermi`) && lower.includes(`energy`)) {
             // Match patterns like "Fermi Energy = 5.123" or "fermi_energy: -0.5"
-            const match = /(?:=|:)\s*([-+]?\d*\.?\d+(?:[eE][-+]?\d+)?)/i.exec(line);
+            const match = /(?:=|:)\s*(?<value>[-+]?\d*\.?\d+(?:[eE][-+]?\d+)?)/i.exec(line);
             if (match) {
-                fermi_energy = parseFloat(match[1]);
+                fermi_energy = Number(match[1]);
                 break;
             }
         }
@@ -149,12 +150,12 @@ function parse_frmsf(content) {
     // 0 = Monkhorst-Pack grid (shifted)
     // 1 = Gamma-centered grid
     // 2 = Gamma + half grid shift
-    const lshift = parseInt(lines[line_idx++], 10);
+    const lshift = Math.trunc(Number(lines[line_idx++]));
     if (isNaN(lshift) || lshift < 0 || lshift > 2) {
         throw new Error(`FRMSF: Invalid lshift value (expected 0, 1, or 2)`);
     }
     // Line 3: number of bands
-    const n_bands = parseInt(lines[line_idx++], 10);
+    const n_bands = Math.trunc(Number(lines[line_idx++]));
     if (isNaN(n_bands) || n_bands <= 0) {
         throw new Error(`FRMSF: Invalid number of bands`);
     }
@@ -183,17 +184,14 @@ function parse_frmsf(content) {
         energies[spin_idx] = [];
         for (let band_idx = 0; band_idx < n_bands; band_idx++) {
             const energy_values = [];
-            // Read energy values (first value per line only, ignore auxiliary columns)
+            // Read energy values (first value per line only, ignore auxiliary columns;
+            // normalize Fortran D-exponents like parse_floats does)
             while (energy_values.length < total_points && line_idx < lines.length) {
-                const line = lines[line_idx];
-                const first_token = line?.split(/\s+/)[0];
-                if (first_token && !isNaN(parseFloat(first_token))) {
-                    energy_values.push(parseFloat(first_token));
-                    line_idx++;
-                }
-                else {
+                const energy = parse_leading_num(normalize_scientific_notation(lines[line_idx] ?? ``));
+                if (isNaN(energy))
                     break;
-                }
+                energy_values.push(energy);
+                line_idx++;
             }
             if (energy_values.length < total_points) {
                 throw new Error(`FRMSF spin ${spin_idx} band ${band_idx}: expected ${total_points} values, got ${energy_values.length}`);
@@ -347,7 +345,7 @@ function parse_ifermi_surface(data) {
     const band_indices = new Set();
     let has_spin = false;
     for (const [band_key, iso_list] of Object.entries(isosurfaces_obj)) {
-        const band_index = parseInt(band_key, 10);
+        const band_index = Math.trunc(Number(band_key));
         // spin is determined by sign: positive = up, negative = down
         const spin = band_index < 0 ? `down` : `up`;
         const abs_band_idx = Math.abs(band_index);
@@ -373,13 +371,13 @@ function parse_ifermi_surface(data) {
             let dimensionality;
             if (ifermi_iso.dimensionality) {
                 const dim = ifermi_iso.dimensionality.toLowerCase();
-                if (dim.includes(`1d`) || dim === `1d`)
+                if (dim.includes(`1d`))
                     dimensionality = `1D`;
-                else if (dim.includes(`quasi`) || dim === `quasi-2d`)
+                else if (dim.includes(`quasi`))
                     dimensionality = `quasi-2D`;
-                else if (dim.includes(`2d`) || dim === `2d`)
+                else if (dim.includes(`2d`))
                     dimensionality = `2D`;
-                else if (dim.includes(`3d`) || dim === `3d`)
+                else if (dim.includes(`3d`))
                     dimensionality = `3D`;
             }
             // Compute area for this isosurface using fan triangulation for N-gons
@@ -438,7 +436,7 @@ export function parse_fermi_file(content, filename) {
             return parse();
         }
         catch (error) {
-            errors.push(`${format}: ${error instanceof Error ? error.message : String(error)}`);
+            errors.push(`${format}: ${to_error(error).message}`);
             console.error(`${format} parse error:`, error);
             return null;
         }