nexaai 1.0.19rc7__cp310-cp310-macosx_14_0_universal2.whl → 1.0.19rc8__cp310-cp310-macosx_14_0_universal2.whl

nexaai/binds/nexa_mlx/py-lib/vlm/modeling/models/aya_vision/vision.py

@@ -1,503 +0,0 @@
-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/binds/nexa_mlx/py-lib/vlm/modeling/models/base.py

@@ -1,202 +0,0 @@
-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