cubecl_core/frontend/

ranges.rs

use alloc::boxed::Box;
use core::ops::{Range, RangeFrom, RangeFull, RangeInclusive, RangeTo, RangeToInclusive};
use cubecl_ir::{ExpandValue, OpInserter, dialect::branch::RangeLoopOp};
use num_traits::NumCast;

use crate as cubecl;
use cubecl::prelude::*;

#[derive_expand(CubeType)]
pub struct Range<Idx: Scalar> {
    start: Idx,
    end: Idx,
}

impl<Idx: Scalar> RangeExpand<Idx> {
    pub fn new(start: NativeExpand<Idx>, end: NativeExpand<Idx>) -> Self {
        Self { start, end }
    }
}

impl<I: Int> RangeExpand<I> {
    pub fn __expand_step_by_method(
        self,
        _: &Scope,
        n: impl Into<NativeExpand<I>>,
    ) -> SteppedRangeExpand<I> {
        SteppedRangeExpand {
            start: self.start,
            end: self.end,
            step: n.into(),
            inclusive: false,
        }
    }
}

impl<I: Int> Iterable for RangeExpand<I> {
    type Item = NativeExpand<I>;

    fn expand_unroll(
        self,
        scope: &Scope,
        body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType),
    ) {
        iter_expand_unroll(scope, self.start, self.end, false, body);
    }

    fn expand(self, scope: &Scope, body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType)) {
        iter_expand(scope, self.start, self.end, false, body);
    }

    fn const_len(&self) -> Option<usize> {
        let start = self.start.expand.as_const()?.as_i64();
        let end = self.end.expand.as_const()?.as_i64();
        Some(start.abs_diff(end) as usize)
    }
}

impl IntoSliceIndices for RangeExpand<usize> {
    fn into_slice_indices<E: CubePrimitive>(
        self,
        _scope: &Scope,
        _list: &impl ListExpand<E>,
    ) -> (NativeExpand<usize>, NativeExpand<usize>) {
        (self.start, self.end)
    }
}

#[derive_expand(CubeType)]
pub struct RangeFrom<Idx: Scalar> {
    start: Idx,
}

impl<Idx: Scalar> RangeFromExpand<Idx> {
    pub fn new(start: NativeExpand<Idx>) -> Self {
        Self { start }
    }
}

impl IntoSliceIndices for RangeFromExpand<usize> {
    fn into_slice_indices<E: CubePrimitive>(
        self,
        scope: &Scope,
        list: &impl ListExpand<E>,
    ) -> (NativeExpand<usize>, NativeExpand<usize>) {
        let end = list.__expand_len_method(scope);
        (self.start, end)
    }
}

#[derive_expand(CubeType)]
pub struct RangeFull;

impl IntoSliceIndices for RangeFullExpand {
    fn into_slice_indices<E: CubePrimitive>(
        self,
        scope: &Scope,
        list: &impl ListExpand<E>,
    ) -> (NativeExpand<usize>, NativeExpand<usize>) {
        let start = NativeExpand::from_lit(scope, 0);
        let end = list.__expand_len_method(scope);
        (start, end)
    }
}

#[derive_expand(CubeType)]
pub struct RangeInclusive<Idx: Scalar> {
    start: Idx,
    last: Idx,
}

impl<Idx: Scalar> RangeInclusiveExpand<Idx> {
    pub fn new(start: NativeExpand<Idx>, last: NativeExpand<Idx>) -> Self {
        Self { start, last }
    }
}

impl<I: Int> RangeInclusiveExpand<I> {
    pub fn __expand_step_by_method(
        self,
        _: &Scope,
        n: impl Into<NativeExpand<I>>,
    ) -> SteppedRangeExpand<I> {
        SteppedRangeExpand {
            start: self.start,
            end: self.last,
            step: n.into(),
            inclusive: true,
        }
    }
}

impl<I: Int> Iterable for RangeInclusiveExpand<I> {
    type Item = NativeExpand<I>;

