def METHOD_NAME(coords1: np.ndarray, coords2: np.ndarray) -> np.ndarray:
    """Returns matrix of pairwise Euclidean distances.
    Parameters
    ----------
    coords1: np.ndarray
        A numpy array of shape `(N, 3)`
    coords2: np.ndarray
        A numpy array of shape `(M, 3)`
    Returns
    -------
    np.ndarray
        A `(N,M)` array with pairwise distances.
    """
    return np.sum((coords1[None, :] - coords2[:, None])**2, -1)**0.5