nexaai 1.0.29__cp310-cp310-macosx_14_0_universal2.whl

nexaai/mlx_backend/vlm/modeling/models/aya_vision/vision.py

@@ -0,0 +1,503 @@
+import inspect
+from dataclasses import dataclass
+from typing import Optional
+
+import mlx.core as mx
+import mlx.nn as nn
+import numpy as np
+
+
+@dataclass
+class VisionConfig:
+    model_type: str
+    hidden_size: int
+    num_attention_heads: int
+    patch_size: int
+    num_hidden_layers: int = 12
+    intermediate_size: int = 3072
+    image_size: int = 224
+    num_channels: int = 3
+    layer_norm_eps: float = 1e-6
+
+    @classmethod
+    def from_dict(cls, params):
+        return cls(
+            **{
+                k: v
+                for k, v in params.items()
+                if k in inspect.signature(cls).parameters
+            }
+        )
+
+
+def check_array_shape(arr):
+    shape = arr.shape
+
+    # Check if the shape has 4 dimensions
+    if len(shape) != 4:
+        return False
+
+    out_channels, kH, KW, _ = shape
+
+    # Check if out_channels is the largest, and kH and KW are the same
+    if (out_channels >= kH) and (out_channels >= KW) and (kH == KW):
+        return True
+    else:
+        return False
+
+
+class Attention(nn.Module):
+    def __init__(
+        self,
+        dims: int,
+        num_heads: int,
+        query_input_dims: Optional[int] = None,
+        key_input_dims: Optional[int] = None,
+        value_input_dims: Optional[int] = None,
+        value_dims: Optional[int] = None,
+        value_output_dims: Optional[int] = None,
+        bias: bool = True,
+    ):
+        super().__init__()
+
+        if (dims % num_heads) != 0:
+            raise ValueError(
+                "The input feature dimensions should be divisible by the "
+                f"number of heads ({dims} % {num_heads}) != 0"
+            )
+
+        query_input_dims = query_input_dims or dims
+        key_input_dims = key_input_dims or dims
+        value_input_dims = value_input_dims or key_input_dims
+        value_dims = value_dims or dims
+        value_output_dims = value_output_dims or dims
+
+        self.num_heads = num_heads
+        head_dim = dims // num_heads
+        self.scale = head_dim**-0.5
+
+        self.q_proj = nn.Linear(query_input_dims, dims, bias=bias)
+        self.k_proj = nn.Linear(key_input_dims, dims, bias=bias)
+        self.v_proj = nn.Linear(value_input_dims, value_dims, bias=bias)
+        self.out_proj = nn.Linear(value_dims, value_output_dims, bias=bias)
+
+    def __call__(self, x, mask=None):
+        queries = self.q_proj(x)
+        keys = self.k_proj(x)
+        values = self.v_proj(x)
+
+        num_heads = self.num_heads
+        B, L, D = queries.shape
+        _, S, _ = keys.shape
+        queries = queries.reshape(B, L, num_heads, -1).transpose(0, 2, 1, 3)
+        keys = keys.reshape(B, S, num_heads, -1).transpose(0, 2, 1, 3)
+        values = values.reshape(B, S, num_heads, -1).transpose(0, 2, 1, 3)
+
+        output = mx.fast.scaled_dot_product_attention(
+            queries, keys, values, scale=self.scale, mask=mask
+        )
+        output = output.transpose(0, 2, 1, 3).reshape(B, L, -1)
+        return self.out_proj(output)
+
+
+class MLP(nn.Module):
+    def __init__(self, config: VisionConfig):
+        super().__init__()
+        self.activation_fn = nn.GELU(approx="precise")
+        self.fc1 = nn.Linear(config.hidden_size, config.intermediate_size, bias=True)
+        self.fc2 = nn.Linear(config.intermediate_size, config.hidden_size, bias=True)
+
+    def __call__(self, x: mx.array) -> mx.array:
+        x = self.fc1(x)
+        x = self.activation_fn(x)
+        x = self.fc2(x)
+        return x
+
+
+class EncoderLayer(nn.Module):
+    def __init__(self, config: VisionConfig):
+        super().__init__()
+        self.embed_dim = config.hidden_size
+        self.self_attn = Attention(
+            config.hidden_size, config.num_attention_heads, bias=True
+        )
+        self.layer_norm1 = nn.LayerNorm(self.embed_dim, eps=config.layer_norm_eps)
+        self.mlp = MLP(config)
+        self.layer_norm2 = nn.LayerNorm(self.embed_dim, eps=config.layer_norm_eps)
+
+    def __call__(self, x: mx.array, mask: Optional[mx.array] = None) -> mx.array:
+        r = self.self_attn(self.layer_norm1(x), mask)
+        h = x + r
+        r = self.mlp(self.layer_norm2(h))
+        return h + r
+
+
+class Encoder(nn.Module):
+    def __init__(self, config: VisionConfig):
+        super().__init__()
+        self.layers = [EncoderLayer(config) for _ in range(config.num_hidden_layers)]
+
+    def __call__(
+        self,
+        x: mx.array,
+        output_hidden_states: Optional[bool] = None,
+        mask: Optional[mx.array] = None,
+    ) -> mx.array:
+        encoder_states = (x,) if output_hidden_states else None
+        h = x
+        for l in self.layers:
+            x = l(x, mask=mask)
+            if output_hidden_states:
+                encoder_states = encoder_states + (x,)
+
+            h = x
+
+        return (h, encoder_states)
+
+
+def gaussian_blur_axis(image, sigma, axis):
+    """
+    Applies a 1D Gaussian blur along the given axis.
+    This version works for arrays with any number of dimensions.
+    """
+    radius = int(3 * sigma)
+    if radius < 1:
+        return image
+    x = mx.arange(-radius, radius + 1)
+    kernel = mx.exp(-(x**2) / (2 * sigma**2))
+    kernel = kernel / mx.sum(kernel)
+
+    # MLX doesn't have a direct apply_along_axis equivalent,
+    # so we'll implement the convolution differently based on the axis
+
+    # Helper function to apply 1D convolution along specific axis
+    def conv_1d(array, kernel, axis):
+        # Reshape kernel to broadcast along the right dimensions
+        kernel_shape = [1] * image.ndim
+        kernel_shape[axis] = len(kernel)
+        kernel_reshaped = kernel.reshape(kernel_shape)
+
+        # Pad the array
+        pad_width = [(0, 0)] * image.ndim
+        pad_width[axis] = (radius, radius)
+        padded = mx.pad(array, pad_width, mode="edge")
+
+        # Perform convolution via sliding window sum
+        result = mx.zeros_like(array)
+        slices = [slice(None)] * padded.ndim
+
+        for i in range(2 * radius + 1):
+            slices[axis] = slice(i, i + array.shape[axis])
+            result = result + padded[tuple(slices)] * kernel_reshaped
+
+        return result
+
+    return conv_1d(image, kernel, axis)
+
+
+def bilinear_interpolate(image, new_height, new_width, align_corners=False):
+    """
+    Performs bilinear interpolation on an array whose spatial dimensions are the first two.
+    It supports extra dimensions (e.g. channels or batch dimensions that have been moved to the trailing axes).
+    """
+    # image is assumed to have shape (H, W, ...) where H and W are spatial dimensions.
+    H_in, W_in = image.shape[0], image.shape[1]
+
+    # Compute sampling positions in the input image.
+    if new_height == 1:
+        row_positions = mx.array([0.0])
+    else:
+        if align_corners:
+            row_positions = mx.linspace(0, H_in - 1, new_height)
+        else:
+            row_positions = (mx.arange(new_height) + 0.5) * H_in / new_height - 0.5
+
+    if new_width == 1:
+        col_positions = mx.array([0.0])
+    else:
+        if align_corners:
+            col_positions = mx.linspace(0, W_in - 1, new_width)
+        else:
+            col_positions = (mx.arange(new_width) + 0.5) * W_in / new_width - 0.5
+
+    # Compute floor and ceil indices.
+    row_floor = mx.floor(row_positions).astype(mx.int32)
+    col_floor = mx.floor(col_positions).astype(mx.int32)
+    row_ceil = row_floor + 1
+    col_ceil = col_floor + 1
+
+    row_floor = mx.clip(row_floor, 0, H_in - 1)
+    row_ceil = mx.clip(row_ceil, 0, H_in - 1)
+    col_floor = mx.clip(col_floor, 0, W_in - 1)
+    col_ceil = mx.clip(col_ceil, 0, W_in - 1)
+
+    row_weight = row_positions - row_floor  # shape (new_height,)
+    col_weight = col_positions - col_floor  # shape (new_width,)
+
+    # Use advanced indexing for gather operations
+    # Create meshgrid for coordinates
+    row_floor_grid, col_floor_grid = mx.meshgrid(row_floor, col_floor, indexing="ij")
+    row_ceil_grid, col_floor_grid = mx.meshgrid(row_ceil, col_floor, indexing="ij")
+    row_floor_grid, col_ceil_grid = mx.meshgrid(row_floor, col_ceil, indexing="ij")
+    row_ceil_grid, col_ceil_grid = mx.meshgrid(row_ceil, col_ceil, indexing="ij")
+
+    # Gather the four surrounding pixels using take_along_axis
+    # For higher dimensional arrays, we'll need to reshape and broadcast
+    extra_dims = image.ndim - 2
+
+    def gather_pixels(row_indices, col_indices):
+        # Flatten the spatial dimensions for gathering
+        flat_indices = row_indices * W_in + col_indices
+        flat_image = mx.reshape(image, (-1,) + image.shape[2:])
+        # Gather and reshape back
+        gathered = mx.take(flat_image, flat_indices.reshape(-1), axis=0)
+        return mx.reshape(gathered, (new_height, new_width) + image.shape[2:])
+
+    top_left = gather_pixels(row_floor_grid, col_floor_grid)
+    top_right = gather_pixels(row_floor_grid, col_ceil_grid)
+    bottom_left = gather_pixels(row_ceil_grid, col_floor_grid)
+    bottom_right = gather_pixels(row_ceil_grid, col_ceil_grid)
+
+    # Expand the weights to have shape (new_height, new_width, *[1]*extra_dims)
+    r_weight = row_weight.reshape(new_height, 1, *([1] * extra_dims))
+    c_weight = col_weight.reshape(1, new_width, *([1] * extra_dims))
+
+    # Perform bilinear interpolation.
+    result = (
+        (1 - r_weight) * (1 - c_weight) * top_left
+        + (1 - r_weight) * c_weight * top_right
+        + r_weight * (1 - c_weight) * bottom_left
+        + r_weight * c_weight * bottom_right
+    )
+    return result
+
+
+def resize_bilinear(image, new_size, align_corners=False, antialias=True):
+    """
+    Resizes an image (or embedding tensor) to new_size=(new_height, new_width)
+    using bilinear interpolation with MLX.
+
+    Supports:
+      - 2D: (H, W)
+      - 3D: (H, W, C)
+      - 4D: (B, C, H, W)  (assumed for typical image batches)
+    """
+    new_height, new_width = new_size
+
+    # Convert numpy arrays to MLX arrays if needed
+    if isinstance(image, np.ndarray):
+        image = mx.array(image)
+
+    if image.ndim == 2 or image.ndim == 3:
+        # Assume spatial dims are the first two.
+        resized = image
+        H_in, W_in = image.shape[:2]
+        if antialias:
+            if new_height < H_in:
+                scale_y = new_height / H_in
+                sigma_y = (1 / scale_y - 1) / 2.0  # heuristic
+                if sigma_y > 0:
+                    resized = gaussian_blur_axis(resized, sigma_y, axis=0)
+            if new_width < W_in:
+                scale_x = new_width / W_in
+                sigma_x = (1 / scale_x - 1) / 2.0
+                if sigma_x > 0:
+                    resized = gaussian_blur_axis(resized, sigma_x, axis=1)
+        resized = bilinear_interpolate(
+            resized, new_height, new_width, align_corners=align_corners
+        )
+        return resized
+
+    elif image.ndim == 4:
+        # Assume shape is (B, C, H, W) (typical PyTorch/MLX format).
+        B, C, H_in, W_in = image.shape
+        # Permute to bring spatial dims to the front: (H, W, B, C)
+        image_perm = mx.transpose(image, (2, 3, 0, 1))
+        resized = image_perm
+        if antialias:
+            if new_height < H_in:
+                scale_y = new_height / H_in
+                sigma_y = (1 / scale_y - 1) / 2.0
+                if sigma_y > 0:
+                    resized = gaussian_blur_axis(resized, sigma_y, axis=0)
+            if new_width < W_in:
+                scale_x = new_width / W_in
+                sigma_x = (1 / scale_x - 1) / 2.0
+                if sigma_x > 0:
+                    resized = gaussian_blur_axis(resized, sigma_x, axis=1)
+        resized = bilinear_interpolate(
+            resized, new_height, new_width, align_corners=align_corners
+        )
+        # Permute back to (B, C, new_height, new_width)
+        resized = mx.transpose(resized, (2, 3, 0, 1))
+        return resized
+
+    else:
+        raise ValueError("Unsupported image dimensions.")
+
+
+class VisionEmbeddings(nn.Module):
+    def __init__(self, config: VisionConfig):
+        super().__init__()
+        self.config = config
+        self.embed_dim = config.hidden_size
+        self.image_size = config.image_size
+        self.patch_size = config.patch_size
+
+        self.patch_embedding = nn.Conv2d(
+            config.num_channels,
+            config.hidden_size,
+            kernel_size=self.patch_size,
+            stride=self.patch_size,
+        )
+
+        self.num_patches = (self.image_size // self.patch_size) ** 2
+        self.num_positions = self.num_patches
+        self.position_embedding = nn.Embedding(self.num_positions, self.embed_dim)
+
+    @staticmethod
+    def resize_positional_embeddings(
+        positional_embeddings: mx.array,
+        spatial_shapes: mx.array,
+        max_length: int,
+    ) -> mx.array:
+        """
+        Resize positional embeddings to image-specific size and pad to a fixed size.
+
+        Args:
+            positional_embeddings (`torch.Tensor`):
+                Position embeddings of shape (height, width, embed_dim)
+            spatial_shapes (`torch.LongTensor`):
+                Spatial shapes of shape (batch_size, 2) to resize the positional embeddings to
+            max_length (`int`):
+                Maximum length of the positional embeddings to pad resized positional embeddings to
+
+        Returns:
+            `torch.Tensor`: Embeddings of shape (batch_size, max_length, embed_dim)
+        """
+        batch_size = spatial_shapes.shape[0]
+        embed_dim = positional_embeddings.shape[-1]
+        source_dtype = positional_embeddings.dtype
+
+        resulted_positional_embeddings = mx.zeros(
+            (batch_size, max_length, embed_dim)
+        ).astype(source_dtype)
+
+        # (height, width, embed_dim) -> (1, embed_dim, height, width) for interpolation
+        positional_embeddings = positional_embeddings.transpose(2, 0, 1).reshape(
+            1, embed_dim, -1
+        )
+
+        # Upcast to float32 on CPU because antialias is not supported for bfloat16/float16 on CPU
+        if positional_embeddings.device.type == "cpu":
+            positional_embeddings = positional_embeddings.astype(mx.float32)
+
+        for i in range(batch_size):
+            # (1, dim, height, width) -> (1, dim, target_height, target_width)
+            height, width = spatial_shapes[i]
+            # Then upsample width dimension
+            resized_embeddings = resize_bilinear(
+                positional_embeddings,
+                (height, width),
+                align_corners=False,
+                antialias=True,
+            )
+
+            # (1, dim, target_height, target_width) -> (target_height * target_width, dim)
+            resized_embeddings = resized_embeddings.reshape(
+                embed_dim, height * width
+            ).transpose(0, 1)
+
+            # Cast to original dtype
+            resized_embeddings = resized_embeddings.astype(source_dtype)
+
+            resulted_positional_embeddings[i, : height * width] = resized_embeddings
+            resulted_positional_embeddings[i, height * width :] = resized_embeddings[0]
+
+        return resulted_positional_embeddings
+
+    def __call__(
+        self, x: mx.array, spatial_shapes: Optional[mx.array] = None
+    ) -> mx.array:
+        batch_size = x.shape[0]
+        patch_embeddings = self.patch_embedding(x)
+        patch_embeddings = mx.flatten(patch_embeddings, start_axis=1, end_axis=2)
+        if spatial_shapes is None:
+            position_ids = mx.array(np.arange(self.num_positions)[None, :])
+            embeddings = patch_embeddings
+            embeddings += self.position_embedding(position_ids)
+
+        else:
+            # Get positional resized and padded positional embeddings
+            positional_embeddings = self.position_embedding.weight.reshape(
+                self.position_embedding_size, self.position_embedding_size, -1
+            )
+
+            resized_positional_embeddings = self.resize_positional_embeddings(
+                positional_embeddings, spatial_shapes, max_length=x.shape[1]
+            )
+
+            # Add positional embeddings to patch embeddings
+            embeddings = patch_embeds + resized_positional_embeddings
+        return embeddings
+
+
+class SigLipVisionModel(nn.Module):
+    def __init__(self, config: VisionConfig):
+        super().__init__()
+
+        self.embeddings = VisionEmbeddings(config)
+        self.encoder = Encoder(config)
+        self.post_layernorm = nn.LayerNorm(config.hidden_size)
+
+    def __call__(
+        self,
+        x: mx.array,
+        spatial_shapes: mx.array,
+        output_hidden_states: Optional[bool] = None,
+    ) -> mx.array:
+        x = self.embeddings(x, spatial_shapes)
+        x = x.astype(self.embeddings.patch_embedding.weight.dtype)
+        encoder_outputs = self.encoder(
+            x=x, output_hidden_states=output_hidden_states, mask=None
+        )
+        pooler_output = self.post_layernorm(encoder_outputs[0])
+        return pooler_output, x, encoder_outputs[-1]
+
+
+class VisionModel(nn.Module):
+    def __init__(self, config: VisionConfig):
+        super().__init__()
+        self.model_type = config.model_type
+        if self.model_type not in ["siglip_vision_model"]:
+            raise ValueError(f"Unsupported model type: {self.model_type}")
+
+        self.vision_model = SigLipVisionModel(config)
+
+    def __call__(
+        self,
+        x: mx.array,
+        spatial_shapes: Optional[mx.array] = None,
+        output_hidden_states: Optional[bool] = None,
+    ) -> mx.array:
+        return self.vision_model(x, spatial_shapes, output_hidden_states)
+
+    def sanitize(self, weights):
+        sanitized_weights = {}
+        for k, v in weights.items():
+            if "position_ids" in k:
+                # Remove unused position_ids
+                continue
+            elif "patch_embedding.weight" in k:
+                # PyTorch conv2d weight tensors have shape:
+                #   [out_channels, in_channels, kH, KW]
+                # MLX conv2d expects the weight be of shape:
+                #   [out_channels, kH, KW, in_channels]
+                if check_array_shape(v):
+                    sanitized_weights[k] = v
+                else:
+                    sanitized_weights[k] = v.transpose(0, 2, 3, 1)
+            else:
+                sanitized_weights[k] = v
+
+        return sanitized_weights

nexaai/mlx_backend/vlm/modeling/models/base.py

@@ -0,0 +1,202 @@
+import math
+from abc import ABC, abstractmethod
+from dataclasses import dataclass
+from typing import Any, Dict, List, Optional
+
+import mlx.core as mx
+import mlx.nn as nn
+from mlx_lm.models.base import create_attention_mask, scaled_dot_product_attention
+from mlx_lm.models.cache import RotatingKVCache
+from PIL import Image
+from transformers.image_processing_utils import BaseImageProcessor as ImageProcessor
+from transformers.image_processing_utils import get_size_dict
+from transformers.image_utils import ChannelDimension, PILImageResampling
+
+
+@dataclass
+class LanguageModelOutput:
+    logits: mx.array
+    cross_attention_states: Optional[List[mx.array]] = None
+    encoder_outputs: Optional[List[mx.array]] = None
+
+
+def expand2square(pil_img, background_color):
+    width, height = pil_img.size
+    if width == height:
+        return pil_img
+    elif width > height:
+        result = Image.new(pil_img.mode, (width, width), background_color)
+        result.paste(pil_img, (0, (width - height) // 2))
+        return result
+    else:
+        result = Image.new(pil_img.mode, (height, height), background_color)
+        result.paste(pil_img, ((height - width) // 2, 0))
+        return result
+
+
+def check_array_shape(arr):
+    shape = arr.shape
+
+    # Check if the shape has 4 dimensions
+    if len(shape) == 4:
+        out_channels, kH, KW, _ = shape
+        # Check if out_channels is the largest, and kH and KW are the same
+        if (out_channels >= kH) and (out_channels >= KW) and (kH == KW):
+            return True
+        else:
+            return False
+    # Check if the shape has 3 dimensions
+    elif len(shape) == 3:
+        _, kW, out_channels = shape
+        # Check if out_channels is the largest
+        if kW >= out_channels:
+            return True
+        else:
+            return False
+    else:
+        return False
+
+
+class BaseImageProcessor(ImageProcessor):
+    def __init__(
+        self,
+        image_mean=(0.5, 0.5, 0.5),
+        image_std=(0.5, 0.5, 0.5),
+        size=(384, 384),
+        crop_size: Dict[str, int] = None,
+        resample=PILImageResampling.BICUBIC,
+        rescale_factor=1 / 255,
+        data_format=ChannelDimension.FIRST,
+    ):
+        crop_size = (
+            crop_size if crop_size is not None else {"height": 384, "width": 384}
+        )
+        crop_size = get_size_dict(
+            crop_size, default_to_square=True, param_name="crop_size"
+        )
+
+        self.image_mean = image_mean
+        self.image_std = image_std
+        self.size = size
+        self.resample = resample
+        self.rescale_factor = rescale_factor
+        self.data_format = data_format
+        self.crop_size = crop_size
+
+    @abstractmethod
+    def preprocess(self, images):
+        pass
+
+
+# Add this code to visualize the chunked attention mask
+def visualize_attention_mask(mask):
+    """Visualize attention mask with symbols for better readability."""
+    if mask is None:
+        print("No mask")
+        return
+
+    seq_len = mask.shape[0]
+
+    print("        ", end="")
+    for i in range(seq_len):
+        print(f"{i:2d} ", end="")
+    print()
+
+    for i in range(seq_len):
+        print(f"Token {i:2d}: ", end="")
+        for j in range(seq_len):
+            if mask[i, j]:
+                print(" ■ ", end="")
+            else:
+                print(" ⬚ ", end="")
+        print()
+
+
+def check_activation_stats(name, tensor):
+    """Helper function to check for anomalies and log stats."""
+
+    print(f"--- Activation Stats: {name} ---")
+    # Check for NaNs/Infs
+    has_nan = mx.isnan(tensor).any()
+    has_inf = mx.isinf(tensor).any()
+    if has_nan:
+        print(f"WARNING: Found NaN in {name}")
+    if has_inf:
+        print(f"WARNING: Found Inf in {name}")
+
+    # Calculate and print stats (ensure computation happens)
+    min_val = mx.min(tensor).item()
+    max_val = mx.max(tensor).item()
+    mean_val = mx.mean(tensor).item()
+    std_val = mx.std(tensor).item()
+    print(f"  Shape: {tensor.shape}")
+    print(f"  Min: {min_val:.4f}, Max: {max_val:.4f}")
+    print(f"  Mean: {mean_val:.4f}, Std: {std_val:.4f}")
+    print("-" * (len(name) + 24))
+
+
+def pixel_shuffle(input_tensor, shuffle_ratio):
+    # input_tensor: [batch_size, num_patches, channels]
+    batch_size, num_patches, channels = input_tensor.shape
+    patch_size = int(math.sqrt(num_patches))
+
+    input_tensor = input_tensor.reshape(batch_size, patch_size, patch_size, -1)
+    batch_size, height, width, channels = input_tensor.shape
+
+    reshaped_tensor = input_tensor.reshape(
+        batch_size, height, int(width * shuffle_ratio), int(channels / shuffle_ratio)
+    )
+    reshaped_tensor = reshaped_tensor.transpose(0, 2, 1, 3)
+
+    reshaped_tensor = reshaped_tensor.reshape(
+        batch_size,
+        int(height * shuffle_ratio),
+        int(width * shuffle_ratio),
+        int(channels / (shuffle_ratio**2)),
+    )
+    reshaped_tensor = reshaped_tensor.transpose(0, 2, 1, 3)
+
+    output_tensor = reshaped_tensor.reshape(batch_size, -1, reshaped_tensor.shape[-1])
+    return output_tensor
+
+
+def interpolate(pos_embed, size, mode="cubic", align_corners=False):
+    """
+    MLX implementation of PyTorch's F.interpolate with bicubic mode
+
+    Args:
+        pos_embed: MLX array with shape [B, C, H_src, W_src] or [C, H_src, W_src]
+        size: Tuple (H_dst, W_dst) - target size
+        align_corners: Boolean - whether to align corners
+
+    Returns:
+        Interpolated array with shape [B, C, H_dst, W_dst] or [C, H_dst, W_dst]
+    """
+    # Handle different input shapes
+    input_dim = pos_embed.ndim
+    original_shape = pos_embed.shape
+
+    if input_dim == 3:
+        # [C, H, W] -> [1, C, H, W]
+        pos_embed = pos_embed.reshape(1, *original_shape)
+
+    # Get source dimensions
+    h_src, w_src = pos_embed.shape[-2:]
+    h_dst, w_dst = size
+
+    # Calculate scale factors
+    scale_h = h_dst / h_src
+    scale_w = w_dst / w_src
+
+    # Create upsampler
+    upsampler = nn.Upsample(
+        scale_factor=(scale_h, scale_w), mode=mode, align_corners=align_corners
+    )
+
+    # Apply upsampling
+    result = upsampler(pos_embed)
+
+    # Return in the original dimension format
+    if input_dim == 3:
+        return result.reshape(original_shape[0], *size)
+    return result