    fn expand_unroll(
        self,
        scope: &Scope,
        body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType),
    ) {
        iter_expand_unroll(scope, self.start, self.last, true, body);
    }

    fn expand(self, scope: &Scope, body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType)) {
        iter_expand(scope, self.start, self.last, true, body);
    }

    fn const_len(&self) -> Option<usize> {
        let start = self.start.expand.as_const()?.as_i64();
        let end = self.last.expand.as_const()?.as_i64();
        Some(start.abs_diff(end) as usize + 1)
    }
}

impl IntoSliceIndices for RangeInclusiveExpand<usize> {
    fn into_slice_indices<E: CubePrimitive>(
        self,
        scope: &Scope,
        _list: &impl ListExpand<E>,
    ) -> (NativeExpand<usize>, NativeExpand<usize>) {
        let end = self
            .last
            .__expand_add_method(scope, NativeExpand::from_lit(scope, 1));
        (self.start, end)
    }
}

#[derive_expand(CubeType)]
pub struct RangeTo<Idx: Scalar> {
    end: Idx,
}

impl<Idx: Scalar> RangeToExpand<Idx> {
    pub fn new(end: NativeExpand<Idx>) -> Self {
        Self { end }
    }
}

impl<I: Int> RangeToExpand<I> {
    pub fn __expand_step_by_method(
        self,
        scope: &Scope,
        n: impl Into<NativeExpand<I>>,
    ) -> SteppedRangeExpand<I> {
        SteppedRangeExpand {
            start: NativeExpand::from_lit(scope, I::new(0)),
            end: self.end,
            step: n.into(),
            inclusive: false,
        }
    }
}

impl<I: Int> Iterable for RangeToExpand<I> {
    type Item = NativeExpand<I>;

    fn expand_unroll(
        self,
        scope: &Scope,
        body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType),
    ) {
        let start = NativeExpand::from_lit(scope, I::new(0));
        iter_expand_unroll(scope, start, self.end, false, body);
    }

    fn expand(self, scope: &Scope, body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType)) {
        let start = NativeExpand::from_lit(scope, I::new(0));
        iter_expand(scope, start, self.end, false, body);
    }

    fn const_len(&self) -> Option<usize> {
        Some(self.end.expand.as_const()?.as_usize())
    }
}

impl IntoSliceIndices for RangeToExpand<usize> {
    fn into_slice_indices<E: CubePrimitive>(
        self,
        scope: &Scope,
        _list: &impl ListExpand<E>,
    ) -> (NativeExpand<usize>, NativeExpand<usize>) {
        let start = NativeExpand::from_lit(scope, 0);
        (start, self.end)
    }
}

#[derive_expand(CubeType)]
pub struct RangeToInclusive<Idx: Scalar> {
    last: Idx,
}

impl<Idx: Scalar> RangeToInclusiveExpand<Idx> {
    pub fn new(last: NativeExpand<Idx>) -> Self {
        Self { last }
    }
}

impl<I: Int> RangeToInclusiveExpand<I> {
    pub fn __expand_step_by_method(
        self,
        scope: &Scope,
        n: impl Into<NativeExpand<I>>,
    ) -> SteppedRangeExpand<I> {
        SteppedRangeExpand {
            start: NativeExpand::from_lit(scope, I::new(0)),
            end: self.last,
            step: n.into(),
            inclusive: true,
        }
    }
}

impl<I: Int> Iterable for RangeToInclusiveExpand<I> {
    type Item = NativeExpand<I>;

    fn expand_unroll(
        self,
        scope: &Scope,
        body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType),
    ) {
        let start = NativeExpand::from_lit(scope, I::new(0));
        iter_expand_unroll(scope, start, self.last, true, body);
    }

    fn expand(self, scope: &Scope, body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType)) {
        let start = NativeExpand::from_lit(scope, I::new(0));
        iter_expand(scope, start, self.last, true, body);
    }

    fn const_len(&self) -> Option<usize> {
        Some(self.last.expand.as_const()?.as_usize() + 1)
    }
}

