matterviz 0.3.7 → 0.4.1

package/dist/symmetry/SymmetryElementControls.svelte

@@ -0,0 +1,69 @@
+<!-- Legend + per-kind visibility toggles for symmetry-element overlays. Renders one
+checkbox per element kind present in `elements`, with a color swatch matching the
+overlay render colors and the element count. Bind `show_kinds` and pass it through to
+SymmetryElements (via symmetry_elements_props on the structure viewer). -->
+<script lang="ts">
+  import type { ShowSymmetryKinds, SymmetryElement } from './symmetry-elements'
+  import {
+    count_symmetry_elements,
+    DEFAULT_SHOW_SYM_KINDS,
+    SYM_ELEM_KIND_INFO,
+    SYM_ELEM_KINDS,
+  } from './symmetry-elements'
+
+  let {
+    elements = [],
+    show_kinds = $bindable({ ...DEFAULT_SHOW_SYM_KINDS }),
+    ...rest
+  }: {
+    elements?: SymmetryElement[]
+    show_kinds?: ShowSymmetryKinds
+    [key: string]: unknown
+  } = $props()
+
+  const counts = $derived(count_symmetry_elements(elements))
+  const present_kinds = $derived(SYM_ELEM_KINDS.filter((kind) => counts[kind]))
+
+  function toggle_kind(kind: (typeof SYM_ELEM_KINDS)[number], checked: boolean) {
+    // Reassign (not mutate) so bound parents always see the change
+    show_kinds = { ...show_kinds, [kind]: checked }
+  }
+</script>
+
+{#if present_kinds.length > 0}
+  <div class="sym-elem-controls" {...rest}>
+    {#each present_kinds as kind (kind)}
+      <label>
+        <input
+          type="checkbox"
+          checked={show_kinds[kind] ?? false}
+          onchange={(evt) => toggle_kind(kind, evt.currentTarget.checked)}
+        />
+        <span class="swatch" style:background={SYM_ELEM_KIND_INFO[kind].color}></span>
+        {SYM_ELEM_KIND_INFO[kind].label} ({counts[kind]})
+      </label>
+    {/each}
+  </div>
+{/if}
+
+<style>
+  .sym-elem-controls {
+    display: flex;
+    flex-direction: column;
+    gap: 3pt;
+  }
+  label {
+    display: flex;
+    align-items: center;
+    gap: 6pt;
+    cursor: pointer;
+    font-size: 0.95em;
+  }
+  .swatch {
+    width: 0.9em;
+    height: 0.9em;
+    border-radius: 2px;
+    border: 1px solid var(--border-color, #ccc);
+    flex-shrink: 0;
+  }
+</style>

package/dist/symmetry/SymmetryElementControls.svelte.d.ts

@@ -0,0 +1,9 @@
+import type { ShowSymmetryKinds, SymmetryElement } from './symmetry-elements';
+type $$ComponentProps = {
+    elements?: SymmetryElement[];
+    show_kinds?: ShowSymmetryKinds;
+    [key: string]: unknown;
+};
+declare const SymmetryElementControls: import("svelte").Component<$$ComponentProps, {}, "show_kinds">;
+type SymmetryElementControls = ReturnType<typeof SymmetryElementControls>;
+export default SymmetryElementControls;

package/dist/symmetry/SymmetryElements.svelte

@@ -0,0 +1,354 @@
+<!-- Render crystallographic symmetry elements (rotation/screw axes, mirror/glide
+planes, inversion centers, rotoinversion axes) inside a Threlte scene. Elements come
+from symmetry_elements_from_ops and are expressed in fractional coordinates of the cell
+described by `lattice` — make sure both refer to the SAME cell (moyo operations are in
+the input-cell frame, so pass the original structure's lattice).
+
+Visual conventions (loosely following ITA diagram conventions: translation-carrying
+elements are dashed/striped):
+- rotation axes: solid cylinders, colored by order
+- screw axes: DASHED cylinders (short dashes, gaps narrower than dashes; same order
+  colors, slightly thinner)
+- mirror planes: solid translucent fills with opaque outlines
+- glide planes: STRIPED translucent fills (stripes run along the glide-translation
+  direction) with opaque outlines
+- inversion centers / rotoinversion markers: small faceted octahedra — themselves
+  centrosymmetric, and clearly distinct from the smooth spheres used for atoms
+
+For performance, geometries are merged per material group (one draw call per distinct
+color/opacity instead of one mesh per element) and disposed on change/unmount. -->
+<script lang="ts">
+  import type { Matrix3x3, Vec2, Vec3 } from '../math'
+  import * as math from '../math'
+  import type { ShowSymmetryKinds, SymmetryElement } from './symmetry-elements'
+  import {
+    clip_line_to_cell,
+    clip_plane_to_cell,
+    dash_segments,
+    DEFAULT_SHOW_SYM_KINDS,
+    frac_to_cart_direction,
+  } from './symmetry-elements'
+  import { T } from '@threlte/core'
+  import {
+    BufferAttribute,
+    BufferGeometry,
+    ClampToEdgeWrapping,
+    CylinderGeometry,
+    DataTexture,
+    DoubleSide,
+    LinearFilter,
+    Matrix4,
+    OctahedronGeometry,
+    Quaternion,
+    RepeatWrapping,
+    RGBAFormat,
+    Vector3,
+  } from 'three'
+  import { mergeGeometries } from 'three/examples/jsm/utils/BufferGeometryUtils.js'
+
+  let {
+    elements = [],
+    lattice,
+    // Per-kind visibility. Defaults to rotation axes ONLY: drawing every kind at once
+    // buries the structure under overlays for high-symmetry cells. Toggle additional
+    // kinds individually (e.g. via SymmetryElementControls).
+    show_kinds = DEFAULT_SHOW_SYM_KINDS,
+    // hide 2-fold axes that are sub-elements of higher-order axes on the same line
+    hide_redundant_axes = true,
+    axis_radius = 0.04,
+    screw_radius = 0.03,
+    // [dash, gap] in Å for dashed screw axes — gap narrower than dash so the line
+    // reads as continuous-but-broken rather than sparse
+    screw_dash = [0.25, 0.1] as Vec2,
+    inversion_radius = 0.12,
+    plane_opacity = 0.2,
+    glide_opacity = 0.15,
+    // opaque polygon outlines make overlapping translucent planes legible
+    show_plane_edges = true,
+    plane_edge_opacity = 0.9,
+    // stripe period in Å for glide-plane fills (stripes run along the glide direction)
+    glide_stripe_period = 0.7,
+    axis_colors = { 2: `#e63946`, 3: `#2a9d8f`, 4: `#3a6fb0`, 6: `#9c27b0` },
+    mirror_color = `#ffb703`,
+    glide_color = `#8ecae6`,
+    inversion_color = `#555555`,
+  }: {
+    elements?: SymmetryElement[]
+    lattice: Matrix3x3
+    show_kinds?: ShowSymmetryKinds
+    hide_redundant_axes?: boolean
+    axis_radius?: number
+    screw_radius?: number
+    screw_dash?: Vec2
+    inversion_radius?: number
+    plane_opacity?: number
+    glide_opacity?: number
+    show_plane_edges?: boolean
+    plane_edge_opacity?: number
+    glide_stripe_period?: number
+    axis_colors?: Record<number, string>
+    mirror_color?: string
+    glide_color?: string
+    inversion_color?: string
+  } = $props()
+
+  const UP = new Vector3(0, 1, 0)
+  const UNIT_SCALE = new Vector3(1, 1, 1)
+
+  // 1D stripe alpha texture for glide fills (three.js alphaMap samples the GREEN
+  // channel): ~55% full-alpha stripe, ~45% faint background so the plane stays
+  // contiguous between stripes. Repeats along U; V is constant.
+  const stripe_texture = (() => {
+    const width = 16
+    const data = new Uint8Array(width * 4)
+    for (let px = 0; px < width; px++) {
+      const val = px < 9 ? 255 : 56
+      data.set([val, val, val, 255], px * 4)
+    }
+    const tex = new DataTexture(data, width, 1, RGBAFormat)
+    tex.wrapS = RepeatWrapping
+    tex.wrapT = ClampToEdgeWrapping
+    tex.magFilter = LinearFilter
+    tex.minFilter = LinearFilter
+    tex.needsUpdate = true
+    return tex
+  })()
+
+  // Key identifying the geometric line of an axis element (direction + intercept)
+  const line_key = (elem: SymmetryElement): string => {
+    const axis = elem.axis
+    if (!axis) return ``
+    const lambda = math.dot(elem.point, axis) / math.dot(axis, axis)
+    const intercept = elem.point.map((val, idx) => val - lambda * axis[idx])
+    return `${axis.join(`,`)}|${intercept.map((val) => val.toFixed(4)).join(`,`)}`
+  }
+
+  type MaterialGroup = {
+    geometry: BufferGeometry
+    color: string
+    opacity?: number
+    striped?: boolean
+  }
+
+  // Cylinder of given radius/length centered at `center` pointing along unit `dir`
+  const oriented_cylinder = (
+    center: Vector3,
+    dir_unit: Vector3,
+    radius: number,
+    length: number,
+  ): CylinderGeometry =>
+    new CylinderGeometry(radius, radius, length, 12).applyMatrix4(
+      new Matrix4().compose(
+        center,
+        new Quaternion().setFromUnitVectors(UP, dir_unit),
+        UNIT_SCALE,
+      ),
+    )
+
+  // Axes (cylinders) + rotoinversion center markers (spheres), merged per color/radius.
+  // Pure rotations render solid; screw axes render DASHED so the two are
+  // distinguishable at a glance (translation-carrying elements are dashed, as in ITA
+  // diagrams).
+  const axis_groups: MaterialGroup[] = $derived.by(() => {
+    const axis_elements = elements.filter(
+      (elem) =>
+        (elem.kind === `rotation` || elem.kind === `screw` ||
+          elem.kind === `rotoinversion`) &&
+        show_kinds[elem.kind] && elem.axis,
+    )
+    // Drop 2-fold axes coincident with a higher-order axis (4 contains 2, 6 contains
+    // 2 and 3, -4 contains 2, …) to reduce visual clutter. Computed over the VISIBLE
+    // elements only, so 2-folds reappear when their enclosing higher-order kind is
+    // toggled off.
+    const max_order_by_line = new Map<string, number>()
+    for (const elem of axis_elements) {
+      const key = line_key(elem)
+      max_order_by_line.set(key, Math.max(max_order_by_line.get(key) ?? 0, elem.order))
+    }
+
+    const parts_by_group = new Map<string, BufferGeometry[]>()
+    for (const elem of axis_elements) {
+      if (
+        hide_redundant_axes && elem.order < (max_order_by_line.get(line_key(elem)) ?? 0)
+      ) continue
+      const clipped = clip_line_to_cell(elem.point, elem.axis as Vec3, lattice)
+      if (!clipped) continue
+      const [start, end] = clipped
+      const span = new Vector3(...math.subtract(end, start))
+      const length = span.length()
+      if (length < 1e-6) continue
+      const dir_unit = span.clone().normalize()
+      const start_vec = new Vector3(...start)
+
+      // Radius is baked into each geometry, so one merged group per color suffices
+      const color = axis_colors[elem.order] ?? `#777777`
+      const group = parts_by_group.get(color) ?? []
+
+      if (elem.kind === `screw`) {
+        // Dashed cylinder: segments along the axis, touching both cell faces
+        for (const dash of dash_segments(length, screw_dash[0], screw_dash[1])) {
+          const center = start_vec.clone().addScaledVector(dir_unit, dash.center)
+          group.push(oriented_cylinder(center, dir_unit, screw_radius, dash.length))
+        }
+      } else {
+        const center = start_vec.clone().addScaledVector(dir_unit, length / 2)
+        group.push(oriented_cylinder(center, dir_unit, axis_radius, length))
+      }
+      if (elem.kind === `rotoinversion`) {
+        const [cx, cy, cz] = frac_to_cart_direction(elem.point, lattice)
+        group.push(new OctahedronGeometry(inversion_radius * 0.8).translate(cx, cy, cz))
+      }
+      parts_by_group.set(color, group)
+    }
+
+    return [...parts_by_group.entries()].flatMap(([color, geometries]) => {
+      const merged = mergeGeometries(geometries)
+      geometries.forEach((geo) => geo.dispose())
+      return merged ? [{ geometry: merged, color }] : []
+    })
+  })
+
+  // Visible mirror/glide planes clipped to the cell — computed once, shared by the
+  // fill and outline groups below. stripe_dir is the unit Cartesian glide direction
+  // (null for mirrors and translation-less entries).
+  const visible_planes = $derived.by(() => {
+    const planes: {
+      polygon: Vec3[]
+      color: string
+      opacity: number
+      stripe_dir: Vec3 | null
+    }[] = []
+    for (const elem of elements) {
+      if ((elem.kind !== `mirror` && elem.kind !== `glide`) || !elem.axis) continue
+      if (!show_kinds[elem.kind]) continue
+      const polygon = clip_plane_to_cell(elem.point, elem.axis, lattice)
+      if (polygon.length < 3) continue
+      const is_mirror = elem.kind === `mirror`
+      planes.push({
+        polygon,
+        color: is_mirror ? mirror_color : glide_color,
+        opacity: is_mirror ? plane_opacity : glide_opacity,
+        stripe_dir: elem.translation
+          ? math.normalize_vec(frac_to_cart_direction(elem.translation, lattice))
+          : null,
+      })
+    }
+    return planes
+  })
+
+  // Plane FILLS: triangles concatenated per color+opacity into one geometry. Glide
+  // fills carry per-vertex UVs whose U coordinate measures Cartesian distance along
+  // the glide-translation direction, so the stripe alphaMap renders stripes running
+  // along the glide direction — the pattern shows both that the plane glides and
+  // where it translates.
+  const plane_groups: MaterialGroup[] = $derived.by(() => {
+    const groups = new Map<string, { positions: number[]; uvs: number[] }>()
+    for (const { polygon, color, opacity, stripe_dir } of visible_planes) {
+      const group_key = `${color}|${opacity}|${stripe_dir ? 1 : 0}`
+      const group = groups.get(group_key) ?? { positions: [], uvs: [] }
+      // Stripe coordinate: Cartesian distance along the glide direction / period
+      const stripe_u = (vert: Vec3): number =>
+        stripe_dir ? math.dot(vert, stripe_dir) / glide_stripe_period : 0
+      // Fan triangulation of the convex polygon
+      for (let idx = 1; idx < polygon.length - 1; idx++) {
+        for (const vert of [polygon[0], polygon[idx], polygon[idx + 1]]) {
+          group.positions.push(...vert)
+          group.uvs.push(stripe_u(vert), 0.5)
+        }
+      }
+      groups.set(group_key, group)
+    }
+    return [...groups.entries()].map(([group_key, { positions, uvs }]) => {
+      const geometry = new BufferGeometry()
+      geometry.setAttribute(`position`, new BufferAttribute(new Float32Array(positions), 3))
+      geometry.setAttribute(`uv`, new BufferAttribute(new Float32Array(uvs), 2))
+      geometry.computeVertexNormals()
+      const [color, opacity, striped] = group_key.split(`|`)
+      return { geometry, color, opacity: Number(opacity), striped: striped === `1` }
+    })
+  })
+
+  // Opaque plane OUTLINES (line segments per color): crisp borders keep overlapping
+  // translucent planes individually legible instead of blending into a single wash.
+  const plane_edge_groups: MaterialGroup[] = $derived.by(() => {
+    if (!show_plane_edges) return []
+    const segments_by_color = new Map<string, number[]>()
+    for (const { polygon, color } of visible_planes) {
+      const positions = segments_by_color.get(color) ?? []
+      for (let idx = 0; idx < polygon.length; idx++) {
+        positions.push(...polygon[idx], ...polygon[(idx + 1) % polygon.length])
+      }
+      segments_by_color.set(color, positions)
+    }
+    return [...segments_by_color.entries()].map(([color, positions]) => {
+      const geometry = new BufferGeometry()
+      geometry.setAttribute(`position`, new BufferAttribute(new Float32Array(positions), 3))
+      return { geometry, color }
+    })
+  })
+
+  // Inversion centers: faceted octahedra (centrosymmetric, unlike the smooth spheres
+  // used for atoms) merged into a single geometry
+  const inversion_group: MaterialGroup | null = $derived.by(() => {
+    if (!show_kinds.inversion) return null
+    const markers = elements
+      .filter((elem) => elem.kind === `inversion`)
+      .map((elem) => {
+        const [cx, cy, cz] = frac_to_cart_direction(elem.point, lattice)
+        return new OctahedronGeometry(inversion_radius).translate(cx, cy, cz)
+      })
+    if (markers.length === 0) return null
+    const merged = mergeGeometries(markers)
+    markers.forEach((geo) => geo.dispose())
+    return merged ? { geometry: merged, color: inversion_color } : null
+  })
+
+  // Dispose each group's merged geometries when that group recomputes or on unmount. One
+  // $effect per group (rather than a single combined one): the deriveds recompute
+  // independently — e.g. tweaking a plane-only prop rebuilds plane_groups but not
+  // axis_groups — and a combined effect would dispose the still-mounted axis geometry.
+  const dispose_geometries = (groups: { geometry: BufferGeometry }[]) => () =>
+    groups.forEach((group) => group.geometry.dispose())
+  $effect(() => dispose_geometries(axis_groups))
+  $effect(() => dispose_geometries(plane_groups))
+  $effect(() => dispose_geometries(plane_edge_groups))
+  $effect(() => dispose_geometries(inversion_group ? [inversion_group] : []))
+
+  // Dispose the (non-reactive) stripe texture on unmount
+  $effect(() => () => stripe_texture.dispose())
+</script>
+
+{#each axis_groups as group, idx (idx)}
+  <T.Mesh geometry={group.geometry}>
+    <T.MeshStandardMaterial color={group.color} />
+  </T.Mesh>
+{/each}
+
+{#each plane_groups as group, idx (idx)}
+  <T.Mesh geometry={group.geometry}>
+    <T.MeshStandardMaterial
+      color={group.color}
+      transparent
+      opacity={group.opacity}
+      alphaMap={group.striped ? stripe_texture : null}
+      side={DoubleSide}
+      depthWrite={false}
+    />
+  </T.Mesh>
+{/each}
+
+{#each plane_edge_groups as group, idx (idx)}
+  <T.LineSegments geometry={group.geometry}>
+    <T.LineBasicMaterial
+      color={group.color}
+      transparent={plane_edge_opacity < 1}
+      opacity={plane_edge_opacity}
+    />
+  </T.LineSegments>
+{/each}
+
+{#if inversion_group}
+  <T.Mesh geometry={inversion_group.geometry}>
+    <T.MeshStandardMaterial color={inversion_group.color} />
+  </T.Mesh>
+{/if}

package/dist/symmetry/SymmetryElements.svelte.d.ts

@@ -0,0 +1,24 @@
+import type { Matrix3x3, Vec2 } from '../math';
+import type { ShowSymmetryKinds, SymmetryElement } from './symmetry-elements';
+type $$ComponentProps = {
+    elements?: SymmetryElement[];
+    lattice: Matrix3x3;
+    show_kinds?: ShowSymmetryKinds;
+    hide_redundant_axes?: boolean;
+    axis_radius?: number;
+    screw_radius?: number;
+    screw_dash?: Vec2;
+    inversion_radius?: number;
+    plane_opacity?: number;
+    glide_opacity?: number;
+    show_plane_edges?: boolean;
+    plane_edge_opacity?: number;
+    glide_stripe_period?: number;
+    axis_colors?: Record<number, string>;
+    mirror_color?: string;
+    glide_color?: string;
+    inversion_color?: string;
+};
+declare const SymmetryElements: import("svelte").Component<$$ComponentProps, {}, "">;
+type SymmetryElements = ReturnType<typeof SymmetryElements>;
+export default SymmetryElements;

package/dist/symmetry/SymmetryStats.svelte

@@ -5,7 +5,15 @@
   import type { HTMLAttributes } from 'svelte/elements'
   import { SETTINGS_CONFIG } from '../settings'
   import type { SymmetrySettings } from './index'
-  import { default_sym_settings, wyckoff_positions_from_moyo } from './index'
+  import {
+    count_structure_free_params,
+    default_sym_settings,
+    enrich_wyckoff_rows,
+    spacegroup_settings,
+    spacegroup_wyckoff_positions,
+    wyckoff_positions_from_moyo,
+    wyckoff_sequence,
+  } from './index'
   import * as spg from './spacegroups'
 
   type SymmetrySnippet = Snippet<
@@ -29,10 +37,35 @@
     header?: SymmetrySnippet
   } = $props()
 
-  const wyckoff_count = $derived(
-    sym_data ? wyckoff_positions_from_moyo(sym_data).length : 0,
+  // Space-group Wyckoff database + all settings (empty when WASM is not initialized,
+  // e.g. SSR or unit tests — the stats below degrade gracefully)
+  const wyckoff_db = $derived(
+    sym_data ? spacegroup_wyckoff_positions(sym_data.hall_number) : [],
   )
-  const display_hm_symbol = $derived(sym_data?.hm_symbol?.replace(/\s+/g, ``) ?? `?`)
+  const settings_entries = $derived(sym_data ? spacegroup_settings(sym_data.number) : [])
+  const current_setting = $derived(
+    settings_entries.find((entry) => entry.hall_number === sym_data?.hall_number),
+  )
+  const occupied_rows = $derived(
+    sym_data ? enrich_wyckoff_rows(wyckoff_positions_from_moyo(sym_data), wyckoff_db) : [],
+  )
+  const wyckoff_count = $derived(occupied_rows.length)
+  const wyckoff_seq = $derived(wyckoff_sequence(occupied_rows))
+  // Internal degrees of freedom (null when ITA coordinates are unavailable)
+  const free_params = $derived(count_structure_free_params(occupied_rows))
+  const display_hm_symbol = $derived(sym_data?.hm_symbol?.replaceAll(/\s+/g, ``) ?? `?`)
+  // Crystal system, plus the lattice system in parens when it differs (e.g. trigonal
+  // space groups split into rhombohedral R-centered and hexagonal P lattices)
+  const crystal_system_label = $derived.by(() => {
+    if (!sym_data) return `?`
+    const crystal_sys = spg.spacegroup_to_crystal_sys(sym_data.number)
+    const lattice_sys = spg.spacegroup_num_to_lattice_system(sym_data.number)
+    if (!crystal_sys) return `?`
+    const suffix = lattice_sys && lattice_sys !== crystal_sys
+      ? ` (${lattice_sys} lattice)`
+      : ``
+    return `${crystal_sys}${suffix}`
+  })
 
   const sym_ops_counts = $derived.by(() => {
     const EPS = 1e-10
@@ -63,7 +96,7 @@
     space_group:
       `International Tables Space group number (1-230) - unique identifier for each space group. Higher numbers indicate more symmetries in the crystal.`,
     crystal_system:
-      `Crystal system classification based on the unit cell symmetry. Seven systems: triclinic, monoclinic, orthorhombic, tetragonal, trigonal, hexagonal, and cubic.`,
+      `Crystal system classification based on the unit cell symmetry. Seven systems: triclinic, monoclinic, orthorhombic, tetragonal, trigonal, hexagonal, and cubic. For trigonal space groups, the lattice system (rhombohedral for R-centered groups, hexagonal otherwise) is shown in parentheses when it differs from the crystal system.`,
     hermann_mauguin:
       `Hermann-Mauguin symbol describes symmetry operations. Format: Lattice type + Point group symmetry. Example: P4/mmm = Primitive + 4-fold rotation + mirror planes`,
     hall_number:
@@ -74,6 +107,12 @@
       `Total symmetry operations that map the crystal structure onto itself. Includes rotations, translations, and combinations.`,
     distinct_orbits:
       `Number of unique Wyckoff positions (symmetry-equivalent atomic sites) in the crystal structure.`,
+    wyckoff_sequence:
+      `Wyckoff sequence: letters of all occupied Wyckoff positions in descending alphabetical order, with superscript counts for letters occupied by multiple orbits. A standard structure-type fingerprint (complements the Pearson symbol).`,
+    free_params:
+      `Internal degrees of freedom: number of free fractional-coordinate parameters (x, y, z in the ITA representative coordinates) summed over occupied Wyckoff orbits. 0 means all atomic positions are fully fixed by symmetry.`,
+    settings:
+      `All settings of this space group in the International Tables (origin choices, unique axes, cell choices, hexagonal vs rhombohedral axes). The setting detected for this structure is highlighted.`,
     translations: `Number of translations in the crystal structure.`,
     rotations: `Number of rotations in the crystal structure.`,
     roto_translations: `Number of roto-translations in the crystal structure.`,
@@ -83,7 +122,7 @@
   function get_step_from_order_of_magnitude(value: number): number {
     if (!Number.isFinite(value) || value <= 0) return 1e-5
     const exponent = Math.floor(Math.log10(value))
-    return Math.pow(10, exponent)
+    return 10 ** exponent
   }
 
   const symprec_step = $derived(get_step_from_order_of_magnitude(settings.symprec))
@@ -151,10 +190,14 @@
         Space Group <strong>{sym_data.number} ({display_hm_symbol})</strong>
       </div>
       <div title={tooltips?.crystal_system} {@attach tooltip()}>
-        Crystal System <strong>{spg.spacegroup_to_crystal_sys(sym_data.number)}</strong>
+        Crystal System <strong>{crystal_system_label}</strong>
       </div>
       <div title={tooltips?.hall_number} {@attach tooltip()}>
-        Hall Number <strong>{sym_data.hall_number}</strong>
+        Hall Number <strong>
+          {sym_data.hall_number}{
+            current_setting ? ` (${current_setting.hall_symbol})` : ``
+          }
+        </strong>
       </div>
       <div title={tooltips?.pearson_symbol} {@attach tooltip()}>
         Pearson <strong>{sym_data.pearson_symbol}</strong>
@@ -162,6 +205,16 @@
       <div title={tooltips?.distinct_orbits} {@attach tooltip()}>
         Wyckoff Positions <strong>{wyckoff_count}</strong>
       </div>
+      {#if wyckoff_seq}
+        <div title={tooltips?.wyckoff_sequence} {@attach tooltip()}>
+          Wyckoff Sequence <strong>{wyckoff_seq}</strong>
+        </div>
+      {/if}
+      {#if free_params !== null}
+        <div title={tooltips?.free_params} {@attach tooltip()}>
+          Free Parameters <strong>{free_params}</strong>
+        </div>
+      {/if}
       <div
         class="sym-ops-summary"
         title="{sym_ops_counts.translations} translations + {sym_ops_counts.rotations} rotations + {sym_ops_counts.roto_translations} roto-translations"
@@ -173,6 +226,39 @@
         }RT)
       </div>
     </div>
+    {#if settings_entries.length > 1}
+      <details class="settings-explorer">
+        <summary title={tooltips?.settings} {@attach tooltip()}>
+          {settings_entries.length} settings of space group {sym_data.number}
+        </summary>
+        <table>
+          <thead>
+            <tr>
+              <th>Hall #</th>
+              <th>H-M full</th>
+              <th>Hall symbol</th>
+              <th>Setting</th>
+              <th>Centering</th>
+            </tr>
+          </thead>
+          <tbody>
+            {#each settings_entries as entry (entry.hall_number)}
+              {@const is_current = entry.hall_number === sym_data.hall_number}
+              <tr
+                class:current={is_current}
+                title={is_current ? `Setting detected for this structure` : null}
+              >
+                <td>{entry.hall_number}</td>
+                <td>{entry.hm_full}</td>
+                <td>{entry.hall_symbol}</td>
+                <td>{entry.setting || `—`}</td>
+                <td>{entry.centering}</td>
+              </tr>
+            {/each}
+          </tbody>
+        </table>
+      </details>
+    {/if}
   {:else}
     <div class="no-data">
       <p>No symmetry data available</p>
@@ -210,6 +296,25 @@
   .stats-grid strong {
     margin: var(--sym-stats-strong-margin, 0 0 0 3pt);
   }
+  .settings-explorer {
+    margin-block: var(--sym-stats-grid-margin-block, 1ex);
+  }
+  .settings-explorer summary {
+    cursor: pointer;
+    color: var(--text-muted, #666);
+  }
+  .settings-explorer table {
+    margin-top: 1ex;
+    border-collapse: collapse;
+  }
+  .settings-explorer :is(th, td) {
+    padding: 1px 8px;
+    text-align: left;
+  }
+  .settings-explorer tr.current {
+    background: color-mix(in srgb, var(--accent-color, #0066cc) 18%, transparent);
+    font-weight: bold;
+  }
   .no-data {
     display: flex;
     align-items: center;