diff --git a/gnm/shape/data/versions/v3_0/gnm_head.npz b/gnm/shape/data/versions/v3_0/gnm_head.npz
index cd3acf68..198711f1 100644
Binary files a/gnm/shape/data/versions/v3_0/gnm_head.npz and b/gnm/shape/data/versions/v3_0/gnm_head.npz differ
diff --git a/gnm/shape/gnm_common.py b/gnm/shape/gnm_common.py
index dd5a3109..3a5f9d67 100644
--- a/gnm/shape/gnm_common.py
+++ b/gnm/shape/gnm_common.py
@@ -14,10 +14,41 @@
 
 """Backend-agnostic core math functions of the GNM model using etils.enp."""
 
+import typing
 from typing import Any, Sequence
 from etils import enp
 
-enpt = enp.typing
+if typing.TYPE_CHECKING:
+
+  class FloatArray:
+
+    @classmethod
+    def __class_getitem__(cls, _):
+      return typing.Any
+
+  class IntArray:
+
+    @classmethod
+    def __class_getitem__(cls, _):
+      return typing.Any
+
+  class BoolArray:
+
+    @classmethod
+    def __class_getitem__(cls, _):
+      return typing.Any
+
+  class _EnpTypingMock:
+    # pylint: disable=invalid-name
+    FloatArray = FloatArray
+    IntArray = IntArray
+    BoolArray = BoolArray
+    # pylint: enable=invalid-name
+
+  enpt = _EnpTypingMock()
+else:
+  enpt = enp.typing
+
 
 _EPSILON = 1e-8
 
diff --git a/gnm/shape/gnm_numpy.py b/gnm/shape/gnm_numpy.py
index b64563ce..3ebe5691 100644
--- a/gnm/shape/gnm_numpy.py
+++ b/gnm/shape/gnm_numpy.py
@@ -147,7 +147,7 @@ def compute_vertex_normals(
     vertex_normals = np.zeros_like(vertices_flat)
     np.add.at(
         vertex_normals,
-        (slice(None), self.triangles, slice(None)),
+        (slice(None), self.triangles, slice(None)),  # pyrefly: ignore[bad-argument-type]
         face_normals_area[:, :, None, :],
     )
 
diff --git a/gnm/shape/gnm_xnp.py b/gnm/shape/gnm_xnp.py
index 396384e9..5f524165 100644
--- a/gnm/shape/gnm_xnp.py
+++ b/gnm/shape/gnm_xnp.py
@@ -21,6 +21,7 @@
 from collections.abc import Mapping, Sequence
 import dataclasses
 import functools
+import typing
 from typing import Any, Self
 
 from absl import logging
@@ -32,7 +33,41 @@
 import numpy as np
 import numpy.typing as npt
 
-enpt = enp.typing
+if typing.TYPE_CHECKING:
+
+  class FloatArray:
+
+    @classmethod
+    def __class_getitem__(cls, _):
+      return typing.Any
+
+  class IntArray:
+
+    @classmethod
+    def __class_getitem__(cls, _):
+      return typing.Any
+
+  class BoolArray:
+
+    @classmethod
+    def __class_getitem__(cls, _):
+      return typing.Any
+
+  class _EnpTypingMock:
+    # pylint: disable=invalid-name
+    FloatArray = FloatArray
+    IntArray = IntArray
+    BoolArray = BoolArray
+    # pylint: enable=invalid-name
+
+  enpt = _EnpTypingMock()
+
+  V = typing.Any
+  VPruned = typing.Any
+  L = typing.Any
+else:
+  enpt = enp.typing
+
 
 _NONZERO_THRESHOLD = 1e-4
 _EPSILON = 1e-8
@@ -107,26 +142,26 @@ class GNM(gnm_base.GNMBase):
 
   version: gnm_specs.GNMVersion
   variant: gnm_specs.GNMVariant
-  template_vertex_positions: enpt.FloatArray
-  template_joint_positions: enpt.FloatArray
-  vertex_identity_basis: enpt.FloatArray
-  joint_identity_basis: enpt.FloatArray
-  expression_basis: enpt.FloatArray
+  template_vertex_positions: enpt.FloatArray['V 3']
+  template_joint_positions: enpt.FloatArray['J 3']
+  vertex_identity_basis: enpt.FloatArray['I V 3']
+  joint_identity_basis: enpt.FloatArray['I J 3']
+  expression_basis: enpt.FloatArray['E V 3']
   identity_names: Sequence[str]
   joint_names: Sequence[str]
   expression_names: Sequence[str]
   joint_parent_indices: Sequence[int]
-  skinning_weights: enpt.FloatArray
-  quads: enpt.IntArray
-  triangles: enpt.IntArray
-  quad_uvs: enpt.FloatArray
-  triangle_uvs: enpt.FloatArray
+  skinning_weights: enpt.FloatArray['J V']
+  quads: enpt.IntArray['Q 4']
+  triangles: enpt.IntArray['T 3']
+  quad_uvs: enpt.FloatArray['Q 4 2']
+  triangle_uvs: enpt.FloatArray['T 3 2']
   mesh_component_names: Sequence[str]
-  mirror_indices: enpt.IntArray
-  joint_regressor: enpt.FloatArray
-  pose_correctives_regressor: enpt.FloatArray
-  bone_aligned_template_joint_orientations: enpt.FloatArray
-  vertex_groups: enpt.FloatArray
+  mirror_indices: enpt.IntArray['V']
+  joint_regressor: enpt.FloatArray['J V']
+  pose_correctives_regressor: enpt.FloatArray['9*J 3*V']
+  bone_aligned_template_joint_orientations: enpt.FloatArray['J 3 3']
+  vertex_groups: enpt.FloatArray['G V']
   vertex_group_names: Sequence[str]
 
   def __init__(self, *args: Any, **kwargs: Any) -> None:
@@ -143,7 +178,8 @@ def __post_init__(self):
     }
     self._xnp = enp.get_np_module(self.template_vertex_positions)
     self._landmarks: dict[
-        gnm_landmarks.GNMLandmarksType, tuple[enpt.IntArray, enpt.FloatArray]
+        gnm_landmarks.GNMLandmarksType,
+        tuple[enpt.IntArray['L'], enpt.FloatArray['L']],
     ] = {}
 
   @property
@@ -571,7 +607,7 @@ def vertex_group(self, name: str) -> npt.NDArray[np.floating]:
 
   def vertex_group_mask(
       self, *names: str, threshold: float = _NONZERO_THRESHOLD
-  ) -> npt.NDArray[bool]:
+  ) -> npt.NDArray[np.bool_]:
     result_mask = np.zeros(self.num_vertices, dtype=bool)
     for name in names:
       operator, inverse = '|', False
@@ -629,7 +665,7 @@ def compute_vertex_normals(
         'Subclasses must implement compute_vertex_normals.'
     )
 
-  def prune_vertices(self, keep_vertices: enpt.IntArray['V_pruned']) -> None:
+  def prune_vertices(self, keep_vertices: enpt.IntArray['VPruned']) -> None:
     """Prunes model vertices in-place."""
     xnp = self.xnp
     num_vertices = self.num_vertices