impl IntoSliceIndices for RangeToInclusiveExpand<usize> {
    fn into_slice_indices<E: CubePrimitive>(
        self,
        scope: &Scope,
        _list: &impl ListExpand<E>,
    ) -> (NativeExpand<usize>, NativeExpand<usize>) {
        let start = NativeExpand::from_lit(scope, 0);
        let end = self
            .last
            .__expand_add_method(scope, NativeExpand::from_lit(scope, 1));
        (start, end)
    }
}

pub(crate) trait IntoSliceIndices {
    fn into_slice_indices<E: CubePrimitive>(
        self,
        scope: &Scope,
        list: &impl ListExpand<E>,
    ) -> (NativeExpand<usize>, NativeExpand<usize>);
}

macro_rules! impl_slice_ranges {
    ($ty: ty, $range: ty) => {
        impl<E: CubePrimitive> IndexExpand<$range> for $ty {
            type Output = SliceExpand<E>;

            fn __expand_index_method(&self, scope: &Scope, index: $range) -> &Self::Output {
                let (start, end) = index.into_slice_indices(scope, self);
                self.__expand_slice_method(scope, start, end)
            }
        }

        impl<E: CubePrimitive> IndexMutExpand<$range> for $ty {
            fn __expand_index_mut_method(
                &mut self,
                scope: &Scope,
                index: $range,
            ) -> &mut <Self as IndexExpand<$range>>::Output {
                let (start, end) = index.into_slice_indices(scope, self);
                self.__expand_slice_mut_method(scope, start, end)
            }
        }
    };
    ($ty: ty) => {
        impl_slice_ranges!($ty, RangeExpand<usize>);
        impl_slice_ranges!($ty, RangeFromExpand<usize>);
        impl_slice_ranges!($ty, RangeFullExpand);
        impl_slice_ranges!($ty, RangeInclusiveExpand<usize>);
        impl_slice_ranges!($ty, RangeToExpand<usize>);
        impl_slice_ranges!($ty, RangeToInclusiveExpand<usize>);
    };
}
pub(crate) use impl_slice_ranges;

fn iter_expand_unroll<I: Int>(
    scope: &Scope,
    start: NativeExpand<I>,
    end: NativeExpand<I>,
    inclusive: bool,
    body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType),
) {
    let start = start
        .expand
        .as_const()
        .expect("Only constant start can be unrolled.")
        .as_i64();
    let end = end
        .expand
        .as_const()
        .expect("Only constant end can be unrolled.")
        .as_i64();

    if inclusive {
        for i in start..=end {
            let val = I::from_int(i);
            body(scope, val.into())
        }
    } else {
        for i in start..end {
            let val = I::from_int(i);
            body(scope, val.into())
        }
    }
}

fn iter_expand<I: Int>(
    scope: &Scope,
    start: NativeExpand<I>,
    end: NativeExpand<I>,
    inclusive: bool,
    body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType),
) {
    let start = I::__expand_cast_from(scope, start).expand;
    let mut end = I::__expand_cast_from(scope, end);
    let step: ExpandValue = I::new(1).into();

    if inclusive {
        end = end.__expand_add_method(scope, I::new(1).into());
    }

    let start = start.read_value(scope);
    let end = end.read_value(scope);
    let step = step.read_value(scope);

    let range_loop = RangeLoopOp::new(scope.ctx_mut(), start, end, step);
    let i = range_loop.iter_var(scope.ctx());
    let body_block = range_loop.loop_body(scope.ctx());
    let child = scope.loop_child(OpInserter::new_at_block_end(body_block));

    body(&child, i.into());
    child.terminate_yield();

    register_range_loop::<I>(scope, &range_loop, &child);
    scope.set_may_return(&[child]);
}

pub struct SteppedRangeExpand<I: Int> {
    start: NativeExpand<I>,
    end: NativeExpand<I>,
    step: NativeExpand<I>,
    inclusive: bool,
}

impl<I: Int + Into<ExpandValue>> Iterable for SteppedRangeExpand<I> {
    type Item = NativeExpand<I>;

    fn expand(self, scope: &Scope, body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType)) {
        let mut end = self.end;
        if self.inclusive {
            end = end.__expand_add_method(scope, I::new(1).into());
        }

        let start = self.start.read_value(scope);
        let end = end.read_value(scope);
        let step = self.step.read_value(scope);

        let range_loop = RangeLoopOp::new(scope.ctx_mut(), start, end, step);
        let i = range_loop.iter_var(scope.ctx());
        let body_block = range_loop.loop_body(scope.ctx());
        let child = scope.loop_child(OpInserter::new_at_block_end(body_block));

        body(&child, i.into());

        register_range_loop::<I>(scope, &range_loop, &child);
        scope.set_may_return(&[child]);
    }

    fn expand_unroll(
        self,
        scope: &Scope,
        body: &mut dyn FnMut(&Scope, <I as CubeType>::ExpandType),
    ) {
        let start = self
            .start
            .expand
            .as_const()
            .expect("Only constant start can be unrolled.")
            .as_i128();
        let end = self
            .end
            .expand
            .as_const()
            .expect("Only constant end can be unrolled.")
            .as_i128();
        let step = self
            .step
            .expand
            .as_const()
            .expect("Only constant step can be unrolled.")
            .as_i128();

        match (self.inclusive, step.is_negative()) {
            (true, true) => {
                for i in (end..=start).rev().step_by(step.unsigned_abs() as usize) {
                    let val = I::from_int_128(i);
                    body(scope, val.into())
                }
            }
            (true, false) => {
                for i in (start..=end).step_by(step.unsigned_abs() as usize) {
                    let val = I::from_int_128(i);
                    body(scope, val.into())
                }
            }
            (false, true) => {
                for i in (end..start).rev().step_by(step.unsigned_abs() as usize) {
                    let val = I::from_int_128(i);
                    body(scope, val.into())
                }
            }
            (false, false) => {
                for i in (start..end).step_by(step.unsigned_abs() as usize) {
                    let val = I::from_int_128(i);
                    body(scope, val.into())
                }
            }
        }
    }

    fn const_len(&self) -> Option<usize> {
        let start = self.start.constant()?.as_i128();
        let end = self.end.constant()?.as_i128();
        let step = self.step.constant()?.as_i128().unsigned_abs();
        Some((start.abs_diff(end) / step) as usize)
    }
}

/// integer range. Equivalent to:
///
/// ```ignore
/// start..end
/// ```
pub fn range<T: Int>(start: T, end: T) -> impl Iterator<Item = T> {
    let start: i64 = start.to_i64().unwrap();
    let end: i64 = end.to_i64().unwrap();
    (start..end).map(<T as NumCast>::from).map(Option::unwrap)
}

pub mod range {
    use cubecl_ir::Scope;

    use crate::prelude::{Int, NativeExpand};

    use super::RangeExpand;

    pub fn expand<I: Int>(
        _scope: &Scope,
        start: NativeExpand<I>,
        end: NativeExpand<I>,
    ) -> RangeExpand<I> {
        RangeExpand { start, end }
    }
}

/// Stepped range. Equivalent to:
///
/// ```ignore
/// (start..end).step_by(step)
/// ```
///
/// Allows using any integer for the step, instead of just usize
pub fn range_stepped<I: Int>(start: I, end: I, step: I) -> Box<dyn Iterator<Item = I>> {
    let start = start.to_i128().unwrap();
    let end = end.to_i128().unwrap();
    let step = step.to_i128().unwrap();

    if step < 0 {
        Box::new(
            (end..start)
                .rev()
                .step_by(step.unsigned_abs() as usize)
                .map(<I as NumCast>::from)
                .map(Option::unwrap),
        )
    } else {
        Box::new(
            (start..end)
                .step_by(step.unsigned_abs() as usize)
                .map(<I as NumCast>::from)
                .map(Option::unwrap),
        )
    }
}

pub mod range_stepped {
    use cubecl_ir::Scope;

    use crate::prelude::{Int, NativeExpand};

    use super::SteppedRangeExpand;

    pub fn expand<I: Int>(
        _scope: &Scope,
        start: NativeExpand<I>,
        end: NativeExpand<I>,
        step: NativeExpand<I>,
    ) -> SteppedRangeExpand<I> {
        SteppedRangeExpand {
            start,
            end,
            step,
            inclusive: false,
        }
    }
}